[petsc-dev] strangness in Chebyshev estimate of eigenvalues

Jed Brown jed at jedbrown.org
Tue May 17 13:18:37 CDT 2016


Barry Smith <bsmith at mcs.anl.gov> writes:
>   Of course, as always, with different number of processes one will
>   get different results, even with Jacobi + Chebyshev because the
>   random numbers generated will be different with different number of
>   processes, is this important?

I can't think of a way to get a "noisy" vector deterministically and
independent of the distribution.  We could use a hash of the global row
index, which would help with debugging in the sense that we could read
in a stored matrix on a different number of cores and get the same
Chebyshev/Jacobi estimates.  I have a minor preference for this instead
of random initialization, but it wouldn't help with most discretizations
because the distribution usually changes the ordering.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20160517/bd2f53bb/attachment.sig>


More information about the petsc-dev mailing list