@itwin/core-geometry 5.0.0-dev.9 → 5.0.0-dev.90

package/lib/esm/bspline/BSpline1dNd.d.ts

@@ -3,108 +3,144 @@
  */
 import { Point3d } from "../geometry3d/Point3dVector3d";
 import { BSplineWrapMode, KnotVector } from "./KnotVector";
-/** Bspline knots and poles for 1d-to-Nd.
+/**
+ * Knots and poles for a B-spline function mapping R to R^n.
  * * The "pole" (aka control point) of this class is a block of `poleLength` numbers.
  * * Derived classes (not this class) assign meaning such as x,y,z,w.
- * * for instance, an instance of this class with `poleLength===3` does not know if its poles are x,y,z or weighted 2D x,y,w
+ * * For instance, an instance of this class with `poleLength===3` does not know if its poles are x,y,z or
+ * weighted 2D x,y,w.
  * @public
  */
 export declare class BSpline1dNd {
-    /** knots of the bspline */
+    /** Knots of the bspline. */
     knots: KnotVector;
-    /** poles, packed in blocks of `poleLength` doubles. */
+    /** Poles packed in blocks of `poleLength` doubles. */
     packedData: Float64Array;
-    /** (property accessor) Return the number of numeric values per pole. */
+    /** The number of numeric values per pole. */
     poleLength: number;
     /** (property accessor) Return the degree of the polynomials. */
     get degree(): number;
-    /** (property accessor) Return the number of order (one more than degree) of the polynomials */
+    /** (property accessor) Return the order (one more than degree) of the polynomials. */
     get order(): number;
-    /** (property accessor) Return the number of bezier spans (including null spans at multiple knots)*/
+    /** (property accessor) Return the number of bezier spans (including null spans at multiple knots). */
     get numSpan(): number;
-    /** (property accessor)  Return the number of poles*/
+    /** (property accessor) Return the number of poles. */
     get numPoles(): number;
-    /** copy 3 values of pole `i` into a point.
+    /**
+     * Copy 3 values of pole `i` into a point.
      * * The calling class is responsible for knowing if this is an appropriate access to the blocked data.
      */
     getPoint3dPole(i: number, result?: Point3d): Point3d | undefined;
-    /** preallocated array (length === `order`) used as temporary in evaluations */
+    /**
+     * Values of the `order` relevant B-spline basis functions at a parameter.
+     * * Preallocated to length `order` in the constructor and used as a temporary in evaluations.
+     */
     basisBuffer: Float64Array;
-    /** preallocated array (length === `poleLength`) used as temporary in evaluations */
-    poleBuffer: Float64Array;
-    /** preallocated array (length === `order`) used as temporary in evaluations */
+    /**
+     * Derivatives of the `order` relevant B-spline basis functions at a parameter.
+     * * Preallocated to length `order` in the constructor and used as a temporary in evaluations.
+     */
     basisBuffer1: Float64Array;
-    /** preallocated array (length === `order`) used as temporary in evaluations */
+    /**
+     * Second derivatives of the `order` relevant B-spline basis functions at a parameter.
+     * * Preallocated to length `order` in the constructor and used as a temporary in evaluations.
+     */
     basisBuffer2: Float64Array;
-    /** preallocated array (length === `poleLength`) used as temporary in evaluations */
+    /**
+     * Temporary to hold a single point.
+     * * Preallocated to length `poleLength` in the constructor.
+     */
+    poleBuffer: Float64Array;
+    /**
+     * Temporary to hold a single derivative vector.
+     * * Preallocated to length `poleLength` in the constructor.
+     */
     poleBuffer1: Float64Array;
-    /** preallocated array (length === `poleLength`) used as temporary in evaluations */
+    /**
+     * Temporary to hold a single second derivative vector.
+     * * Preallocated to length `poleLength` in the constructor.
+     */
     poleBuffer2: Float64Array;
     /**
-     * initialize arrays for given spline dimensions.
-     * @param numPoles number of poles
-     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted, 3 for 2d weighted)
-     * @param order number of poles in support for a section of the bspline
-     * @param knots KnotVector.  This is captured, not cloned.
+     * Initialize arrays for given spline dimensions.
+     * @param numPoles number of poles.
+     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted,
+     * 3 for 2d weighted).
+     * @param order number of poles defining a Bezier segment of the B-spline function.
+     * @param knots the KnotVector. This is captured, not cloned.
      */
     protected constructor(numPoles: number, poleLength: number, order: number, knots: KnotVector);
     /**
-     * create a 1Bspline1dNd`
-     * @param numPoles number of poles
-     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted, 3 for 2d weighted)
-     * @param order number of poles in support for a section of the bspline
-     * @param knots KnotVector.  This is captured, not cloned.
+     * Create a `BSpline1dNd`.
+     * @param numPoles number of poles.
+     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted,
+     * 3 for 2d weighted).
+     * @param order number of poles defining a Bezier segment of the B-spline function.
+     * @param knots the KnotVector. This is captured, not cloned.
      */
-    static create(numPoles: number, poleLength: number, order: number, knots: KnotVector): BSpline1dNd | undefined;
-    /** Map a span index and local fraction to knot value. */
-    spanFractionToKnot(span: number, localFraction: number): number;
+    static create(numPoles: number, poleLength: number, order: number, knots: KnotVector): BSpline1dNd;
+    /** Map a span index and span fraction to knot value. */
+    spanFractionToKnot(spanIndex: number, spanFraction: number): number;
     /**
-     * Evaluate the `order` basis functions (and optionally one or two derivatives) at a given fractional position within indexed span.
-     * @returns true if and only if output arrays are sufficiently sized
+     * Evaluate the `order` basis functions (and optionally one or two derivatives) at a given fractional position within
+     * indexed span.
+     * @returns true if and only if output arrays are sufficiently sized.
     */
     evaluateBasisFunctionsInSpan(spanIndex: number, spanFraction: number, f: Float64Array, df?: Float64Array, ddf?: Float64Array): boolean;
+    /**
+     * Compute the linear combination of the given span's poles and the weights in `basisBuffer`, and store this point
+     * in `poleBuffer`.
+     */
+    sumPoleBufferForSpan(spanIndex: number): void;
+    /**
+     * Compute the linear combination of the given span's poles and the weights in `basisBuffer1`, and store this
+     * derivative vector in `poleBuffer1`.
+     */
+    sumPoleBuffer1ForSpan(spanIndex: number): void;
+    /**
+     * Compute the linear combination of the given span's poles and the weights in `basisBuffer2`, and store this
+     * second derivative vector in `poleBuffer2`.
+     */
+    sumPoleBuffer2ForSpan(spanIndex: number): void;
     /**
      * * Evaluate the basis functions at spanIndex and fraction.
-     *   * Evaluations are stored in the preallocated `this.basisBuffer`
-     * * Immediately do the summations of the basis values times the respective control points
+     *   * Evaluations are stored in the preallocated `this.basisBuffer`.
+     * * Immediately do the summations of the basis values times the respective poles.
      *   * Summations are stored in the preallocated `this.poleBuffer`
      * */
     evaluateBuffersInSpan(spanIndex: number, spanFraction: number): void;
     /**
      * * Evaluate the basis functions and one derivative at spanIndex and fraction.
-     *   * Evaluations are stored in the preallocated `this.basisBuffer`
-     * * Immediately do the summations of the basis values times the respective control points
+     *   * Function evaluations are stored in the preallocated `this.basisBuffer`; derivative evaluations in `this.basisBuffer1`.
+     * * Immediately do the summations of the basis values times the respective poles.
      *   * Summations are stored in the preallocated `this.poleBuffer` and `this.poleBuffer1`
      * */
     evaluateBuffersInSpan1(spanIndex: number, spanFraction: number): void;
-    /** sum poles in `poleBuffer` at span `spanIndex` by the weights in the `basisBuffer` */
-    sumPoleBufferForSpan(spanIndex: number): void;
-    /** sum poles in `poleBuffer1` at span `spanIndex` by the weights in the `basisBuffer1`, i.e. form first derivatives */
-    sumPoleBuffer1ForSpan(spanIndex: number): void;
-    /** sum poles in `poleBuffer2` at span `spanIndex` by the weights in the `basisBuffer2`, i.e. form second derivatives */
-    sumPoleBuffer2ForSpan(spanIndex: number): void;
-    /** Evaluate the function values and 1 or 2 derivatives into `this.poleBuffer`, `this.poleBuffer1` and `this.poleBuffer2` */
-    evaluateBuffersAtKnot(u: number, numDerivative?: number): void;
     /**
-     * Reverse the (blocked) poles (in `this.packedData` in place.
+     * Evaluate the B-spline function and optional derivatives at the given parameter in knot space.
+     * * Function value is stored in `poleBuffer`; optional derivative vectors in `poleBuffer1` and `poleBuffer2`.
      */
+    evaluateBuffersAtKnot(u: number, numDerivative?: number): void;
+    /** Reverse the instance poles and knots in place. */
     reverseInPlace(): void;
     /**
-     * Test if the leading and trailing polygon coordinates are replicated in the manner of a "closed" bspline polygon which has been expanded
-     * to act as a normal bspline.
-     * @returns true if `degree` leading and trailing polygon blocks match
-     * @deprecated in 4.x. Use testClosablePolygon instead.
+     * Test if the leading and trailing polygon coordinates are replicated in the manner of a "closed" bspline polygon
+     * which has been expanded to act as a normal bspline.
+     * @returns true if `degree` leading and trailing polygon blocks match.
+     * @deprecated in 4.x. Use `testClosablePolygon` instead.
      */
     testCloseablePolygon(mode?: BSplineWrapMode): boolean;
     /**
-     * Test if the leading and trailing polygon coordinates are replicated in the manner of a "closed" bspline polygon which has been expanded
-     * to act as a normal bspline.
-     * @returns true if `degree` leading and trailing polygon blocks match
+     * Test if the leading and trailing poles are replicated in the manner of a "closed" B-spline function with wraparound
+     * control polygon.
+     * @param mode wrap mode, indicating the number of wraparound poles to check. If undefined, `knots.wrappable` is used.
+     * @returns true if the expected leading and trailing poles match, according to `mode`.
      */
     testClosablePolygon(mode?: BSplineWrapMode): boolean;
-    /** Insert knot and resulting pole into the instance, optionally multiple times.
-     * @param knot the knot to be inserted (may already exist in the KnotVector)
-     * @param totalMultiplicity the total multiplicity of the knot on return
+    /**
+     * Insert the knot and resulting pole into the instance, optionally multiple times.
+     * @param knot the knot to be inserted (may already exist in the KnotVector).
+     * @param totalMultiplicity the total multiplicity of the knot on return.
      */
     addKnot(knot: number, totalMultiplicity: number): boolean;
 }

