June 2015 - mpiwg-fortran - mailman3.cels.anl.gov

Re: [MPIWG Fortran] What compilers support TS 29113?
by Jeff Squyres (jsquyres) 23 Jun '15

23 Jun '15

On Jun 23, 2015, at 11:51 AM, Bill Long <longb(a)cray.com> wrote: > > gfortran - yes (4.9 works on my Mac) I thought fortran 4.9 did not have the standardized TS 29113 dope vector...? > Cray - yes (8.3 for sure, or later). I knew about Cray, but I thought it was the only one. That may be incorrect and/or outdated information. > Not sure about others. People who build libraries (mpich, OpenMPI, ...) would know based on which ones can compiler their mpi_f08 module. Open MPI currently only has the non-TS-29113 interfaces for mpi_f08. We're looking into whether it's time to implement the TS 29113 interfaces; that's the origin of my question. :-) -- Jeff Squyres jsquyres(a)cisco.com For corporate legal information go to: http://www.cisco.com/web/about/doing_business/legal/cri/

1 0

Re: [MPIWG Fortran] What compilers support TS 29113?
by Bill Long 23 Jun '15

23 Jun '15

On Jun 23, 2015, at 1:44 PM, Jeff Squyres (jsquyres) <jsquyres(a)cisco.com> wrote: > I have not been following the development of modern Fortran compilers recently... > > Have any new Fortran compiler releases added support for TS 29113? > > E.g., does gfortran 5.0 support it? Ifort 15? …? gfortran - yes (4.9 works on my Mac) Cray - yes (8.3 for sure, or later). Not sure about others. People who build libraries (mpich, OpenMPI, ...) would know based on which ones can compiler their mpi_f08 module. Cheers, Bill > > -- > Jeff Squyres > jsquyres(a)cisco.com > For corporate legal information go to: http://www.cisco.com/web/about/doing_business/legal/cri/ > > _______________________________________________ > mpiwg-fortran mailing list > mpiwg-fortran(a)lists.mpi-forum.org > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran Bill Long longb(a)cray.com Fortran Technical Support & voice: 651-605-9024 Bioinformatics Software Development fax: 651-605-9142 Cray Inc./ Cray Plaza, Suite 210/ 380 Jackson St./ St. Paul, MN 55101

1 0

What compilers support TS 29113?
by Jeff Squyres (jsquyres) 23 Jun '15

23 Jun '15

I have not been following the development of modern Fortran compilers recently... Have any new Fortran compiler releases added support for TS 29113? E.g., does gfortran 5.0 support it? Ifort 15? ...? -- Jeff Squyres jsquyres(a)cisco.com For corporate legal information go to: http://www.cisco.com/web/about/doing_business/legal/cri/

1 0

Re: [MPIWG Fortran] Question about named constants
by William Gropp 19 Jun '15

19 Jun '15

