VTK  9.2.6
vtkAMRResampleFilter.h
Go to the documentation of this file.
1/*=========================================================================
2
3 Program: Visualization Toolkit
4 Module: vtkAMRResampleFilter.h
5
6 Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
7 All rights reserved.
8 See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
9
10 This software is distributed WITHOUT ANY WARRANTY; without even
11 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
12 PURPOSE. See the above copyright notice for more information.
13
14 =========================================================================*/
36#ifndef vtkAMRResampleFilter_h
37#define vtkAMRResampleFilter_h
38
39#include "vtkFiltersAMRModule.h" // For export macro
41#include <vector> // For STL vector
42
43class vtkInformation;
45class vtkUniformGrid;
49class vtkFieldData;
50class vtkCellData;
51class vtkPointData;
52class vtkIndent;
53
54class vtkAMRBox;
55class VTKFILTERSAMR_EXPORT vtkAMRResampleFilter : public vtkMultiBlockDataSetAlgorithm
56{
57public:
60 void PrintSelf(ostream& oss, vtkIndent indent) override;
61
63
67 vtkSetVector3Macro(NumberOfSamples, int);
68 vtkGetVector3Macro(NumberOfSamples, int);
70
72
75 vtkSetMacro(TransferToNodes, int);
76 vtkGetMacro(TransferToNodes, int);
78
80
84 vtkSetMacro(DemandDrivenMode, int);
85 vtkGetMacro(DemandDrivenMode, int);
87
89
92 vtkSetMacro(NumberOfPartitions, int);
93 vtkGetMacro(NumberOfPartitions, int);
95
97
100 vtkSetVector3Macro(Min, double);
101 vtkGetVector3Macro(Min, double);
103
105
108 vtkSetVector3Macro(Max, double);
109 vtkGetVector3Macro(Max, double);
111
113
116 vtkSetMacro(UseBiasVector, bool);
117 vtkGetMacro(UseBiasVector, bool);
119
121
126 vtkSetVector3Macro(BiasVector, double);
127 vtkGetVector3Macro(BiasVector, double);
129
131
135 vtkGetObjectMacro(Controller, vtkMultiProcessController);
137
138 // Standard pipeline routines
139
145 vtkInformationVector* outputVector) override;
146
148 int FillInputPortInformation(int port, vtkInformation* info) override;
149 int FillOutputPortInformation(int port, vtkInformation* info) override;
150
155
156protected:
159
161 vtkMultiBlockDataSet* ROI; // Pointer to the region of interest.
162 int NumberOfSamples[3];
163 int GridNumberOfSamples[3];
164 double Min[3];
165 double Max[3];
166 double GridMin[3];
167 double GridMax[3];
174 double BiasVector[3];
175
176 // Debugging Stuff
185
186 std::vector<int> BlocksToLoad; // Holds the ids of the blocks to load.
187
192
197 bool IsRegionMine(const int regionIdx);
198
203 int GetRegionProcessId(const int regionIdx);
204
208 void ComputeCellCentroid(vtkUniformGrid* g, const vtkIdType cellIdx, double c[3]);
209
217
222
227 bool FoundDonor(double q[3], vtkUniformGrid*& donorGrid, int& cellIdx);
228
234 bool SearchForDonorGridAtLevel(double q[3], vtkOverlappingAMR* amrds, unsigned int level,
235 unsigned int& gridId, int& donorCellIdx);
236
243 int ProbeGridPointInAMR(double q[3], unsigned int& donorLevel, unsigned int& donorGridId,
244 vtkOverlappingAMR* amrds, unsigned int maxLevel, bool hadDonorGrid);
245
252 int ProbeGridPointInAMRGraph(double q[3], unsigned int& donorLevel, unsigned int& donorGridId,
253 vtkOverlappingAMR* amrds, unsigned int maxLevel, bool useCached);
254
260
266
271
277
282 bool IsBlockWithinBounds(double* grd);
283
290
295 vtkOverlappingAMR* amrds, int N[3], double min[3], double max[3], double h[3]);
296
300 void GetDomainParameters(vtkOverlappingAMR* amr, double domainMin[3], double domainMax[3],
301 double h[3], int dims[3], double& rf);
302
307 double domainMin[3], double domainMax[3], double regionMin[3], double regionMax[3]);
308
313 void AdjustNumberOfSamplesInRegion(const double Rh[3], const bool outside[6], int N[3]);
314
321 const int N[3], const double h0[3], const double L[3], const double rf);
322
329 void SnapBounds(const double h0[3], const double domainMin[3], const double domainMax[3],
330 const int dims[3], bool outside[6]);
331
338
342 void GetRegion(double h[3]);
343
347 bool GridsIntersect(double* g1, double* g2);
348
353
367 void SearchGridDecendants(double q[3], vtkOverlappingAMR* amrds, unsigned int maxLevel,
368 unsigned int& level, unsigned int& gridId, int& id);
369
375 double q[3], vtkOverlappingAMR* amrds, unsigned int& level, unsigned int& gridId, int& id);
376
377private:
379 void operator=(const vtkAMRResampleFilter&) = delete;
380};
381
382#endif /* vtkAMRResampleFilter_h */
Encloses a rectangular region of voxel like cells.
Definition vtkAMRBox.h:40
This filter is a concrete instance of vtkMultiBlockDataSetAlgorithm and provides functionality for ex...
bool IsBlockWithinBounds(double *grd)
Checks if the AMR block, described by a uniform grid, is within the bounds of the ROI perscribed by t...
int RequestData(vtkInformation *, vtkInformationVector **, vtkInformationVector *) override
This is called by the superclass.
int FillOutputPortInformation(int port, vtkInformation *info) override
Fill the output port information objects for this algorithm.
void SearchGridDecendants(double q[3], vtkOverlappingAMR *amrds, unsigned int maxLevel, unsigned int &level, unsigned int &gridId, int &id)
Writes a uniform grid to a file.
void TransferToGridNodes(vtkUniformGrid *g, vtkOverlappingAMR *amrds)
Transfer the solution from the AMR dataset to the nodes of the given uniform grid.
static vtkAMRResampleFilter * New()
void ExtractRegion(vtkOverlappingAMR *amrds, vtkMultiBlockDataSet *mbds, vtkOverlappingAMR *metadata)
Extract the region (as a multiblock) from the given AMR dataset.
vtkOverlappingAMR * AMRMetaData
bool IsParallel()
Checks if this filter instance is running on more than one processes.
bool GridsIntersect(double *g1, double *g2)
Checks if two uniform grids intersect.
void TransferToCellCenters(vtkUniformGrid *g, vtkOverlappingAMR *amrds)
Transfers the solution from the AMR dataset to the cell-centers of the given uniform grid.
int ProbeGridPointInAMR(double q[3], unsigned int &donorLevel, unsigned int &donorGridId, vtkOverlappingAMR *amrds, unsigned int maxLevel, bool hadDonorGrid)
Finds the AMR grid that contains the point q.
bool RegionIntersectsWithAMR(double domainMin[3], double domainMax[3], double regionMin[3], double regionMax[3])
Checks if the domain and requested region intersect.
vtkMultiBlockDataSet * ROI
void CopyData(vtkFieldData *target, vtkIdType targetIdx, vtkCellData *src, vtkIdType srcIdx)
Copies the data to the target from the given source.
vtkUniformGrid * GetReferenceGrid(vtkOverlappingAMR *amrds)
Returns a reference grid from the amrdataset.
void AdjustNumberOfSamplesInRegion(const double Rh[3], const bool outside[6], int N[3])
This method adjust the numbers of samples in the region, N, if the requested region falls outside,...
void ComputeAMRBlocksToLoad(vtkOverlappingAMR *metadata)
Given a user-supplied region of interest and the metadata by a module upstream, this method generates...
void ComputeRegionParameters(vtkOverlappingAMR *amrds, int N[3], double min[3], double max[3], double h[3])
Computes the region parameters.
void PrintSelf(ostream &oss, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
void InitializeFields(vtkFieldData *f, vtkIdType size, vtkCellData *src)
Given the source cell data of an AMR grid, this method initializes the field values,...
void SnapBounds(const double h0[3], const double domainMin[3], const double domainMax[3], const int dims[3], bool outside[6])
This method snaps the bounds s.t.
bool SearchForDonorGridAtLevel(double q[3], vtkOverlappingAMR *amrds, unsigned int level, unsigned int &gridId, int &donorCellIdx)
Given a query point q and a target level, this method finds a suitable grid at the given level that c...
int RequestUpdateExtent(vtkInformation *, vtkInformationVector **, vtkInformationVector *) override
Performs upstream requests to the reader.
void GetDomainParameters(vtkOverlappingAMR *amr, double domainMin[3], double domainMax[3], double h[3], int dims[3], double &rf)
This method accesses the domain boundaries.
virtual void SetController(vtkMultiProcessController *)
Set & Get macro for the multi-process controller.
~vtkAMRResampleFilter() override
void ComputeCellCentroid(vtkUniformGrid *g, const vtkIdType cellIdx, double c[3])
Given a cell index and a grid, this method computes the cell centroid.
void TransferSolution(vtkUniformGrid *g, vtkOverlappingAMR *amrds)
Transfers the solution.
vtkMultiProcessController * Controller
bool IsRegionMine(const int regionIdx)
Given the Region ID this function returns whether or not the region belongs to this process or not.
int GetRegionProcessId(const int regionIdx)
Given the Region ID, this method computes the corresponding process ID that owns the region based on ...
int FillInputPortInformation(int port, vtkInformation *info) override
Fill the input port information objects for this algorithm.
bool FoundDonor(double q[3], vtkUniformGrid *&donorGrid, int &cellIdx)
Given a query point q and a candidate donor grid, this method checks for the corresponding donor cell...
void GetRegion(double h[3])
This method gets the region of interest as perscribed by the user.
int RequestInformation(vtkInformation *rqst, vtkInformationVector **inputVector, vtkInformationVector *outputVector) override
Gets the metadata from upstream module and determines which blocks should be loaded by this instance.
int ProbeGridPointInAMRGraph(double q[3], unsigned int &donorLevel, unsigned int &donorGridId, vtkOverlappingAMR *amrds, unsigned int maxLevel, bool useCached)
Finds the AMR grid that contains the point q.
std::vector< int > BlocksToLoad
void ComputeAndAdjustRegionParameters(vtkOverlappingAMR *amrds, double h[3])
This method computes and adjusts the region parameters s.t.
void ComputeLevelOfResolution(const int N[3], const double h0[3], const double L[3], const double rf)
This method computes the level of resolution based on the number of samples requested,...
bool SearchGridAncestors(double q[3], vtkOverlappingAMR *amrds, unsigned int &level, unsigned int &gridId, int &id)
Find an ancestor of the specified grid that contains the point.
represent and manipulate cell attribute data
Definition vtkCellData.h:42
represent and manipulate fields of data
a simple class to control print indentation
Definition vtkIndent.h:40
Store zero or more vtkInformation instances.
Store vtkAlgorithm input/output information.
Superclass for algorithms that produce only vtkMultiBlockDataSet as output.
Composite dataset that organizes datasets into blocks.
Multiprocessing communication superclass.
hierarchical dataset of vtkUniformGrids
represent and manipulate point attribute data
Computes the portion of a dataset which is inside a selection.
image data with blanking
boost::graph_traits< vtkGraph * >::vertex_descriptor target(boost::graph_traits< vtkGraph * >::edge_descriptor e, vtkGraph *)
int vtkIdType
Definition vtkType.h:332
#define max(a, b)