CarpetLib: use map instead of set for dh::gfs

some routines assume that elements in gfs are sorted by variable index
when traversing the container (eg recompose_allocate does)

1 files changed, 7 insertions, 1 deletions

diff --git a/Carpet/CarpetLib/src/dh.hh b/Carpet/CarpetLib/src/dh.hh
index b558b013f..55b75218c 100644
--- a/Carpet/CarpetLib/src/dh.hh
+++ b/Carpet/CarpetLib/src/dh.hh
@@ -5,6 +5,7 @@
 #include <cstdlib>
 #include <iostream>
 #include <set>
+#include <map>
 #include <string>
 #include <vector>
 
@@ -248,7 +249,12 @@ public:                         // should be readonly
   level_mboxes level_boxes;     // grid hierarchy [ml][rl]
   fast_mboxes fast_boxes;       // grid hierarchy [ml][rl][p]
   
-  set<ggf*> gfs;                // all grid functions
+private:
+
+  // this needs to be sorted by varindex so that when iteration through the
+  // container in order with a forward iterator, vector leaders are processed
+  // first
+  map<int,ggf*> gfs;            // all grid functions
   
 public: