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? 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
