NAN : After reset matrix A, you must call KSPSetOperators(ksp,...) again. See petsc-dev/src/ksp/ksp/examples/tutorials/ex5.c on how to solve modified matrices.
Hong > > I need to solve Ax=b in a iteratively, and the A is b is ajusted in each > iteration. I am using matzeroentries to zeroout the values in A (A is in > MPIBAIJ format). Somehow, I found this cause me some unexpected wrong > solutions. I do not want to retain the non-zero structures in A when I try > to zero it. Does anyone have a idea to do that. Or I need to destroy the A > matrix each time I give new values to it (the non-zero structure is > changing). > > Thanks, > Nan > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121030/344f1b5a/attachment.html>
