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

[Michael Povolotskyi]

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.
Michael.