[petsc-dev] Getting Richardson output

Smith, Barry F. bsmith at mcs.anl.gov
Mon Mar 11 20:59:52 CDT 2019

   Richardson has this "special case where it calls PCApplyRichardson() under certain circumstances (for efficiency) but KSPSolve_Richardson() has the check

  if (exists && maxit > 0 && richardsonP->scale == 1.0 && !ksp->numbermonitors && !ksp->transpose_solve && !nullsp) {

so if monitoring is turned on it should run "the regular way" and call the monitor. I don't think there is any special code associated with its use in multigrid so I would expect it to work. I guess you need to check in the debugger why it is not following the appropriate route.


  Also check that it is not running with preonly which has its own special case.

> On Mar 11, 2019, at 8:45 PM, Matthew Knepley via petsc-dev <petsc-dev at mcs.anl.gov> wrote:
> Does Richardson not respond to -mg_levels_ksp_monitor by default? If so, how do you turn it on?
>   Thanks,
>      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
> https://www.cse.buffalo.edu/~knepley/

More information about the petsc-dev mailing list