Refactor the code that substitutes unary ops in model equations

In particular, make it more visible that the substitution also occurs in PAC
equations (and not only in VAR and TCM equations).