lapack f90 routines or not?

Open discussion regarding features, bugs, issues, vendors, etc.

lapack f90 routines or not?

Postby alexatlapack » Tue Jul 04, 2006 6:04 am

Hi all,
we actually revise our dft simulation package code and try to get rid of some warnings and pitfalls.
our code is written in f90.
for some compilers it has been discussed, that they eventually make a copy of some matrices when calling lapack.
my question now is: does there exist some f90 interface which prevents from that? does a stable f90 interface exist at all? (this could be interesting in order to get a consistent code....)

thanks for any hint
alex
alexatlapack
 
Posts: 4
Joined: Sun Feb 19, 2006 12:59 pm

Postby Piotr Luszczek » Tue Jul 04, 2006 11:49 am

You can have a look at LAPACK95 and LAPACK3E. But from your
perspective they will be just thin layers on top of old LAPACK.

The problem with any kind of solution is the lack of BLAS interface
that would support Fortran 90. Imagine that the whole LAPACK was
rewritten in Fortran 90. As unlikely as it may be, it still doesn't
solve your problem. You potentially avoid a copy when calling
LAPACK but LAPACK will be making copies when calling BLAS (old
BLAS).

Until there is BLAS in Fortran 90 there is no point in having fancy
LAPACK interface. It's just better to make one copy when you call
LAPACK then to make multiple copies when you call BLAS inside LAPACK.

And considering the fact that it took some 30 years to make Fotran 77
BLAS a widely accepted standard, I would not bet much money on
having Fortran 90 BLAS any time soon.
Piotr Luszczek
 
Posts: 8
Joined: Wed Dec 08, 2004 8:23 pm
Location: Knoxville, TN, USA

Postby alexatlapack » Tue Jul 04, 2006 11:53 am

... that sounds reasonable...
alexatlapack
 
Posts: 4
Joined: Sun Feb 19, 2006 12:59 pm


Return to User Discussion

Who is online

Users browsing this forum: No registered users and 2 guests