Alasdair McAndrew wrote:
> Thank you for the solution of first compiling lapack as root - it worked
> fine. I had no trouble with this, as I've had to compile code as root for
> other software, so it can then be accessible by all users.
>
> My next problem is - are there any Maxima interfaces to the lapack
> routines? Suppose, for example, I wanted to calculate the eigenvalues
> numerically of a 50x50 matrix. Then:
>
> load(lapack);
> M:genmatrix(lambda([i,j],random(1.0)),50,50)$
> es:dgeev(M);
>
> But I have to know about the lapack function dgeev to do this. It
> seems to
> me that the lapack routines should be hidden from the average user, so
> that,
> for example
>
> eigenvalues(M),numer;
>
> (or whatever) should call dgeev automatically. In other words, if the user
I think the intent is to keep the LAPACK routine names so that users can
use LAPACK documentation to figure out what routine to use.
And yes, I think some higher level routines should take advantage of the
LAPACK routines, but no one I know of has tried to do that yet or even
what the interface should really be.
Ray