On Sat, 1 Feb 2014, Jens Lohne Eftang wrote: > I used these arguments for 3D elasticity and they seemed to work well > also in combination with DirichletBoundary (why does that break symmetry?)
When you constrain your element matrix and vector (any constraints, but including the heterogenous constraints from DirichletBoundary), the default argument "asymmetric_constraint_rows=true" causes your matrix to be left in an asymmetric state. The advantage is that in this case the linear solver can handle your constraint equations for you; the disadvantage is an asymmetric state. --- Roy ------------------------------------------------------------------------------ Managing the Performance of Cloud-Based Applications Take advantage of what the Cloud has to offer - Avoid Common Pitfalls. Read the Whitepaper. http://pubads.g.doubleclick.net/gampad/clk?id=121051231&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
