1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
|
#ifndef TH_HH
#define TH_HH
#include <cctk.h>
#include <cassert>
#include <cmath>
#include <iostream>
#include <vector>
#include "defs.hh"
#include "gh.hh"
using namespace std;
// Forward declaration
class th;
// Input
istream& operator>> (istream& is, th& t);
// Output
ostream& operator<< (ostream& os, const th& t);
// The time hierarchy (information about the current time)
class th {
static list<th*> allth;
list<th*>::iterator allthi;
public: // should be readonly
// Fields
gh& h; // hierarchy
gh::th_handle gh_handle;
bool const time_interpolation_during_regridding;
int timelevels; // const
private:
vector<int> reffacts; // const
vector<vector<vector<CCTK_REAL> > > times; // current times [ml][rl][tl]
vector<vector<CCTK_REAL> > deltas; // time steps [ml][rl]
public:
// Constructors
th (gh& h, bool time_interpolation_during_regridding);
th (gh& h, int timelevels, vector<int> const& reffacts,
bool time_interpolation_during_regridding);
// Destructors
~th ();
// Modifiers
void regrid ();
void regrid_free ();
// Time management
void set_time (int const ml, int const rl, int const tl, CCTK_REAL const& t)
{
assert (ml>=0 and ml<h.mglevels());
assert (rl>=0 and rl<h.reflevels());
assert (tl>=0 and tl<timelevels);
// assert (isfinite(t));
times.AT(ml).AT(rl).AT(tl) = t;
}
CCTK_REAL get_time (int const ml, int const rl, int const tl)
const
{
assert (ml>=0 and ml<h.mglevels());
assert (rl>=0 and rl<h.reflevels());
assert (tl>=0 and tl<timelevels);
CCTK_REAL const t = times.AT(ml).AT(rl).AT(tl);
// assert (isfinite(t));
return t;
}
void set_delta (int const ml, int const rl, CCTK_REAL const& dt)
{
assert (ml>=0 and ml<h.mglevels());
assert (rl>=0 and rl<h.reflevels());
// assert (isfinite(dt));
deltas.AT(ml).AT(rl) = dt;
}
CCTK_REAL get_delta (int const ml, int const rl) const
{
assert (ml>=0 and ml<h.mglevels());
assert (rl>=0 and rl<h.reflevels());
CCTK_REAL const dt = deltas.AT(ml).AT(rl);
// assert (isfinite(dt));
return dt;
}
void advance_time (int const ml, int const rl);
void retreat_time (int const ml, int const rl);
void flip_timelevels (int const ml, int const rl);
// Output
size_t memory () const CCTK_MEMBER_ATTRIBUTE_PURE;
static size_t allmemory () CCTK_MEMBER_ATTRIBUTE_PURE;
istream& input (istream& is);
ostream& output (ostream& os) const;
};
inline size_t memoryof (th const & t)
{
return t.memory ();
}
inline istream& operator>> (istream& is, th& t) {
return t.input(is);
}
inline ostream& operator<< (ostream& os, const th& t) {
return t.output(os);
}
#endif // TH_HH
|
