[petsc-dev] I do not think this is the right solution

[Jed Brown]

On Sun, Mar 18, 2012 at 18:45, Matthew Knepley <knepley at gmail.com> wrote:

> No, I don't want the nonzeros removed, I want the values set to zero,
> exactly like I said. This is for robustly checking
> that a matrix created in parallel is the same as one made on a different
> number of processors. Something I do for ex62,
> but we obviously do not give a shit about doing.
>

What? Why not do a MatGetRedundantMatrix() on the parallel thing and
compute a norm of the difference to the serial thing?

And if you really want to stay parallel and print all the non-small values,
MatDuplicate and fill up the new matrix with the thresholded values.


Note, there was a request a while back to be able to make an ASCII viewer
use a user-defined floating-point format string.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120318/fcff745a/attachment.html>