On Thu, 26 Mar 2009, Jed Brown wrote: > 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.
If only I was confident that none of our bugs were too subtle to show up on the first timestep. Unfortunately the last ghosted vector bug we fixed only showed up after half a dozen more more refinement steps, and IIRC our outstanding ParallelMesh bug takes at least three refinement and coarsening steps to trigger. > 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). Right, but for regression testing we don't care as much about scalability as repeatability. This looks like it will be a big help, thanks! --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
