36 #ifndef G4ErrorSymMatrix_hh
37 #define G4ErrorSymMatrix_hh
218 std::vector<G4double >
m;
299 #include "G4ErrorSymMatrix.icc"
void house_with_update2(G4ErrorSymMatrix *a, G4ErrorMatrix *v, G4int row=1, G4int col=1)
const G4double & fast(G4int row, G4int col) const
void invert4(G4int &ifail)
std::vector< ExP01TrackerHit * > a
friend void tridiagonal(G4ErrorSymMatrix *a, G4ErrorMatrix *hsm)
static const G4double CHOLESKY_THRESHOLD_5x5
std::ostream & operator<<(std::ostream &, const BasicVector3D< float > &)
friend void diag_step(G4ErrorSymMatrix *t, G4int begin, G4int end)
void invertCholesky6(G4int &ifail)
const G4ErrorSymMatrix & _a
void invertBunchKaufman(G4int &ifail)
G4ErrorSymMatrix_row(G4ErrorSymMatrix &, G4int)
static constexpr double m2
const G4double & operator()(G4int row, G4int col) const
friend G4ErrorMatrix operator*(const G4ErrorSymMatrix &m1, const G4ErrorSymMatrix &m2)
friend G4ErrorSymMatrix operator+(const G4ErrorSymMatrix &m1, const G4ErrorSymMatrix &m2)
G4ErrorSymMatrix apply(G4double(*f)(G4double, G4int, G4int)) const
virtual ~G4ErrorSymMatrix()
void invert(G4int &ifail)
friend G4double condition(const G4ErrorSymMatrix &m)
G4double condition(const G4ErrorSymMatrix &m)
void invert6(G4int &ifail)
static G4ThreadLocal G4double posDefFraction5x5
G4ErrorSymMatrix operator-() const
void invertCholesky5(G4int &ifail)
BasicVector3D< float > operator*(const BasicVector3D< float > &v, double a)
static G4ThreadLocal G4double posDefFraction6x6
static constexpr double m
void tridiagonal(G4ErrorSymMatrix *a, G4ErrorMatrix *hsm)
friend G4ErrorMatrix diagonalize(G4ErrorSymMatrix *s)
G4ErrorMatrix dsum(const G4ErrorMatrix &, const G4ErrorMatrix &)
BasicVector3D< float > operator+(const BasicVector3D< float > &v)
static G4ThreadLocal G4double adjustment6x6
G4ErrorSymMatrix similarity(const G4ErrorMatrix &m1) const
static G4ThreadLocal G4double adjustment5x5
G4double & operator[](G4int)
G4ErrorSymMatrix & operator+=(const G4ErrorSymMatrix &m2)
const G4double & operator[](G4int) const
std::vector< G4double > m
G4ErrorSymMatrix similarityT(const G4ErrorMatrix &m1) const
friend void house_with_update2(G4ErrorSymMatrix *a, G4ErrorMatrix *v, G4int row, G4int col)
void assign(const G4ErrorMatrix &m2)
static const G4double CHOLESKY_CREEP_6x6
G4ErrorSymMatrix & operator=(const G4ErrorSymMatrix &m2)
void invert5(G4int &ifail)
G4ErrorSymMatrix sub(G4int min_row, G4int max_row) const
G4ErrorSymMatrix_row_const(const G4ErrorSymMatrix &, G4int)
G4ErrorMatrix diagonalize(G4ErrorSymMatrix *s)
BasicVector3D< float > operator/(const BasicVector3D< float > &v, double a)
void invertHaywood4(G4int &ifail)
static const G4double CHOLESKY_THRESHOLD_6x6
G4double determinant() const
void diag_step(G4ErrorSymMatrix *t, G4int begin, G4int end)
static const G4double CHOLESKY_CREEP_5x5
G4ErrorSymMatrix T() const
void invertHaywood6(G4int &ifail)
void invertHaywood5(G4int &ifail)
G4ErrorSymMatrix & operator-=(const G4ErrorSymMatrix &m2)
G4ErrorSymMatrix & operator*=(G4double t)
G4ErrorSymMatrix inverse(G4int &ifail) const
BasicVector3D< float > operator-(const BasicVector3D< float > &v)
G4ErrorSymMatrix_row operator[](G4int)
G4ErrorSymMatrix & operator/=(G4double t)