Dear all, 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/e618c3eb/attachment.html>
