comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Lanczos diagonalisation with IDL
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
Re: Lanczos diagonalisation with IDL [message #51811] Thu, 07 December 2006 19:07 Go to next message
Kenneth P. Bowman is currently offline  Kenneth P. Bowman
Messages: 585
Registered: May 2000
Senior Member
In article <1165497198.963840.266160@l12g2000cwl.googlegroups.com>,
philip.tregenna@psi.ch wrote:

> I would like to diagonalise a large, sparse, symmetric matrix, and find
> a limited set of eigenvalues. The best way to do this is to use an
> algorithm based upon the Lanczos method. Has a Lanczos diagonalisation
> routine been written in the IDL language? The procedure with a
> functionality close to that which I need is "LA_EIGENQL". However,
> this routine is not ideal for large sparse matrices as it requires as
> input all the elements of the matrix, including those with a value of
> zero.

Haven't used it, but this looks like a job for IDL Analyst (aka IMSL)

http://www.ittvis.com/idl/pdfs/IDL63_Analyst_RoutinesFunc.pd f

It includes sparse matrix routines. It costs extra, but consider the
cost of your time.

Ken Bowman
Re: Lanczos diagonalisation with IDL [message #51812 is a reply to message #51811] Thu, 07 December 2006 20:13 Go to previous message
Richard French is currently offline  Richard French
Messages: 173
Registered: December 2000
Senior Member
You might check out the IMSL routines in IDL Analyst - extra bucks for this
software, I know, but it might be worth it. For example, IMSL_SP_LUSOL
solves a sparse system of linear equations without having to specify the
full sparse matrix. IMSL_SP_LUFAC computes an LU factorization of a sparse
matrix without having to store the full matrix. The documentation for IDL
Analyst is available on the RSI web site someplace.

Dick French

On 12/7/06 8:13 AM, in article
1165497198.963840.266160@l12g2000cwl.googlegroups.com,
"philip.tregenna@psi.ch" <philip.tregenna@psi.ch> wrote:

> I would like to diagonalise a large, sparse, symmetric matrix, and find
> a limited set of eigenvalues. The best way to do this is to use an
> algorithm based upon the Lanczos method. Has a Lanczos diagonalisation
> routine been written in the IDL language? The procedure with a
> functionality close to that which I need is "LA_EIGENQL". However,
> this routine is not ideal for large sparse matrices as it requires as
> input all the elements of the matrix, including those with a value of
> zero.
>
> Regards,
>
> Philip.
>
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Re: false color display
Next Topic: histogram along y-axis

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 14:32:36 PDT 2025

Total time taken to generate the page: 1.52134 seconds