Confusing phase numbering
In e.g. the 2p2c model, the phase numbering is a bit confusing. The formulation allows for p0s1 or p1s0, but the setState() function for the priVars allows firstPhaseOnly and secondPhaseOnly. In addition, it should be checked, if the privarsswitchnewton could print the phase's names instead of something like "First phase dissapears .... s0 < 0".