sincos_poly.cxx File Reference

#include "GeoSpecialShapes/LArWheelCalculator.h"
#include "CxxUtils/sincos.h"
#include "TMath.h"
#include "TMatrixD.h"
#include "TVectorD.h"
#include "TMatrixDLazy.h"
#include "TDecompLU.h"
#include "TDecompSVD.h"
#include <sys/time.h>
#include <iostream>
#include <iomanip>
#include <math.h>
#include <mutex>
#include "GeoModelKernel/Units.h"

Go to the source code of this file.

Macros
#define	DEBUGPRINT   0

Typedefs
typedef TVectorD	VECTOR
typedef TMatrixD	MATRIX
typedef TMatrixDSym	SYMMATRIX

Functions
template<typename T, typename Q>
std::ostream &	operator<< (std::ostream &ostr, const TVectorT< T > &v)
static VECTOR	findLinearApproximation (const Int_t dataLen, const Int_t nBasisFuntions, const VECTOR &y, const MATRIX &bf)

Macro Definition Documentation

DEBUGPRINT

#define DEBUGPRINT   0

Definition at line 38 of file sincos_poly.cxx.

Typedef Documentation

MATRIX

typedef TMatrixD MATRIX

Definition at line 21 of file sincos_poly.cxx.

SYMMATRIX

typedef TMatrixDSym SYMMATRIX

Definition at line 22 of file sincos_poly.cxx.

VECTOR

typedef TVectorD VECTOR

Definition at line 20 of file sincos_poly.cxx.

Function Documentation

findLinearApproximation()

VECTOR findLinearApproximation ( const Int_t dataLen, const Int_t nBasisFuntions, const VECTOR & y, const MATRIX & bf )

static

Definition at line 60 of file sincos_poly.cxx.

{
#ifndef SINCOSPOLY_USE_EIGEN
  SYMMATRIX A(nBasisFuntions);
  VECTOR    vY(nBasisFuntions);
#else
  SYMMATRIX A=SYMMATRIX::Zero(nBasisFuntions,nBasisFuntions);
  VECTOR    vY=VECTOR::Zero(nBasisFuntions);
#endif
  for(Int_t j = 0; j < nBasisFuntions; ++ j){
    for(Int_t k = 0; k < nBasisFuntions; ++ k){
      Double_t Ajk = 0.0;
      for(Int_t i = 0; i < dataLen; ++ i){
        Ajk += bf(j, i) * bf(k, i);
      }
      A(j, k) = Ajk;
    }
  }
 
  for(Int_t k = 0; k < nBasisFuntions; ++ k){
    Double_t vYk = 0.0;
    for(Int_t i = 0; i < dataLen; ++ i){
      vYk += y[i]*bf(k,i);
    }
    vY[k] = vYk;
  }
#ifndef SINCOSPOLY_USE_EIGEN 
  SYMMATRIX Ainv(A);
  Ainv.Invert();
  return Ainv*vY;
#else
  return A.inverse()*vY;
#endif
}

operator<<()

template<typename T, typename Q>

std::ostream & operator<<	(	std::ostream &	ostr,
		const TVectorT< T > &	v )

Definition at line 44 of file sincos_poly.cxx.

{
  std::ios_base::fmtflags save_flags(ostr.flags());
  ostr << '[';
  ostr << std::scientific;
  for(Int_t idx=v.GetLwb();idx<v.GetUpb();idx++) {
    ostr << v[idx] << ", ";
  }
  ostr << v[v.GetUpb()];
  ostr << ']';
  ostr.flags(save_flags);
  return ostr;
}