Hi folks, Has anybody tried AuxiliaryOperatorPC with coefficients dependent on the Newton state i.e. using the variable being solved for? I’m getting some quite odd behaviour (sometimes NaNs or badly solved direct solve on the PC) which I haven’t quite pinned down yet. All the best --cjc