[Bug-gsl] [bug #25413] [patch] monte test failure under mingw

[Brian Gough]

Follow-up Comment #5, bug #25413 (project gsl):

The first lines of your test.log show

The limits of integation are:

xl[0]=0.000000    xu[0]=1.000000

num_dim=1, calls=299, it_num=1, max_it_num=5 veb=1, alph=1.50,

The correct value of calls should be 300.  

There are only a few steps in the calculation of calls, if you step through
them you should find where it is wrong.  Starting from line 166 of vegas.c I
get

boxes=166 
box_per_bin=max(166/50,1)=3
bins=min(166/3,50)=50
boxes=50*3=150
tot_boxes=pow(150,1)=150
calls_per_box=max(333/150,2)=2
calls=2*150=300


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?25413>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/