Remove alltrans_startup

artis-mcrt · Sep 8, 2024 · 843b6bb · 843b6bb
1 parent 6821e72
commit 843b6bb
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 14 deletions.
diff --git a/atomic.h b/atomic.h
@@ -373,7 +373,8 @@ inline auto get_includedlevels() -> int { return includedlevels; }
 }
 
 [[nodiscard]] inline auto get_uptranslist(const int element, const int ion, const int level) -> LevelTransition * {
-  return globals::alltrans + globals::elements[element].ions[ion].levels[level].alltrans_startup;
+  const auto &levelref = globals::elements[element].ions[ion].levels[level];
+  return globals::alltrans + levelref.alltrans_startdown + levelref.ndowntrans;
 }
 
 [[nodiscard]] inline auto get_uptransspan(const int element, const int ion, const int level) {

diff --git a/globals.h b/globals.h
@@ -81,11 +81,10 @@ struct PhotoionTarget {
 };
 
 struct EnergyLevel {
-  double epsilon{-1};  // Excitation energy of this level relative to the neutral ground level.
-  int nuptrans{0};
-  int alltrans_startup{};  // Allowed upward transitions from this level
+  double epsilon{-1};        // Excitation energy of this level relative to the neutral ground level.
+  int alltrans_startdown{};  // Allowed downward transitions from this level
   int ndowntrans{0};
-  int alltrans_startdown{};        // Allowed downward transitions from this level
+  int nuptrans{0};
   PhotoionTarget *phixstargets{};  // pointer to table of target states and probabilities
   int phixsstart{-1};              // index to start of photoionisation cross-sections table in global::allphixs
   int nphixstargets{0};            // length of phixstargets array:

diff --git a/input.cc b/input.cc
@@ -491,8 +491,6 @@ void add_transitions_to_unsorted_linelist(const int element, const int ion, cons
       for (int level = 0; level < nlevelsmax; level++) {
         globals::elements[element].ions[ion].levels[level].alltrans_startdown = alltransindex;
         alltransindex += get_ndowntrans(element, ion, level);
-
-        globals::elements[element].ions[ion].levels[level].alltrans_startup = alltransindex;
         alltransindex += get_nuptrans(element, ion, level);
 
         set_ndowntrans(element, ion, level, 0);
@@ -559,13 +557,13 @@ void add_transitions_to_unsorted_linelist(const int element, const int ion, cons
                                   .coll_str = transition.coll_str,
                                   .osc_strength = f_ul,
                                   .forbidden = transition.forbidden};
-          temp_alltranslist[globals::elements[element].ions[ion].levels[lowerlevel].alltrans_startup + nloweruptrans -
-                            1] = {.lineindex = -1,
-                                  .targetlevelindex = level,
-                                  .einstein_A = transition.A,
-                                  .coll_str = transition.coll_str,
-                                  .osc_strength = f_ul,
-                                  .forbidden = transition.forbidden};
+          temp_alltranslist[globals::elements[element].ions[ion].levels[lowerlevel].alltrans_startdown +
+                            nupperdowntrans + nloweruptrans - 1] = {.lineindex = -1,
+                                                                    .targetlevelindex = level,
+                                                                    .einstein_A = transition.A,
+                                                                    .coll_str = transition.coll_str,
+                                                                    .osc_strength = f_ul,
+                                                                    .forbidden = transition.forbidden};
         }
 
       } else if (pass == 1 && globals::rank_in_node == 0) {