Dear Wolfgang Another issue related to this question is that, ApproximateMassSchur has a vmult member function defined, so it can be used in CG solver to get the product of its inverse with some vector, just like step-20 does. However unlike the dealii SolverCG class, PETScWrappers::SolverCG only takes MatrixBase type. In this case, is deriving my ApproximateMassSchur from PETScWrappers::MatrixBase with an overridden vmult function a correct solution? Or I simply should compute and store B*diag(M)^{-1}*B^T as a PETScWrappers::MPI::SparseMatrix?
Thanks Jie -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.