Am Montag, 29. April 2019 22:36:55 UTC+2 schrieb Wolfgang Bangerth: > > On 4/29/19 7:35 AM, 'Maxi Miller' via deal.II User Group wrote: > > Could get it to work, now the code works equally fine for one or several > > MPI threads. > > What did you have to do? I was going to reproduce your problem today by > installing a Trilinos version that has NOX enabled. Is this moot now, > i.e., was it a bug in your code or did you just work around the issue in > some way that doesn't expose it? > > NOX expects vectors containing only the local elements, but no ghost elements. Thus I had to initialize all vectors going in or out from any NOX-related function using locally_owned_dofs, and copy accordingly if external vectors contained ghost elements.
> > > Still, the next open question is how to apply non-zero > > Dirichlet boundary conditions (as in the file, for cos(pi*x)*cos(pi*y). > > When applying them as I do in the file, I get the warning that the > > solver is not converged. When using Neumann-Conditions, the solver > > converges and I get the correct output. What would be the solution for > that? > > The question to ask is who is right. When you get the message that the > solver is not converged, is this correct? That is, does it converge if > you allow more iterations? What if you use a case where you know the > exact solution -- is the computed solution actively wrong? > The solver does not converge, and the output looks as if it is using Dirichlet-conditions with u = 0, independently of the "real" boundary conditions. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > <javascript:> > www: http://www.math.colostate.edu/~bangerth/ > -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
