Hi all,
I was reading through the dof constraints code recently, and two questions
occurred to me that I was hoping to get some input on:
1. In the last part of DofMap::build_constraint_matrix_and_vector() we have:
if ((C.n() == Cnew.m()) && // If the constraint matrix
(Cnew.n() == elem_dofs.size())) // is constrained...
{
C.right_multiply(Cnew);
// If x = Cy + h and y = Dz + g
// Then x = (CD)z + (Cg + h)
C.vector_mult_add(H, 1, Hnew);
}
It seems to me that according to the comment we should do
C.vector_mult_add(H, 1, Hnew); before C.right_multiply(Cnew), no? Isn't the
code currently using "CDg + h" instead of "Cg + h" ? (Most likely I've just
misunderstood the code... but anyway, I thought I'd ask.)
2. What happens if we have a "cycle" of dof constraints? Would the
recursive constraints code handle that properly?
Thanks,
David
------------------------------------------------------------------------------
Don't Limit Your Business. Reach for the Cloud.
GigeNET's Cloud Solutions provide you with the tools and support that
you need to offload your IT needs and focus on growing your business.
Configured For All Businesses. Start Your Cloud Today.
https://www.gigenetcloud.com/
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
