Re: [Libmesh-users] Set an Initial Guess of a Nonlinear System

2014-01-15 Thread Dafang Wang
Hi Paul, Thanks for your tips. I extracted the current_solution in the follow way, but I got an all-zero result when running with multiple processors. Only the single processor yielded the right result. Do I miss something? NonlinearSystem system; PetscVector v; v.init( *(system.current_local_s

Re: [Libmesh-users] Set an Initial Guess of a Nonlinear System

2014-01-15 Thread Paul T. Bauman
I've done something similar to this. You'll need to catch the error the solver throws (ConvergenceFailure, I think, double check), then reset the solution to your cached, good solution. This can be done by grabbing the solution vector from the system. Don't forget to close() the solution vector onc

[Libmesh-users] Set an Initial Guess of a Nonlinear System

2014-01-14 Thread Dafang Wang
Hi folks, Does anyone know how to set the initial guess solution of a nonlinear system from a known vector? My code is in the following manner: NonlinearImplicitSystem system. Vector v1 = system.current_solution(); //Step1 system.solve(); if the nonlinear solve diverges Reset system.current_