> This generic works on the basis of both rank and the extra
> arguments. (Either would be sufficient by itself to make the generic
> unambiguous.) Array ranks higher than 1 have to go through an extra level
> of procedure reference, but there should be no copying of the array, so
> that cost should be small and the subroutine dummy1 should be just as
> efficient as Aleksandar wants. The nuisance is the necessity of creating
> all those pass through routines for each dimension above 1.
Thanks,
This makes sense, and the overhead should be low, since the compiler knows
the array is already contiguous and should not generate additional copy-in,
copy-out overheads.
I will add these interfaces to my module (I use a powerful preprocessor
which makes making those dummy wrapper routines easy--more on this after the
weekend!).
Aleksandar
_____________________________________________
Aleksandar Donev
http://www.pa.msu.edu/~donev/
[log in to unmask]
(517) 432-6770
Department of Physics and Astronomy
Michigan State University
East Lansing, MI 48824-1116
_____________________________________________
|