[Trilinos-Users] Amesos and complex numbers

[Heroux, Michael A]

Nik,

As long as your original matrix is well-conditioned, the answer you get back from the equivalent real form will be unique.

On a related note, I have completed a class that forms an equivalent real system using two Epetra matrices in much the same way that the Komplex package works.  I am testing it now and should have it in Trilinos 9.0.

I don't think we will have time to make Amesos work directly with Tpetra complex matrices before the Trilinos 9.0 release.

Mike


On 7/23/08 1:56 PM, "Nikhil Kriplani" <nmkripla at ncsu.edu> wrote:

Hi

I am trying to use Amesos LU factorization for solving a matrix that
consists of complex numbers. I've set up 2 real matrices using Epetra,
Areal and Aimag, and built up an augmented system as described on the
komplex web page.
http://trilinos.sandia.gov/packages/docs/dev/packages/komplex/doc/html/index.html#example
So essentially, I am doing an LU factorization on a real matrix of
twice the dimension.

My question is, after the LU factors are generated and the solution
generated, is this solution unique? I guess this is more a
mathematical question that a software one, but do I have to some
compensation/scaling to the solution with this scheme?

Also, when Tpetra complex number support comes out in Sep, will Amesos
LU solver be compatible with factorizing a complex matrix?

Thanks,
Nik

_______________________________________________
Trilinos-Users mailing list
Trilinos-Users at software.sandia.gov
http://software.sandia.gov/mailman/listinfo/trilinos-users


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20080724/815a69d5/attachment-0001.html