VTK
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 =========================================================================*/
37 #ifndef vtkSPHQuarticKernel_h
38 #define vtkSPHQuarticKernel_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 
48 class VTKFILTERSPOINTS_EXPORT vtkSPHQuarticKernel : public vtkSPHKernel
49 {
50 public:
52 
55  static vtkSPHQuarticKernel *New();
57  void PrintSelf(ostream& os, vtkIndent indent) override;
59 
65  vtkPointData *pd) override;
66 
68 
71  double ComputeFunctionWeight(const double d) override
72  {
73  double tmp1 = 2.5 - std::min(d,2.5);
74  double tmp2 = 1.5 - std::min(d,1.5);
75  double tmp3 = 0.5 - std::min(d,0.5);
76  return (tmp1*tmp1*tmp1*tmp1 - 5.0*tmp2*tmp2*tmp2*tmp2 +
77  10.0*tmp3*tmp3*tmp3*tmp3);
78  }
80 
82 
86  double ComputeDerivWeight(const double d) override
87  {
88  double tmp1 = 2.5 - std::min(d,2.5);
89  double tmp2 = 1.5 - std::min(d,1.5);
90  double tmp3 = 0.5 - std::min(d,0.5);
91  return (-4.0*tmp1*tmp1*tmp1 + 20.0*tmp2*tmp2*tmp2 - 40.0*tmp3*tmp3*tmp3);
92  }
94 
95 protected:
97  ~vtkSPHQuarticKernel() override;
98 
99 private:
100  vtkSPHQuarticKernel(const vtkSPHQuarticKernel&) = delete;
101  void operator=(const vtkSPHQuarticKernel&) = delete;
102 };
103 
104 #endif
represent and manipulate point attribute data
Definition: vtkPointData.h:31
abstract class to specify dataset behavior
Definition: vtkDataSet.h:56
double ComputeFunctionWeight(const double d) override
Compute weighting factor given a normalized distance from a sample point.
dynamic, self-adjusting array of double
a simple class to control print indentation
Definition: vtkIndent.h:33
double ComputeDerivWeight(const double d) override
Compute weighting factor for derivative quantities given a normalized distance from a sample point.
abstract class to quickly locate points in 3-space
list of point or cell ids
Definition: vtkIdList.h:30
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.
a quartic SPH interpolation kernel
static vtkObject * New()
Create an object with Debug turned off, modified time initialized to zero, and reference counting on.
a family of SPH interpolation kernels
Definition: vtkSPHKernel.h:68