Hello, It seems that by default FEMSystem always assembles the matrix and ignores the flag System::assemble_before_solve. The DiffSolvers, NewtonSolver and PetscDiffSolver always call System::assembly(false,true). The workarounds I see would be to create my own TimeSolver to override TimeSolver::solve () and call some solver that can call System::assembly(false, System::assemble_before_false). Either that or rewrite the DIffSolvers to call System::assembly() like that.
I think the question should be: Is FEMSystem designed to only handle nonlinear equations? If not, are there other alternatives to recycle the assembly process within FEMSystem? Thanks Miguel ------------------------------------------------------------------------------ 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 [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
