On Fri, 11 Jul 2008, Maxime Debon wrote: > Back to Lagrange, the new file uses the penalty method.
This one gives me a perfect straight line answer, with the SVN libMesh and my default options, in parallel or serial... with the exception that the lower boundary is at a value of 1e-10 instead of 0, of course. Ah, but it looks like if I get rid of my default PETSc options, I get that nasty curve instead! It's still an incomplete solve; you're just hitting the default max iterations count of 5000 rather than a user-imposed count. I'm going to copy this to libmesh-users, even though it's just the same problem again; we ought to be able to come up with some behavior better than "quit silently with the wrong solution" when a max iteration count is reached. It feels to me like the "right" thing to do would be to throw an exception, so that an outer nonlinear or time solver would have a chance to recover but a user of the linear solver alone could get some feedback as to what went wrong. However, we've still got exception throwing as an optional libMesh feature (one that's not even in a non-SVN release yet, I think), and this would be the first time we used it for anything other than an unexpected error. --- Roy ------------------------------------------------------------------------- Sponsored by: SourceForge.net Community Choice Awards: VOTE NOW! Studies have shown that voting for your favorite open source project, along with a healthy diet, reduces your potential for chronic lameness and boredom. Vote Now at http://www.sourceforge.net/community/cca08 _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
