/usr/include/vtk-7.1/vtkSPHQuarticKernel.h is in libvtk7-dev 7.1.1+dfsg1-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | /*=========================================================================
Program: Visualization Toolkit
Module: vtkSPHQuarticKernel.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 vtkSPHQuarticKernel
* @brief a quartic SPH interpolation kernel
*
*
* vtkSPHQuarticKernel is an smooth particle hydrodynamics interpolation kernel as
* described by D.J. Price. This is a quartic 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 vtkSPHQuarticKernel_h
#define vtkSPHQuarticKernel_h
#include "vtkFiltersPointsModule.h" // For export macro
#include "vtkSPHKernel.h"
#include <algorithm> // For std::min()
class vtkIdList;
class vtkDoubleArray;
class VTKFILTERSPOINTS_EXPORT vtkSPHQuarticKernel : public vtkSPHKernel
{
public:
//@{
/**
* Standard methods for instantiation, obtaining type information, and printing.
*/
static vtkSPHQuarticKernel *New();
vtkTypeMacro(vtkSPHQuarticKernel,vtkSPHKernel);
void PrintSelf(ostream& os, vtkIndent indent);
//@}
/**
* Produce the computational parameters for the kernel. Invoke this method
* after setting initial values like SpatialStep.
*/
virtual void Initialize(vtkAbstractPointLocator *loc, vtkDataSet *ds,
vtkPointData *pd);
//@{
/**
* Compute weighting factor given a normalized distance from a sample point.
*/
virtual double ComputeFunctionWeight(const double d)
{
double tmp1 = 2.5 - std::min(d,2.5);
double tmp2 = 1.5 - std::min(d,1.5);
double tmp3 = 0.5 - std::min(d,0.5);
return (tmp1*tmp1*tmp1*tmp1 - 5.0*tmp2*tmp2*tmp2*tmp2 +
10.0*tmp3*tmp3*tmp3*tmp3);
}
//@}
//@{
/**
* Compute weighting factor for derivative quantities given a normalized
* distance from a sample point.
*/
virtual double ComputeDerivWeight(const double d)
{
double tmp1 = 2.5 - std::min(d,2.5);
double tmp2 = 1.5 - std::min(d,1.5);
double tmp3 = 0.5 - std::min(d,0.5);
return (-4.0*tmp1*tmp1*tmp1 + 20.0*tmp2*tmp2*tmp2 - 40.0*tmp3*tmp3*tmp3);
}
//@}
protected:
vtkSPHQuarticKernel();
~vtkSPHQuarticKernel();
private:
vtkSPHQuarticKernel(const vtkSPHQuarticKernel&) VTK_DELETE_FUNCTION;
void operator=(const vtkSPHQuarticKernel&) VTK_DELETE_FUNCTION;
};
#endif
|