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