Hi, On Thursday, September 7, 2017 at 5:01:10 PM UTC-4, Maxi Miller wrote: > > > for (unsigned int i=0; i<dof_handler.n_locally_owned_dofs(); ++i) > if (boundary_dofs[i] == true) > residual(i) = 0; > > This looks wrong. It shouldn't be residual(i) = 0; because there is only one residual(0) This code should throw in debug mode, are you running in release mode? If you do, I strongly suggest that you run in debug mode, you will catch early a lot of error.
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.
