Geant4 Cross Reference |
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 // G4Physics2DVector inline implementation 27 // 28 // Author: Vladimir Ivanchenko, 25.09.2011 29 // -------------------------------------------------------------------- 30 31 inline G4double G4Physics2DVector::Value(G4double x, G4double y) const 32 { 33 std::size_t idx = 0; 34 std::size_t idy = 0; 35 return Value(x, y, idx, idy); 36 } 37 38 inline void G4Physics2DVector::PutX(std::size_t idx, G4double val) 39 { 40 xVector[idx] = val; 41 } 42 43 inline void G4Physics2DVector::PutY(std::size_t idy, G4double val) 44 { 45 yVector[idy] = val; 46 } 47 48 inline void G4Physics2DVector::PutValue(std::size_t idx, std::size_t idy, 49 G4double val) 50 { 51 (*(value[idy]))[idx] = val; 52 } 53 54 inline G4double G4Physics2DVector::GetX(std::size_t index) const 55 { 56 return xVector[index]; 57 } 58 59 inline G4double G4Physics2DVector::GetY(std::size_t index) const 60 { 61 return yVector[index]; 62 } 63 64 inline G4double G4Physics2DVector::GetValue(std::size_t idx, 65 std::size_t idy) const 66 { 67 return (*(value[idy]))[idx]; 68 } 69 70 inline G4double G4Physics2DVector::FindLinearX(G4double rand, G4double y) const 71 { 72 std::size_t idy = 0; 73 return FindLinearX(rand, y, idy); 74 } 75 76 inline std::size_t G4Physics2DVector::GetLengthX() const 77 { 78 return numberOfXNodes; 79 } 80 81 inline std::size_t G4Physics2DVector::GetLengthY() const 82 { 83 return numberOfYNodes; 84 } 85 86 inline G4PhysicsVectorType G4Physics2DVector::GetType() const { return type; } 87 88 inline void G4Physics2DVector::SetBicubicInterpolation(G4bool val) 89 { 90 useBicubic = val; 91 } 92 93 inline std::size_t G4Physics2DVector::FindBin(const G4double z, 94 const G4PV2DDataVector& v, 95 const std::size_t idx, 96 const std::size_t idxmax) const 97 { 98 std::size_t id = idx; 99 if(z <= v[1]) 100 { 101 id = 0; 102 } 103 else if(z >= v[idxmax]) 104 { 105 id = idxmax; 106 } 107 else if(idx > idxmax || z < v[idx] || z > v[idx + 1]) 108 { 109 id = std::lower_bound(v.begin(), v.end(), z) - v.begin() - 1; 110 } 111 return id; 112 } 113 114 inline std::size_t G4Physics2DVector::FindBinLocationX( 115 const G4double x, const std::size_t idx) const 116 { 117 return FindBin(x, xVector, idx, numberOfXNodes - 2); 118 } 119 120 inline std::size_t G4Physics2DVector::FindBinLocationY( 121 const G4double y, const std::size_t idy) const 122 { 123 return FindBin(y, yVector, idy, numberOfYNodes - 2); 124 } 125 126 inline void G4Physics2DVector::SetVerboseLevel(G4int val) 127 { 128 verboseLevel = val; 129 } 130 131 inline G4double G4Physics2DVector::DerivativeX(std::size_t idx, std::size_t idy, 132 G4double fac) const 133 { 134 std::size_t i1 = idx; 135 if(i1 > 0) 136 { 137 --i1; 138 } 139 std::size_t i2 = idx; 140 if(i2 + 1 < numberOfXNodes) 141 { 142 ++i2; 143 } 144 return fac * (GetValue(i2, idy) - GetValue(i1, idy)) / (GetX(i2) - GetX(i1)); 145 } 146 147 inline G4double G4Physics2DVector::DerivativeY(std::size_t idx, std::size_t idy, 148 G4double fac) const 149 { 150 std::size_t i1 = idy; 151 if(i1 > 0) 152 { 153 --i1; 154 } 155 std::size_t i2 = idy; 156 if(i2 + 1 < numberOfYNodes) 157 { 158 ++i2; 159 } 160 return fac * (GetValue(idx, i2) - GetValue(idx, i1)) / (GetY(i2) - GetY(i1)); 161 } 162 163 inline G4double G4Physics2DVector::DerivativeXY(std::size_t idx, 164 std::size_t idy, 165 G4double fac) const 166 { 167 std::size_t i1 = idx; 168 if(i1 > 0) 169 { 170 --i1; 171 } 172 std::size_t i2 = idx; 173 if(i2 + 1 < numberOfXNodes) 174 { 175 ++i2; 176 } 177 std::size_t j1 = idy; 178 if(j1 > 0) 179 { 180 --j1; 181 } 182 std::size_t j2 = idy; 183 if(j2 + 1 < numberOfYNodes) 184 { 185 ++j2; 186 } 187 return fac * 188 (GetValue(i2, j2) - GetValue(i1, j2) - GetValue(i2, j1) + 189 GetValue(i1, j1)) / 190 ((GetX(i2) - GetX(i1)) * (GetY(j2) - GetY(j1))); 191 } 192