So I am wanting to compare the performance of various FEM discretization with their respective "best" possible solver/pre conditioner. There are saddle-point systems which HDiv formulations like RT0 work, but then there are others like LSFEM that are naturally SPD and so the CG solver can be used (though finding a good preconditioner is still an open problem). <div><br></div><div>I have read and learned that the advantage of LSFEM is that it will always give you an SPD system, even for non-linear problems (because what you do is linearize the problem first and then minimize/take the Gateaux derivative to get the weak form). But after talking to some people and reading some stuff online, it seems one could also make non SPD systems SPD (hence eliminating what may be the only advantage of LSFEM).<br><br>Two of said people happen to be PETSc developers but I forgot to ask them how one would achieve that. Or if one really only can achieve S or PD and not both :)</div><div><br>On Saturday, November 28, 2015, Patrick Sanan <<a href="mailto:patrick.sanan@gmail.com">patrick.sanan@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On Sat, Nov 28, 2015 at 06:31:31AM -0600, Matthew Knepley wrote:<br>
> On Sat, Nov 28, 2015 at 12:10 AM, Justin Chang <<a href="javascript:;" onclick="_e(event, 'cvml', 'jychang48@gmail.com')">jychang48@gmail.com</a>> wrote:<br>
><br>
> > Hi all,<br>
> ><br>
> > Say I have a saddle-point system for the mixed-poisson equation:<br>
> ><br>
> > [I  -grad] [u]  = [0]<br>
> > [-div  0  ] [p]     [-f]<br>
> ><br>
> > The above is symmetric but indefinite. I have heard that one could make<br>
> > the above symmetric and positive definite (SPD). How would I do that? And<br>
> > if that's the case, would this allow me to use CG instead of GMRES?<br>
> ><br>
><br>
> I believe you just multiply the bottom row by -1. You can use CG for an SPD<br>
> system, but you can<br>
> use MINRES for symmetric indefinite.<br>
If I'm remembering correctly, flipping that sign lets you make your system alternately P.D. or<br>
symmetric, but not both. Maybe you were hearing about the Bramble-Pasciak preconditioner or a related approach?<br>
><br>
>    Matt<br>
><br>
><br>
> > Thanks,<br>
> > Justin<br>
> ><br>
><br>
><br>
><br>
> --<br>
> What most experimenters take for granted before they begin their<br>
> experiments is infinitely more interesting than any results to which their<br>
> experiments lead.<br>
> -- Norbert Wiener<br>
</blockquote></div>