Dear Developers, I'm in the process of upgrading my code to libmesh 1.0.0
In doing so I found a new behavior of of the constrain_element_matrix_and_vector. My problem is as follows: I solve a system and I have to define some periodic boundary conditions and some Dirichlet boundary conditions. My application is an old one, it was developed before DirichletBoundaries and PeriodicBeoundaries existed in libmesh. So at that time (around 2005) I developed the following solution. I added constraints to the dof_map according to my boundary conditions. This was not a very easy piece of code, since I had to add constraints in terms of unconstrained DOFs only, but it was working. When I upgraded to libmesh 1.0.0, I found that in order to get correct results, I had to call constrain_element_matrix and constrain_element_vector with asymmetric_constraint_rows equal to false. Does it make sense to you? Michael. ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
