#include <NumericalIntegration.h>

Public Member Functions
double	error () const

	GaussIntegrator ()

	GaussIntegrator (double epsilon)

template<typename F >
double	operator() (const F &f, double a, double b) const

Private Attributes
double	aa

double	aconst

double	bb

double	c1

double	c2

double	epsilon_

double	error_

double	f1

double	f2

double	h

unsigned int	i

double	s16

double	s8

double	u

double	xx

Static Private Attributes
static const double	kCST = 5. / 1000

static const double	kHF = 0.5

static const double	w [12]

static const double	x [12]

Detailed Description

Definition at line 65 of file NumericalIntegration.h.

Constructor & Destructor Documentation

funct::GaussIntegrator::GaussIntegrator ( )

inline

Definition at line 67 of file NumericalIntegration.h.

67 {}

funct::GaussIntegrator::GaussIntegrator ( double epsilon )

inline

Definition at line 68 of file NumericalIntegration.h.

68 : epsilon_(epsilon) {}

geometryDiff.epsilon

int epsilon

Definition: geometryDiff.py:26

funct::GaussIntegrator::epsilon_

double epsilon_

Definition: NumericalIntegration.h:119

Member Function Documentation

double funct::GaussIntegrator::error ( ) const

inline

Definition at line 115 of file NumericalIntegration.h.

Referenced by Page1Parser.Page1Parser::check_for_whole_start_tag().

115 { return error_; }

funct::GaussIntegrator::error_

double error_

Definition: NumericalIntegration.h:118

template<typename F >

double funct::GaussIntegrator::operator()	(	const F &	f,
		double	a,
		double	b
	)		const

inline

Definition at line 70 of file NumericalIntegration.h.

References a, funct::abs(), b, alignmentValidation::c1, f, DeadROC_duringRun::f1, DeadROC_duringRun::f2, h, mps_fire::i, Gflash::kHF, w, and geometryCSVtoXML::xx.

                                                             {
       h = 0;
       if (b == a)
         return h;
       aconst = kCST / std::abs(b - a);
       bb = a;
     CASE1:
       aa = bb;
       bb = b;
     CASE2:
       c1 = kHF * (bb + aa);
       c2 = kHF * (bb - aa);
       s8 = 0;
       for (i = 0; i < 4; ++i) {
         u = c2 * x[i];
         xx = c1 + u;
         f1 = f(xx);
         xx = c1 - u;
         f2 = f(xx);
         s8 += w[i] * (f1 + f2);
       }
       s16 = 0;
       for (i = 4; i < 12; ++i) {
         u = c2 * x[i];
         xx = c1 + u;
         f1 = f(xx);
         xx = c1 - u;
         f2 = f(xx);
         s16 += w[i] * (f1 + f2);
       }
       s16 = c2 * s16;
       if (std::abs(s16 - c2 * s8) <= epsilon_ * (1. + std::abs(s16))) {
         h += s16;
         if (bb != b)
           goto CASE1;
       } else {
         bb = c1;
         if (1. + aconst * std::abs(c2) != 1)
           goto CASE2;
         h = s8;
       }
 
       error_ = std::abs(s16 - c2 * s8);
       return h;
     }

Member Data Documentation

double funct::GaussIntegrator::aa

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

Referenced by geometryXMLparser.Alignable::covariance().

double funct::GaussIntegrator::aconst

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

double funct::GaussIntegrator::bb

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

Referenced by geometryXMLparser.Alignable::covariance().

double funct::GaussIntegrator::c1

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

double funct::GaussIntegrator::c2

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

double funct::GaussIntegrator::epsilon_

private

Definition at line 119 of file NumericalIntegration.h.

double funct::GaussIntegrator::error_

mutableprivate

Definition at line 118 of file NumericalIntegration.h.

double funct::GaussIntegrator::f1

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

double funct::GaussIntegrator::f2

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

double funct::GaussIntegrator::h

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

unsigned int funct::GaussIntegrator::i

mutableprivate

Definition at line 124 of file NumericalIntegration.h.

const double funct::GaussIntegrator::kCST = 5. / 1000

staticprivate

Definition at line 122 of file NumericalIntegration.h.

Referenced by funct::GaussLegendreIntegrator::GaussLegendreIntegrator().

const double funct::GaussIntegrator::kHF = 0.5

staticprivate

Definition at line 121 of file NumericalIntegration.h.

Referenced by funct::GaussLegendreIntegrator::GaussLegendreIntegrator().

double funct::GaussIntegrator::s16

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

double funct::GaussIntegrator::s8

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

double funct::GaussIntegrator::u

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

const double funct::GaussIntegrator::w

staticprivate

Initial value:

= {0.10122853629037626,
                                              0.22238103445337447,
                                              0.31370664587788729,
                                              0.36268378337836198,
                                              0.02715245941175409,
                                              0.06225352393864789,
                                              0.09515851168249278,
                                              0.12462897125553387,
                                              0.14959598881657673,
                                              0.16915651939500254,
                                              0.18260341504492359,
                                              0.18945061045506850}

Definition at line 120 of file NumericalIntegration.h.

Referenced by funct::GaussLegendreIntegrator::GaussLegendreIntegrator().

const double funct::GaussIntegrator::x

staticprivate

Initial value:

= {0.96028985649753623,
                                              0.79666647741362674,
                                              0.52553240991632899,
                                              0.18343464249564980,
                                              0.98940093499164993,
                                              0.94457502307323258,
                                              0.86563120238783174,
                                              0.75540440835500303,
                                              0.61787624440264375,
                                              0.45801677765722739,
                                              0.28160355077925891,
                                              0.09501250983763744}

Definition at line 120 of file NumericalIntegration.h.

Referenced by svgfig.Curve.Sample::__repr__(), svgfig.Ellipse::__repr__(), Vispa.Gui.WidgetContainer.WidgetContainer::autosize(), Vispa.Gui.VispaWidget.VispaWidget::boundingRect(), funct::GaussLegendreIntegrator::GaussLegendreIntegrator(), geometryXMLparser.Alignable::pos(), Vispa.Gui.ConnectableWidget.ConnectableWidget::positionizeMenuWidget(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().

double funct::GaussIntegrator::xx

mutableprivate

Definition at line 123 of file NumericalIntegration.h.

Referenced by geometryXMLparser.Alignable::covariance().

Public Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation