From: Jed B. <je...@59...> - 2009-03-26 17:54:34
|
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 |