66
77#include < algorithm>
88#include < cmath>
9+ #include < ranges>
910
1011
1112void Nuclide::setSpinParity () const
@@ -233,51 +234,51 @@ void Nuclide::setSeparationEnergies(const std::vector<Nuclide>& nuc) const
233234 // Loop from the penultimate isotope towards the beginning.
234235 // As the vector is ordered by A (low to high), this will
235236 // remove a large number of checks as the vector get bigger.
236- for (auto previous = nuc. crbegin (); previous != nuc. crend (); ++previous )
237+ for (const auto & previous : std::ranges::reverse_view ( nuc) )
237238 {
238239 // Single particle energies.
239- if (A - previous-> A == 1 )
240+ if (A - previous. A == 1 )
240241 {
241242 // S_p(Z,N) = M(Z-1,N) - M(Z,N) + M(1,0)
242- if (N - previous-> N == 0 && Z - previous-> Z == 1 )
243+ if (N - previous. N == 0 && Z - previous. Z == 1 )
243244 {
244- s_p = previous-> NUBASE_ME - NUBASE_ME + p_ME;
245- ds_p = errorQuadrature (previous-> NUBASE_dME , NUBASE_dME, p_dME);
245+ s_p = previous. NUBASE_ME - NUBASE_ME + p_ME;
246+ ds_p = errorQuadrature (previous. NUBASE_dME , NUBASE_dME, p_dME);
246247 numSeparationEnergiesRead++;
247248 }
248249 // S_n(Z,N) = M(Z,N-1) - M(Z,N) + M(0,1)
249- else if (Z - previous-> Z == 0 && N - previous-> N == 1 )
250+ else if (Z - previous. Z == 0 && N - previous. N == 1 )
250251 {
251- s_n = previous-> NUBASE_ME - NUBASE_ME + n_ME;
252- ds_n = errorQuadrature (previous-> NUBASE_dME , NUBASE_dME, n_dME);
252+ s_n = previous. NUBASE_ME - NUBASE_ME + n_ME;
253+ ds_n = errorQuadrature (previous. NUBASE_dME , NUBASE_dME, n_dME);
253254 numSeparationEnergiesRead++;
254255 }
255256 }
256257 // Two particle energies.
257- else if (A - previous-> A == 2 )
258+ else if (A - previous. A == 2 )
258259 {
259260 // S_2p(Z,N) = M(Z-2,N) - M(Z,N) + 2*M(1,0)
260- if (N - previous-> N == 0 && Z - previous-> Z == 2 )
261+ if (N - previous. N == 0 && Z - previous. Z == 2 )
261262 {
262- s_2p = previous-> NUBASE_ME - NUBASE_ME + 2 * p_ME;
263- ds_2p = errorQuadrature (previous-> NUBASE_dME , NUBASE_dME, p_dME, p_dME);
263+ s_2p = previous. NUBASE_ME - NUBASE_ME + 2 * p_ME;
264+ ds_2p = errorQuadrature (previous. NUBASE_dME , NUBASE_dME, p_dME, p_dME);
264265 numSeparationEnergiesRead++;
265266 }
266267 // S_2n(Z,N) = M(Z,N-2) - M(Z,N) + 2*M(0,1)
267- else if (Z - previous-> Z == 0 && N - previous-> N == 2 )
268+ else if (Z - previous. Z == 0 && N - previous. N == 2 )
268269 {
269- s_2n = previous-> NUBASE_ME - NUBASE_ME + 2 * n_ME;
270- ds_2n = errorQuadrature (previous-> NUBASE_dME , NUBASE_dME, n_dME, n_dME);
270+ s_2n = previous. NUBASE_ME - NUBASE_ME + 2 * n_ME;
271+ ds_2n = errorQuadrature (previous. NUBASE_dME , NUBASE_dME, n_dME, n_dME);
271272
272273 // |dV_pn(Z,N)| = 1/4*[S_2p(Z,N) - S_2p(Z,N-2)]
273- dV_pn = fabs (0.25 * (s_2p - previous-> s_2p ));
274- ddV_pn = 0.25 * errorQuadrature (previous-> ds_2p , ds_2p);
274+ dV_pn = fabs (0.25 * (s_2p - previous. s_2p ));
275+ ddV_pn = 0.25 * errorQuadrature (previous. ds_2p , ds_2p);
275276 numSeparationEnergiesRead++;
276277 }
277278 }
278279 // Once the difference in A is greater than 2 we wont get any more useful data
279280 // so set the 'exit variable' to the get out value.
280- else if (A - previous-> A >= 3 )
281+ else if (A - previous. A >= 3 )
281282 {
282283 numSeparationEnergiesRead = 4 ;
283284 }
@@ -296,9 +297,9 @@ void Nuclide::setIsomerData(std::vector<Nuclide>& nuc, const int state) const
296297 // Loop from the penultimate isotope towards the beginning.
297298 // Original order is ground state followed by ascending states,
298299 // theoretically we could just modify nuc.back(), but that's not safe
299- for (auto previous = nuc. rbegin (); previous != nuc. rend (); ++previous )
300+ for (auto & previous : std::ranges::reverse_view ( nuc) )
300301 {
301- if (A == previous-> A && Z == previous-> Z )
302+ if (A == previous. A && Z == previous. Z )
302303 {
303304 double energy{ 0.0 };
304305 double error{ 0.0 };
@@ -313,7 +314,7 @@ void Nuclide::setIsomerData(std::vector<Nuclide>& nuc, const int state) const
313314
314315 setIsomerEnergyError (error);
315316
316- previous-> energy_levels .emplace_back (State ( state, energy, error) );
317+ previous. energy_levels .emplace_back (state, energy, error);
317318 break ;
318319 }
319320 }
0 commit comments