<div class="gmail_quote">On Wed, May 23, 2012 at 2:26 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
  This looks good:<br>
<div class="im"><br>
$ $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<br>
lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
  Linear solve converged due to CONVERGED_RTOL iterations 10<br>
  Linear solve converged due to CONVERGED_RTOL iterations 7<br>
  Linear solve converged due to CONVERGED_RTOL iterations 8<br>
  Linear solve converged due to CONVERGED_RTOL iterations 8<br>
  Linear solve converged due to CONVERGED_RTOL iterations 9<br>
  Linear solve converged due to CONVERGED_RTOL iterations 9<br>
  Linear solve converged due to CONVERGED_RTOL iterations 8<br>
  Linear solve converged due to CONVERGED_RTOL iterations 8<br>
  Linear solve converged due to CONVERGED_RTOL iterations 7<br>
Number of SNES iterations = 9<br>
<br>
</div>   I assume it is independent of the number of levels of refinement?<br></blockquote><div><br></div><div><div>$ $base -da_refine 0 -snes_grid_sequence 5 -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</div>
<div>lid velocity = 100, prandtl # = 1, grashof # = 10000</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 13</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 7</div>
<div>            Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 9</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 7</div>
<div>            Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 10</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 7</div>
<div>            Linear solve converged due to CONVERGED_RTOL iterations 12</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>            Linear solve converged due to CONVERGED_RTOL iterations 11</div>
<div>          Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>          Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>          Linear solve converged due to CONVERGED_RTOL iterations 11</div>
<div>          Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>          Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>        Linear solve converged due to CONVERGED_RTOL iterations 14</div>
<div>        Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>        Linear solve converged due to CONVERGED_RTOL iterations 11</div><div>        Linear solve converged due to CONVERGED_RTOL iterations 10</div>
<div>      Linear solve converged due to CONVERGED_RTOL iterations 9</div><div>      Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>      Linear solve converged due to CONVERGED_RTOL iterations 8</div>
<div>      Linear solve converged due to CONVERGED_RTOL iterations 7</div><div>      Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>    Linear solve converged due to CONVERGED_RTOL iterations 8</div><div>
    Linear solve converged due to CONVERGED_RTOL iterations 6</div><div>    Linear solve converged due to CONVERGED_RTOL iterations 7</div><div>    Linear solve converged due to CONVERGED_RTOL iterations 6</div><div>  Linear solve converged due to CONVERGED_RTOL iterations 11</div>
<div>  Linear solve converged due to CONVERGED_RTOL iterations 7</div><div>  Linear solve converged due to CONVERGED_RTOL iterations 5</div><div>  Linear solve converged due to CONVERGED_RTOL iterations 5</div><div>Number of SNES iterations = 4</div>
</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
  Note that you could use -pc_type eisenstat perhaps in this case instead. Might save lots of flops?  I've often wondered about doing Mark's favorite chebyshev smoother with Eisenstat, seems like it should be a good match.<br>
</blockquote><div><br></div><div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------------</div><div>[0]PETSC ERROR: No support for this operation for this object type!</div><div>[0]PETSC ERROR: Cannot have different mat and pmat!</div>
</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span class="HOEnZb"><font color="#888888"><br>
<br>
<br>
   Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
On May 23, 2012, at 11:48 AM, Jed Brown wrote:<br>
<br>
> If ex19 is what matters...<br>
><br>
> $ base='mpiexec.hydra -n 4 ./ex19 -da_refine 3 -lidvelocity 100 -grashof 1e4 -pc_type mg -ksp_type fgmres -pc_mg_type full -ksp_converged_reason'<br>
> $ $base<br>
> lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 23<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 12<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 12<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 17<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 15<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 13<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 19<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 21<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 18<br>
> Number of SNES iterations = 9<br>
> $ $base -mg_levels_ksp_type chebyshev -mg_levels_ksp_chebyshev_estimate_eigenvalues 0.1,1.1 -mg_levels_ksp_max_it 1 -mg_levels_pc_type pbjacobi<br>
> lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 25<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 27<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 27<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 30<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 28<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 42<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 49<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 52<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 48<br>
> Number of SNES iterations = 9<br>
> $ $base -mg_levels_ksp_type chebyshev -mg_levels_ksp_chebyshev_estimate_eigenvalues 0.1,1.1 -mg_levels_ksp_max_it 1 -mg_levels_pc_type sor<br>
> lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 22<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 16<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 17<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 19<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 19<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 17<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 20<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 18<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 17<br>
> Number of SNES iterations = 9<br>
> $ $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 pbjacobi<br>
> lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 15<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 13<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 16<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 15<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 16<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 16<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 18<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 17<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 16<br>
> Number of SNES iterations = 9<br>
> $ $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<br>
> lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 10<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 7<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 8<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 8<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 9<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 9<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 8<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 8<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 7<br>
> Number of SNES iterations = 9<br>
> $ $base -mg_levels_ksp_type chebyshev -mg_levels_ksp_chebyshev_estimate_eigenvalues 0.1,1.1 -mg_levels_ksp_max_it 3 -mg_levels_pc_type pbjacobi<br>
> lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 12<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 9<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 11<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 11<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 11<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 12<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 12<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 11<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 11<br>
> Number of SNES iterations = 9<br>
> $ $base -mg_levels_ksp_type chebyshev -mg_levels_ksp_chebyshev_estimate_eigenvalues 0.1,1.1 -mg_levels_ksp_max_it 3 -mg_levels_pc_type sor<br>
> lid velocity = 100, prandtl # = 1, grashof # = 10000<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 9<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 6<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 6<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 7<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 7<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 6<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 7<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 6<br>
>   Linear solve converged due to CONVERGED_RTOL iterations 6<br>
> Number of SNES iterations = 9<br>
><br>
<br>
</div></div></blockquote></div><br>