The problem was the newly added check that the preconditioner must be set to "none" when a direct method is used. This wasn't checked before.
The default values have now been changed (in LinearSolver, LinearVariationalSoler and NonlinearVariationalSolver) to "default" "none" This will result in LU being chosen in LinearSolver.cpp. All the logic for this is in the constructor of LinearSolver.cpp if more tuning is needed. -- Anders On Thu, Oct 20, 2011 at 10:22:45AM +0100, Garth N. Wells wrote: > What was the problem that has been 'fixed'? > > Garth > > ---------- Forwarded message ---------- > From: <[email protected]> > Date: 20 October 2011 10:11 > Subject: [Branch ~dolfin-core/dolfin/main] Rev 6347: Fix problem with > default values for linear solver and preconditioner, > To: Garth Wells <[email protected]> > > > ------------------------------------------------------------ > revno: 6347 > committer: Anders Logg <[email protected]> > branch nick: work > timestamp: Thu 2011-10-20 11:09:10 +0200 > message: > Fix problem with default values for linear solver and preconditioner, > now set to "default" and "none". > modified: > dolfin/fem/LinearVariationalSolver.h > dolfin/fem/NonlinearVariationalSolver.h > dolfin/la/LinearSolver.h > > > === modified file 'dolfin/fem/LinearVariationalSolver.h' > --- dolfin/fem/LinearVariationalSolver.h 2011-06-22 10:01:17 +0000 > +++ dolfin/fem/LinearVariationalSolver.h 2011-10-20 09:09:10 +0000 > @@ -18,7 +18,7 @@ > // Modified by Marie E. Rognes, 2011. > // > // First added: 2011-01-14 (2008-12-26 as VariationalProblem) > -// Last changed: 2011-06-22 > +// Last changed: 2011-10-20 > > #ifndef __LINEAR_VARIATIONAL_SOLVER_H > #define __LINEAR_VARIATIONAL_SOLVER_H > @@ -53,8 +53,8 @@ > { > Parameters p("linear_variational_solver"); > > - p.add("linear_solver", "lu"); > - p.add("preconditioner", "default"); > + p.add("linear_solver", "default"); > + p.add("preconditioner", "none"); > p.add("symmetric", false); > p.add("reset_jacobian", true); > > > === modified file 'dolfin/fem/NonlinearVariationalSolver.h' > --- dolfin/fem/NonlinearVariationalSolver.h 2011-08-24 22:15:11 +0000 > +++ dolfin/fem/NonlinearVariationalSolver.h 2011-10-20 09:09:10 +0000 > @@ -18,7 +18,7 @@ > // Modified by Marie E. Rognes, 2011. > // > // First added: 2011-01-14 (2008-12-26 as VariationalProblem.h) > -// Last changed: 2011-06-22 > +// Last changed: 2011-10-20 > > #ifndef __NONLINEAR_VARIATIONAL_SOLVER_H > #define __NONLINEAR_VARIATIONAL_SOLVER_H > @@ -57,8 +57,8 @@ > { > Parameters p("nonlinear_variational_solver"); > > - p.add("linear_solver", "lu"); > - p.add("preconditioner", "default"); > + p.add("linear_solver", "default"); > + p.add("preconditioner", "none"); > p.add("symmetric", false); > p.add("reset_jacobian", true); > > > === modified file 'dolfin/la/LinearSolver.h' > --- dolfin/la/LinearSolver.h 2011-10-07 12:24:30 +0000 > +++ dolfin/la/LinearSolver.h 2011-10-20 09:09:10 +0000 > @@ -44,8 +44,8 @@ > public: > > /// Create linear solver > - LinearSolver(std::string method = "lu", > - std::string preconditioner = "ilu"); > + LinearSolver(std::string method = "default", > + std::string preconditioner = "none"); > > /// Destructor > ~LinearSolver(); > > _______________________________________________ > Mailing list: https://launchpad.net/~dolfin > Post to : [email protected] > Unsubscribe : https://launchpad.net/~dolfin > More help : https://help.launchpad.net/ListHelp _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
