[petsc-users] MUMPS with symmetric matrices
Gianluca Meneghello
gianmail at gmail.com
Wed Nov 4 18:18:33 CST 2015
I have just read that there is no special algorithm for Hermitian matrices
in MUMPS (sorry, I meant Hermitian, not symmetric... the matrix is complex).
Sorry for this. In any case, if there is any suggestion it is more than
welcome!
Thanks for your help and your work,
Gianluca
On Wed, Nov 4, 2015 at 3:37 PM, Gianluca Meneghello <gianmail at gmail.com>
wrote:
> It is a discretization of the differential operator, of which I would need
> the inverse (or LU decomposition). My goal is frequency response
> (resolvant) analysis of the linearized Navier-Stokes operator.
>
> There was a reason I was not using MatTransposeMatMult, that is the
> matrix is complex and I would need MatTransposeHemitianMatMult (or
> something like that). It seems to me that is not available (or does
> MatTransposeMatMult compute the Hermitian transpose?)
>
> Any suggestion is of course welcome!
>
> Thanks
>
> Gianluca
>
>
>
>
>
> On Wed, Nov 4, 2015 at 1:16 PM, Jed Brown <jed at jedbrown.org> wrote:
>
>> Gianluca Meneghello <gianmail at gmail.com> writes:
>>
>> > That is correct... I will try with -pc_type cholesky and use
>> > MatTransposeMatMult.
>> >
>> > Using cholesky I do not need to specify mumps as a solver, am I right?
>>
>> Of course you do.
>>
>> > A is a linearization of the Navier Stokes equation.
>>
>> Of the differential operator, its inverse, or a map from some parameters
>> to observations?
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151104/87324448/attachment-0001.html>
More information about the petsc-users
mailing list