Given the intent of the wording with respect to C, I don’t think the test is valid (in C, the clear intent was that handles could be pointers resolved only at link time). Bill On Jun 18, 2015, at 10:19 PM, Bill Long <longb(a)cray.com> wrote: > > On Jun 18, 2015, at 5:52 PM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: > >> In OpenMPI F08, those predefined MPI datatypes are not "parameter variables". Thus OMPI failed to compile an MPICH F08 test, which has >> >> INTEGER, PARAMETER :: nb=2 >> TYPE(MPI_DATATYPE) :: types(nb)=(/MPI_DOUBLE_PRECISION, MPI_CHAR/) >> >> Jeff Squyres and I discussed in emails. I don't know whether this is correct according to the MPI spec. > > Given the wording in the spec, I’m not sure either. BTW, there is no such thing as a “parameter variable”. Something with the parameter attribute is a constant, and hence not a variable. It is true that an initialization expression in a type declaration is required to be a constant expression. If OMPI is declaring these datatypes as variables instead of constants, then there would be a problem compiling the above code using their module. > > Cheers, > Bill > > > >> >> >> >> --Junchao Zhang >> >> On Thu, Jun 18, 2015 at 12:48 PM, Bill Long <longb(a)cray.com> wrote: >> >> On Jun 18, 2015, at 11:05 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: >> >>> Hello, Fortran experts, >>> I am a little puzzled by MPI-3.1 Chapter 2.5.4, >> >> Not surprising. The (mis)use of terminology is confusing at best. >> >>> >>> All named constants, with the exceptions noted below for Fortran, can be used in >>> initialization expressions or assignments, but not necessarily in array declarations or as >>> labels in C switch or Fortran select/case statements. This implies named constants >>> to be link-time but not necessarily compile-time constants. >> >> In Fortran “named constant” is a defined term. They are declared with the PARAMETER attribute. And can be used in most places where literal constants can be used, including array declarations. They cannot be used as statement labels, though I think the reference to “Fortran select/case statements” intends to talk about case-values, and not labels. Named constants can be used as case-values. >> >>> >>> Then in F08, should MPI_INTEGER be declared as a parameter variable? >> >> Yes. That’s what I did in my prototype mpi_f08 modules. >> >>> If no, then I get an error when compiling this initialization expression with gfortran >>> type(MPI_Datatype) :: dtype = MPI_INTEGER >>> >>> Error: Parameter 'mpi_integer' at (1) has not been declared or is a variable, which does not reduce to a constant expression >>> >>> If yes, then of course I can use MPI_INTEGER in Fortran select/case. >> >> Right. >> >>> >>> Another question is what is the meaning of "link-time constants". I could not find in MPI. Does it mean I can use MPI_INTEGER after main before MPI_Init? >> >> I suspect this is a historical artifact where some special head standing was needed to correctly define MPI_INTEGER based on whether the code was compiled with -i4 or -i8 (or equivalent, depending on the compiler). For the modern case of using the mpi_f08 module, you just compile the mpi_f08 module twice (once with -i4 and once with -i8) and then have the compile command select the right actual module to reference based on how the using program unit is compiled. Generally, values of objects are not set at link time; addresses of objects are set at link time. MPI objects that represent addresses (MPI_BOTTOM ??) are examples. If MPI_INTEGER is accessible, you can reference it. If should not have anything to do with whether MPI_Init has been called. >> >> Perhaps some of the MPI historians can shed more light on this terminology. I suspect it would be completely moot if we got rid of mpif.h. >> >> Cheers, >> Bill >> >> >>> >>> Thank you! >>> >>> --Junchao Zhang >>> _______________________________________________ >>> mpiwg-fortran mailing list >>> mpiwg-fortran(a)lists.mpi-forum.org >>> http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran >> >> Bill Long longb(a)cray.com >> Fortran Technical Support & voice: 651-605-9024 >> Bioinformatics Software Development fax: 651-605-9142 >> Cray Inc./ Cray Plaza, Suite 210/ 380 Jackson St./ St. Paul, MN 55101 >> >> >> _______________________________________________ >> mpiwg-fortran mailing list >> mpiwg-fortran(a)lists.mpi-forum.org >> http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran >> >> _______________________________________________ >> mpiwg-fortran mailing list >> mpiwg-fortran(a)lists.mpi-forum.org >> http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran > > Bill Long longb(a)cray.com > Fortran Technical Support & voice: 651-605-9024 > Bioinformatics Software Development fax: 651-605-9142 > Cray Inc./ Cray Plaza, Suite 210/ 380 Jackson St./ St. Paul, MN 55101 > > > _______________________________________________ > mpiwg-fortran mailing list > mpiwg-fortran(a)lists.mpi-forum.org > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran

1 0

Re: [MPIWG Fortran] Question about named constants
by Bill Long 18 Jun '15

18 Jun '15

