/usr/include/sc_bspline.h is in libp4est-dev 1.1-4.
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 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 | /*
This file is part of the SC Library.
The SC Library provides support for parallel scientific applications.
Copyright (C) 2010 The University of Texas System
The SC Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The SC Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the SC Library; if not, write to the Free Software
Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
02110-1301, USA.
*/
#ifndef SC_BSPLINE_H
#define SC_BSPLINE_H
#include <sc_dmatrix.h>
SC_EXTERN_C_BEGIN;
typedef struct
{
int d; /** Dimensionality of control points */
int p; /** Number of control points is p + 1 */
int n; /** Polynomial degree is n >= 0 */
int m; /** Number of knots is m + 1 = n + p + 2 */
int l; /** Number of internal intervals l = m - 2 * n > 0 */
int cacheknot; /* previously evaluated knot interval */
sc_dmatrix_t *points; /* (p + 1) x d array of points, not owned */
sc_dmatrix_t *knots; /* m + 1 array of knots */
int knots_owned;
sc_dmatrix_t *works; /* Workspace ((n + 1) * (n + 1)) x d */
int works_owned;
}
sc_bspline_t;
/** Compute the minimum required number of points for a certain degree.
* \param [in] n Polynomial degree of the spline functions, n >= 0.
* \return Return minimum point number = p + 1 >= n + 1.
*/
int sc_bspline_min_number_points (int n);
/** Compute the minimum required number of knots for a certain degree.
* \param [in] n Polynomial degree of the spline functions, n >= 0.
* \return Return minimum knot number = m + 1 >= 2 * n + 2.
*/
int sc_bspline_min_number_knots (int n);
/** Create a uniform B-spline knot vector.
* \param [in] n Polynomial degree of the spline functions, n >= 0.
* \param [in] points (p + 1) x d array of points in R^d, p >= 0, d >= 1.
* \return (n + p + 2) x 1 array of knots.
*/
sc_dmatrix_t *sc_bspline_knots_new (int n, sc_dmatrix_t * points);
/** Create a B-spline knots array roughly proportional to the arc length.
* This works only for at least linear B-splines, n >= 1.
* \param [in] n Polynomial degree of the spline functions, n >= 1.
* \param [in] points (p + 1) x d array of points in R^d, p >= 0, d >= 1.
* \return (n + p + 2) x 1 array of knots.
*/
sc_dmatrix_t *sc_bspline_knots_new_length (int n,
sc_dmatrix_t * points);
/** Create a uniform B-spline knot vector for a periodic B-spline.
* \param [in] n Polynomial degree of the spline functions, n >= 0.
* \param [in] points (p + 1) x d array of points in R^d, p >= 0, d >= 1.
* \return (n + p + 2) x 1 array of knots.
*/
sc_dmatrix_t *sc_bspline_knots_new_periodic (int n,
sc_dmatrix_t * points);
/** Create a B-spline knots array roughly proportional to the arc length for a
* periodic B-spline.
* This works only for at least linear B-splines, n >= 1.
* \param [in] n Polynomial degree of the spline functions, n >= 1.
* \param [in] points (p + 1) x d array of points in R^d, p >= 0, d >= 1.
* \return (n + p + 2) x 1 array of knots.
*/
sc_dmatrix_t *sc_bspline_knots_new_length_periodic (int n,
sc_dmatrix_t *
points);
/** Take a vector of points and make them appropriate for a periodic B-spine.
* \param [in] n Polynomial degree of the spline functions, n >= 0.
* \param [in,out] points On input, an (l x d) array of points in R^d, l > 0,
* d >= 1. If n is odd, these points are associated
* with the left endpoints of the intervals over
* which the B-spline is fully defined; if n is odd,
* they are associated with the midpoints.
* On output, points is a (l + n) x d, array of points
* appropriate for use as control points for a periodic
* B-spline.
*/
void sc_bspline_make_points_periodic (int n, sc_dmatrix_t *
points);
/** Create workspace for B-spline evaluation.
* \param [in] n Polynomial degree of the spline functions, n >= 0.
* \param [in] d Dimension of the control points in R^d, d >= 1.
* \return Workspace ((n + 1) * (n + 1)) x d.
*/
sc_dmatrix_t *sc_bspline_workspace_new (int n, int d);
/** Create a new B-spline structure.
* \param [in] n Polynomial degree of the spline functions, n >= 0.
* \param [in] points (p + 1) x d array of points in R^d, p >= 0, d >= 1.
* Borrowed, matrix is not copied so it must not be
* destroyed while the B-spline structure is in use.
* \param [in] knots (n + p + 2) x 1 array of knots. Borrowed.
* If NULL the knots are computed equidistantly.
* \param [in] works Workspace ((n + 1) * (n + 1)) x d. Borrowed.
* If NULL the workspace is allocated internally.
*/
sc_bspline_t *sc_bspline_new (int n, sc_dmatrix_t * points,
sc_dmatrix_t * knots,
sc_dmatrix_t * works);
/** Destroy a B-spline structure.
*/
void sc_bspline_destroy (sc_bspline_t * bs);
/** Evaluate a B-spline at a certain point.
* \param [in] bs B-spline structure.
* \param [in] t Value that must be within the range of the knots.
* \param [out] result The computed point in R^d is placed here.
*/
void sc_bspline_evaluate (sc_bspline_t * bs,
double t, double *result);
/** Evaluate a B-spline derivative at a certain point.
* \param [in] bs B-spline structure.
* \param [in] t Value that must be within the range of the knots.
* \param [out] result The computed derivative in R^d is placed here.
*/
void sc_bspline_derivative (sc_bspline_t * bs,
double t, double *result);
/** Evaluate any order B-spline derivative at a certain point.
* \param [in] bs B-spline structure.
* \param [in] order Order of the derivative >= 0.
* \param [in] t Value that must be within the range of the knots.
* \param [out] result The computed derivative in R^d is placed here.
*/
void sc_bspline_derivative_n (sc_bspline_t * bs, int order,
double t, double *result);
/** Evaluate a B-spline derivative at a certain point. Obsolete.
* \param [in] bs B-spline structure.
* \param [in] t Value that must be within the range of the knots.
* \param [out] result The computed derivative in R^d is placed here.
*/
void sc_bspline_derivative2 (sc_bspline_t * bs,
double t, double *result);
SC_EXTERN_C_END;
#endif /* !SC_BSPLINE_H */
|