Hi, I'm using PETScWrapper to parallelize my code. In my preconditioner for the GMRES solver, there is one step that requires a matrix copied from the system matrix, and set all the elements to be the absolute value. It was fine in serial because I can iterator over all the entries simply using the iterator, but it does not work in parallel since each process only owns local entries. When I run the code using iteration the error message will tell me "argument out of range; Not local rows". It also tells me "Object is in wrong state; Not for unassembled matrix" despite the fact that I called compress(VectorOperation:add) right before the iteration.
Is there any way to bypass the iteration to get an absolute value matrix or to iterate locally and not trigger the "unassembled matrix" error? Thanks, Feimi -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.