[petsc-users] question about LU decomposition of a conjugated matrix

Michael Povolotskyi mpovolot at purdue.edu
Fri Nov 22 07:18:18 CST 2013

On 11/22/2013 06:07 PM, Jed Brown wrote:
> Michael Povolotskyi <mpovolot at purdue.edu> 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.

More information about the petsc-users mailing list