Svetlana Tkachenko <[email protected]> writes: > The example has a line to compute matrix A from Ax = b, and a line to compute > source term (vector b). > >> call KSPSetComputeOperators(ksp,ComputeMatrix,ctx,ierr) >> call KSPSetComputeRHS(ksp,ComputeRHS,ctx,ierr) > > What do I use to give it the initial guess on x?
Those lines are registering callbacks, in order to have a uniform interface that can rediscretize the original operator and on coarse levels of multigrid. Look at the implementation of ComputeMatrix and ComputeRHS for more.
pgpbzRlXhFfpf.pgp
Description: PGP signature
