Changeset 392 for library/bdm/stat/merger.h
- Timestamp:
- 06/19/09 17:57:53 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/stat/merger.h
r388 r392 81 81 void set_sources (const Array<mpdf*> &Sources, bool own) { 82 82 compositepdf::set_elements (Sources,own); 83 Nsources=mpdfs.length(); 83 84 //set sizes 84 85 dls.set_size (Sources.length()); … … 128 129 smpi (j) = XYZ (j) (0); /* beginning of the interval*/ 129 130 it_assert(gridsize(j)!=0.0,"Zeros in gridsize!"); 130 steps (j) = ( smpi(j) - XYZ(j)(1))/gridsize(j);131 steps (j) = ( XYZ(j)(1)-smpi(j) )/gridsize(j); 131 132 } 132 133 // fill samples … … 138 139 for (int j = 0;j < dim;j++) { 139 140 if (ind (j) == gridsize (j) - 1) { //j-th index is full 140 ind (j) = 0; //shift back141 // ind (j) = 0; //shift back 141 142 smpi(j) = XYZ(j)(0); 142 143 143 ind (j + 1) ++; //increase the next dimension;144 // ind (j + 1) ++; //increase the next dimension; 144 145 smpi(j+1) += steps(j+1); 145 146 146 147 if (ind (j + 1) < gridsize (j + 1) - 1) break; 147 148 } else { 148 ind (j) ++;149 // ind (j) ++; 149 150 smpi(j) +=steps(j); 150 151 break; … … 272 273 set_method (GEOMETRIC); 273 274 else if (!strcmp (meth_str.c_str(), "lognormal")) { 274 set_method ( GEOMETRIC);275 set_method (LOGNORMAL); 275 276 set.lookupValue( "beta",beta); 276 277 } 277 278 } 278 validate();279 279 } 280 280