def dpsi2_dmuS(self, dL_dpsi2, Z, mu, S, target_mu, target_S):
self._psi_computations(Z, mu, S)<EOL>tmp = (self.inv_lengthscale2 * self._psi2[:, :, :, None]) / self._psi2_denom<EOL>target_mu += -<NUM_LIT>*(dL_dpsi2[:, :, :, None] * tmp * self._psi2_mudist).sum(<NUM_LIT:1>).sum(<NUM_LIT:1>)<EOL>target_S += (dL_dpsi2[:, :, :, None] * tmp * (<NUM_LIT>*self._psi2_mudist_sq - <NUM_LIT:1>)).sum(<NUM_LIT:1>).sum(<NUM_LIT:1>)<EOL>
Think N,num_inducing,num_inducing,input_dim