Hello,
I have two issues with the lapack implementation in maxima.
First of all, it goes without saying that many thanks are due to Ray
Toy for his porting of the package to maxima.
Issue 1: (not terribly important)
In maxima 5.13.0, I can't load lapack with the command
load(lapack);
Here is what happens: (%i1) build_info();
Maxima version: 5.13.0
Maxima build date: 21:33 9/7/2007
host type: i686-pc-linux-gnu
lisp-implementation-type: SBCL
lisp-implementation-version: 1.0.7
(%o1)
(%i2) load(lapack);
Maxima version: 5.13.0
Maxima build date: 21:33 9/7/2007
host type: i686-pc-linux-gnu
lisp-implementation-type: SBCL
lisp-implementation-version: 1.0.7
stdin:16:Incorrect syntax: : is not a prefix operator
As root, this problem does not happen. Also, if I put
load("/usr/local/share/maxima/5.13.0/share/lapack/lapack.mac");
in my maxima-init.mac file, then the problem disappears and lapack does
load.
Issue #2: documentation.
Is there anyway to get documentation for the maxima version of
'lapack' (other than reading the sources of course)?
For instance, I looked up the fortran version and one calls
dgeev with too many variable to remember. Maybe it works with fewer
variables, I did not try it since I don't use fortran that much anymore.
I noticed in the archives that
'dgeev(M)'
produces the eigenvalues of the matrix 'M'. So, that worked. But the
fortran version docs state that 'dgeevx(M)' gives the eigenvalues and
eigenvectors.
Maxima does not seem to know about 'dgeevx'--at least I don't know
how to call it.
Playing around, I found that
'dgeev(M,1) or 'dgeev(M,M)' --or any second variable character I
guess, gives the eigenvectors also.
I suspect that Ray is too busy to deal with the documentation of so
many packages.
Perhaps, is some users find things out (as I apparently did with
'dgeev(M,1)',
Is there anyway (except reading through the sources, of course) to find
out how to implement the lapack?
Any ideas?
TIA,
-sen
P.S. Is it hard to set up a 'maxima-user-list' at sourceforge? Maybe
it is time for two lists: one for development and bugs, and another for
users.