88 #ifndef vtkLagrangianParticleTracker_h
89 #define vtkLagrangianParticleTracker_h
93 #include "vtkFiltersFlowPathsModule.h"
114 struct IntegratingFunctor;
125 STEP_LAST_CELL_LENGTH = 0,
126 STEP_CUR_CELL_LENGTH = 1,
127 STEP_LAST_CELL_VEL_DIR = 2,
128 STEP_CUR_CELL_VEL_DIR = 3,
129 STEP_LAST_CELL_DIV_THEO = 4,
130 STEP_CUR_CELL_DIV_THEO = 5
131 } CellLengthComputation;
157 vtkSetMacro(GeneratePolyVertexInteractionOutput,
bool);
158 vtkGetMacro(GeneratePolyVertexInteractionOutput,
bool);
185 vtkSetMacro(CellLengthComputationMode,
int);
186 vtkGetMacro(CellLengthComputationMode,
int);
193 vtkSetMacro(StepFactor,
double);
194 vtkGetMacro(StepFactor,
double);
201 vtkSetMacro(StepFactorMin,
double);
202 vtkGetMacro(StepFactorMin,
double);
209 vtkSetMacro(StepFactorMax,
double);
210 vtkGetMacro(StepFactorMax,
double);
217 vtkSetMacro(MaximumNumberOfSteps,
int);
218 vtkGetMacro(MaximumNumberOfSteps,
int);
226 vtkSetMacro(MaximumIntegrationTime,
double);
227 vtkGetMacro(MaximumIntegrationTime,
double);
237 vtkSetMacro(AdaptiveStepReintegration,
bool);
238 vtkGetMacro(AdaptiveStepReintegration,
bool);
239 vtkBooleanMacro(AdaptiveStepReintegration,
bool);
247 vtkSetMacro(GenerateParticlePathsOutput,
bool);
248 vtkGetMacro(GenerateParticlePathsOutput,
bool);
249 vtkBooleanMacro(GenerateParticlePathsOutput,
bool);
323 std::queue<vtkLagrangianParticle*>& particles,
vtkPointData* seedData);
326 int nVar, std::queue<vtkLagrangianParticle*>& particles);
353 std::queue<vtkLagrangianParticle*>&,
vtkPolyData* particlePathsOutput,
360 vtkIdList* particlePathPointId,
bool prev =
false);
366 unsigned int interactedSurfaceFlatIndex,
vtkDataObject* interactionOutput);
374 double t,
double& delT,
double& delTActual,
double minStep,
double maxStep,
double cellLength,
387 bool GenerateParticlePathsOutput =
true;
406 friend struct IntegratingFunctor;
Proxy object to connect input/output ports.
Fast, simple class for dealing with 3D bounds.
object to represent cell connectivity
abstract superclass for arrays of numeric data
Superclass for algorithms that produce only data object as output.
general representation of visualization data
abstract class to specify dataset behavior
dynamic, self-adjusting array of double
list of point or cell ids
a simple class to control print indentation
Integrate a set of ordinary differential equations (initial value problem) in time.
vtkFunctionSet abstract implementation to be used in the vtkLagrangianParticleTracker integrator.
Filter to inject and track particles in a flow.
vtkDataObject * GetSurface()
vtkIdType IntegratedParticleCounterIncrement
double MinimumVelocityMagnitude
void SetSourceConnection(vtkAlgorithmOutput *algOutput)
Specify the source object used to generate particle initial position (seeds).
vtkLagrangianBasicIntegrationModel * IntegrationModel
virtual bool UpdateSurfaceCacheIfNeeded(vtkDataObject *&surfaces)
vtkDataObject * SurfacesCache
virtual vtkIdType GetNewParticleId()
Get an unique id for a particle This method is thread safe.
double MaximumIntegrationTime
void InsertPathOutputPoint(vtkLagrangianParticle *particle, vtkPolyData *particlePathsOutput, vtkIdList *particlePathPointId, bool prev=false)
This method is thread safe.
vtkInitialValueProblemSolver * Integrator
int FillInputPortInformation(int port, vtkInformation *info) override
Declare input port type.
virtual void GetParticleFeed(std::queue< vtkLagrangianParticle * > &particleQueue)
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
bool GeneratePolyVertexInteractionOutput
int RequestData(vtkInformation *request, vtkInformationVector **inputVector, vtkInformationVector *outputVector) override
Process inputs to integrate particle and generate output.
vtkBoundingBox FlowBoundsCache
vtkMTimeType SurfacesTime
bool ComputeNextStep(vtkInitialValueProblemSolver *integrator, double *xprev, double *xnext, double t, double &delT, double &delTActual, double minStep, double maxStep, double cellLength, int &integrationRes, vtkLagrangianParticle *particle)
This method is thread safe.
vtkLagrangianParticleTracker()
virtual void InitializeSurface(vtkDataObject *&surfaces)
virtual void GenerateParticles(const vtkBoundingBox *bounds, vtkDataSet *seeds, vtkDataArray *initialVelocities, vtkDataArray *initialIntegrationTimes, vtkPointData *seedData, int nVar, std::queue< vtkLagrangianParticle * > &particles)
static vtkLagrangianParticleTracker * New()
virtual bool InitializeParticles(const vtkBoundingBox *bounds, vtkDataSet *seeds, std::queue< vtkLagrangianParticle * > &particles, vtkPointData *seedData)
static void InsertPolyVertexCell(vtkPolyData *polydata)
double ComputeCellLength(vtkLagrangianParticle *particle)
void SetIntegrator(vtkInitialValueProblemSolver *integrator)
Set/Get the integrator.
void SetSurfaceData(vtkDataObject *source)
Specify the source object used to compute surface interaction with Note that this method does not con...
static void InsertVertexCells(vtkPolyData *polydata)
vtkMTimeType GetMTime() override
Get the tracker modified time taking into account the integration model and the integrator.
std::atomic< vtkIdType > IntegratedParticleCounter
void SetSurfaceConnection(vtkAlgorithmOutput *algOutput)
Specify the object used to compute surface interaction with.
void SetIntegrationModel(vtkLagrangianBasicIntegrationModel *integrationModel)
Set/Get the integration model.
virtual int Integrate(vtkInitialValueProblemSolver *integrator, vtkLagrangianParticle *, std::queue< vtkLagrangianParticle * > &, vtkPolyData *particlePathsOutput, vtkPolyLine *particlePath, vtkDataObject *interactionOutput)
This method is thread safe.
std::atomic< vtkIdType > ParticleCounter
virtual bool InitializeFlow(vtkDataObject *flow, vtkBoundingBox *bounds)
vtkDataObject * FlowCache
int FillOutputPortInformation(int port, vtkInformation *info) override
Declare output port type.
bool AdaptiveStepReintegration
int RequestDataObject(vtkInformation *, vtkInformationVector **, vtkInformationVector *) override
Create outputs objects.
virtual bool InitializePathsOutput(vtkPointData *seedData, vtkIdType numberOfSeeds, vtkPolyData *&particlePathsOutput)
This method is thread safe.
virtual bool InitializeInteractionOutput(vtkPointData *seedData, vtkDataObject *surfaces, vtkDataObject *&interractionOutput)
This method is thread safe.
vtkDataObject * GetSource()
double MinimumReductionFactor
void SetSourceData(vtkDataObject *source)
Specify the source object used to generate particle initial position (seeds).
void InsertInteractionOutputPoint(vtkLagrangianParticle *particle, unsigned int interactedSurfaceFlatIndex, vtkDataObject *interactionOutput)
This method is thread safe.
int CellLengthComputationMode
virtual bool FinalizeOutputs(vtkPolyData *particlePathsOutput, vtkDataObject *interactionOutput)
~vtkLagrangianParticleTracker() override
Basis class for Lagrangian particles.
Composite dataset that organizes datasets into blocks.
composite dataset to encapsulates pieces of dataset.
represent and manipulate point attribute data
represent and manipulate 3D points
concrete dataset represents vertices, lines, polygons, and triangle strips
cell represents a set of 1D lines
boost::graph_traits< vtkGraph * >::vertex_descriptor source(boost::graph_traits< vtkGraph * >::edge_descriptor e, vtkGraph *)
vtkTypeUInt32 vtkMTimeType