Dear Isabel, > so I mean I would have to use "TrilinosWrappers::PreconditionStokes > Class", wouldn't I? > > And this class is in the developing part of the library DealII, isn't it?
Yes, you can give it a try (but I'm not sure whether this class will remain in the library in the current form also in the future, since it is very specific, whereas the library is rather meant to be general). As the author of that piece of code I have to admit that it is not really tested (actually I only tested the Stokes part rigorously, but not the necessary changes for Navier-Stokes, so I'm pretty sure there still are some bugs). You will need to have the matrix in TrilinosWrappers::BlockSparseMatrix format and the vector as TrilinosWrappers::BlockVector. And then there need to be some more matrices around, like a pressure mass, Laplace and convection-diffusion matrix (the last one you build just as the Navier-Stokes velocity matrix, but on the scalar pressure space)... So just try it and report me any problems! Best, Martin _______________________________________________
