def _update_dPrxy(self):
if '<STR_LIT>' in self.freeparams:<EOL><INDENT>scipy.copyto(self.dPrxy['<STR_LIT>'], self.Prxy / self.kappa,<EOL>where=CODON_TRANSITION)<EOL>_fill_diagonals(self.dPrxy['<STR_LIT>'], self._diag_indices)<EOL><DEDENT>if '<STR_LIT>' in self.freeparams:<EOL><INDENT>scipy.copyto(self.dPrxy['<STR_LIT>'], self.Frxy_no_omega * self.Qxy,<EOL>where=CODON_NONSYN)<EOL>_fill_diagonals(self.dPrxy['<STR_LIT>'], self._diag_indices)<EOL><DEDENT>if '<STR_LIT>' in self.freeparams:<EOL><INDENT>self.dPrxy['<STR_LIT>'].fill(<NUM_LIT:0>)<EOL>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.Prxy *<EOL>(<NUM_LIT:1> / self.beta + (self.piAx_piAy_beta *<EOL>(self.ln_piAx_piAy_beta / self.beta) /<EOL>(<NUM_LIT:1> - self.piAx_piAy_beta))), where=CODON_NONSYN)<EOL><DEDENT>scipy.copyto(self.dPrxy['<STR_LIT>'], self.Prxy/self.beta *<EOL>(<NUM_LIT:1> - self.piAx_piAy_beta), where=scipy.logical_and(<EOL>CODON_NONSYN, scipy.fabs(<NUM_LIT:1> - self.piAx_piAy_beta)<EOL>< ALMOST_ZERO))<EOL>_fill_diagonals(self.dPrxy['<STR_LIT>'], self._diag_indices)<EOL><DEDENT>if '<STR_LIT>' in self.freeparams:<EOL><INDENT>for i in range(N_NT - <NUM_LIT:1>):<EOL><INDENT>for w in range(i, N_NT):<EOL><INDENT>scipy.copyto(self.dPrxy['<STR_LIT>'][i], self.Prxy / (self.eta[i]<EOL>- int(i == w)), where=CODON_NT_MUT[w])<EOL><DEDENT>_fill_diagonals(self.dPrxy['<STR_LIT>'][i], self._diag_indices)<EOL><DEDENT><DEDENT>
Update `dPrxy`.