I ended up creating a function similar to the one in LInearImplicitSystem to get around this issue.
Working out fine now. Thanks! Manav > On Jan 29, 2018, at 1:59 PM, John Peterson <jwpeter...@gmail.com> wrote: > > > > On Mon, Jan 29, 2018 at 11:42 AM, Manav Bhatia <bhatiama...@gmail.com > <mailto:bhatiama...@gmail.com>> wrote: > Except for a sensitivity solve. > > Hmm... so you are referring to the get_linear_solver() calls in > > ImplicitSystem::adjoint_solve() > ImplicitSystem::sensitivity_solve() > ImplicitSystem::weighted_sensitivity_adjoint_solve() > ImplicitSystem::weighted_sensitivity_solve() > > ? > > If I understand correctly, those should all be calling the derived class > LinearImplicitSystem::get_linear_solver() method in practice. > > In the class you are writing that is derived from NonlinearImplicitSystem, if > you hang on to your own LinearSolver object, you could write a function > similar to the one above to provide access to it. > > -- > John ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users