Forgot to copy the list:
On Wed, 7 Jan 2015, Paul T. Bauman wrote:
On Wed, Jan 7, 2015 at 5:08 PM, Roy Stogner <[email protected]>
wrote:
You're sure your system is invertible? I get a zero pivot if I try
"-pc_type lu", but that may just be because PETSc only does partial
pivoting.
Throw -pc_factor_mat_solver_package superlu (or mumps or superlu_dist) in
addition on your command line. If that throws a zero pivot, you can be
confident that there's a zero pivot.
Was using PETSc without any good solver packages built in, but
"-mat_view ::ascii_matlab" and octave were enough to catch the red
herring; the matrix is fine.
---
Roy
------------------------------------------------------------------------------
Dive into the World of Parallel Programming! The Go Parallel Website,
sponsored by Intel and developed in partnership with Slashdot Media, is your
hub for all things parallel software development, from weekly thought
leadership blogs to news, videos, case studies, tutorials and more. Take a
look and join the conversation now. http://goparallel.sourceforge.net
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
