Now it's working. I added the vectors I mentioned before to the system using the routine System::add_vector() and the cg works after the refinement. It's weird because it was working before after refinement when using gmres.
Could you help me out with the preconditioners please? Just to mention how you set up the preconditioners you usually employ. Like you said: AMG on each processor wins on problems where it works, ILU (usually with ~4 levels) wins otherwise, tweaking Block Jacobi vs Additive Schwarz to tie it together can speed things up slightly. Do you set them up from the command line? Or they require more settings inside the code? Thanks Miguel On Fri, Sep 5, 2014 at 9:51 AM, Roy Stogner <[email protected]> wrote: > > On Thu, 4 Sep 2014, Miguel Angel Salazar de Troya wrote: > > I'm using a FEMSystem with a NewtonSolver. I realized that in the >> FEMSystem, in the assembly routines, DofMap::constrain_element_ >> matrix_and_vector >> is called with asymmetric_constraint_rows = false. Also, because it's a >> nonlinear solver. DofMap::enforce_constraints_exactly(system) is called >> as >> well. Would this ensure me a symmetric matrix even for the refinement >> case? >> > > Assuming your element matrices are symmetric before refinement, this > should keep them symmetric after refinement, yes. > --- > Roy > -- *Miguel Angel Salazar de Troya* Graduate Research Assistant Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign (217) 550-2360 [email protected] ------------------------------------------------------------------------------ Slashdot TV. Video for Nerds. Stuff that matters. http://tv.slashdot.org/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
