Hi, I am doing solid mechanics calculations with FEMSystem and a moving mesh and need periodic boundary conditions. As the PR-BC code in libmesh only works for unmodified meshes, I wrote some thing similar which saves a mapping of the nodes on the positive side to the nodes of the negative side in the first timestep (on the unmodified mesh) and uses this mapping in all following timesteps to impose the periodic boundary conditions.
As especially the PR-BC causes entries in "unusual" regions of the system matrix, I have currently a rather large performance issue due to the resorting of the matrix while adding the PR-BCs. So how can I influence the default sparsity pattern used to create the system matrix to avoid these performance issues? Robert ------------------------------------------------------------------------------ All of the data generated in your IT infrastructure is seriously valuable. Why? It contains a definitive record of application performance, security threats, fraudulent activity, and more. Splunk takes this data and makes sense of it. IT sense. And common sense. http://p.sf.net/sfu/splunk-d2dcopy2 _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
