Just a little warning, PETSc MUMPS Cholesky support was broken until the latest 3.5.3 release. See https://bitbucket.org/petsc/petsc/issue/81
Jan On Sun, 01 Mar 2015 15:21:58 +0100 Miro Kuchta <[email protected]> wrote: > Hi, > > please consider the attached script. Following this > <https://bitbucket.org/fenics-project/dolfin/pull-request/2/use-cholesky-rather-than-lu-decomposition/diff#chg-dolfin/fem/LinearVariationalSolver.cpp> > > discussion, if method is mumps, petsc or pastix and > we have symmetric=True the linear system is solved with Cholesky > factorization (it this so?). While testing different method/symmetry > combinations I noticed that PETSc's own symmetric solver is easily > 10 times slower then mumps (I don't have pastix to compare against). > Can anyone else reproduce > this? Thanks. > > Regards, Miro _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
