/usr/include/libwildmagic/Wm5IntpAkimaUniform2.h is in libwildmagic-dev 5.13-1+b2.
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 | // Geometric Tools, LLC
// Copyright (c) 1998-2014
// Distributed under the Boost Software License, Version 1.0.
// http://www.boost.org/LICENSE_1_0.txt
// http://www.geometrictools.com/License/Boost/LICENSE_1_0.txt
//
// File Version: 5.0.1 (2010/10/01)
#ifndef WM5INTPAKIMAUNIFORM2_H
#define WM5INTPAKIMAUNIFORM2_H
#include "Wm5MathematicsLIB.h"
namespace Wm5
{
template <typename Real>
class WM5_MATHEMATICS_ITEM IntpAkimaUniform2
{
public:
// Construction and destruction. IntpAkimaUniform2 does not accept
// responsibility for deleting the input array. The application must do
// so. The interpolator is for uniformly spaced (x,y)-values. The
// function values are assumed to be organized as f(x,y) = F[y][x].
IntpAkimaUniform2 (int xBound, int yBound, Real xMin, Real xSpacing,
Real yMin, Real ySpacing, Real** F);
~IntpAkimaUniform2 ();
class WM5_MATHEMATICS_ITEM Polynomial
{
public:
Polynomial ();
// P(x,y) = (1,x,x^2,x^3)*A*(1,y,y^2,y^3). The matrix term A[ix][iy]
// corresponds to the polynomial term x^{ix} y^{iy}.
Real& A (int ix, int iy);
Real operator() (Real x, Real y) const;
Real operator() (int xOrder, int yOrder, Real x, Real y) const;
protected:
Real mCoeff[4][4];
};
int GetXBound () const;
int GetYBound () const;
int GetQuantity () const;
Real** GetF () const;
Polynomial** GetPolynomials () const;
const Polynomial& GetPolynomial (int ix, int iy) const;
Real GetXMin () const;
Real GetXMax () const;
Real GetXSpacing () const;
Real GetYMin () const;
Real GetYMax () const;
Real GetYSpacing () const;
// Evaluate the function and its derivatives. The application is
// responsible for ensuring that xmin <= x <= xmax and ymin <= y <= ymax.
// If (x,y) is outside the extremes, the function returns MAXREAL. The
// first operator is for function evaluation. The second operator is for
// function or derivative evaluations. The uiXOrder argument is the order
// of the x-derivative and the uiYOrder argument is the order of the
// y-derivative. Both orders are zero to get the function value itself.
Real operator() (Real x, Real y) const;
Real operator() (int xOrder, int yOrder, Real x, Real y) const;
private:
Real ComputeDerivative (Real* slope) const;
void Construct (Polynomial& poly, Real F[2][2], Real FX[2][2],
Real FY[2][2], Real FXY[2][2]);
bool XLookup (Real x, int& xIndex, Real& dx) const;
bool YLookup (Real y, int& yIndex, Real& dy) const;
int mXBound, mYBound, mQuantity;
Real** mF;
Polynomial** mPoly;
Real mXMin, mXMax, mXSpacing;
Real mYMin, mYMax, mYSpacing;
};
typedef IntpAkimaUniform2<float> IntpAkimaUniform2f;
typedef IntpAkimaUniform2<double> IntpAkimaUniform2d;
}
#endif
|