Hi there,
I’m using libMesh to perform topology optimization on a mesh with a spherical surface. I only have a single LinearImplicitSystem. The normal displacement of nodes on the spherical surface needs to be constrained to force the surface to remain spherical. I am using Lagrange multipliers to do this, but have found that adding a large number of scalar variables results in EquationSystems::init() taking an extraordinary amount of time to build the system. For a mesh with 9585 nodes, 45728 elements, and first-order interpolation, EquationSystems::init() takes practically no time at all without scalar variables for the MPCs. However, when I add 1332 scalar variables (one for each node on the spherical surface, EquationSystems::init() runs for over an hour. Is this to be expected with libMesh? Is there a more efficient way to apply MPCs to numerous nodes? Bailey C ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users