Malcolm Cohen wrote:
>
> Pierre Hugonnet said:
> > dot product was just an example. Consider routines which do
> > not have intrinsic equivalent (or even dot product on >1 dimension
> > arrays, or NRM2(),...).
>
UNI-C in Denmark (sorry, I don't have a URL handy) have written an F90
interface to lapack. I use their interface with the Digital Fortran dxml
library and it gives you proper F90 access to all of lapack. No more
temporary arrays that have to be passed or any of that junk. And it
takes care of a lot of the higher dimensional LA routines you should
need.
It must have been real drudgery to write it, but I think it is great. I
don't know whether it adds much runtime over-head, but the reduction in
programming effort makes it well worth it.
Adriaan
Oh, I found it. Its on the netlib mirror at www.hensa.ac.uk
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|