On 11/22/2013 06:07 PM, Jed Brown wrote:
Michael Povolotskyi <[email protected]> writes:
Dear developers of Petsc and users,
I have the following numerical problem.
I need to solve the system Ax = B, where A is a square matrix, x and B
are rectangular matrix
The system is solved by doing LU decomposition.
You can use MatLUFactor(), MatMatSolve(). There is no
MatMatSolveTranspose() at this time.
In the next step I need to solve a system A'y = C, where A' is a complex
conjugate of A.
My question: is it possible to reuse the LU decomposition of A for doing
this?
Not presently; you'd have to use MatHermitianTranspose() and refactor.
Support for reuse would not be very difficult to implement.
Thank you.
I'm actually interested in implementing it.
Could you, please, tell what has to be done in MatMatSolveTranspose()?
By transpose I mean transpose and complex conjugation.
Michael.
