better sparse logdet()

lbl-camera · Dec 15, 2022 · 051e24a · 051e24a
1 parent ba5f9ed
commit 051e24a
Show file tree

Hide file tree

Showing 3 changed files with 529 additions and 5 deletions.
diff --git a/fvgp/gp.py b/fvgp/gp.py
@@ -565,8 +565,9 @@ def log_likelihood(self,hyperparameters):
         y = self.y_data - mean
         sign, logdet = self.slogdet(K)
         n = len(y)
-        if sign == 0.0: return (0.5 * (y.T @ x)) + (0.5 * n * np.log(2.0*np.pi))
-        return (0.5 * (y.T @ x)) + (0.5 * sign * logdet) + (0.5 * n * np.log(2.0*np.pi))
+        #if sign == 0.0: return (0.5 * (y.T @ x)) + (0.5 * n * np.log(2.0*np.pi))
+        #return (0.5 * (y.T @ x)) + (0.5 * sign * logdet) + (0.5 * n * np.log(2.0*np.pi))
+        return (0.5 * (y.T @ x)) + (0.5 * logdet) + (0.5 * n * np.log(2.0*np.pi))
     ##################################################################################
     def log_likelihood_gradient(self, hyperparameters):
         """