2017-09-08 8:35 GMT-04:00 'Maxi Miller' via deal.II User Group <[email protected]>: > Nevertheless, it did not fix the problem... This won't fix any thing. I am just surprise that your code goes as far as the l2_norm(), I think that it should crash earlier. What happens if you comment these line
for (unsigned int i=0; i<dof_handler.n_locally_owned_dofs(); ++i) if (boundary_dofs[i] == true) residual(i) = 0; residual.compress(VectorOperation::insert);//Have to call it after setting the boundary elements Best, Bruno -- 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.
