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 // Author: Ivana Hrivnacova, IJCLab IN2P3/CNRS 28 29 #include <algorithm> 30 31 // 32 // public functions 33 // 34 35 //____________________________________________ 36 template <class T, std::size_t N> 37 G4AccArray<T, N>::G4AccArray( 38 const G4String& name, const T& value, G4Merg 39 : G4VAccumulable(name, mergeMode), 40 fInitValue(value), 41 fMergeFunction(GetMergeFunction<T>(mergeMode 42 { 43 if (G4Accumulables::VerboseLevel > 1 ) { 44 G4cout << "G4AccArray ctor1" << G4endl; 45 } 46 47 for (auto& element : fArray ) { 48 element = fInitValue; 49 } 50 } 51 52 //____________________________________________ 53 template <class T, std::size_t N> 54 template <typename... Args> 55 G4AccArray<T, N>::G4AccArray(Args&&... args) 56 : G4VAccumulable(), 57 fArray{{std::forward<Args>(args)...}}, fMerg 58 { 59 if (G4Accumulables::VerboseLevel > 1 ) { 60 G4cout << "G4AccArray ctor2" << G4endl; 61 } 62 } 63 64 //____________________________________________ 65 template <class T, std::size_t N> 66 template <typename First, typename... Args> 67 G4AccArray<T, N>::G4AccArray(const First& firs 68 : G4VAccumulable(firstArg), 69 fArray{{std::forward<Args>(args)...}}, 70 fMergeFunction(GetMergeFunction<T>(fMergeMod 71 { 72 if (G4Accumulables::VerboseLevel > 1 ) { 73 G4cout << "G4AccArray ctor3" << G4endl; 74 } 75 } 76 77 //____________________________________________ 78 template <class T, std::size_t N> 79 void G4AccArray<T, N>::Merge(const G4VAccumula 80 { 81 const auto& otherArrayAcc = static_cast<cons 82 const auto& otherArray = otherArrayAcc.GetAr 83 84 if (G4Accumulables::VerboseLevel > 2 ) { 85 G4cout << "G4AccArray<T, N>::Merge: " << G 86 G4cout << "destination: "; 87 for (auto v : this->fArray) { 88 G4cout << v << ", "; 89 } 90 G4cout << G4endl; 91 G4cout << "merged data: "; 92 for (auto v : otherArray) { 93 G4cout << v << ", "; 94 } 95 G4cout << G4endl; 96 } 97 98 std::transform(fArray.begin(), fArray.end(), 99 fArray.begin(), fMergeFunctio 100 101 if (G4Accumulables::VerboseLevel > 1 ) { 102 G4cout << "G4AccArray<T, N>::Merge: done" 103 } 104 } 105 106 //____________________________________________ 107 template <class T, std::size_t N> 108 void G4AccArray<T, N>::Reset() 109 { 110 for (auto& value : fArray) { 111 value = fInitValue; 112 } 113 } 114 115 //____________________________________________ 116 template <class T, std::size_t N> 117 void G4AccArray<T, N>::Print(G4PrintOptions op 118 { 119 if (options.Has(G4PrintOptions::kType)) { 120 G4cout << "array<" << typeid(fInitValue).n 121 } 122 123 PrintBase(options); 124 125 bool first = true; 126 for (auto& value : fArray) { 127 if (! first) { G4cout << ", "; } 128 G4cout << value; 129 first = false; 130 } 131 G4cout << G4endl; 132 } 133 134 //____________________________________________ 135 template <class T, std::size_t N> 136 void G4AccArray<T, N>::SetMergeMode(G4MergeMod 137 { 138 G4VAccumulable::SetMergeMode(value); 139 fMergeFunction = GetMergeFunction<T>(value); 140 } 141 142 //____________________________________________ 143 template <class T, std::size_t N> 144 std::array<T,N>& G4AccArray<T, N>::GetArray() 145 { 146 return fArray; 147 } 148 149 //____________________________________________ 150 template <class T, std::size_t N> 151 const std::array<T,N>& G4AccArray<T, N>::GetA 152 { 153 return fArray; 154 } 155