On Jun 18, 2015, at 5:52 PM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: > In OpenMPI F08, those predefined MPI datatypes are not "parameter variables". Thus OMPI failed to compile an MPICH F08 test, which has > > INTEGER, PARAMETER :: nb=2 > TYPE(MPI_DATATYPE) :: types(nb)=(/MPI_DOUBLE_PRECISION, MPI_CHAR/) > > Jeff Squyres and I discussed in emails. I don't know whether this is correct according to the MPI spec. Given the wording in the spec, I’m not sure either. BTW, there is no such thing as a “parameter variable”. Something with the parameter attribute is a constant, and hence not a variable. It is true that an initialization expression in a type declaration is required to be a constant expression. If OMPI is declaring these datatypes as variables instead of constants, then there would be a problem compiling the above code using their module. Cheers, Bill > > > > --Junchao Zhang > > On Thu, Jun 18, 2015 at 12:48 PM, Bill Long <longb(a)cray.com> wrote: > > On Jun 18, 2015, at 11:05 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: > > > Hello, Fortran experts, > > I am a little puzzled by MPI-3.1 Chapter 2.5.4, > > Not surprising. The (mis)use of terminology is confusing at best. > > > > > All named constants, with the exceptions noted below for Fortran, can be used in > > initialization expressions or assignments, but not necessarily in array declarations or as > > labels in C switch or Fortran select/case statements. This implies named constants > > to be link-time but not necessarily compile-time constants. > > In Fortran “named constant” is a defined term. They are declared with the PARAMETER attribute. And can be used in most places where literal constants can be used, including array declarations. They cannot be used as statement labels, though I think the reference to “Fortran select/case statements” intends to talk about case-values, and not labels. Named constants can be used as case-values. > > > > > Then in F08, should MPI_INTEGER be declared as a parameter variable? > > Yes. That’s what I did in my prototype mpi_f08 modules. > > > If no, then I get an error when compiling this initialization expression with gfortran > > type(MPI_Datatype) :: dtype = MPI_INTEGER > > > > Error: Parameter 'mpi_integer' at (1) has not been declared or is a variable, which does not reduce to a constant expression > > > > If yes, then of course I can use MPI_INTEGER in Fortran select/case. > > Right. > > > > > Another question is what is the meaning of "link-time constants". I could not find in MPI. Does it mean I can use MPI_INTEGER after main before MPI_Init? > > I suspect this is a historical artifact where some special head standing was needed to correctly define MPI_INTEGER based on whether the code was compiled with -i4 or -i8 (or equivalent, depending on the compiler). For the modern case of using the mpi_f08 module, you just compile the mpi_f08 module twice (once with -i4 and once with -i8) and then have the compile command select the right actual module to reference based on how the using program unit is compiled. Generally, values of objects are not set at link time; addresses of objects are set at link time. MPI objects that represent addresses (MPI_BOTTOM ??) are examples. If MPI_INTEGER is accessible, you can reference it. If should not have anything to do with whether MPI_Init has been called. > > Perhaps some of the MPI historians can shed more light on this terminology. I suspect it would be completely moot if we got rid of mpif.h. > > Cheers, > Bill > > > > > > Thank you! > > > > --Junchao Zhang > > _______________________________________________ > > mpiwg-fortran mailing list > > mpiwg-fortran(a)lists.mpi-forum.org > > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran > > Bill Long longb(a)cray.com > Fortran Technical Support & voice: 651-605-9024 > Bioinformatics Software Development fax: 651-605-9142 > Cray Inc./ Cray Plaza, Suite 210/ 380 Jackson St./ St. Paul, MN 55101 > > > _______________________________________________ > mpiwg-fortran mailing list > mpiwg-fortran(a)lists.mpi-forum.org > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran > > _______________________________________________ > mpiwg-fortran mailing list > mpiwg-fortran(a)lists.mpi-forum.org > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran Bill Long longb(a)cray.com Fortran Technical Support & voice: 651-605-9024 Bioinformatics Software Development fax: 651-605-9142 Cray Inc./ Cray Plaza, Suite 210/ 380 Jackson St./ St. Paul, MN 55101