package/lib/esm/bspline/BSpline1dNd.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"BSpline1dNd.d.ts","sourceRoot":"","sources":["../../../src/bspline/BSpline1dNd.ts"],"names":[],"mappings":"AAIA;;GAEG;AAIH,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,eAAe,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE3D;;;;;GAKG;AACH,qBAAa,WAAW;IACtB,2BAA2B;IACpB,KAAK,EAAE,UAAU,CAAC;IACzB,uDAAuD;IAChD,UAAU,EAAE,YAAY,CAAC;IAChC,wEAAwE;IACjE,UAAU,EAAE,MAAM,CAAC;IAC1B,gEAAgE;IAChE,IAAW,MAAM,IAAI,MAAM,CAA8B;IACzD,+FAA+F;IAC/F,IAAW,KAAK,IAAI,MAAM,CAAkC;IAC5D,oGAAoG;IACpG,IAAW,OAAO,IAAI,MAAM,CAA8C;IAC1E,qDAAqD;IACrD,IAAW,QAAQ,IAAI,MAAM,CAAqD;IAClF;;OAEG;IACI,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IACvE,+EAA+E;IACxE,WAAW,EAAE,YAAY,CAAC;IACjC,oFAAoF;IAC7E,UAAU,EAAE,YAAY,CAAC;IAChC,+EAA+E;IACxE,YAAY,EAAE,YAAY,CAAC;IAClC,+EAA+E;IACxE,YAAY,EAAE,YAAY,CAAC;IAClC,oFAAoF;IAC7E,WAAW,EAAE,YAAY,CAAC;IACjC,oFAAoF;IAC7E,WAAW,EAAE,YAAY,CAAC;IAEjC;;;;;;OAMG;IACH,SAAS,aAAa,QAAQ,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,UAAU;IAY5F;;;;;;OAMG;WACW,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,UAAU,GAAG,WAAW,GAAG,SAAS;IAGrH,yDAAyD;IAClD,kBAAkB,CAAC,IAAI,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,GAAG,MAAM;IAItE;;;MAGE;IACK,4BAA4B,CAAC,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EAAE,CAAC,EAAE,YAAY,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,GAAG,CAAC,EAAE,YAAY,GAAG,OAAO;IAS7I;;;;;SAKK;IACE,qBAAqB,CAAC,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM;IAIpE;;;;;SAKK;IACE,sBAAsB,CAAC,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM;IAKrE,wFAAwF;IACjF,oBAAoB,CAAC,SAAS,EAAE,MAAM;IAO7C,uHAAuH;IAChH,qBAAqB,CAAC,SAAS,EAAE,MAAM;IAS9C,wHAAwH;IACjH,qBAAqB,CAAC,SAAS,EAAE,MAAM;IAS9C,4HAA4H;IACrH,qBAAqB,CAAC,CAAC,EAAE,MAAM,EAAE,aAAa,GAAE,MAAU;IAgBjE;;OAEG;IACI,cAAc,IAAI,IAAI;IAgB7B;;;;;OAKG;IACI,oBAAoB,CAAC,IAAI,CAAC,EAAE,eAAe,GAAG,OAAO;IAG5D;;;;OAIG;IACI,mBAAmB,CAAC,IAAI,CAAC,EAAE,eAAe,GAAG,OAAO;IAqB3D;;;OAGG;IACI,OAAO,CAAC,IAAI,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,GAAG,OAAO;CA2DjE"}
+{"version":3,"file":"BSpline1dNd.d.ts","sourceRoot":"","sources":["../../../src/bspline/BSpline1dNd.ts"],"names":[],"mappings":"AAIA;;GAEG;AAGH,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,eAAe,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE3D;;;;;;;GAOG;AACH,qBAAa,WAAW;IACtB,4BAA4B;IACrB,KAAK,EAAE,UAAU,CAAC;IACzB,sDAAsD;IAC/C,UAAU,EAAE,YAAY,CAAC;IAChC,6CAA6C;IACtC,UAAU,EAAE,MAAM,CAAC;IAC1B,gEAAgE;IAChE,IAAW,MAAM,IAAI,MAAM,CAE1B;IACD,sFAAsF;IACtF,IAAW,KAAK,IAAI,MAAM,CAEzB;IACD,sGAAsG;IACtG,IAAW,OAAO,IAAI,MAAM,CAE3B;IACD,sDAAsD;IACtD,IAAW,QAAQ,IAAI,MAAM,CAE5B;IACD;;;OAGG;IACI,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAGvE;;;OAGG;IACI,WAAW,EAAE,YAAY,CAAC;IACjC;;;OAGG;IACI,YAAY,EAAE,YAAY,CAAC;IAClC;;;OAGG;IACI,YAAY,EAAE,YAAY,CAAC;IAClC;;;OAGG;IACI,UAAU,EAAE,YAAY,CAAC;IAChC;;;OAGG;IACI,WAAW,EAAE,YAAY,CAAC;IACjC;;;OAGG;IACI,WAAW,EAAE,YAAY,CAAC;IACjC;;;;;;;OAOG;IACH,SAAS,aAAa,QAAQ,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,UAAU;IAW5F;;;;;;;OAOG;WACW,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,UAAU,GAAG,WAAW;IAGzG,wDAAwD;IACjD,kBAAkB,CAAC,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,GAAG,MAAM;IAG1E;;;;MAIE;IACK,4BAA4B,CACjC,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EAAE,CAAC,EAAE,YAAY,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,GAAG,CAAC,EAAE,YAAY,GAC9F,OAAO;IAWV;;;OAGG;IACI,oBAAoB,CAAC,SAAS,EAAE,MAAM;IAO7C;;;OAGG;IACI,qBAAqB,CAAC,SAAS,EAAE,MAAM;IAO9C;;;OAGG;IACI,qBAAqB,CAAC,SAAS,EAAE,MAAM;IAQ9C;;;;;SAKK;IACE,qBAAqB,CAAC,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM;IAIpE;;;;;SAKK;IACE,sBAAsB,CAAC,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM;IAKrE;;;OAGG;IACI,qBAAqB,CAAC,CAAC,EAAE,MAAM,EAAE,aAAa,GAAE,MAAU;IAgBjE,qDAAqD;IAC9C,cAAc,IAAI,IAAI;IAS7B;;;;;OAKG;IACI,oBAAoB,CAAC,IAAI,CAAC,EAAE,eAAe,GAAG,OAAO;IAG5D;;;;;OAKG;IACI,mBAAmB,CAAC,IAAI,CAAC,EAAE,eAAe,GAAG,OAAO;IAoB3D;;;;OAIG;IACI,OAAO,CAAC,IAAI,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,GAAG,OAAO;CAqDjE"}

