Hi, I have a code that I am writing based on Miscellaneous ex 7. I set the tolerances in the constructor (that derives from TransientNonlinearImplicitSystem) through
nonlinear_solver->absolute_residual_tolerance = _solbub._tol; nonlinear_solver->relative_residual_tolerance = _solbub._rel_tol; However when the residual and jacobian functions are called - I don't get the same tolerances that were set in the constructor. I am running the code with the following, ./main-opt --use-petsc-dm -snes_type virs --snes_linesearch_type basic so that it uses petscdmlibmesh.C and petsc_dm_nonlinear_solver.C . I looked inside petsc_dm_nonlinear_solver.C and all seems to be right, and the values are taken from the parent classes. Do I need to do something else for this? As an aside - how do I pass the petsc options from a constructor instead of passing the options at the command line. Subramanya Sadasiva "But memory eventually fades. Turbulences damp out, internal strains yield to plastic flow, concentration inhomogeneities diffuse to uniformity. Systems tend to subside to very simple states,independent of their specific history." Herbert Callen . Thermodynamics and an Introduction to Thermostatics. ------------------------------------------------------------------------------ The Go Parallel Website, sponsored by Intel - in partnership with Geeknet, is your hub for all things parallel software development, from weekly thought leadership blogs to news, videos, case studies, tutorials, tech docs, whitepapers, evaluation guides, and opinion stories. Check out the most recent posts - join the conversation now. http://goparallel.sourceforge.net/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
