Mercurial > hg > qm-dsp

diff hmm/hmm.c @ 487:5998ee1042d3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use M_PI
author Chris Cannam <cannam@all-day-breakfast.com>
date Fri, 31 May 2019 16:33:55 +0100
parents d48276a3ae24
children
line wrap: on
line diff
--- a/hmm/hmm.c	Fri May 31 16:33:46 2019 +0100
+++ b/hmm/hmm.c	Fri May 31 16:33:55 2019 +0100
@@ -627,7 +627,7 @@
 
     s = cblas_ddot(L, z, 1, y, 1);
 	
-    return exp(-s/2.0) / (pow(2*PI, L/2.0) * sqrt(detcov));
+    return exp(-s/2.0) / (pow(2 * M_PI, L/2.0) * sqrt(detcov));
 }
 
 /* log probability of multivariate Gaussian given mean, inverse and determinant of covariance */
@@ -647,7 +647,7 @@
 
     s = cblas_ddot(L, z, 1, y, 1);
 	
-    ret = -0.5 * (s + L * log(2*PI) + log(detcov));
+    ret = -0.5 * (s + L * log(2 * M_PI) + log(detcov));
 	
     return ret;
 }