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

Matthew Knepley knepley at gmail.com
Sun Mar 18 18:54:59 CDT 2012


On Sun, Mar 18, 2012 at 6:52 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> 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?
>

That has nothing to do with anything we are talking about.


> 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.
>

Neither does this.


> Note, there was a request a while back to be able to make an ASCII viewer
> use a user-defined floating-point format string.
>

That would not help. For a small number, you do not want ANY significant
digits.

    Matt

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120318/6a86c49c/attachment.html>


More information about the petsc-dev mailing list