pub const GMRES: KrylovMethod;Expand description
The dual solve by GMRES, for a nonsymmetric dual operator, restarting only after a generous 100 iterations since the solve converges in tens.
Not the default: conjugate gradients also refuses a singular subdomain through its positive-definiteness check, which GMRES has no way to do.