Dynamic topology phase IV

packages/seacas/libraries/ioss/src/Ioss_DatabaseIO.h

@@ -303,6 +303,8 @@ namespace Ioss {
       return groups_describe_nl(return_full_names);
     }
 
+    IOSS_NODISCARD virtual std::string get_group_name() const { return "/"; }
+
     /** \brief Set the database to the given State.
      *
      *  All transitions must begin from the 'STATE_CLOSED' state or be to
@@ -356,7 +358,13 @@ namespace Ioss {
     void get_step_times()
     {
       IOSS_FUNC_ENTER(m_);
-      return get_step_times_nl();
+      get_step_times_nl();
+    }
+
+    std::vector<double> get_db_step_times()
+    {
+      IOSS_FUNC_ENTER(m_);
+      return get_db_step_times_nl();
     }
 
     IOSS_NODISCARD virtual bool internal_edges_available() const { return false; }
@@ -805,6 +813,7 @@ namespace Ioss {
 
     virtual void read_meta_data_nl() = 0;
     virtual void get_step_times_nl() {}
+    virtual std::vector<double> get_db_step_times_nl() { return std::vector<double>(); }
 
     virtual bool begin_state_nl(int state, double time);
     virtual bool end_state_nl(int state, double time);

packages/seacas/libraries/ioss/src/Ioss_Region.C

@@ -308,6 +308,56 @@ namespace {
       entity->field_erase(role);
     }
   }
+
+  struct NearestGroupState
+  {
+    std::string group{"/"};
+    int nearestState{-1};
+    double nearestTime{-std::numeric_limits<double>::max()};
+  };
+
+  void get_state_time_round_down(Ioss::DatabaseIO *db, const double targetTime, NearestGroupState& loc)
+  {
+    std::vector<double> timesteps = db->get_db_step_times();
+    int stepCount = timesteps.size();
+
+    double minTimeDiff = loc.nearestState < 0 ? std::numeric_limits<double>::max() : std::abs(loc.nearestTime - targetTime);
+
+    for(int istep = 1; istep <= stepCount; istep++) {
+      double stateTime = timesteps[istep-1];
+      double stepTimeDiff = std::abs(stateTime - targetTime);
+      if(stepTimeDiff <= minTimeDiff) {
+        minTimeDiff = stepTimeDiff;
+        loc.nearestTime = stateTime;
+        loc.nearestState = istep;
+        loc.group = db->get_group_name();
+      }
+    }
+  }
+
+  void locate_state(Ioss::DatabaseIO *db, const double targetTime, NearestGroupState& loc)
+  {
+    // Get state count and all states...
+    std::vector<double> timesteps = db->get_db_step_times();
+    int stepCount = timesteps.size();
+
+    if(targetTime < 0.0) {
+      get_state_time_round_down(db, targetTime, loc);
+    }
+    else {
+      double minTimeDiff = loc.nearestState < 0 ? std::numeric_limits<double>::max() : std::fabs(targetTime - loc.nearestTime);
+      for(int istep = 1; istep <= stepCount; istep++) {
+        double stateTime = timesteps[istep-1];
+        double stepTimeDiff = std::fabs(targetTime - stateTime);
+        if(stepTimeDiff < minTimeDiff) {
+          minTimeDiff = stepTimeDiff;
+          loc.nearestTime = stateTime;
+          loc.nearestState = istep;
+          loc.group = db->get_group_name();
+        }
+      }
+    }
+  }
 } // namespace
 
 namespace Ioss {
@@ -3055,4 +3105,28 @@ namespace Ioss {
       topologyObserver->reset_topology_modification();
     }
   }
+
+  std::string Region::get_group_name() const
+  {
+    return get_database()->get_group_name();
+  }
+
+  std::tuple<std::string, int, double> Region::locate_db_state(const double targetTime)
+  {
+    NearestGroupState loc;
+
+    auto db = get_database();
+    std::string currentGroup = db->get_group_name();
+
+    for(int i=0; i<db->num_child_group(); i++) {
+      db->open_root_group();
+      db->open_child_group(i);
+      locate_state(db, targetTime, loc);
+    }
+
+    db->open_root_group();
+    db->open_group(currentGroup);
+
+    return std::make_tuple(loc.group, loc.nearestState, loc.nearestTime);
+  }
 } // namespace Ioss

packages/seacas/libraries/ioss/src/Ioss_Region.h

@@ -33,6 +33,7 @@
 #include <memory>
 #include <sstream>
 #include <string>  // for string, operator<
+#include <tuple>
 #include <utility> // for pair
 #include <vector>  // for vector
 
@@ -313,6 +314,10 @@ namespace Ioss {
     bool load_group_mesh(const std::string &child_group_name);
     bool load_group_mesh(const int child_group_index);
 
+    std::string get_group_name() const;
+
+    IOSS_NODISCARD std::tuple<std::string, int, double> locate_db_state(const double targetTime);
+
   protected:
     void update_dynamic_topology();
     void clone_and_replace_output_database(int steps = 0);

packages/seacas/libraries/ioss/src/catalyst/Iocatalyst_DatabaseIO.C

@@ -357,6 +357,16 @@ namespace Iocatalyst {
       return true;
     }
 
+    std::vector<double> getTime()
+    {
+      std::vector<double> times;
+      auto &node = this->DBNode;
+      if (node.has_path(getTimePath())) {
+        times.push_back(node[getTimePath()].as_float64());
+      }
+      return times;
+    }
+
     int64_t putField(const std::string &containerName, const Ioss::GroupingEntity *entityGroup,
                      const Ioss::Field &field, void *data, size_t data_size, bool deep_copy)
     {
@@ -1426,6 +1436,12 @@ namespace Iocatalyst {
     impl.readTime(region);
   }
 
+  std::vector<double> DatabaseIO::get_db_step_times_nl()
+  {
+    auto &impl = (*this->Impl.get());
+    return impl.getTime();
+  }
+
   void *DatabaseIO::get_catalyst_conduit_node()
   {
     auto &impl = (*this->Impl.get());

packages/seacas/libraries/ioss/src/catalyst/Iocatalyst_DatabaseIO.h

@@ -97,6 +97,7 @@ namespace Iocatalyst {
 
     void read_meta_data_nl() override;
     void get_step_times_nl() override;
+    std::vector<double> get_db_step_times_nl() override;
 
     bool begin_state_nl(int state, double time) override;
     bool end_state_nl(int state, double time) override;

packages/seacas/libraries/ioss/src/cgns/Iocgns_DatabaseIO.C

@@ -1640,6 +1640,16 @@ namespace Iocgns {
                           myProcessor);
   }
 
+  std::vector<double> DatabaseIO::get_db_step_times_nl()
+  {
+    std::vector<double> timesteps;
+
+    Utils::get_step_times(get_file_pointer(), timesteps, nullptr, timeScaleFactor,
+                          myProcessor);
+
+    return timesteps;
+  }
+
   void DatabaseIO::write_adjacency_data()
   {
     // Determine adjacency information between unstructured blocks.

packages/seacas/libraries/ioss/src/cgns/Iocgns_DatabaseIO.h

@@ -108,6 +108,7 @@ namespace Iocgns {
     size_t finalize_structured_blocks();
     void   finalize_database() const override;
     void   get_step_times_nl() override;
+    std::vector<double> get_db_step_times_nl() override;
 
     void create_unstructured_block(int base, int zone, size_t &num_node);
     void write_adjacency_data();

packages/seacas/libraries/ioss/src/cgns/Iocgns_ParallelDatabaseIO.C

@@ -806,6 +806,16 @@ namespace Iocgns {
                           myProcessor);
   }
 
+  std::vector<double> ParallelDatabaseIO::get_db_step_times_nl()
+  {
+    std::vector<double> timesteps;
+
+    Utils::get_step_times(get_file_pointer(), timesteps, nullptr, timeScaleFactor,
+                          myProcessor);
+
+    return timesteps;
+  }
+
   void ParallelDatabaseIO::write_adjacency_data()
   {
     // Determine adjacency information between unstructured blocks.

packages/seacas/libraries/ioss/src/cgns/Iocgns_ParallelDatabaseIO.h

@@ -113,6 +113,7 @@ namespace Iocgns {
     int64_t handle_node_ids(void *ids, int64_t num_to_get) const;
     void    finalize_database() const override;
     void    get_step_times_nl() override;
+    std::vector<double> get_db_step_times_nl() override;
     void    write_adjacency_data();
 
     int64_t get_field_internal(const Ioss::Region *reg, const Ioss::Field &field, void *data,

packages/seacas/libraries/ioss/src/cgns/Iocgns_Utils.C

@@ -2396,7 +2396,9 @@ int Iocgns::Utils::get_step_times(int cgns_file_ptr, std::vector<double> &timest
 
   timesteps.reserve(num_timesteps);
   for (int i = 0; i < num_timesteps; i++) {
-    region->add_state(times[i] * timeScaleFactor);
+    if(nullptr != region) {
+      region->add_state(times[i] * timeScaleFactor);
+    }
     timesteps.push_back(times[i]);
   }
   return num_timesteps;

packages/seacas/libraries/ioss/src/exodus/Ioex_BaseDatabaseIO.h

@@ -90,6 +90,8 @@ namespace Ioex {
     // database supports that type (e.g. return_value & Ioss::FACESET)
     IOSS_NODISCARD unsigned entity_field_support() const override;
 
+    IOSS_NODISCARD std::string get_group_name() const override { return m_groupName; }
+
   protected:
     // Check to see if database state is ok...
     // If 'write_message' true, then output a warning message indicating the problem.