1 0

Re: [MPIWG Fortran] Question about named constants
by William Gropp 18 Jun '15

18 Jun '15

This has more to do with whether things like MPI_INTEGER are compile-time constants (as they are in MPICH) or pointers to memory that are determined when the application is linked with the library. Bill On Jun 18, 2015, at 12:48 PM, Bill Long <longb(a)cray.com> wrote: > > > Another question is what is the meaning of "link-time constants". I could not find in MPI. Does it mean I can use MPI_INTEGER after main before MPI_Init? >

1 0

Re: [MPIWG Fortran] Question about named constants
by Junchao Zhang 18 Jun '15

18 Jun '15

In OpenMPI F08, those predefined MPI datatypes are not "parameter variables". Thus OMPI failed to compile an MPICH F08 test, which has INTEGER, PARAMETER :: nb=2 TYPE(MPI_DATATYPE) :: types(nb)=(/MPI_DOUBLE_PRECISION, MPI_CHAR/) Jeff Squyres and I discussed in emails. I don't know whether this is correct according to the MPI spec. --Junchao Zhang On Thu, Jun 18, 2015 at 12:48 PM, Bill Long <longb(a)cray.com> wrote: > > On Jun 18, 2015, at 11:05 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: > > > Hello, Fortran experts, > > I am a little puzzled by MPI-3.1 Chapter 2.5.4, > > Not surprising. The (mis)use of terminology is confusing at best. > > > > > All named constants, with the exceptions noted below for Fortran, can be > used in > > initialization expressions or assignments, but not necessarily in array > declarations or as > > labels in C switch or Fortran select/case statements. This implies > named constants > > to be link-time but not necessarily compile-time constants. > > In Fortran “named constant” is a defined term. They are declared with the > PARAMETER attribute. And can be used in most places where literal > constants can be used, including array declarations. They cannot be used > as statement labels, though I think the reference to “Fortran select/case > statements” intends to talk about case-values, and not labels. Named > constants can be used as case-values. > > > > > Then in F08, should MPI_INTEGER be declared as a parameter variable? > > Yes. That’s what I did in my prototype mpi_f08 modules. > > > If no, then I get an error when compiling this initialization > expression with gfortran > > type(MPI_Datatype) :: dtype = MPI_INTEGER > > > > Error: Parameter 'mpi_integer' at (1) has not been declared or is a > variable, which does not reduce to a constant expression > > > > If yes, then of course I can use MPI_INTEGER in Fortran select/case. > > Right. > > > > > Another question is what is the meaning of "link-time constants". I > could not find in MPI. Does it mean I can use MPI_INTEGER after main before > MPI_Init? > > I suspect this is a historical artifact where some special head standing > was needed to correctly define MPI_INTEGER based on whether the code was > compiled with -i4 or -i8 (or equivalent, depending on the compiler). For > the modern case of using the mpi_f08 module, you just compile the mpi_f08 > module twice (once with -i4 and once with -i8) and then have the compile > command select the right actual module to reference based on how the using > program unit is compiled. Generally, values of objects are not set at > link time; addresses of objects are set at link time. MPI objects that > represent addresses (MPI_BOTTOM ??) are examples. If MPI_INTEGER is > accessible, you can reference it. If should not have anything to do with > whether MPI_Init has been called. > > Perhaps some of the MPI historians can shed more light on this > terminology. I suspect it would be completely moot if we got rid of mpif.h. > > Cheers, > Bill > > > > > > Thank you! > > > > --Junchao Zhang > > _______________________________________________ > > mpiwg-fortran mailing list > > mpiwg-fortran(a)lists.mpi-forum.org > > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran > > Bill Long > longb(a)cray.com > Fortran Technical Support & voice: > 651-605-9024 > Bioinformatics Software Development fax: 651-605-9142 > Cray Inc./ Cray Plaza, Suite 210/ 380 Jackson St./ St. Paul, MN 55101 > > > _______________________________________________ > mpiwg-fortran mailing list > mpiwg-fortran(a)lists.mpi-forum.org > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran >

