[Trilinos-Users] Help, Epetra and Complex Equivalent Problem
Thornquist, Heidi K
hkthorn at sandia.gov
Thu Dec 18 10:49:38 MST 2008
Hi Davood,
If you look at the Komplex_LinearProblem class in the komplex package, it seems to have what you want.
It is a general expression for a complex matrix that takes Epetra_CrsMatrix and Epetra_MultiVector objects as
arguments. You can then extract an Epetra_LinearProblem from this class and use it in another
Trilinos solver/preconditioner.
Thanks,
Heidi
On 12/18/08 10:24 AM, "Davood Ansari" <david.ansari at gmail.com> wrote:
Hi all
I have a complex matrix problem to solve (FEM matrices) :
(Ar + j Ai) x = b
where x = xr + j xi
and b=br + j bi
I have got the real and imaginary parts of the matrix A and the RHS b assembled in separate arrays (i.e. Ar Ai br bi).
I know that this is equal to a real matrix problem of twice the dimensions:
Ar -Ai
Ai Ar
xr
xi
=
br
bi
So I would like to construct this equivalent Epetra_LinearProblem without having to store the Ar and Ai (-Ai) twice for each.
I know that the komplex package would probably do this but I would rather use some thing based on Epetra since I already have
the matrices ready to go in Epetra_CrsMatrix objects.
Please Advise
Davood
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20081218/519da872/attachment.html
More information about the Trilinos-Users
mailing list