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
|
/*=========================================================================
Program: Visualization Toolkit
Module: vtkAMRUtilities.h
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
// .NAME vtkAMRUtilities -- Support for serial AMR operations
//
// .SECTION Description
// A concrete instance of vtkObject that employs a singleton design
// pattern and implements functionality for AMR specific operations.
//
// .SECTION See Also
// vtkOverlappingAMR, vtkAMRBox
#ifndef VTKAMRUTILITIES_H_
#define VTKAMRUTILITIES_H_
#include "vtkCommonDataModelModule.h" // For export macro
#include "vtkObject.h"
#include <vector> // For C++ vector
// Forward declarations
class vtkFieldData;
class vtkOverlappingAMR;
class vtkUniformGrid;
class VTKCOMMONDATAMODEL_EXPORT vtkAMRUtilities : public vtkObject
{
public:
// Standard Routines
vtkTypeMacro(vtkAMRUtilities,vtkObject);
void PrintSelf(ostream& os, vtkIndent indent );
// Description:
// This method detects and strips partially overlapping cells from a
// given AMR dataset. If ghost layers are detected, they are removed and
// new grid instances are created to represent the stripped
// data-set otherwise, each block is shallow-copied.
//
// .SECTION Assumptions
// 1) The ghosted AMR data must have complete metadata information.
static void StripGhostLayers(
vtkOverlappingAMR *ghostedAMRData,
vtkOverlappingAMR *strippedAMRData);
// Description:
// A quick test of whether partially overlapping ghost cells exist. This test
// starts from the highest-res boxes and checks if they have partially
// overlapping cells. The code returns with true once partially overlapping
// cells are detected. Otherwise, false is returned.
static bool HasPartiallyOverlappingGhostCells(vtkOverlappingAMR *amr);
// Description:
// Blank cells in overlapping AMR
static void BlankCells(vtkOverlappingAMR* amr);
protected:
vtkAMRUtilities() {}
~vtkAMRUtilities() {}
// Description:
// Given the real-extent w.r.t. the ghosted grid, this method copies the
// field data (point/cell) data on the stripped grid.
static void CopyFieldsWithinRealExtent(
int realExtent[6],
vtkUniformGrid *ghostedGrid,
vtkUniformGrid *strippedGrid);
// Description:
// Copies the fields from the given source to the given target.
static void CopyFieldData(
vtkFieldData *target, vtkIdType targetIdx,
vtkFieldData *source, vtkIdType sourceIdx );
// Description:
// Strips ghost layers from the given grid according to the given ghost
// vector which encodes the number of cells to remote from each of the
// 6 sides {imin,imax,jmin,jmax,kmin,kmax}. For example, a ghost vector
// of {0,2,0,2,0,0} would indicate that there exist 2 ghost cells on the
// imax and jmax side.
static vtkUniformGrid* StripGhostLayersFromGrid(
vtkUniformGrid* grid, int ghost[6]);
static void BlankGridsAtLevel(vtkOverlappingAMR* amr, int levelIdx,
std::vector<std::vector<unsigned int> >& children,
const std::vector<int>& processMap);
private:
vtkAMRUtilities(const vtkAMRUtilities&); // Not implemented
void operator=(const vtkAMRUtilities&); // Not implemented
};
#endif /* VTKAMRUTILITIES_H_ */
|
