On 05/10/2012 02:54 PM, Vijay S. Mahadevan wrote:
With PETSc, you will have to specify that you are using a non-zero
initial guess using
http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/KSP/KSPSetInitialGuessNonzero.html
KSPSetInitialGuessNonzero(ksp, PETSC_TRUE)
If this is not set, it will reset the vector to zero and sure, then
there is nothing to be gained. I haven't looked if a similar behavior
exists for the Deal.II internal solvers.
The deal.II solvers simply start with whatever values are in the
solution vector whenever you call
cg.solve (system_matrix, solution, system_rhs, preconditioner);
Best
W.
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
