Updates will be applied April 15th at 12pm EDT (UTC-0400). GitLab could be a little slow between 12 - 12:45pm EDT.

Commit 58e8bc6d authored by Mathieu Westphal's avatar Mathieu Westphal

Renable Divergence Theorem for vtkVoxel cells.

parent fce6c638
......@@ -1406,7 +1406,7 @@ double vtkLagrangianParticleTracker::ComputeCellLength(vtkLagrangianParticle* pa
}
else if ((this->CellLengthComputationMode == STEP_CUR_CELL_DIV_THEO ||
this->CellLengthComputationMode == STEP_LAST_CELL_DIV_THEO) &&
vtkMath::Norm(vel) > 0.0 && cell->GetCellType() != VTK_VOXEL) //TODO ?
vtkMath::Norm(vel) > 0.0)
{
double velHat[3] = { vel[0], vel[1], vel[2] };
vtkMath::Normalize(velHat);
......
......@@ -173,12 +173,10 @@ public:
* and the edges of the last cell the particle was in.
* - STEP_LAST_CELL_DIV_THEO :
* Compute cell length using the particle velocity
* and the divergence theorem, not supported
* with vtkVoxel, fallback to STEP_LAST_CELL_LENGTH
* and the divergence theorem.
* - STEP_CUR_CELL_DIV_THEO :
* Compute cell length using the particle velocity
* and the divergence theorem, not supported
* with vtkVoxel, fallback to STEP_CUR_CELL_LENGTH
* and the divergence theorem.
* Default is STEP_LAST_CELL_LENGTH.
*/
vtkSetMacro(CellLengthComputationMode, int);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment