26 lines
380 B
C++
26 lines
380 B
C++
// $Id$
|
|
|
|
#include <halton.h>
|
|
|
|
HaltonSequenceGenerator_t::HaltonSequenceGenerator_t(int b)
|
|
{
|
|
base=b;
|
|
fbase=(float) b;
|
|
seed=1;
|
|
|
|
}
|
|
|
|
float HaltonSequenceGenerator_t::GetElement(int elem)
|
|
{
|
|
int tmpseed=seed;
|
|
float ret=0.0;
|
|
float base_inv=1.0/fbase;
|
|
while(tmpseed)
|
|
{
|
|
int dig=tmpseed % base;
|
|
ret+=((float) dig)*base_inv;
|
|
base_inv/=fbase;
|
|
tmpseed/=base;
|
|
}
|
|
return ret;
|
|
}
|