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 <[email protected]> wrote: > Gianluca Meneghello <[email protected]> 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? >
