Geant4  10.00.p03
G4ErrorSymMatrix.icc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4ErrorSymMatrix.icc 66892 2013-01-17 10:57:59Z gunter $
27 //
28 // ------------------------------------------------------------
29 // GEANT 4 class inline implementation
30 // ------------------------------------------------------------
31 
32 
33 inline G4ErrorSymMatrix::G4ErrorSymMatrix()
34  : m(0), nrow(0), size(0)
35 {
36 }
37 
38 inline G4int G4ErrorSymMatrix::num_row() const
39 {
40  return nrow;
41 }
42 
43 inline G4int G4ErrorSymMatrix::num_col() const
44 {
45  return nrow;
46 }
47 
48 inline G4int G4ErrorSymMatrix::num_size() const
49 {
50  return size;
51 }
52 
53 inline G4double & G4ErrorSymMatrix::fast(G4int row,G4int col)
54 {
55  return *(m.begin()+(row*(row-1))/2+(col-1));
56 }
57 
58 inline const G4double & G4ErrorSymMatrix::fast(G4int row,G4int col) const
59 {
60  return *(m.begin()+(row*(row-1))/2+(col-1));
61 }
62 
63 inline G4double &
64 G4ErrorSymMatrix::operator()(G4int row, G4int col)
65 {
66  return (row>=col? fast(row,col) : fast(col,row));
67 }
68 
69 inline const G4double &
70 G4ErrorSymMatrix::operator()(G4int row, G4int col) const
71 {
72  return (row>=col? fast(row,col) : fast(col,row));
73 }
74 
75 inline void G4ErrorSymMatrix::assign(const G4ErrorSymMatrix &mat)
76 {
77  (*this)=mat;
78 }
79 
80 inline G4ErrorSymMatrix G4ErrorSymMatrix::T() const
81 {
82  return G4ErrorSymMatrix(*this);
83 }
84 
85 inline G4ErrorSymMatrix::G4ErrorSymMatrix_row
86 G4ErrorSymMatrix::operator[] (G4int r)
87 {
88  G4ErrorSymMatrix_row b(*this,r);
89  return b;
90 }
91 
92 inline G4ErrorSymMatrix::G4ErrorSymMatrix_row_const
93 G4ErrorSymMatrix::operator[] (G4int r) const
94 {
95  const G4ErrorSymMatrix_row_const b(*this,r);
96  return b;
97 }
98 
99 inline G4double &G4ErrorSymMatrix::G4ErrorSymMatrix_row::operator[](G4int c)
100 {
101  if (_r >= c )
102  {
103  return *(_a.m.begin() + (_r+1)*_r/2 + c);
104  }
105  else
106  {
107  return *(_a.m.begin() + (c+1)*c/2 + _r);
108  }
109 }
110 
111 inline const G4double &
112 G4ErrorSymMatrix::G4ErrorSymMatrix_row_const::operator[](G4int c) const
113 {
114  if (_r >= c )
115  {
116  return *(_a.m.begin() + (_r+1)*_r/2 + c);
117  }
118  else
119  {
120  return *(_a.m.begin() + (c+1)*c/2 + _r);
121  }
122 }
123 
124 inline G4ErrorSymMatrix::
125 G4ErrorSymMatrix_row::G4ErrorSymMatrix_row(G4ErrorSymMatrix &a, G4int r)
126  : _a(a), _r(r)
127 {}
128 
129 inline G4ErrorSymMatrix::
130 G4ErrorSymMatrix_row_const::G4ErrorSymMatrix_row_const
131 (const G4ErrorSymMatrix&a,G4int r)
132  : _a(a), _r(r)
133 {}
134 
135 inline G4ErrorSymMatrix G4ErrorSymMatrix::inverse(G4int &ifail) const
136 {
137  G4ErrorSymMatrix mTmp(*this);
138  mTmp.invert(ifail);
139  return mTmp;
140 }