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
131
|
/***************************************************************************
data.hh - Data storage
-------------------
begin : Sun Jun 11 2000
copyright : (C) 2000 by Erik Schnetter
email : schnetter@astro.psu.edu
$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/data.hh,v 1.2 2001/03/05 14:31:03 eschnett Exp $
***************************************************************************/
/***************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
* *
***************************************************************************/
#ifndef DATA_HH
#define DATA_HH
#include <cassert>
#include <string>
#include "defs.hh"
#include "dist.hh"
#include "bbox.hh"
#include "bboxset.hh"
#include "gdata.hh"
#include "vect.hh"
// Forward definition
template<class T,int D> class data;
// Output
template<class T,int D>
ostream& operator<< (ostream& os, const data<T,D>& d);
// A real data storage
template<class T,int D>
class data: public generic_data<D> {
// Types
typedef vect<int,D> ivect;
typedef bbox<int,D> ibbox;
typedef bboxset<int,D> ibset;
// Fields
T* restrict _storage; // the data (if located on this processor)
public:
// Constructors
data ();
data (const ibbox& extent, const int proc);
// Destructors
virtual ~data ();
// Pseudo constructors
virtual data* make_typed (const ibbox& extent, const int proc) const;
// Storage management
virtual void allocate (const ibbox& extent, const int proc);
virtual void free ();
virtual void transfer_from (generic_data<D>* gsrc);
// Processor management
virtual void change_processor (const int newproc);
// Accessors
virtual const T* storage () const {
assert (_has_storage);
return _storage;
}
virtual T* storage () {
assert (_has_storage);
return _storage;
}
// Data accessors
const T& operator[] (const ivect& index) const {
assert (_storage);
return _storage[offset(index)];
}
T& operator[] (const ivect& index) {
assert (_storage);
return _storage[offset(index)];
}
// Data manipulators
virtual void copy_from (const generic_data<D>* gsrc, const ibbox& b);
virtual void interpolate_from (const generic_data<D>* gsrc,
const ibbox& box);
virtual void interpolate_from (const generic_data<D>* gsrc,
const double sfact,
const generic_data<D>* gtrc,
const double tfact,
const ibbox& box);
void write_ascii_output_element (ofstream& file, const ivect& index) const;
// void write_ieee (const string name, const int time,
// const int tl, const int rl, const int c, const int ml)
// const;
// void write_hdf (const string name, const int time,
// const int tl, const int rl, const int c, const int ml)
// const;
// void write_h5 (const string name, const int time,
// const int tl, const int rl, const int c, const int ml)
// const;
public:
// Output
ostream& out (ostream& os) const;
};
#if defined(TMPL_IMPLICIT)
# include "data.cc"
#endif
#endif // DATA_HH
|
