I'm honestly not sure why I did that, however what I thought would be the obvious fix ( 'Mat h = *A' ) isn't, I'm still getting the error.
Is there a reason I can't do that? -Andrew On May 16, 2012, at 7:17 AM, Jed Brown wrote: > On Wed, May 16, 2012 at 7:08 AM, Andrew Spott <andrew.spott at gmail.com> > wrote: > Not until you pointed it out? :). Though it probably should have been. It > was late. > > Though I'm not sure how to fix that? should I do a MatDestroy(A) at the > beginning of that function every time?. > > Is *A a pointer to the last used matrix, should I just do the mat duplicate > into A? > > Do you really mean to only change A (the operator), but not B (the > preconditioning matrix)? > > Why do you want to create a new matrix instead of reusing those? -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120516/b52b1c4c/attachment.htm>
