nojhan/eodev
Archived
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

BUGFIX: should use the standard deviation when sampling a mono-normal law, because it have the same scale than the mean

Browse source
This commit is contained in:
nojhan 2011-11-09 11:38:40 +01:00
commit 1af6715ddc
1 changed files with 8 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

16
edo/src/edoSamplerNormalMono.h
View file

@ -28,6 +28,8 @@ Authors:
#ifndef _edoSamplerNormalMono_h
#define _edoSamplerNormalMono_h
#include <cmath>
#include <utils/eoRNG.h>
#include "edoSampler.h"
@ -47,27 +49,25 @@ public:
edoSamplerNormalMono( edoRepairer<EOT> & repairer ) : edoSampler< D >( repairer) {}
EOT sample( edoNormalMono< EOT >& distrib )
EOT sample( edoNormalMono<EOT>& distrib )
{
unsigned int size = distrib.size();
assert(size > 0);
// Point we want to sample to get higher a set of points
// The point we want to draw
// (coordinates in n dimension)
// x = {x1, x2, ..., xn}
EOT solution;
// Sampling all dimensions
for (unsigned int i = 0; i < size; ++i)
{
for (unsigned int i = 0; i < size; ++i) {
AtomType mean = distrib.mean()[i];
AtomType variance = distrib.variance()[i];
AtomType random = rng.normal(mean, variance);
assert(variance >= 0);
// should use the standard deviation, which have the same scale than the mean
AtomType random = rng.normal(mean, sqrt(variance) );
solution.push_back(random);
}
}
return solution;
}