Beefy Boxes and Bandwidth Generously Provided by pair Networks
Don't ask to ask, just ask
 
PerlMonks  

Re: Lapack, inline and Perl

by rg0now (Chaplain)
on Mar 26, 2005 at 19:06 UTC ( [id://442538]=note: print w/replies, xml ) Need Help??


in reply to Lapack, inline and Perl

If your matrix is not sooo large, you could use a pure Perl solution, Math::MatrixReal for diagonalization instead of messing with Lapack.

Math::MatrixReal has a nice ($l, $V) = $matrix->sym_diagonalize() method that lets you to compute the eigensystem of a quadratic and symmetric matrix. Note that the implementation is really cool (it worked for me even for some astonishingly large matrices I threw at it just for fun): it uses Householder reduction followed by the QL algorithm, so it is in fact O(N^3)!!!

rg0now

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://442538]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others avoiding work at the Monastery: (4)
As of 2024-04-18 19:27 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found