Initial revision

darcs-hash:20010301114010-f6438-12fb8a9ffcc80e86c0a97e37b5b0dae0dbc59b79.gz

1 files changed, 104 insertions, 0 deletions

diff --git a/Carpet/CarpetLib/src/th.hh b/Carpet/CarpetLib/src/th.hh
new file mode 100644
index 000000000..12c77d782
--- /dev/null
+++ b/Carpet/CarpetLib/src/th.hh
@@ -0,0 +1,104 @@
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/th.hh,v 1.11 2004/03/23 19:30:14 schnetter Exp $
+
+#ifndef TH_HH
+#define TH_HH
+
+#include <assert.h>
+
+#include <iostream>
+#include <vector>
+
+#include "cctk.h"
+
+#include "defs.hh"
+#include "gh.hh"
+
+using namespace std;
+
+
+
+// Forward declaration
+template<int D> class th;
+
+// Output
+template<int D>
+ostream& operator<< (ostream& os, const th<D>& t);
+
+
+
+// The time hierarchy (information about the current time)
+template<int D>
+class th {
+  
+public:				// should be readonly
+  
+  // Fields
+  gh<D>& h;                     // hierarchy
+  
+private:
+  
+  CCTK_REAL delta;		// time step
+  vector<vector<CCTK_REAL> > times; // current times
+  vector<vector<CCTK_REAL> > deltas; // time steps
+  
+public:
+  
+  // Constructors
+  th (gh<D>& h, const CCTK_REAL basedelta);
+  
+  // Destructors
+  ~th ();
+  
+  // Modifiers
+  void recompose ();
+  
+  // Time management
+  CCTK_REAL get_time (const int rl, const int ml) const {
+    assert (rl>=0 && rl<h.reflevels());
+    assert (ml>=0 && ml<h.mglevels(rl,0));
+    return times.at(rl).at(ml);
+  }
+  
+  void set_time (const int rl, const int ml, const CCTK_REAL t) {
+    assert (rl>=0 && rl<h.reflevels());
+    assert (ml>=0 && ml<h.mglevels(rl,0));
+    times.at(rl).at(ml) = t;
+  }
+  
+  void advance_time (const int rl, const int ml) {
+    set_time(rl,ml, get_time(rl,ml) + get_delta(rl,ml));
+  }
+  
+  CCTK_REAL get_delta (const int rl, const int ml) const {
+    assert (rl>=0 && rl<h.reflevels());
+    assert (ml>=0 && ml<h.mglevels(rl,0));
+    return deltas.at(rl).at(ml);
+  }
+  
+  void set_delta (const int rl, const int ml, const CCTK_REAL dt) {
+    assert (rl>=0 && rl<h.reflevels());
+    assert (ml>=0 && ml<h.mglevels(rl,0));
+    deltas.at(rl).at(ml) = dt;
+  }
+  
+  CCTK_REAL time (const int tl, const int rl, const int ml) const {
+    assert (rl>=0 && rl<h.reflevels());
+    assert (ml>=0 && ml<h.mglevels(rl,0));
+    return get_time(rl, ml) + tl * get_delta(rl, ml);
+  }
+  
+  // Output
+  void output (ostream& os) const;
+};
+
+
+
+template<int D>
+inline ostream& operator<< (ostream& os, const th<D>& t) {
+  t.output(os);
+  return os;
+}
+
+
+
+#endif // TH_HH