From b092df6902aab2a23498b4aa85a35b757be72017 Mon Sep 17 00:00:00 2001
From: Olga Pustovalova <162949+olp-cs@users.noreply.github.com>
Date: Sat, 7 Jun 2025 09:50:22 +0000
Subject: [PATCH] Fix the formula for logp computation

---
 pymc/distributions/multivariate.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/pymc/distributions/multivariate.py b/pymc/distributions/multivariate.py
index ac2a13431a..a9f99d79c1 100644
--- a/pymc/distributions/multivariate.py
+++ b/pymc/distributions/multivariate.py
@@ -347,7 +347,7 @@ def precision_mv_normal_logp(op: PrecisionMvNormalRV, value, rng, size, mean, ta
     delta = value - mean
     quadratic_form = delta.T @ tau @ delta
     logdet, posdef = _logdet_from_cholesky(nan_lower_cholesky(tau))
-    logp = -0.5 * (k * pt.log(2 * np.pi) + quadratic_form) + logdet
+    logp = -0.5 * (k * pt.log(2 * np.pi) + quadratic_form + logdet)
 
     return check_parameters(
         logp,