> > - When using the LinearOperator, the function vmult() is used when > solving the system with a GMRES-solver (or similar). Does the same hold up > for the MatrixFree-framework? > > Yes.
> > - If the above is correct: In each GMRES-iteration I have to calculate > the residual of my equation, once for the current solution and once for the > current solution + a small step. At the moment I calculate the residual by > summing the terms into my residual vector, and then distributing the values > while setting the boundary condition to 0. If I understand it correctly, > that happens as soon as I call phi.distribute_local_to_global(dst), with > phi a FEEvaluation-element. > > Yes. > > - This means that I first have to assemble the residual vector, then > distribute it, and then can use the vmult-function. Is that correct, or do > I have to use another approach? > > Isn't assembling and distributing already what `vmult` is supposed to do? Can you elaborate on what you are doing in each GMRES iteration? Best, Daniel -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAOYDWbLm%3D8p%3DY33AonJE837Q%2BxKibmBBRiN2e-6iexVu7r1UMw%40mail.gmail.com.
