On 26 Mar 2014, at 16:26, Jan Blechta <[email protected]> wrote:
> On Wed, 26 Mar 2014 16:16:25 +0100 > Johannes Ring <[email protected]> wrote: > >> On Wed, Mar 26, 2014 at 1:39 PM, Jan Blechta >> <[email protected]> wrote: >>> As a follow-up of 'Broken PETSc wrappers?' thread on this list, can >>> anyone reproduce incorrect (orders of magnitude) norm using >>> superlu_dist on following example? Both in serial and parallel. >>> Thanks, >> >> This is the result I got: >> >> Serial: >> >> L2 norm mumps 0.611356580181 >> L2 norm superlu_dist 92.4733890983 >> >> Parallel (2 processes): >> >> L2 norm mumps 0.611356580181 >> L2 norm superlu_dist 220.027905995 >> L2 norm mumps 0.611356580181 >> L2 norm superlu_dist 220.027905995 > > superlu_dist results are obviously wrong. Do we have broken > installations or is there something wrong with the library? > > In the latter case I would suggest switching the default back to MUMPS. > (Additionally, MUMPS has Cholesky factorization!) What was your > motivation for switching to superlu_dist, Garth? > MUMPS often fails in parallel with global dofs, and there is no indication that MUMPS developers are willing to fix bugs. Garth > Jan > >> >> Johannes >> _______________________________________________ >> fenics-support mailing list >> [email protected] >> http://fenicsproject.org/mailman/listinfo/fenics-support _______________________________________________ fenics-support mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics-support
