- //avoid a disaster value - c'est du bricolage mais ca marche
-/* if(x == 0) x = 10e-8; */
- return var->df/x;
-}
-
-/*
- * For Vegas fpip: $-\frac{\alpha D_f}{x^2}$
- */
-double func_vegas_fpip(lmm_variable_t var, double x){
- //avoid a disaster value - c'est du bricolage mais ca marche
-/* if(x == 0) x = 10e-8; */
- return -( var->df/(x*x) ) ;
-}
-
-
-/*
- * For Reno f: $\frac{\sqrt{\frac{3}{2}}}{D_f} \arctan\left(\sqrt{\frac{3}{2}}x_f D_f\right)$
- */
-double func_reno_f(lmm_variable_t var, double x){
- xbt_assert0(var->df>0.0,"Don't call me with stupid values!");
- // \sqrt{3/2} = 0.8164965808
- return (0.8164965808 / var->df) * atan( (0.8164965808 / var->df)*x );
-}
-
-/*
- * For Reno fp: $\frac{3}{3 {D_f}^2 x^2 + 2}$
- */
-double func_reno_fp(lmm_variable_t var, double x){
- return 3 / (3*var->df*var->df*x*x + 2);
+ xbt_assert0(x>0.0,"Don't call me with stupid values!");
+ return VEGAS_SCALING*var->df/x;