Open
Description
Describe the issue:
At this line,
logp = -0.5 * (k * pt.log(2 * np.pi) + quadratic_form) + logdet
should be
logp = -0.5 * (k * pt.log(2 * np.pi) + quadratic_form + logdet)
Reproduceable code example:
N/A
Error message:
PyMC version information:
5.21.1
Context for the issue:
No response
Activity
ricardoV94 commentedon Apr 9, 2025
@ckrapu wanna open a PR by any chance?
shanegladson commentedon Apr 28, 2025
This issue seems inactive so I'm happy to take it over if available
PrecisionMvNormalRV
missing factor of 1/2 for log determinant inlogp
function #7811