With this other options it actually works but I suspect the solver is not really changing since the residual of all iterations of Newton's method are the same between all solvers (direct or not).
My issue is that the solution is worse whenever I increase the degrees of freedom, maybe there is some singularities that I am not aware of... Thanks On Sun, Oct 19, 2014 at 3:24 PM, Jed Brown <[email protected]> wrote: > Rodrigo Broggi <[email protected]> writes: > > > Result: > > > > [0]PETSC ERROR: --------------------- Error Message > > ------------------------------------ > > [0]PETSC ERROR: No support for this operation for this object type! > > [0]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc LU! > > [0]PETSC ERROR: > > ------------------------------------------------------------------------ > > [0]PETSC ERROR: Petsc Release Version 3.4.2, Jul, 02, 2013 > > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [0]PETSC ERROR: See docs/index.html for manual pages. > > > > Does that mean that I don't have the library? > > PETSc does not have a "native" distributed LU factorization. I > recommend debugging in serial, but if you want a parallel factorization, > use -pc_type lu -pc_factor_mat_solver_package superlu_dist (or > -pc_factor_mat_solver_package mumps). This requires configuring PETSc > with the external package SuperLU_DIST, e.g., as > --download-superlu_dist (or similar for MUMPS). > ------------------------------------------------------------------------------ Comprehensive Server Monitoring with Site24x7. Monitor 10 servers for $9/Month. Get alerted through email, SMS, voice calls or mobile push notifications. Take corrective actions from your mobile device. http://p.sf.net/sfu/Zoho _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
