Legend:
- Unmodified
- Added
- Removed
-
bdm/stat/libEF.cpp
r102 r129 59 59 for ( i=0; i<rv.count(); i++ ) { 60 60 GamRNG.setup ( alpha ( i ),beta ( i ) ); 61 #pragma omp critical 61 62 smp ( i ) = GamRNG(); 62 63 } -
bdm/stat/libEF.h
r124 r129 205 205 double evalpdflog ( const vec &val ) const {return lnk;} 206 206 vec sample() const { 207 vec smp ( rv.count() ); UniRNG.sample_vector ( rv.count(),smp ); 207 vec smp ( rv.count() ); 208 #pragma omp critical 209 UniRNG.sample_vector ( rv.count(),smp ); 208 210 return low+elem_mult(distance,smp); 209 211 }