From fc7e9967fc443e5b740629df1dac3f22828f1179 Mon Sep 17 00:00:00 2001
From: Mark Woolrich <woolrich@fmrib.ox.ac.uk>
Date: Wed, 29 Aug 2007 14:48:52 +0000
Subject: [PATCH] *** empty log message ***
---
miscprob.cc | 4 ++++
miscprob.h | 1 +
2 files changed, 5 insertions(+)
diff --git a/miscprob.cc b/miscprob.cc
index b88eaa0..633267a 100644
--- a/miscprob.cc
+++ b/miscprob.cc
@@ -172,6 +172,10 @@ ReturnMatrix gammapdf(const RowVector& vals, const float mu, const float var)
return res;
}
+ float normpdf(const float val, const float mu, const float var)
+ {
+ return std::exp(-0.5*(std::pow(val-mu,2)/var))*std::pow(2*M_PI*var,-0.5);
+ }
ReturnMatrix normpdf(const RowVector& vals, const RowVector& mu, const RowVector& var)
{
diff --git a/miscprob.h b/miscprob.h
index a801e91..100eda4 100644
--- a/miscprob.h
+++ b/miscprob.h
@@ -35,6 +35,7 @@ namespace MISCMATHS {
// returns nsamps*nparams matrix:
ReturnMatrix mvnrnd(const RowVector& mu, const SymmetricMatrix& covar, int nsamp = 1);
+ float normpdf(const float val, const float mu = 0, const float var = 1);
ReturnMatrix normpdf(const RowVector& vals, const float mu = 0, const float var = 1);
--
GitLab