Classical solution methods in fast Fourier transform-based computational micromechanics operate on, either, compatible strain fields or equilibrated stress fields. By contrast, polarization schemes are primal-dual whose iterates neither nor equilibrated. Recently, it was demonstrated that may outperform the classical methods. Unfortunately, their power critically depends on a judicious choice o...