On 22 February 2016 at 23:16, Timothée Nicolas <timothee.nico...@gmail.com> wrote:
> Hi all, > > It sounds it should be obvious but I can't find it somehow. I would like > to use weighted jacobi (more precisely point block jacobi) as a smoother > for multigrid, but I don't find the options to set a weight omega, after I > set -mg_levels_pc_type pbjacobi > Why not use -ksp_type richardson -ksp_richardson_scale XXX ? Cheers, Dave So I was thinking I could just artificially scale the second matrix in > KSPSetOperators with a weight omega. Since it is the one used in retrieving > the diagonal, it should be equivalent to have a weight and should not > affect the residuals, right ? However, I notice no change of behaviour > whatsoever when I do this. > > Is it normal ? Have I grossly misunderstood something ? > > Thanks > > Timothee >
