LAPACK interfaces

[PierU]

@RonShepard

The extra copy is required because ultimately there is a dummy argument that is assumed size or explicit shape. That is because of the limitation to f77 argument conventions in the working routines. If those routines in LAPACK had been updated to use, for example, assumed shape arrays, then many (or all) of the copies that are currently required would be eliminated.

How does a code based on assumed shape arguments compare to a code based on INC* and LD* along with assumed size arguments (current LAPACK), in terms of performance?