def _K_compute_ode_eq(self, transpose=False):
if self._t2 is not None:<EOL><INDENT>if transpose:<EOL><INDENT>t_eq = self._t[self._index==<NUM_LIT:0>]<EOL>t_ode = self._t2[self._index2><NUM_LIT:0>]<EOL>index_ode = self._index2[self._index2><NUM_LIT:0>]-<NUM_LIT:1><EOL><DEDENT>else:<EOL><INDENT>t_eq = self._t2[self._index2==<NUM_LIT:0>]<EOL>t_ode = self._t[self._index><NUM_LIT:0>]<EOL>index_ode = self._index[self._index><NUM_LIT:0>]-<NUM_LIT:1><EOL><DEDENT><DEDENT>else:<EOL><INDENT>t_eq = self._t[self._index==<NUM_LIT:0>]<EOL>t_ode = self._t[self._index><NUM_LIT:0>]<EOL>index_ode = self._index[self._index><NUM_LIT:0>]-<NUM_LIT:1><EOL><DEDENT>if t_ode.size==<NUM_LIT:0> or t_eq.size==<NUM_LIT:0>:<EOL><INDENT>if transpose:<EOL><INDENT>self._K_eq_ode = np.zeros((t_eq.shape[<NUM_LIT:0>], t_ode.shape[<NUM_LIT:0>]))<EOL><DEDENT>else:<EOL><INDENT>self._K_ode_eq = np.zeros((t_ode.shape[<NUM_LIT:0>], t_eq.shape[<NUM_LIT:0>]))<EOL><DEDENT>return<EOL><DEDENT>t_ode_mat = t_ode[:, None]<EOL>t_eq_mat = t_eq[None, :]<EOL>if self.delay is not None:<EOL><INDENT>t_ode_mat -= self.delay[index_ode, None]<EOL><DEDENT>diff_t = (t_ode_mat - t_eq_mat)<EOL>inv_sigma_diff_t = <NUM_LIT:1.>/self.sigma*diff_t<EOL>decay_vals = self.decay[index_ode][:, None]<EOL>half_sigma_d_i = <NUM_LIT:0.5>*self.sigma*decay_vals<EOL>if self.is_stationary:<EOL><INDENT>ln_part, signs = ln_diff_erfs(inf, half_sigma_d_i - inv_sigma_diff_t, return_sign=True)<EOL><DEDENT>else:<EOL><INDENT>ln_part, signs = ln_diff_erfs(half_sigma_d_i + t_eq_mat/self.sigma, half_sigma_d_i - inv_sigma_diff_t, return_sign=True)<EOL><DEDENT>sK = signs*np.exp(half_sigma_d_i*half_sigma_d_i - decay_vals*diff_t + ln_part)<EOL>sK *= <NUM_LIT:0.5><EOL>if not self.is_normalized:<EOL><INDENT>sK *= np.sqrt(np.pi)*self.sigma<EOL><DEDENT>if transpose:<EOL><INDENT>self._K_eq_ode = sK.T<EOL><DEDENT>else:<EOL><INDENT>self._K_ode_eq = sK<EOL><DEDENT>
Compute the cross covariances between latent exponentiated quadratic and observed ordinary differential equations.

        :param transpose: if set to false the exponentiated quadratic is on the rows of the matrix and is computed according to self._t, if set to true it is on the columns and is computed according to self._t2 (default=False).
        :type transpose: bool