| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 // -------------------------------------------
28 // GEANT 4 class inline implementation
29 // -------------------------------------------
30
31 inline G4ErrorSymMatrix::G4ErrorSymMatrix()
32 : m(0)
33 , nrow(0)
34 , size(0)
35 {}
36
37 inline G4int G4ErrorSymMatrix::num_row() const
38
39 inline G4int G4ErrorSymMatrix::num_col() const
40
41 inline G4int G4ErrorSymMatrix::num_size() cons
42
43 inline G4double& G4ErrorSymMatrix::fast(G4int
44 {
45 return *(m.begin() + (row * (row - 1)) / 2 +
46 }
47
48 inline const G4double& G4ErrorSymMatrix::fast(
49 {
50 return *(m.begin() + (row * (row - 1)) / 2 +
51 }
52
53 inline G4double& G4ErrorSymMatrix::operator()(
54 {
55 return (row >= col ? fast(row, col) : fast(c
56 }
57
58 inline const G4double& G4ErrorSymMatrix::opera
59 {
60 return (row >= col ? fast(row, col) : fast(c
61 }
62
63 inline void G4ErrorSymMatrix::assign(const G4E
64 {
65 (*this) = mat;
66 }
67
68 inline G4ErrorSymMatrix G4ErrorSymMatrix::T()
69 {
70 return G4ErrorSymMatrix(*this);
71 }
72
73 inline G4ErrorSymMatrix::G4ErrorSymMatrix_row
74 G4int r)
75 {
76 G4ErrorSymMatrix_row b(*this, r);
77 return b;
78 }
79
80 inline G4ErrorSymMatrix::G4ErrorSymMatrix_row_
81 operator[](G4int r) const
82 {
83 const G4ErrorSymMatrix_row_const b(*this, r)
84 return b;
85 }
86
87 inline G4double& G4ErrorSymMatrix::G4ErrorSymM
88 {
89 if(_r >= c)
90 {
91 return *(_a.m.begin() + (_r + 1) * _r / 2
92 }
93 else
94 {
95 return *(_a.m.begin() + (c + 1) * c / 2 +
96 }
97 }
98
99 inline const G4double& G4ErrorSymMatrix::G4Err
100 G4int c) const
101 {
102 if(_r >= c)
103 {
104 return *(_a.m.begin() + (_r + 1) * _r / 2
105 }
106 else
107 {
108 return *(_a.m.begin() + (c + 1) * c / 2 +
109 }
110 }
111
112 inline G4ErrorSymMatrix::G4ErrorSymMatrix_row:
113 G4ErrorSymMatrix& a, G4int r)
114 : _a(a)
115 , _r(r)
116 {}
117
118 inline G4ErrorSymMatrix::G4ErrorSymMatrix_row_
119 const G4ErrorSymMatrix& a, G4int r)
120 : _a(a)
121 , _r(r)
122 {}
123
124 inline G4ErrorSymMatrix G4ErrorSymMatrix::inve
125 {
126 G4ErrorSymMatrix mTmp(*this);
127 mTmp.invert(ifail);
128 return mTmp;
129 }
130