Hello, when starting to use FEMSystem I discovered some strange behaviour which looks like a bug: The FEMContext has a attribute request_jacobian which I assumed to have the same meaning as the parameter to e.g. FEMSystem::side_constraint(..). But it the value in FEMContext never is updated to reflect the current state of request_jacobian. So either remove request_jacobian from FEMContext or set it to the same value as the parameter.
Another thing I like to know is whether there is some mechanism to store additional information for nodes or quadrature points so it gets distributed and reassembled when using MPI. I intend to use this e.g. to for store the progress of damage or other internal state of the material. Robert ------------------------------------------------------------------------------ The ultimate all-in-one performance toolkit: Intel(R) Parallel Studio XE: Pinpoint memory and threading errors before they happen. Find and fix more than 250 security defects in the development cycle. Locate bottlenecks in serial and parallel code that limit performance. http://p.sf.net/sfu/intel-dev2devfeb _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
