/usr/include/paraview/vtkSPHQuarticKernel.h is in paraview-dev 5.1.2+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 | /*=========================================================================
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.
=========================================================================*/
// .NAME vtkSPHQuarticKernel - a quartic SPH interpolation kernel
// .SECTION Description
// vtkSPHQuarticKernel is an smooth particle hydrodynamics interpolation kernel as
// described by D.J. Price. This is a quartic formulation.
//
// .SECTION Caveats
// For more information see D.J. Price, Smoothed particle hydrodynamics and
// magnetohydrodynamics, J. Comput. Phys. 231:759-794, 2012. Especially
// equation 49.
// .SECTION Acknowledgments
// The following work has been generously supported by Altair Engineering
// and FluiDyna GmbH, Please contact Steve Cosgrove or Milos Stanic for
// more information.
// .SECTION See Also
// 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:
// Description:
// Standard methods for instantiation, obtaining type information, and printing.
static vtkSPHQuarticKernel *New();
vtkTypeMacro(vtkSPHQuarticKernel,vtkSPHKernel);
void PrintSelf(ostream& os, vtkIndent indent);
// Description:
// 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);
}
// Description:
// Compute weighting factor for derivative quantities given a normalized
// distance from a sample point.
virtual double ComputeGradientWeight(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 - 5.0*tmp2*tmp2*tmp2 + 10.0*tmp3*tmp3*tmp3);
}
protected:
vtkSPHQuarticKernel();
~vtkSPHQuarticKernel();
private:
vtkSPHQuarticKernel(const vtkSPHQuarticKernel&); // Not implemented.
void operator=(const vtkSPHQuarticKernel&); // Not implemented.
};
#endif
|