forked from Kitware/VTK
-
Notifications
You must be signed in to change notification settings - Fork 0
/
vtkSPHQuinticKernel.h
105 lines (90 loc) · 3.03 KB
/
vtkSPHQuinticKernel.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
/*=========================================================================
Program: Visualization Toolkit
Module: vtkSPHQuinticKernel.h
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
/**
* @class vtkSPHQuinticKernel
* @brief a quintic SPH interpolation kernel
*
*
* vtkSPHQuinticKernel is an smooth particle hydrodynamics interpolation kernel as
* described by D.J. Price. This is a quintic formulation.
*
* @warning
* FOr more information see D.J. Price, Smoothed particle hydrodynamics and
* magnetohydrodynamics, J. Comput. Phys. 231:759-794, 2012. Especially
* equation 49.
*
* @par Acknowledgments:
* The following work has been generously supported by Altair Engineering
* and FluiDyna GmbH. Please contact Steve Cosgrove or Milos Stanic for
* more information.
*
* @sa
* vtkSPHKernel vtkSPHInterpolator
*/
#ifndef vtkSPHQuinticKernel_h
#define vtkSPHQuinticKernel_h
#include "vtkFiltersPointsModule.h" // For export macro
#include "vtkSPHKernel.h"
#include <algorithm> // For std::min()
class vtkIdList;
class vtkDoubleArray;
class VTKFILTERSPOINTS_EXPORT vtkSPHQuinticKernel : public vtkSPHKernel
{
public:
//@{
/**
* Standard methods for instantiation, obtaining type information, and printing.
*/
static vtkSPHQuinticKernel *New();
vtkTypeMacro(vtkSPHQuinticKernel,vtkSPHKernel);
void PrintSelf(ostream& os, vtkIndent indent) override;
//@}
/**
* Produce the computational parameters for the kernel. Invoke this method
* after setting initial values like SpatialStep.
*/
void Initialize(vtkAbstractPointLocator *loc, vtkDataSet *ds,
vtkPointData *pd) override;
//@{
/**
* Compute weighting factor given a normalized distance from a sample point.
*/
double ComputeFunctionWeight(const double d) override
{
double tmp1 = 3.0 - std::min(d,3.0);
double tmp2 = 2.0 - std::min(d,2.0);
double tmp3 = 1.0 - std::min(d,1.0);
return (tmp1*tmp1*tmp1*tmp1*tmp1 - 6.0*tmp2*tmp2*tmp2*tmp2*tmp2 +
15.0*tmp3*tmp3*tmp3*tmp3*tmp3);
}
//@}
//@{
/**
* Compute weighting factor for derivative quantities given a normalized
* distance from a sample point.
*/
double ComputeDerivWeight(const double d) override
{
double tmp1 = 3.0 - std::min(d,3.0);
double tmp2 = 2.0 - std::min(d,2.0);
double tmp3 = 1.0 - std::min(d,1.0);
return (-5.0*tmp1*tmp1*tmp1*tmp1 + 30.0*tmp2*tmp2*tmp2*tmp2 +
-75.0*tmp3*tmp3*tmp3*tmp3);
}
//@}
protected:
vtkSPHQuinticKernel();
~vtkSPHQuinticKernel() override;
private:
vtkSPHQuinticKernel(const vtkSPHQuinticKernel&) = delete;
void operator=(const vtkSPHQuinticKernel&) = delete;
};
#endif