Sure...  The constraints are applied to an *element* residual vector or Jacobian matrix much in the same way they are applied for the linear case in ex10.  So, if you assemble the residual and Jacobian element-by-element you simply need to call dof_map.constrain_element_matrix(...) and dof_map.constrain_element_vector(...) before inserting them into the global residual/Jacobian.
As I mentioned, the other missing piece is the ability to take an arbitrary global vector and extract the element solution vector from it, accounting for constraints.  I'll implement that shortly and check it it.
I've got a modified version of example 13 which illustrates the use of the nonlinear solver, but it does not address adaptivity.  I'm leaning instead to creating a new example altogether...
-----Original Message-----
From: Michael Povolotskyi []
Sent: Friday, June 10, 2005 9:02 AM
Subject: Re: [Libmesh-users] Nonlinear hanging nodes

Thank you for the answer!
Could you, please, tell more details on how to apply these constraints to the Jacobian matrix and to the  residual vector.


I plan to add such a method to the DofMap after 0.5.0 is released (which should be tonight).  The interface is something like this:
DofMap::build_elem_solution (global_solution, dof_indices, local_solution);
For non-constrained Dofs this will simply set
local_solution[i] = global_solution[dof_indices[i]]
And for constrained Dofs the constraint matrix will be built and applied to the solution vector.
Thanks for your patience...
-----Original Message-----
From: [] On Behalf Of Michael Povolotskyi
Sent: Friday, June 10, 2005 6:27 AM
Subject: [Libmesh-users] Nonlinear hanging nodes

Dear Libmesh developers,
I solve nonlinear equation systems using the NonlinearImplicitSystem class .
I'd like to use the adapative mesh refinement for the nonlinear system.
Does libmesh have a tool that can handle the hanging nodes problem for a nonlinear system?
Thank you,
------------------------------------------------------- This SF.Net email is sponsored by: NEC IT Guy Games. How far can you shotput a projector? How fast can you ride your desk chair down the office luge track? If you want to score the big prize, get to know the little guy. Play to win an NEC 61" plasma display: _______________________________________________ Libmesh-users mailing list

Michael Povolotskyi, Ph.D.
University of Rome "Tor Vergata"
Department of Electronic Engineering
Viale Politecnico, 1 - 00133 Rome - Italy 

Phone + 39 06 72597367
Fax   + 39 06 2020519