On 28/01/14 19:49, Garth N. Wells wrote:
A plus of this approach is that if your system is symmetric you can use CG or MINRES, which you can't with NonlinearDiscreteProblem.
I don't really understand your comment. Perhaps I've misunderstood. You can use CG or MINRES perfectly well without subclassing NonlinearProblem directly. Try it out in the nonlinear-poisson demo. That's using NonlinearDiscreteProblem under the hood. Am I missing something?
I don't have an objection to it being spun out, but it is of limited practical use.
I think it would be quite useful, as it enables people whose codes currently use the solve(F == 0) interface to easily set up custom SNES solvers without copying and pasting the same code everywhere. On 28/01/14 21:13, Anders Logg wrote:
So in your case, either pick (1) or (2). I don't see that there should be an interface in between.
I wasn't suggesting to create a new interface in between. I was suggesting that the way (3) is turned into (1) should be accessible to users, so that codes that currently use (3) can easily use (1) themselves. Cheerio, Patrick _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
