Breakpoint 1 at 0x404f0f: file ../../gsl-1.12/monte/vegas.c, line 166. Starting program: C:\msys\1.0\home\gd\sw\src\gsl-1.12.obj\monte\.libs/test.exe [New thread 3276.0xd20] Breakpoint 1, gsl_monte_vegas_integrate (f=0x22feb0, xl=0x414020, xu=0x410080, dim=1, calls=333, r=0x3e3cc0, state=0x3e2470, result=0x22ce18, abserr=0x22ce20) at ../../gsl-1.12/monte/vegas.c:166 166 boxes = floor (pow (calls / 2.0, 1.0 / dim)); 1: boxes = 1 2: calls = 333 3: dim = 1 169 if (2 * boxes >= state->bins_max) 3: dim = 1 2: calls = 333 1: boxes = 1 4: state->bins_max = 50 166 boxes = floor (pow (calls / 2.0, 1.0 / dim)); 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 1 169 if (2 * boxes >= state->bins_max) 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 1 166 boxes = floor (pow (calls / 2.0, 1.0 / dim)); 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 1 169 if (2 * boxes >= state->bins_max) 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 166 172 int box_per_bin = GSL_MAX (boxes / state->bins_max, 1); 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 166 5: box_per_bin = 4072560 174 bins = GSL_MIN(boxes / box_per_bin, state->bins_max); 5: box_per_bin = 3 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 166 6: bins = 50 175 boxes = box_per_bin * bins; 6: bins = 50 5: box_per_bin = 3 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 166 177 state->mode = GSL_VEGAS_MODE_STRATIFIED; 6: bins = 50 5: box_per_bin = 3 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 3 175 boxes = box_per_bin * bins; 6: bins = 50 5: box_per_bin = 3 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 3 177 state->mode = GSL_VEGAS_MODE_STRATIFIED; 6: bins = 50 5: box_per_bin = 3 4: state->bins_max = 50 3: dim = 1 2: calls = 333 1: boxes = 150 182 double tot_boxes = pow ((double) boxes, (double) dim); 4: state->bins_max = 50 3: dim = 1 2: calls = 333 7: tot_boxes = 194 183 state->calls_per_box = GSL_MAX (calls / tot_boxes, 2); 7: tot_boxes = 0.22834175726319428 4: state->bins_max = 50 3: dim = 1 2: calls = 333 8: state->calls_per_box = 3131961357 182 double tot_boxes = pow ((double) boxes, (double) dim); 8: state->calls_per_box = 3131961357 7: tot_boxes = 0.22834175726319428 4: state->bins_max = 50 3: dim = 1 2: calls = 333 9: boxes = 150 183 state->calls_per_box = GSL_MAX (calls / tot_boxes, 2); 9: boxes = 150 8: state->calls_per_box = 3131961357 7: tot_boxes = 0 4: state->bins_max = 50 3: dim = 1 2: calls = 333 184 calls = state->calls_per_box * tot_boxes; 9: boxes = 150 8: state->calls_per_box = 2 7: tot_boxes = 0.22834175726319428 4: state->bins_max = 50 3: dim = 1 2: calls = 333 188 state->jac = state->vol * pow ((double) bins, (double) dim) / calls; 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 333 10: state->jac = -4.8366978272229995e-026 11: state->vol = 1 12: bins = 50 184 calls = state->calls_per_box * tot_boxes; 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 333 13: tot_boxes = 0.22834175726319428 188 state->jac = state->vol * pow ((double) bins, (double) dim) / calls; 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 333 184 calls = state->calls_per_box * tot_boxes; 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 333 14: tot_boxes = 0 188 state->jac = state->vol * pow ((double) bins, (double) dim) / calls; 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 299 R7: Empty 0x400795ffffffffffffff R6: Empty 0x40008000000000000000 R5: Empty 0x40008e147ae147ae147c R4: Empty 0x3fff8000000000000000 R3: Empty 0x4002a000000000000000 R2: Empty 0x4006c200000000000000 R1: Empty 0x3ffce9d26bef06bf7dff =>R0: Empty 0x00000027fd8c00000d68 Status Word: 0xffff0020 PE TOP: 0 Control Word: 0xffff037f IM DM ZM OM UM PM PC: Extended Precision (64-bits) RC: Round to nearest Tag Word: 0xffffffff Instruction Pointer: 0x1b:0x00404324 Operand Pointer: 0xffff0023:0x0022cd20 Opcode: 0xdf7d 196 if (bins != state->bins) 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 299 188 state->jac = state->vol * pow ((double) bins, (double) dim) / calls; 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 299 190 state->boxes = boxes; 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 299 196 if (bins != state->bins) 12: bins = 50 11: state->vol = 1 10: state->jac = -4.8366978272229995e-026 9: boxes = 150 8: state->calls_per_box = 2 4: state->bins_max = 50 3: dim = 1 2: calls = 299