36 : fFunction(pFunction),
44 : fFunction(pFunction),
45 fTolerance(pTolerance),
60 G4int iterationNumber )
62 G4double Step = (xFinal - xInitial)/iterationNumber ;
65 for(
G4int i=1;i<iterationNumber;i++)
76 G4int iterationNumber )
78 G4double Step = (xFinal - xInitial)/iterationNumber ;
81 for(
G4int i=1;i<iterationNumber;i++)
92 G4int iterationNumber )
95 static const G4double root = 1.0/std::sqrt(3.0) ;
96 G4double Step = (xFinal - xInitial)/(2.0*iterationNumber) ;
99 for(
G4int i=0;i<iterationNumber;i++)
110 G4int iterationNumber )
112 G4double Step = (xFinal - xInitial)/iterationNumber ;
114 G4double xPlus = xInitial + 0.5*Step ;
117 for(
G4int i=1;i<iterationNumber;i++)
125 return mean*Step/3.0 ;
147 static const G4double root = 1.0/std::sqrt(3.0) ;
149 G4double xMean = (xInitial + xFinal)/2.0 ;
150 G4double Step = (xFinal - xInitial)/2.0 ;
166 G4Exception(
"G4SimpleIntegration::AdaptGauss()",
"Error",
169 G4double xMean = (xInitial + xFinal)/2.0 ;
173 if(std::fabs(leftHalf+rightHalf-full) <
fTolerance)