One further caveat is that known symmetric systems (either because you have set MAT_SYMMETRIC, or the matrix format is SBAIJ) are treated with ICC(0), bjacobi/ICC(0), and Cholesky inside the redundant solve in FMG. Note that this is probably not what you want if you have a symmetric indefinite system. Jed