need a little IMSL help, harwell-boeing sparse matrix features

need a little IMSL help, harwell-boeing sparse matrix features

Can I get help with use of the s_hbc_sparse features of IMSL on this forum?  Or should I go to roguewave forums?

I want to know if it is possible to do an LU factor, followed by many linear equation solves, with a sparse matrix using the s_hbc_sparse derived type and overloaded operators.

Thanks,

Brian in Austin, Texas

6 post / 0 nuovi
Ultimo contenuto
Per informazioni complete sulle ottimizzazioni del compilatore, consultare l'Avviso sull'ottimizzazione
Ritratto di Steve Lionel (Intel)

There are some experienced IMSL users here, but I would recommend asking in the Rogue Wave forums.

Steve

I posted at Roguewave last week.  So far no replies.  I gotta believe that if they call it SuperLU, it ought to be able to do what I want.  Isn't that what LU factoring is all about?  I suppose I'll have to make test program and take a stab at it.  I'll just be guessing as to how to get the LU to be saved.  I don't see anything in the IMSL documentation about that.

Brian

Ritratto di Steve Lionel (Intel)

Let me forward your question to my counterparts at RW and see what they say.

Steve
Ritratto di Steve Lionel (Intel)

Ok, got an answer:

Yes, he is correct. For the sparse operator using SuperLU, the LU  factorization is done on the first “call” to .ix. and saved for any subsequent calls for the same s_hbc_sparse definition.

 

Steve

Thanks, Steve.  I am preparing code right now that will try this.  It's in the iteration loop of an Arnoldi eigensolver.

Cheers,

Brian

Accedere per lasciare un commento.