or Matt should post it or fix
mine :)
A
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120523/a73078fe/attachment.html
lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120523/a40245c6/attachment.html
On May 22, 2012, at 4:17 PM, Hong Zhang wrote:
Barry,
Do you recall the paper or detailed Chebyshev algorithm that you implemented
the Chebyshev in petsc?
No
I added 'hybrid' to it, which is not correct. Jie is helping to fix it.
Having an original
algorithm would help.
The code
On May 22, 2012, at 9:32 PM, Jed Brown wrote:
Well, here's the problem:
MatSetBlockSize() used to call mat-ops-setblocksize() which the MPI
matrices forwarded down to their seq submatrices.
But then this patch removed that behavior so MatSetBlockSize() is processed
immediately and not
Cheby can be tricky, textbook algorithms (eg, Saad's book) have some curve
balls. Prometheus has a Cheby implementation that I know is correct. It was
worked on by several good people and implemented in several codes and we wrote
a paper on it.
Mark
On May 23, 2012, at 9:21 AM, Barry Smith
On May 22, 2012, at 11:09 PM, Barry Smith wrote:
On May 22, 2012, at 9:58 PM, Jed Brown wrote:
On Tue, May 22, 2012 at 9:52 PM, Barry Smith bsmith at mcs.anl.gov wrote:
How about a default nonlinear Krylov accelerator? fgrmes is not much more
expensive than gmres under causal
://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120523/9b079114/attachment.html
On May 23, 2012, at 9:36 AM, Barry Smith wrote:
On May 22, 2012, at 9:32 PM, Jed Brown wrote:
Well, here's the problem:
MatSetBlockSize() used to call mat-ops-setblocksize() which the MPI
matrices forwarded down to their seq submatrices.
But then this patch removed that behavior so
comments in
the form of SetBlockSize statements that are commented out in places were I
had trouble.
The problem is fixed.
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120523/01d524a5/attachment.html
iterations 6
Number of SNES iterations = 9
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120523/183d4aa8/attachment.html
On May 23, 2012, at 9:26 AM, Jed Brown wrote:
On Wed, May 23, 2012 at 9:27 AM, Mark F. Adams mark.adams at columbia.edu
wrote:
I have gone through and done this manual inheritance, but 1) I may have
missed some and 2) there were some that I could not figure out how to do.
I'm back
On May 23, 2012, at 9:23 AM, Mark F. Adams wrote:
On May 22, 2012, at 11:09 PM, Barry Smith wrote:
On May 22, 2012, at 9:58 PM, Jed Brown wrote:
On Tue, May 22, 2012 at 9:52 PM, Barry Smith bsmith at mcs.anl.gov wrote:
How about a default nonlinear Krylov accelerator? fgrmes is not
This looks good:
$ $base -mg_levels_ksp_type chebyshev
-mg_levels_ksp_chebyshev_estimate_eigenvalues 0.1,1.1 -mg_levels_ksp_max_it 2
-mg_levels_pc_type sor
lid velocity = 100, prandtl # = 1, grashof # = 1
Linear solve converged due to CONVERGED_RTOL iterations 10
Linear solve
13 matches
Mail list logo