On Thu 2009-03-26 12:40, Roy Stogner wrote: > I'd be tempted to add that behavior as a very-non-default option to > FEMSystem... except, for it to really be useful you'd need a > deterministic solve() too, and I'm sure the same distributed FP > ordering problems come up there.
Nah, just MatView and confirm that they're identical. Also, you should be able to guarantee a deterministic solve with -ksp_type preonly -pc_type redundant which duplicates the entire matrix on every process and solves redundantly on each process using LU (obviously it's not scalable). Jed
pgpJx0EzqI3kv.pgp
Description: PGP signature
------------------------------------------------------------------------------
_______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
