33 #define INCLXX_IN_GEANT4_MODE 1
53 const G4double IFunction1D::integrationCoefficients[] = {
82 nIntervals =
G4int(interval/step);
85 G4int remainder = nIntervals % 9;
87 nIntervals += 9 - remainder;
89 nIntervals = std::max(nIntervals, 9);
92 const G4double dx = interval/nIntervals;
94 for(
G4int j = 1; j<nIntervals; ++j) {
96 const unsigned index = j%9;
100 return result*dx*
sign;
113 return theFunction->integrate(
xMin,x);
117 } *thePrimitive =
new Primitive(
this);
128 normalisation(1./theFunction->integrate(
xMin,
xMax))
132 return std::min(1., normalisation * theFunction->integrate(
xMin,x));
137 } *theInverseCDF =
new InverseCDF(
this);
140 delete theInverseCDF;