package/lib/esm/bspline/BSpline1dNd.js

@@ -5,64 +5,115 @@
 /** @packageDocumentation
  * @module Bspline
  */
-// import { Point2d } from "../Geometry2d";
 import { Geometry } from "../Geometry";
 import { Point3d } from "../geometry3d/Point3dVector3d";
 import { BSplineWrapMode, KnotVector } from "./KnotVector";
-/** Bspline knots and poles for 1d-to-Nd.
+/**
+ * Knots and poles for a B-spline function mapping R to R^n.
  * * The "pole" (aka control point) of this class is a block of `poleLength` numbers.
  * * Derived classes (not this class) assign meaning such as x,y,z,w.
- * * for instance, an instance of this class with `poleLength===3` does not know if its poles are x,y,z or weighted 2D x,y,w
+ * * For instance, an instance of this class with `poleLength===3` does not know if its poles are x,y,z or
+ * weighted 2D x,y,w.
  * @public
  */
 export class BSpline1dNd {
+    /** Knots of the bspline. */
+    knots;
+    /** Poles packed in blocks of `poleLength` doubles. */
+    packedData;
+    /** The number of numeric values per pole. */
+    poleLength;
     /** (property accessor) Return the degree of the polynomials. */
-    get degree() { return this.knots.degree; }
-    /** (property accessor) Return the number of order (one more than degree) of the polynomials */
-    get order() { return this.knots.degree + 1; }
-    /** (property accessor) Return the number of bezier spans (including null spans at multiple knots)*/
-    get numSpan() { return this.numPoles - this.knots.degree; }
-    /** (property accessor)  Return the number of poles*/
-    get numPoles() { return this.packedData.length / this.poleLength; }
-    /** copy 3 values of pole `i` into a point.
+    get degree() {
+        return this.knots.degree;
+    }
+    /** (property accessor) Return the order (one more than degree) of the polynomials. */
+    get order() {
+        return this.knots.degree + 1;
+    }
+    /** (property accessor) Return the number of bezier spans (including null spans at multiple knots). */
+    get numSpan() {
+        return this.numPoles - this.knots.degree;
+    }
+    /** (property accessor) Return the number of poles. */
+    get numPoles() {
+        return this.packedData.length / this.poleLength;
+    }
+    /**
+     * Copy 3 values of pole `i` into a point.
      * * The calling class is responsible for knowing if this is an appropriate access to the blocked data.
      */
-    getPoint3dPole(i, result) { return Point3d.createFromPacked(this.packedData, i, result); }
+    getPoint3dPole(i, result) {
+        return Point3d.createFromPacked(this.packedData, i, result);
+    }
+    /**
+     * Values of the `order` relevant B-spline basis functions at a parameter.
+     * * Preallocated to length `order` in the constructor and used as a temporary in evaluations.
+     */
+    basisBuffer;
+    /**
+     * Derivatives of the `order` relevant B-spline basis functions at a parameter.
+     * * Preallocated to length `order` in the constructor and used as a temporary in evaluations.
+     */
+    basisBuffer1;
+    /**
+     * Second derivatives of the `order` relevant B-spline basis functions at a parameter.
+     * * Preallocated to length `order` in the constructor and used as a temporary in evaluations.
+     */
+    basisBuffer2;
+    /**
+     * Temporary to hold a single point.
+     * * Preallocated to length `poleLength` in the constructor.
+     */
+    poleBuffer;
+    /**
+     * Temporary to hold a single derivative vector.
+     * * Preallocated to length `poleLength` in the constructor.
+     */
+    poleBuffer1;
+    /**
+     * Temporary to hold a single second derivative vector.
+     * * Preallocated to length `poleLength` in the constructor.
+     */
+    poleBuffer2;
     /**
-     * initialize arrays for given spline dimensions.
-     * @param numPoles number of poles
-     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted, 3 for 2d weighted)
-     * @param order number of poles in support for a section of the bspline
-     * @param knots KnotVector.  This is captured, not cloned.
+     * Initialize arrays for given spline dimensions.
+     * @param numPoles number of poles.
+     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted,
+     * 3 for 2d weighted).
+     * @param order number of poles defining a Bezier segment of the B-spline function.
+     * @param knots the KnotVector. This is captured, not cloned.
      */
     constructor(numPoles, poleLength, order, knots) {
         this.knots = knots;
         this.packedData = new Float64Array(numPoles * poleLength);
         this.poleLength = poleLength;
         this.basisBuffer = new Float64Array(order);
-        this.poleBuffer = new Float64Array(poleLength);
         this.basisBuffer1 = new Float64Array(order);
         this.basisBuffer2 = new Float64Array(order);
+        this.poleBuffer = new Float64Array(poleLength);
         this.poleBuffer1 = new Float64Array(poleLength);
         this.poleBuffer2 = new Float64Array(poleLength);
     }
     /**
-     * create a 1Bspline1dNd`
-     * @param numPoles number of poles
-     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted, 3 for 2d weighted)
-     * @param order number of poles in support for a section of the bspline
-     * @param knots KnotVector.  This is captured, not cloned.
+     * Create a `BSpline1dNd`.
+     * @param numPoles number of poles.
+     * @param poleLength number of coordinates per pole (e.g.. 3 for 3D unweighted, 4 for 3d weighted, 2 for 2d unweighted,
+     * 3 for 2d weighted).
+     * @param order number of poles defining a Bezier segment of the B-spline function.
+     * @param knots the KnotVector. This is captured, not cloned.
      */
     static create(numPoles, poleLength, order, knots) {
         return new BSpline1dNd(numPoles, poleLength, order, knots);
     }
-    /** Map a span index and local fraction to knot value. */
-    spanFractionToKnot(span, localFraction) {
-        return this.knots.spanFractionToKnot(span, localFraction);
+    /** Map a span index and span fraction to knot value. */
+    spanFractionToKnot(spanIndex, spanFraction) {
+        return this.knots.spanFractionToKnot(spanIndex, spanFraction);
     }
     /**
-     * Evaluate the `order` basis functions (and optionally one or two derivatives) at a given fractional position within indexed span.
-     * @returns true if and only if output arrays are sufficiently sized
+     * Evaluate the `order` basis functions (and optionally one or two derivatives) at a given fractional position within
+     * indexed span.
+     * @returns true if and only if output arrays are sufficiently sized.
     */
     evaluateBasisFunctionsInSpan(spanIndex, spanFraction, f, df, ddf) {
         if (spanIndex < 0)
@@ -76,57 +127,64 @@ export class BSpline1dNd {
             this.knots.evaluateBasisFunctions(knotIndex0, globalKnot, f);
     }
     /**
-     * * Evaluate the basis functions at spanIndex and fraction.
-     *   * Evaluations are stored in the preallocated `this.basisBuffer`
-     * * Immediately do the summations of the basis values times the respective control points
-     *   * Summations are stored in the preallocated `this.poleBuffer`
-     * */
-    evaluateBuffersInSpan(spanIndex, spanFraction) {
-        this.evaluateBasisFunctionsInSpan(spanIndex, spanFraction, this.basisBuffer);
-        this.sumPoleBufferForSpan(spanIndex);
-    }
-    /**
-     * * Evaluate the basis functions and one derivative at spanIndex and fraction.
-     *   * Evaluations are stored in the preallocated `this.basisBuffer`
-     * * Immediately do the summations of the basis values times the respective control points
-     *   * Summations are stored in the preallocated `this.poleBuffer` and `this.poleBuffer1`
-     * */
-    evaluateBuffersInSpan1(spanIndex, spanFraction) {
-        this.evaluateBasisFunctionsInSpan(spanIndex, spanFraction, this.basisBuffer, this.basisBuffer1);
-        this.sumPoleBufferForSpan(spanIndex);
-        this.sumPoleBuffer1ForSpan(spanIndex);
-    }
-    /** sum poles in `poleBuffer` at span `spanIndex` by the weights in the `basisBuffer` */
+     * Compute the linear combination of the given span's poles and the weights in `basisBuffer`, and store this point
+     * in `poleBuffer`.
+     */
     sumPoleBufferForSpan(spanIndex) {
         this.poleBuffer.fill(0);
         let k = spanIndex * this.poleLength;
-        for (const f of this.basisBuffer) {
-            for (let j = 0; j < this.poleLength; j++) {
+        for (const f of this.basisBuffer)
+            for (let j = 0; j < this.poleLength; j++)
                 this.poleBuffer[j] += f * this.packedData[k++];
-            }
-        }
     }
-    /** sum poles in `poleBuffer1` at span `spanIndex` by the weights in the `basisBuffer1`, i.e. form first derivatives */
+    /**
+     * Compute the linear combination of the given span's poles and the weights in `basisBuffer1`, and store this
+     * derivative vector in `poleBuffer1`.
+     */
     sumPoleBuffer1ForSpan(spanIndex) {
         this.poleBuffer1.fill(0);
         let k = spanIndex * this.poleLength;
-        for (const f of this.basisBuffer1) {
-            for (let j = 0; j < this.poleLength; j++) {
+        for (const f of this.basisBuffer1)
+            for (let j = 0; j < this.poleLength; j++)
                 this.poleBuffer1[j] += f * this.packedData[k++];
-            }
-        }
     }
-    /** sum poles in `poleBuffer2` at span `spanIndex` by the weights in the `basisBuffer2`, i.e. form second derivatives */
+    /**
+     * Compute the linear combination of the given span's poles and the weights in `basisBuffer2`, and store this
+     * second derivative vector in `poleBuffer2`.
+     */
     sumPoleBuffer2ForSpan(spanIndex) {
         this.poleBuffer2.fill(0);
         let k = spanIndex * this.poleLength;
         for (const f of this.basisBuffer2) {
-            for (let j = 0; j < this.poleLength; j++) {
+            for (let j = 0; j < this.poleLength; j++)
                 this.poleBuffer2[j] += f * this.packedData[k++];
-            }
         }
     }
-    /** Evaluate the function values and 1 or 2 derivatives into `this.poleBuffer`, `this.poleBuffer1` and `this.poleBuffer2` */
+    /**
+     * * Evaluate the basis functions at spanIndex and fraction.
+     *   * Evaluations are stored in the preallocated `this.basisBuffer`.
+     * * Immediately do the summations of the basis values times the respective poles.
+     *   * Summations are stored in the preallocated `this.poleBuffer`
+     * */
+    evaluateBuffersInSpan(spanIndex, spanFraction) {
+        this.evaluateBasisFunctionsInSpan(spanIndex, spanFraction, this.basisBuffer);
+        this.sumPoleBufferForSpan(spanIndex);
+    }
+    /**
+     * * Evaluate the basis functions and one derivative at spanIndex and fraction.
+     *   * Function evaluations are stored in the preallocated `this.basisBuffer`; derivative evaluations in `this.basisBuffer1`.
+     * * Immediately do the summations of the basis values times the respective poles.
+     *   * Summations are stored in the preallocated `this.poleBuffer` and `this.poleBuffer1`
+     * */
+    evaluateBuffersInSpan1(spanIndex, spanFraction) {
+        this.evaluateBasisFunctionsInSpan(spanIndex, spanFraction, this.basisBuffer, this.basisBuffer1);
+        this.sumPoleBufferForSpan(spanIndex);
+        this.sumPoleBuffer1ForSpan(spanIndex);
+    }
+    /**
+     * Evaluate the B-spline function and optional derivatives at the given parameter in knot space.
+     * * Function value is stored in `poleBuffer`; optional derivative vectors in `poleBuffer1` and `poleBuffer2`.
+     */
     evaluateBuffersAtKnot(u, numDerivative = 0) {
         const knotIndex0 = this.knots.knotToLeftKnotIndex(u);
         if (numDerivative < 1) {
@@ -145,36 +203,30 @@ export class BSpline1dNd {
             this.sumPoleBuffer2ForSpan(knotIndex0 - this.degree + 1);
         }
     }
-    /**
-     * Reverse the (blocked) poles (in `this.packedData` in place.
-     */
+    /** Reverse the instance poles and knots in place. */
     reverseInPlace() {
-        // reverse poles in blocks ...
-        const b = this.poleLength;
+        const pLen = this.poleLength;
         const data = this.packedData;
-        for (let i0 = 0, j0 = b * (this.numPoles - 1); i0 < j0; i0 += b, j0 -= b) {
-            let t = 0;
-            for (let i = 0; i < b; i++) {
-                t = data[i0 + i];
-                data[i0 + i] = data[j0 + i];
-                data[j0 + i] = t;
-            }
+        for (let i0 = 0, j0 = pLen * (this.numPoles - 1); i0 < j0; i0 += pLen, j0 -= pLen) {
+            for (let i = 0; i < pLen; i++)
+                [data[i0 + i], data[j0 + i]] = [data[j0 + i], data[i0 + i]];
         }
         this.knots.reflectKnots();
     }
     /**
-     * Test if the leading and trailing polygon coordinates are replicated in the manner of a "closed" bspline polygon which has been expanded
-     * to act as a normal bspline.
-     * @returns true if `degree` leading and trailing polygon blocks match
-     * @deprecated in 4.x. Use testClosablePolygon instead.
+     * Test if the leading and trailing polygon coordinates are replicated in the manner of a "closed" bspline polygon
+     * which has been expanded to act as a normal bspline.
+     * @returns true if `degree` leading and trailing polygon blocks match.
+     * @deprecated in 4.x. Use `testClosablePolygon` instead.
      */
     testCloseablePolygon(mode) {
         return this.testClosablePolygon(mode);
     }
     /**
-     * Test if the leading and trailing polygon coordinates are replicated in the manner of a "closed" bspline polygon which has been expanded
-     * to act as a normal bspline.
-     * @returns true if `degree` leading and trailing polygon blocks match
+     * Test if the leading and trailing poles are replicated in the manner of a "closed" B-spline function with wraparound
+     * control polygon.
+     * @param mode wrap mode, indicating the number of wraparound poles to check. If undefined, `knots.wrappable` is used.
+     * @returns true if the expected leading and trailing poles match, according to `mode`.
      */
     testClosablePolygon(mode) {
         if (mode === undefined)
@@ -196,9 +248,10 @@ export class BSpline1dNd {
         }
         return true;
     }
-    /** Insert knot and resulting pole into the instance, optionally multiple times.
-     * @param knot the knot to be inserted (may already exist in the KnotVector)
-     * @param totalMultiplicity the total multiplicity of the knot on return
+    /**
+     * Insert the knot and resulting pole into the instance, optionally multiple times.
+     * @param knot the knot to be inserted (may already exist in the KnotVector).
+     * @param totalMultiplicity the total multiplicity of the knot on return.
      */
     addKnot(knot, totalMultiplicity) {
         if (knot < this.knots.leftKnot || knot > this.knots.rightKnot)