On Thu, Sep 13, 2012 at 3:06 PM, David Knezevic <[email protected]> wrote: > On a related issue, I wanted to mention 3D linear elasticity example > that I added: systems_of_equations_ex6. It seems very slow with the > default iterative solver in PETSc (ILU + GMRES, right?). The example > takes 26 seconds on my laptop. > > But this is a symmetric problem, so you'd expect CG to be good. Calling > ./systems_of_equations_ex6-opt -ksp_type cg > runs fast, but gives the wrong answer (I should add some code to check > the final residual...)
I'd check to be sure your matrix is really symmetric, too. Perhaps something like boundary conditions made it non-symmetric? -- John ------------------------------------------------------------------------------ Live Security Virtual Conference Exclusive live event will cover all the ways today's security and threat landscape has changed and how IT managers can respond. Discussions will include endpoint security, mobile security and the latest in malware threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
