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 elementbyelement 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...
Ben
Original Message
From: Michael Povolotskyi [mailto:povolotskyi@...]
Sent: Friday, June 10, 2005 9:02 AM
To: KIRK, BENJAMIN (JSCEG) (NASA)
Subject: Re: [Libmeshusers] 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.
Michael.
KIRK, BENJAMIN (JSCEG) (NASA) wrote:
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 nonconstrained 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...
Ben
Original Message
From: libmeshusersadmin@...
<mailto:libmeshusersadmin@...>
[mailto:libmeshusersadmin@...
<mailto:libmeshusersadmin@...> ] On Behalf Of Michael
Povolotskyi
Sent: Friday, June 10, 2005 6:27 AM
To: libmeshusers@...
<mailto:libmeshusers@...>
Subject: [Libmeshusers] 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,
Michael.
 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: http://www.necitguy.com/?r=20
<http://www.necitguy.com/?r=20>;
_______________________________________________ Libmeshusers mailing list
Libmeshusers@...
<mailto:Libmeshusers@...>
https://lists.sourceforge.net/lists/listinfo/libmeshusers
<https://lists.sourceforge.net/lists/listinfo/libmeshusers>;


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
http://www.optolab.uniroma2.it/pages/moshe/moshe.html
<http://www.optolab.uniroma2.it/pages/moshe/moshe.html>;

