Thanks for the patch... I'll take a look at all of this shortly and commit it.
As for the Jacobian in parallel... let's call not_implemented(). I would rather get an error than get completely different behavior and not know why. Thanks again! Derek On Oct 17, 2008, at 6:39 AM, Norbert Stoop wrote: > Derek Gaston wrote: > >>> Ok, I'll do that as soon as I have made some more tests. >> >> As long as you're making tests... can you turn that LaplaceYoung >> example >> into a _real_example (just clean it up a bit). We really need to get >> that thing committed so that people have something to look at for >> how to >> use the NonlinearSolver classes. > > Attached is a cleaned up version of it. The biggest change is that I > removed the compute_matvec function which is not used anyway. Is that > okay?. There is a lot of redundant code for the element and quadrature > setup etc., but I think for the moment this is tolerable... :-) > Feel free to make any changes to it. > > nox.patch is for fixing the matrix-free parallel NOX solver - not > really > worth a patch, just the two lines for closing the vectors in computeF. > > A patch for Jacobian support will follow... How should we deal with > its > current limitation to single process execution? Should we default to > use > the Jacobian when run on a single CPU and fall back to MF when run in > parallel? Or just throw a not_implemented for parallel execution? > > Norbert > <ex19.tar.gz>Index: trilinos_nox_nonlinear_solver.C > =================================================================== > --- trilinos_nox_nonlinear_solver.C (revision 3104) > +++ trilinos_nox_nonlinear_solver.C (working copy) > @@ -96,6 +96,8 @@ > if( _solver->residual == NULL ) > return false; > _solver->residual (*sys.current_local_solution.get(), R); > + R.close(); > + X_global.close(); > return true; > } > > ------------------------------------------------------------------------- > This SF.Net email is sponsored by the Moblin Your Move Developer's > challenge > Build the coolest Linux based applications with Moblin SDK & win > great prizes > Grand prize is a trip for two to an Open Source event anywhere in > the world > http://moblin-contest.org/redirect.php?banner_id=100&url=/_______________________________________________ > Libmesh-devel mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-devel ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
