On 8/31/22 11:09, Lucas Myers wrote:
Thanks for the help! I printed out the number of constraints after each
of your suggestions, and after setting the hanging node constraints the
number of constraints is 0 (as it should be, I'm working with a 1D mesh
that has been globally refined some number
Thanks for the help! I printed out the number of constraints after each of
your suggestions, and after setting the hanging node constraints the number
of constraints is 0 (as it should be, I'm working with a 1D mesh that has
been globally refined some number of times), but after interpolating th
On 8/30/22 13:56, Lucas Myers wrote:
To enforce these boundary conditions, I initialize constraints with:
```
constraints.clear();
dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
dealii::VectorTools::
interpolate_boundary_values(dof_handler,
Hi everyone,
I'm trying to solve a nonlinear ODE with Dirichlet conditions using
Newton's method. For this I give an initial guess for the solution ( y =
(1/2) x ) and then update each Newton step with an update that ought to be
zero at the boundaries (hence preserving the boundary conditions o
