def _update_dPrxy(self):
super(ExpCM_empirical_phi_divpressure, self)._update_dPrxy()<EOL>if '<STR_LIT>' in self.freeparams:<EOL><INDENT>with scipy.errstate(divide='<STR_LIT>', under='<STR_LIT>', over='<STR_LIT>',<EOL>invalid='<STR_LIT:ignore>'):<EOL><INDENT>scipy.copyto(self.dPrxy['<STR_LIT>'], -self.ln_piAx_piAy_beta<EOL>* self.Qxy * self.omega /<EOL>(<NUM_LIT:1> - self.piAx_piAy_beta), where=CODON_NONSYN)<EOL><DEDENT>scipy.copyto(self.dPrxy['<STR_LIT>'], self.Qxy * self.omega,<EOL>where=scipy.logical_and(CODON_NONSYN, scipy.fabs(<NUM_LIT:1> -<EOL>self.piAx_piAy_beta) < ALMOST_ZERO))<EOL>for r in range(self.nsites):<EOL><INDENT>self.dPrxy['<STR_LIT>'][r] *= self.deltar[r]<EOL><DEDENT>_fill_diagonals(self.dPrxy['<STR_LIT>'], self._diag_indices)<EOL><DEDENT>
Update `dPrxy`, accounting for dependence of `Prxy` on `omega2`.