<div class="gmail_quote">On Sun, Mar 18, 2012 at 18:45, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>No, I don't want the nonzeros removed, I want the values set to zero, exactly like I said. This is for robustly checking</div><div>that a matrix created in parallel is the same as one made on a different number of processors. Something I do for ex62,</div>

<div>but we obviously do not give a shit about doing.</div></blockquote></div><br><div>What? Why not do a MatGetRedundantMatrix() on the parallel thing and compute a norm of the difference to the serial thing?</div><div>
<br></div><div>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.</div><div><br></div><div><br></div><div>Note, there was a request a while back to be able to make an ASCII viewer use a user-defined floating-point format string.</div>