Hi Adriano, On Fri, Mar 11, 2016 at 1:49 PM, Adriano Côrtes <adrimacor...@gmail.com> wrote:
> Dear Roy, > > I discovered that the example 1 of FEMSystem (NS equation) uses > NewtonSolver, that is, _diff_solver in libMesh::TimeSolver points to > NewtonSolver. Where should I change the code to use PetscDiffSolver? Is > there an example where I can take a look? > Unfortunately there's not an example (that I know of). But, what you can do is when the TimeSolver is first created (line 144 in the unmodifeed fem_system_ex1.C), you can the the diff_solver() ( http://libmesh.github.io/doxygen/classlibMesh_1_1TimeSolver.html#abaf790fa7f3fc464e6e607a9b914a753) and reset it to PetscDiffSolver: system.time_sover.diff_solver().reset( new PetscDiffSolver(system) ) Let us know if you have any trouble. I don't think anyone has hit PetscDiffSolver very hard, so please do let us know if anything peculiar pops up. Thanks, Paul ------------------------------------------------------------------------------ Transform Data into Opportunity. Accelerate data analysis in your applications with Intel Data Analytics Acceleration Library. Click to learn more. http://pubads.g.doubleclick.net/gampad/clk?id=278785111&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
