Change of the algorithm to calculate gap integrals
This commit is contained in:
parent
aab2c7a717
commit
7f6911d7ac
@ -74,7 +74,7 @@ TGapDWave::TGapDWave() {
|
||||
}
|
||||
|
||||
TGapAnSWave::TGapAnSWave() {
|
||||
TAnSWaveGapIntegralSuave *gapint = new TAnSWaveGapIntegralSuave();
|
||||
TAnSWaveGapIntegralCuhre *gapint = new TAnSWaveGapIntegralCuhre();
|
||||
fGapIntegral = gapint;
|
||||
gapint = 0;
|
||||
delete gapint;
|
||||
|
2
src/external/libGapIntegrals/TGapIntegrals.h
vendored
2
src/external/libGapIntegrals/TGapIntegrals.h
vendored
@ -89,7 +89,7 @@ public:
|
||||
double operator()(double, const vector<double>&) const;
|
||||
|
||||
private:
|
||||
TAnSWaveGapIntegralSuave *fGapIntegral;
|
||||
TAnSWaveGapIntegralCuhre *fGapIntegral;
|
||||
mutable vector<double> fTemp;
|
||||
mutable vector<double>::const_iterator fTempIter;
|
||||
mutable vector<double> fIntegralValues;
|
||||
|
6
src/external/libGapIntegrals/TIntegrator.cpp
vendored
6
src/external/libGapIntegrals/TIntegrator.cpp
vendored
@ -75,7 +75,7 @@ double TAnSWaveGapIntegralCuhre::IntegrateFunc()
|
||||
const double EPSABS (1e-6);
|
||||
const unsigned int VERBOSE (0);
|
||||
const unsigned int LAST (4);
|
||||
const unsigned int MINEVAL (10000);
|
||||
const unsigned int MINEVAL (1000);
|
||||
const unsigned int MAXEVAL (1000000);
|
||||
|
||||
const unsigned int KEY (13);
|
||||
@ -107,7 +107,7 @@ double TAnSWaveGapIntegralDivonne::IntegrateFunc()
|
||||
const double EPSREL (1e-4);
|
||||
const double EPSABS (1e-6);
|
||||
const unsigned int VERBOSE (0);
|
||||
const unsigned int MINEVAL (10000);
|
||||
const unsigned int MINEVAL (1000);
|
||||
const unsigned int MAXEVAL (1000000);
|
||||
const unsigned int KEY1 (47);
|
||||
const unsigned int KEY2 (1);
|
||||
@ -149,7 +149,7 @@ double TAnSWaveGapIntegralSuave::IntegrateFunc()
|
||||
const double EPSABS (1e-6);
|
||||
const unsigned int VERBOSE (0);
|
||||
const unsigned int LAST (4);
|
||||
const unsigned int MINEVAL (10000);
|
||||
const unsigned int MINEVAL (1000);
|
||||
const unsigned int MAXEVAL (1000000);
|
||||
|
||||
const unsigned int NNEW (1000);
|
||||
|
Loading…
x
Reference in New Issue
Block a user