-
Notifications
You must be signed in to change notification settings - Fork 0
/
fieldA.h
80 lines (62 loc) · 1.6 KB
/
fieldA.h
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
#ifndef FIELDA_H
#define FIELDA_H
using uint32 = std::uint32_t;
class Field3DA
{
private:
double* data_;
double ***data3D_;
uint32 nx_, ny_, nz_;
public:
void allocate()
{
data_ = new double[nx_ * ny_ * nz_];
data3D_ = new double**[nx_];
for (uint32 i=0; i < nx_; ++i)
{
data3D_[i] = new double*[ny_];
if (data3D_[i] == nullptr)
{
throw std::runtime_error("Error");
}
for (uint32 j=0; j < ny_; ++j)
{
data3D_[i][j] = data_ + i * nz_ * ny_ + j * nz_;
}
}
}
Field3DA(uint32 nx, uint32 ny, uint32 nz)
:data_{nullptr}, data3D_{nullptr},
nx_{nx},ny_{ny},nz_{nz}
{
allocate();
}
std::string name() {return "multiple allocations";}
Field3DA(Field3DA const& source)
: data_{nullptr}, data3D_{nullptr}
, nx_{source.nx_}, ny_{source.ny_}, nz_{source.nz_}
{
allocate();
for (auto i=0; i < nx_*ny_*nz_; ++i)
{
data_[i] = source.data_[i];
}
}
std::array<std::size_t,3> shape() const {return std::array<std::size_t,3>{nx_, ny_, nz_};}
inline double& operator()(uint32 i, uint32 j, uint32 k)
{
return data3D_[i][j][k];
}
~Field3DA()
{
for (uint32 i=0; i < nx_; ++i)
{
if (data3D_[i])
delete [] data3D_[i];
else
std::cout << "no need " << i << "\n";
}
delete [] data3D_;
}
};
#endif