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/
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
