1 files changed, 91 insertions, 0 deletions
diff --git a/Carpet/CarpetLib/src/gf.cc b/Carpet/CarpetLib/src/gf.cc
new file mode 100644
index 000000000..69537d537
--- /dev/null
+++ b/Carpet/CarpetLib/src/gf.cc
@@ -0,0 +1,91 @@
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/gf.cc,v 1.21 2004/08/07 19:47:11 schnetter Exp $
+
+#include <assert.h>
+
+#include "cctk.h"
+
+#include "defs.hh"
+
+#include "gf.hh"
+
+using namespace std;
+
+
+
+// Constructors
+template<class T,int D>
+gf<T,D>::gf (const int varindex, const operator_type transport_operator,
+             th<D>& t, dh<D>& d,
+	     const int tmin, const int tmax, const int prolongation_order_time,
+             const int vectorlength, const int vectorindex,
+             gf* const vectorleader)
+  : ggf<D>(varindex, transport_operator,
+           t, d, tmin, tmax, prolongation_order_time,
+           vectorlength, vectorindex, vectorleader)
+{
+  // this->recompose ();
+  this->recompose_crop ();
+  for (int rl=0; rl<this->h.reflevels(); ++rl) {
+    this->recompose_allocate (rl);
+#if 0
+    for (comm_state<D> state; !state.done(); state.step()) {
+      this->recompose_fill (state, rl, false);
+    }
+#endif
+    this->recompose_free (rl);
+#if 0
+    for (comm_state<D> state; !state.done(); state.step()) {
+      this->recompose_bnd_prolongate (state, rl, false);
+    }
+    for (comm_state<D> state; !state.done(); state.step()) {
+      this->recompose_sync (state, rl, false);
+    }
+#endif
+  } // for rl
+}
+
+// Destructors
+template<class T,int D>
+gf<T,D>::~gf () { }
+
+
+
+// Access to the data
+template<class T,int D>
+const data<T,D>* gf<T,D>::operator() (int tl, int rl, int c, int ml) const {
+  assert (tl>=this->tmin && tl<=this->tmax);
+  assert (rl>=0 && rl<this->h.reflevels());
+  assert (c>=0 && c<this->h.components(rl));
+  assert (ml>=0 && ml<this->h.mglevels(rl,c));
+  return (const data<T,D>*)this->storage.at(tl-this->tmin).at(rl).at(c).at(ml);
+}
+
+template<class T,int D>
+data<T,D>* gf<T,D>::operator() (int tl, int rl, int c, int ml) {
+  assert (tl>=this->tmin && tl<=this->tmax);
+  assert (rl>=0 && rl<this->h.reflevels());
+  assert (c>=0 && c<this->h.components(rl));
+  assert (ml>=0 && ml<this->h.mglevels(rl,c));
+  return (data<T,D>*)this->storage.at(tl-this->tmin).at(rl).at(c).at(ml);
+}
+
+
+
+// Output
+template<class T,int D>
+ostream& gf<T,D>::output (ostream& os) const {
+  T Tdummy;
+  os << "gf<" << typestring(Tdummy) << "," << D << ">:"
+     << this->varindex << "[" << CCTK_VarName(this->varindex) << "],"
+     << "dt=[" << this->tmin << ":" << this->tmax<< "]";
+  return os;
+}
+
+
+
+#define INSTANTIATE(T)				\
+template class gf<T,3>;
+
+#include "instantiate"
+
+#undef INSTANTIATE