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
|
/*=========================================================================
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.
=========================================================================*/
/**
* @class vtkAMRUtilities
*
*
* A concrete instance of vtkObject that employs a singleton design
* pattern and implements functionality for AMR specific operations.
*
* @sa
* vtkOverlappingAMR, vtkAMRBox
*/
#ifndef vtkAMRUtilities_h
#define vtkAMRUtilities_h
#include "vtkCommonDataModelModule.h" // For export macro
#include "vtkObject.h"
#include <vector> // For C++ vector
// Forward declarations
VTK_ABI_NAMESPACE_BEGIN
class vtkDataArray;
class vtkFieldData;
class vtkOverlappingAMR;
class vtkUniformGrid;
class vtkUnsignedCharArray;
class VTKCOMMONDATAMODEL_EXPORT vtkAMRUtilities : public vtkObject
{
public:
// Standard Routines
vtkTypeMacro(vtkAMRUtilities, vtkObject);
void PrintSelf(ostream& os, vtkIndent indent) override;
/**
* 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);
/**
* 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);
/**
* Blank cells in overlapping AMR
*/
static void BlankCells(vtkOverlappingAMR* amr);
protected:
vtkAMRUtilities() = default;
~vtkAMRUtilities() override = default;
/**
* 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);
/**
* Copies the fields from the given source to the given target.
*/
static void CopyFieldData(
vtkFieldData* target, vtkIdType targetIdx, vtkFieldData* source, vtkIdType sourceIdx);
/**
* 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&) = delete;
void operator=(const vtkAMRUtilities&) = delete;
static void MergeGhostArrays(vtkDataArray* existingArray, vtkUnsignedCharArray* ghosts);
};
VTK_ABI_NAMESPACE_END
#endif /* vtkAMRUtilities_h */
|
