aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetLib/src/th.cc
diff options
context:
space:
mode:
Diffstat (limited to 'Carpet/CarpetLib/src/th.cc')
-rw-r--r--Carpet/CarpetLib/src/th.cc28
1 files changed, 14 insertions, 14 deletions
diff --git a/Carpet/CarpetLib/src/th.cc b/Carpet/CarpetLib/src/th.cc
index d849e3c47..5fd8eabee 100644
--- a/Carpet/CarpetLib/src/th.cc
+++ b/Carpet/CarpetLib/src/th.cc
@@ -1,4 +1,4 @@
-// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/th.cc,v 1.12 2003/08/10 21:58:45 schnetter Exp $
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/th.cc,v 1.13 2004/01/25 14:57:30 schnetter Exp $
#include <assert.h>
#include <math.h>
@@ -18,23 +18,23 @@ using namespace std;
// Constructors
template<int D>
-th<D>::th (gh<D>* h, const CCTK_REAL basedelta)
+th<D>::th (gh<D>& h, const CCTK_REAL basedelta)
: h(h), delta(basedelta) {
- h->add(this);
+ h.add(this);
}
// Destructors
template<int D>
th<D>::~th () {
- h->remove(this);
+ h.remove(this);
}
// Modifiers
template<int D>
void th<D>::recompose () {
- times.resize(h->reflevels());
- deltas.resize(h->reflevels());
- for (int rl=0; rl<h->reflevels(); ++rl) {
+ times.resize(h.reflevels());
+ deltas.resize(h.reflevels());
+ for (int rl=0; rl<h.reflevels(); ++rl) {
const int old_mglevels = times[rl].size();
CCTK_REAL mgtime;
// Select default time
@@ -45,15 +45,15 @@ void th<D>::recompose () {
} else {
mgtime = times[rl][old_mglevels-1];
}
- times[rl].resize(h->mglevels(rl,0), mgtime);
- deltas[rl].resize(h->mglevels(rl,0));
- for (int ml=0; ml<h->mglevels(rl,0); ++ml) {
+ times[rl].resize(h.mglevels(rl,0), mgtime);
+ deltas[rl].resize(h.mglevels(rl,0));
+ for (int ml=0; ml<h.mglevels(rl,0); ++ml) {
if (rl==0 && ml==0) {
deltas[rl][ml] = delta;
} else if (ml==0) {
- deltas[rl][ml] = deltas[rl-1][ml] / h->reffact;
+ deltas[rl][ml] = deltas[rl-1][ml] / h.reffact;
} else {
- deltas[rl][ml] = deltas[rl][ml-1] * h->mgfact;
+ deltas[rl][ml] = deltas[rl][ml-1] * h.mgfact;
}
}
}
@@ -66,8 +66,8 @@ template<int D>
void th<D>::output (ostream& os) const {
os << "th<" << D << ">:"
<< "times={";
- for (int rl=0; rl<h->reflevels(); ++rl) {
- for (int ml=0; ml<h->mglevels(rl,0); ++ml) {
+ for (int rl=0; rl<h.reflevels(); ++rl) {
+ for (int ml=0; ml<h.mglevels(rl,0); ++ml) {
if (!(rl==0 && ml==0)) os << ",";
os << rl << ":" << ml << ":"
<< times[rl][ml] << "(" << deltas[rl][ml] << ")";