Here is a really, really raw example, the comments are not clear right now, but I wanted to keep you informed. This requires the latest svn branch to work.
Unpack it in the ./examples directory and run make. Run it as $ ./ex19-dbg -snes_view -r 4 for a successive approximation which will converge linearly, and $ ./ex19-dbg -snes_view -r 4 -snes_mf_operator for a matrix-free approach in the iterative solver which will converge quadratically. > Anyway, I do have a question regarding using PETSc object with LibMesh. I > have been trying to use Petsc objects Mats, Vecs and SNES solver with > Libmesh but the one thing I cant seem figure out is how to set the > nonlinear_solver public attribute of say a NonlinearImplicitSystem object to > a PETSc SNES object which I have created and initialized separately. Since, > the SNES object used in the wrapper PetscNonlinearSolver is private, I don¹t > understand how this can be done. > > Have I missed something and taken a completely wrong path on this ? I would > very much appreciate any comments that you can provide to help me out here. The user interface is totally up for discussion since I am the only one who as exercised it to date. (I am sure Roy will have some comments!) It seems to me the right approach will be to add a method which gives the user access to the SNES object? From there the KSP, Mat, Vec, PC, etc... can be accessed. This would be similar to the approach used in the PetscLinearSolver. -Ben
ex19.tar.gz
Description: Binary data
------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
