diff --git a/src/DIS/integrandsIC.f b/src/DIS/integrandsIC.f index ef0d51cc..1985c05c 100644 --- a/src/DIS/integrandsIC.f +++ b/src/DIS/integrandsIC.f @@ -51,7 +51,7 @@ function integrandsICm(y) double precision integrandsICm * integrandsICm = 0d0 - if(y.ge.1d0) return + if(y.ge.one) return * * Interpolant functions * diff --git a/src/commons/wrapIC.h b/src/commons/wrapIC.h index e5bda1d0..7e5031a9 100644 --- a/src/commons/wrapIC.h +++ b/src/commons/wrapIC.h @@ -6,7 +6,7 @@ double precision m2strange parameter(m2strange = 0.01d0) double precision one - parameter(one=0.9999999999999999d0) + parameter(one=0.999999999999999d0) * double precision Q2IC,m12,m22,m1,m2 double precision Splus,Sminus,Rplus,Rminus