VTK  9.2.5
vtkSPHQuarticKernel.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkSPHQuarticKernel.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 =========================================================================*/
40 #ifndef vtkSPHQuarticKernel_h
41 #define vtkSPHQuarticKernel_h
42 
43 #include "vtkFiltersPointsModule.h" // For export macro
44 #include "vtkSPHKernel.h"
45 #include <algorithm> // For std::min()
46 
47 class vtkIdList;
48 class vtkDoubleArray;
49 
50 class VTKFILTERSPOINTS_EXPORT vtkSPHQuarticKernel : public vtkSPHKernel
51 {
52 public:
54 
59  void PrintSelf(ostream& os, vtkIndent indent) override;
61 
67 
69 
72  double ComputeFunctionWeight(const double d) override
73  {
74  double tmp1 = 2.5 - std::min(d, 2.5);
75  double tmp2 = 1.5 - std::min(d, 1.5);
76  double tmp3 = 0.5 - std::min(d, 0.5);
77  return (tmp1 * tmp1 * tmp1 * tmp1 - 5.0 * tmp2 * tmp2 * tmp2 * tmp2 +
78  10.0 * tmp3 * tmp3 * tmp3 * tmp3);
79  }
81 
83 
87  double ComputeDerivWeight(const double d) override
88  {
89  double tmp1 = 2.5 - std::min(d, 2.5);
90  double tmp2 = 1.5 - std::min(d, 1.5);
91  double tmp3 = 0.5 - std::min(d, 0.5);
92  return (-4.0 * tmp1 * tmp1 * tmp1 + 20.0 * tmp2 * tmp2 * tmp2 - 40.0 * tmp3 * tmp3 * tmp3);
93  }
95 
96 protected:
99 
100 private:
101  vtkSPHQuarticKernel(const vtkSPHQuarticKernel&) = delete;
102  void operator=(const vtkSPHQuarticKernel&) = delete;
103 };
104 
105 #endif
abstract class to quickly locate points in 3-space
abstract class to specify dataset behavior
Definition: vtkDataSet.h:172
dynamic, self-adjusting array of double
list of point or cell ids
Definition: vtkIdList.h:143
a simple class to control print indentation
Definition: vtkIndent.h:119
represent and manipulate point attribute data
Definition: vtkPointData.h:151
a family of SPH interpolation kernels
Definition: vtkSPHKernel.h:68
a quartic SPH interpolation kernel
void Initialize(vtkAbstractPointLocator *loc, vtkDataSet *ds, vtkPointData *pd) override
Produce the computational parameters for the kernel.
void PrintSelf(ostream &os, vtkIndent indent) override
Standard methods for instantiation, obtaining type information, and printing.
static vtkSPHQuarticKernel * New()
Standard methods for instantiation, obtaining type information, and printing.
double ComputeFunctionWeight(const double d) override
Compute weighting factor given a normalized distance from a sample point.
~vtkSPHQuarticKernel() override
double ComputeDerivWeight(const double d) override
Compute weighting factor for derivative quantities given a normalized distance from a sample point.