1 0

Re: [MPIWG Fortran] Question about named constants
by Bill Long 18 Jun '15

18 Jun '15

On Jun 18, 2015, at 11:05 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: > Hello, Fortran experts, > I am a little puzzled by MPI-3.1 Chapter 2.5.4, Not surprising. The (mis)use of terminology is confusing at best. > > All named constants, with the exceptions noted below for Fortran, can be used in > initialization expressions or assignments, but not necessarily in array declarations or as > labels in C switch or Fortran select/case statements. This implies named constants > to be link-time but not necessarily compile-time constants. In Fortran “named constant” is a defined term. They are declared with the PARAMETER attribute. And can be used in most places where literal constants can be used, including array declarations. They cannot be used as statement labels, though I think the reference to “Fortran select/case statements” intends to talk about case-values, and not labels. Named constants can be used as case-values. > > Then in F08, should MPI_INTEGER be declared as a parameter variable? Yes. That’s what I did in my prototype mpi_f08 modules. > If no, then I get an error when compiling this initialization expression with gfortran > type(MPI_Datatype) :: dtype = MPI_INTEGER > > Error: Parameter 'mpi_integer' at (1) has not been declared or is a variable, which does not reduce to a constant expression > > If yes, then of course I can use MPI_INTEGER in Fortran select/case. Right. > > Another question is what is the meaning of "link-time constants". I could not find in MPI. Does it mean I can use MPI_INTEGER after main before MPI_Init? I suspect this is a historical artifact where some special head standing was needed to correctly define MPI_INTEGER based on whether the code was compiled with -i4 or -i8 (or equivalent, depending on the compiler). For the modern case of using the mpi_f08 module, you just compile the mpi_f08 module twice (once with -i4 and once with -i8) and then have the compile command select the right actual module to reference based on how the using program unit is compiled. Generally, values of objects are not set at link time; addresses of objects are set at link time. MPI objects that represent addresses (MPI_BOTTOM ??) are examples. If MPI_INTEGER is accessible, you can reference it. If should not have anything to do with whether MPI_Init has been called. Perhaps some of the MPI historians can shed more light on this terminology. I suspect it would be completely moot if we got rid of mpif.h. Cheers, Bill > > Thank you! > > --Junchao Zhang > _______________________________________________ > mpiwg-fortran mailing list > mpiwg-fortran(a)lists.mpi-forum.org > http://lists.mpi-forum.org/mailman/listinfo.cgi/mpiwg-fortran Bill Long longb(a)cray.com Fortran Technical Support & voice: 651-605-9024 Bioinformatics Software Development fax: 651-605-9142 Cray Inc./ Cray Plaza, Suite 210/ 380 Jackson St./ St. Paul, MN 55101

1 0

Question about named constants
by Junchao Zhang 18 Jun '15

18 Jun '15

Hello, Fortran experts, I am a little puzzled by MPI-3.1 Chapter 2.5.4, All named constants, with the exceptions noted below for Fortran, can be used in initialization expressions or assignments, but not necessarily in array declarations or as labels in C switch or Fortran select/case statements. This implies named constants to be link-time but not necessarily compile-time constants. Then in F08, should MPI_INTEGER be declared as a parameter variable? If no, then I get an error when compiling this initialization expression with gfortran type(MPI_Datatype) :: dtype = MPI_INTEGER Error: Parameter 'mpi_integer' at (1) has not been declared or is a variable, which does not reduce to a constant expression If yes, then of course I can use MPI_INTEGER in Fortran select/case. Another question is what is the meaning of "link-time constants". I could not find in MPI. Does it mean I can use MPI_INTEGER after main before MPI_Init? Thank you! --Junchao Zhang

1 0