VTK  9.1.0
vtkSPHQuinticKernel.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkSPHQuinticKernel.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 =========================================================================*/
37 #ifndef vtkSPHQuinticKernel_h
38 #define vtkSPHQuinticKernel_h
39 
40 #include "vtkFiltersPointsModule.h" // For export macro
41 #include "vtkSPHKernel.h"
42 #include <algorithm> // For std::min()
43 
44 class vtkIdList;
45 class vtkDoubleArray;
46 
47 class VTKFILTERSPOINTS_EXPORT vtkSPHQuinticKernel : public vtkSPHKernel
48 {
49 public:
51 
56  void PrintSelf(ostream& os, vtkIndent indent) override;
58 
64 
66 
69  double ComputeFunctionWeight(const double d) override
70  {
71  double tmp1 = 3.0 - std::min(d, 3.0);
72  double tmp2 = 2.0 - std::min(d, 2.0);
73  double tmp3 = 1.0 - std::min(d, 1.0);
74  return (tmp1 * tmp1 * tmp1 * tmp1 * tmp1 - 6.0 * tmp2 * tmp2 * tmp2 * tmp2 * tmp2 +
75  15.0 * tmp3 * tmp3 * tmp3 * tmp3 * tmp3);
76  }
78 
80 
84  double ComputeDerivWeight(const double d) override
85  {
86  double tmp1 = 3.0 - std::min(d, 3.0);
87  double tmp2 = 2.0 - std::min(d, 2.0);
88  double tmp3 = 1.0 - std::min(d, 1.0);
89  return (-5.0 * tmp1 * tmp1 * tmp1 * tmp1 + 30.0 * tmp2 * tmp2 * tmp2 * tmp2 +
90  -75.0 * tmp3 * tmp3 * tmp3 * tmp3);
91  }
93 
94 protected:
97 
98 private:
100  void operator=(const vtkSPHQuinticKernel&) = delete;
101 };
102 
103 #endif
abstract class to quickly locate points in 3-space
abstract class to specify dataset behavior
Definition: vtkDataSet.h:57
dynamic, self-adjusting array of double
list of point or cell ids
Definition: vtkIdList.h:31
a simple class to control print indentation
Definition: vtkIndent.h:34
represent and manipulate point attribute data
Definition: vtkPointData.h:33
a family of SPH interpolation kernels
Definition: vtkSPHKernel.h:68
a quintic SPH interpolation kernel
void Initialize(vtkAbstractPointLocator *loc, vtkDataSet *ds, vtkPointData *pd) override
Produce the computational parameters for the kernel.
double ComputeDerivWeight(const double d) override
Compute weighting factor for derivative quantities given a normalized distance from a sample point.
double ComputeFunctionWeight(const double d) override
Compute weighting factor given a normalized distance from a sample point.
void PrintSelf(ostream &os, vtkIndent indent) override
Standard methods for instantiation, obtaining type information, and printing.
~vtkSPHQuinticKernel() override
static vtkSPHQuinticKernel * New()
Standard methods for instantiation, obtaining type information, and printing.