Hi Satyendra, >From what I understand, the solver uses the current solution vector (all zero or otherwise) as the starting point. So, when solving Ax=b and you want the solver starting point to be x0, then just set x=x0 before the solver step "solver.solve (A, x, b, preconditioner)".
Hope this helps. J-P On 25 June 2010 17:16, Satyendra Tomar <[email protected]> wrote: > Hi everyone, > > I presume that the CG solver is taking the initial guess to be a zero > vector. Could someone please let me know how to change it to a random > vector? > > Bests, > Satyendra > PS: I haven't received any response to my previous post (Subject: > Timings of the mixed method for Poisson problem). I'm not sure if it > was successfully delivered (though I received a copy of that). > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii >
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
