-
Notifications
You must be signed in to change notification settings - Fork 0
/
CartesianDataBlock.h
116 lines (96 loc) · 2.38 KB
/
CartesianDataBlock.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
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
/*
____ _ __ ____ __ ____
/ __/___(_) / ___ ____/ __ \__ _____ ___ / /_ / _/__ ____
_\ \/ __/ / _ \/ -_) __/ /_/ / // / -_|_-</ __/ _/ // _ \/ __/
/___/\__/_/_.__/\__/_/ \___\_\_,_/\__/___/\__/ /___/_//_/\__(_)
Copyright 2012 SciberQuest Inc.
*/
#ifndef __CartesianDataBlock_h
#define __CartesianDataBlock_h
#include <iostream>
using std::ostream;
#include "CartesianExtent.h"
#include "CartesianBounds.h"
class vtkDataSet;
/// Data and Meta-data describing data on cartesian grid.
class CartesianDataBlock
{
public:
CartesianDataBlock();
~CartesianDataBlock();
/**
Set the block id, it's coordinates in the cartesian
decomposition.
*/
void SetId(int i, int j, int k, int idx);
void SetId(int *I);
int *GetId(){ return this->Id; }
/**
Get the block's flat index.
*/
int GetIndex(){ return this->Id[3]; }
/**
Set the physical area enclosed by the Block.
*/
void SetBounds(double *bounds);
void SetBounds(
double xlo,
double xhi,
double ylo,
double yhi,
double zlo,
double zhi);
/**
Direct acces to internal bounds data.
*/
CartesianBounds &GetBounds(){ return this->Bounds; }
/**
Return non-zero if the point is contained within the block.
*/
int Inside(double *pt){ return this->Bounds.Inside(pt); }
/**
Set the block's extent.
*/
void SetExtent(int *ext);
void SetExtent(
int ilo,
int ihi,
int jlo,
int jhi,
int klo,
int khi);
/**
Direct access to internal extent data.
*/
CartesianExtent &GetExtent(){ return this->Extent; }
/**
Return the lower left corner of the block in world space.
*/
void GetBlockOrigin(double *x0)
{
x0[0]=this->Bounds[0];
x0[1]=this->Bounds[2];
x0[2]=this->Bounds[4];
}
/**
Set and get the data associated with this block.
*/
void SetData(vtkDataSet *data);
vtkDataSet *GetData();
/**
Return the decomp index of the neighboring
block containing the given point.
*/
// void GetNeighborId(double *pt, int *I);
private:
CartesianDataBlock(CartesianDataBlock &other);
void operator=(CartesianDataBlock &other);
friend ostream &operator<<(ostream &os, CartesianDataBlock &b);
private:
int Id[4];
CartesianExtent Extent;
CartesianBounds Bounds;
vtkDataSet *Data;
};
ostream &operator<<(ostream &os, CartesianDataBlock &b);
#endif