Re: [ublas] sparse iterative eigensolver

From: Sorkin Dmitri (dsorkin_at_[hidden])
Date: 2006-02-24 12:03:11

On Fri, 24 Feb 2006, Sorkin Dmitri wrote:
> Hi.
> I have some wrapper for ARPACK which also uses
> bindings for UMFPACK. I didn't benchmark the
> performance. It is pretty good for me.
>
> I attach the files that I use. There is lots
> of "bad" code.The main class is
> the "LAL_TypeEigDriverNonSym".
"LAL" means that I actually solve the generalized
eigen problem
   M*v*s = K*v
in the form
   (L^-1)*K*(U^-1)*x=x*s

where L*U = M (with umfpack)
x = U*v
"s" is the eigenvalue.

I don't know whether the method is
good or bad, I've just coded it.

Dima.