The current implementation treats the two-electron part of the BP Hamiltonian via the SOMF approximation. We discuss special aspects of the application of SOMF to openshell references and propose practical solutions. In particular, L + and L -reduced matrix elements are slightly different. This issue can be solved by averaging. Contribution of the (unphysical) singlet part of the transition density can be important. This issue does not occur if only the triplet component is used.