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 Beiträge / 0 neu
Letzter Beitrag
Nähere Informationen zur Compiler-Optimierung finden Sie in unserem Optimierungshinweis.

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

Steve - Intel Developer Support

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

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

Steve - Intel Developer Support

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 - Intel Developer Support

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

Cheers,

Brian

Kommentar hinterlassen

Bitte anmelden, um einen Kommentar hinzuzufügen. Sie sind noch nicht Mitglied? Jetzt teilnehmen