8b621a409cd07d0f5610a2648b6413322c39a822,GPy/likelihoods/likelihood.py,Likelihood,ep_gradients,#Likelihood#,230

Before Change

                    Y_metadata_list.append(Y_metadata_i)

                val = self.site_derivatives_ep(Y[index], tau[index], v[index], Y_metadata_i)
                dlik_dtheta[:, index] = val.ravel()

            f = partial(self.integrate)
            quads = zip(*map(f, Y.flatten(), mu.flatten(), np.sqrt(sigma2.flatten())))

After Change

                    Y_metadata_list.append(Y_metadata_i)

            if quad_mode == "gk":
                f = partial(self.integrate_gk)
                quads = zip(*map(f, Y.flatten(), mu.flatten(), np.sqrt(sigma2.flatten()), Y_metadata_list))
                quads = np.vstack(quads)
                quads.reshape(self.size, shape[0], shape[1])
            elif quad_mode == "gh":

Instances