Re: using a lapack routine



I prefer to compile LAPACK and create a library. I then, create a
module with interfaces to LAPACK subroutines. For instance I have
created a function for the eigenvalues and eigenvectors of a real or a
real symmetric matrix. So when I want to calculate, say, the
eigenvectors of a real matrix A I type v=eigVectors(A), or if I want
the eigenvectors of a symmetric matrix A I type
v=eigVectors(A,sym="Y"). In this way the cumbersome argument list of
the LAPACK routines does not appear, or have to be used anywhere.

Michael
.



Relevant Pages

  • Re: Matrix Diagonalization
    ... LAPACK and virtually everything else get that one wrong! ... Run it, extracting all eigenvectors, and see :-) ... then the routine will ignore the lower triangle of the ... input and diagonalize a unit matrix of dimension 2x2. ...
    (comp.lang.fortran)
  • Re: Matrix Diagonalization
    ... LAPACK and virtually everything else get that one wrong! ... Run it, extracting all eigenvectors, and see :-) ... The all zeroes is not invalid, ... bug in the error detection; if the former, ...
    (comp.lang.fortran)
  • Re: minimization of a matrix function
    ... It appears that you want the d lowest eigenvectors of the ... I think you need to study eigenvalue problems and perturbation ... LAPACK documentation would be a good place to start, ... can find the eigenvectors by minimizing ...
    (sci.math.num-analysis)
  • Re: Matrix Diagonalization
    ... LAPACK and virtually everything else get that one wrong! ... Run it, extracting all eigenvectors, and see :-) ... The all zeroes is not invalid, ... bug in the error detection; if the former, ...
    (comp.lang.fortran)
  • Re: eigenstructure from eig
    ... > If A is a symmetric matrix, ... how does it choose the eigenvectors? ... You'd need to look at the properties of your matrix A, find which LAPACK ... you'll need to check the LAPACK routines and/or the ...
    (comp.soft-sys.matlab)

Loading