Argument Passing to F77_LOCAL Procedures

A typical Fortran 77 implementation passes arguments by reference, usually by passing the base address of the location of the first data element, and such arguments may also be assumed to be sequence associated. These facts make it most practical for the default method of passing a distributed data structure from HPF to an F77_LOCAL procedure by passing the base address of that section of local memory that has been allocated to it. To allow for sequence association of actual and dummy arguments, data should also be reordered or compressed or both, if necessary, on all processors. This is the safest method of passing distributed data to an EXTRINSIC(F77_LOCAL) procedure, and hence it should be the default one. However, it tends to have the greatest performance costs.

A second argument passing option is to pass distributed array data ``as is'' from a global HPF procedure to the local F77 ones, not guaranteeing sequence association of the dummy arguments in order to avoid unwanted local data motion that might be required to compress or reorder the elements of an array local to a processor. In other words, it should be possible to do no more data motion than if the same argument were being passed to another HPF procedure. The guarantee of a sequence associated dummy argument is sacrificed for the possible gains in performance available because the local components of the actual argument are not reordered or compressed. The local programmer must be able to use the implementation-dependent ordering created by the global HPF program.

A third option that can be useful to permit HPF_LOCAL-style local programming from an EXTRINSIC(F77_LOCAL) procedure call is to pass an array via a descriptor or handle, as is typically done in HPF implementations or for Fortran 90 assumed shape arrays. The local procedure may not access elements of this dummy argument directly but may only pass it on to special utility routines, perhaps to obtain local or global distribution information.

The following attributes suffice to support the above three alternate form of passing data to an EXTRINSIC(F77_LOCAL)x procedure:

• LAYOUT("0DF77_ARRAY"0D) indicates that the rectangular configuration should be FORTRAN 77 sequence associated in local, serial memory.

  For example, many compilers add border elements for ``stencil'' optimizations or pad array allocations on particular processors so that all processors allocate equal amounts of memory for each array. Local reordering eliminates such padding and provides FORTRAN 77 sequence association for actual data values.

  Any local reordering is in addition to any global remapping that may be dictated by DISTRIBUTION or ALIGN directives in the INTERFACE block.

  If no LAYOUT attribute is specified, then LAYOUT("0DF77_ARRAY"0D) is assumed.

• LAYOUT("0DHPF_ARRAY"0D) indicates that an array argument is passed just as it would be to a global HPF procedure, with no local reordering of the data.

  This option is desirable when the programmer decides that the overhead of local reordering should be eliminated or that certain characteristics of the global HPF compiler's ordering (border cells, equal-size allocations among processors, etc.) should be preserved at the local level. It forces the local programmer to access the local data in whatever implementation-dependent style the global HPF compiler employs.

Furthermore, each argument in the INTERFACE block can also have its PASS_BY attribute specified to indicate whether the data is passed by reference, for Fortran 77-style access, or via a special handle, perhaps a descriptor used for HPF variable passing, that permits the global HPF caller to pass special mapping information for use within the local Fortran 77 procedure.

• PASS_BY("0D*"0D) indicates that the local procedure should be able to access the dummy argument locally as an F77-style variable, passed by reference.
• PASS_BY("0DHPF_HANDLE"0D) indicates that the local procedure should receive a reference to a global descriptor that can be used with special inquiry routines to obtain useful distribution information.

Thus, the default dummy argument attributes are LAYOUT("0DF77_ARRAY"0D), a guarantee of sequence association, and PASS_BY("0D*"0D), an indication that data is being passed via a pointer to its location.

Advice to implementors.In addition to providing argument passing and data reordering options, a good EXTRINSIC(F77_LOCAL) implementation should address the problem of declaring arbitrary sized local subgrids and accessing their elements without being able to describe them as assumed-shape arrays, as in HPF. Dealing with the local results of global data distributions within each local procedure initiated by an extrinsic procedure call can also be difficult without Fortran 90 array inquiry functions and the inquiry subroutines in the HPF Library. Special inquiry routines, callable globally or locally, such as the proposed library of Fortran 77 function interfaces in Annex G are recommended as supplements to the EXTRINSIC(F77_LOCAL) procedure interface in order to permit more flexible and efficient use of a broad range of possible global HPF data distributions.

(End of advice to implementors.)