/usr/include/libwildmagic/Wm5NURBSCurve2.h is in libwildmagic-dev 5.13-1ubuntu1.
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 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 | // 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 WM5NURBSCURVE2_H
#define WM5NURBSCURVE2_H
#include "Wm5MathematicsLIB.h"
#include "Wm5SingleCurve2.h"
#include "Wm5BSplineBasis.h"
namespace Wm5
{
template <typename Real>
class WM5_MATHEMATICS_ITEM NURBSCurve2 : public SingleCurve2<Real>
{
public:
// Construction and destruction. The caller is responsible for deleting
// the input arrays if they were dynamically allocated. Internal copies
// of the arrays are made, so to dynamically change control points,
// control weights, or knots, you must use the 'SetControlPoint',
// 'GetControlPoint', 'SetControlWeight', 'GetControlWeight', and 'Knot'
// member functions.
// The homogeneous input points are (x,y,w) where the (x,y) values are
// stored in the ctrlPoin array and the w values are stored in the
// ctrlWeight array. The output points from curve evaluations are of
// the form (x',y') = (x/w,y/w).
// Uniform spline. The number of control points is n+1 >= 2. The degree
// of the spline is d and must satisfy 1 <= d <= n. The knots are
// implicitly calculated in [0,1]. If open is 'true', the spline is
// open and the knots are
// t[i] = 0, 0 <= i <= d
// (i-d)/(n+1-d), d+1 <= i <= n
// 1, n+1 <= i <= n+d+1
// If open is 'false', the spline is periodic and the knots are
// t[i] = (i-d)/(n+1-d), 0 <= i <= n+d+1
// If loop is 'true', extra control points are added to generate a closed
// curve. For an open spline, the control point array is reallocated and
// one extra control point is added, set to the first control point
// C[n+1] = C[0]. For a periodic spline, the control point array is
// reallocated and the first d points are replicated. In either case the
// knot array is calculated accordingly.
NURBSCurve2 (int numCtrlPoints, const Vector2<Real>* ctrlPoint,
const Real* ctrlWeight, int degree, bool loop, bool open);
// Open, nonuniform spline. The knot array must have n-d elements. The
// elements must be nondecreasing. Each element must be in [0,1].
NURBSCurve2 (int numCtrlPoints, const Vector2<Real>* ctrlPoint,
const Real* ctrlWeight, int degree, bool loop, const Real* knot);
virtual ~NURBSCurve2 ();
int GetNumCtrlPoints () const;
int GetDegree () const;
bool IsOpen () const;
bool IsUniform () const;
bool IsLoop () const;
// Control points and weights may be changed at any time. The input index
// should be valid (0 <= i <= n). If it is invalid, the return value of
// GetControlPoint is a vector whose components are all MAX_REAL, and the
// return value of GetControlWeight is MAX_REAL.
// undefined.
void SetControlPoint (int i, const Vector2<Real>& ctrl);
Vector2<Real> GetControlPoint (int i) const;
void SetControlWeight (int i, Real weight);
Real GetControlWeight (int i) const;
// The knot values can be changed only if the basis function is nonuniform
// and the input index is valid (0 <= i <= n-d-1). If these conditions
// are not satisfied, GetKnot returns MAX_REAL.
void SetKnot (int i, Real knot);
Real GetKnot (int i) const;
// The spline is defined for 0 <= t <= 1. If a t-value is outside [0,1],
// an open spline clamps t to [0,1]. That is, if t > 1, t is set to 1;
// if t < 0, t is set to 0. A periodic spline wraps to to [0,1]. That
// is, if t is outside [0,1], then t is set to t-floor(t).
virtual Vector2<Real> GetPosition (Real t) const;
virtual Vector2<Real> GetFirstDerivative (Real t) const;
virtual Vector2<Real> GetSecondDerivative (Real t) const;
virtual Vector2<Real> GetThirdDerivative (Real t) const;
// If you need position and derivatives at the same time, it is more
// efficient to call these functions. Pass the addresses of those
// quantities whose values you want. You may pass 0 in any argument
// whose value you do not want.
void Get (Real t, Vector2<Real>* pos, Vector2<Real>* der1,
Vector2<Real>* der2, Vector2<Real>* der3) const;
// Access the basis function to compute it without control points. This
// is useful for least squares fitting of curves.
BSplineBasis<Real>& GetBasis ();
protected:
// Replicate the necessary number of control points when the Create
// function has loop equal to true, in which case the spline curve must
// be a closed curve.
void CreateControl (const Vector2<Real>* ctrlPoint,
const Real* ctrlWeight);
int mNumCtrlPoints;
Vector2<Real>* mCtrlPoint; // ctrl[n+1]
Real* mCtrlWeight; // weight[n+1]
bool mLoop;
BSplineBasis<Real> mBasis;
int mReplicate; // the number of replicated control points
};
typedef NURBSCurve2<float> NURBSCurve2f;
typedef NURBSCurve2<double> NURBSCurve2d;
}
#endif
|