[dune-functions] Stokes/Taylor-Hood example

Oliver Sander sander at igpm.rwth-aachen.de
Tue Apr 14 17:06:11 CEST 2015


Hi all,

I have just extended the stokes-taylorhood example.  The linear system is now solved
with a GMRes with no preconditioner.  There is no preconditioner because there is no
preconditioner:  all existing dune-istl preconditioners either choke on Matrix<BCRSMatrix>
or expect the diagonal elements to be nonzero.  OTOH, currently the example uses a 4x4 grid.
That should work even without preconditioning...

Unfortunately, the program does not appear to converge against the correct solution.
Is there an existing Stokes/TH implementation in pdelab that I can use for comparison tests?

Cheers,
Oliver

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 473 bytes
Desc: OpenPGP digital signature
URL: <https://lists.dune-project.org/pipermail/dune-functions/attachments/20150414/39c37346/attachment.sig>


More information about the dune-functions mailing list