@itwin/core-geometry 5.0.0-dev.4 → 5.0.0-dev.40

package/lib/esm/numerics/Polynomials.js

@@ -507,6 +507,8 @@ export class SphereImplicit {
  *
  */
 export class AnalyticRoots {
+    static { this._EQN_EPS = 1.0e-9; }
+    static { this._safeDivideFactor = 1.0e-14; }
     /** Absolute zero test with a tolerance that has worked well for the analytic root use case . . . */
     static isZero(x) {
         return Math.abs(x) < this._EQN_EPS;
@@ -921,11 +923,11 @@ export class AnalyticRoots {
      *   * -1 -- beta, gamma are zero, alpha is not.There is no line defined.There are no solutions.
      *   * 0 -- the line is well defined, but passes completely outside the unit circle.
      *     * In this case, (c1, s1) is the circle point closest to the line and(c2, s2) is the line point closest to the circle.
-     * * 1 -- the line is tangent to the unit circle.
-     *   * Tangency is determined by tolerances, which calls a "close approach" point a tangency.
-     *    * (c1, s1) is the closest circle point
-     *    * (c2, s2) is the line point.
-     * * 2 -- two simple intersections.
+     *   * 1 -- the line is tangent to the unit circle.
+     *     * Tangency is determined by tolerances, which calls a "close approach" point a tangency.
+     *       * (c1, s1) is the closest circle point
+     *       * (c2, s2) is the line point.
+     *   * 2 -- two simple intersections.
      * @param alpha constant coefficient on line
      * @param beta x cosine coefficient on line
      * @param gamma y sine coefficient on line
@@ -979,8 +981,6 @@ export class AnalyticRoots {
         return solutionType;
     }
 }
-AnalyticRoots._EQN_EPS = 1.0e-9;
-AnalyticRoots._safeDivideFactor = 1.0e-14;
 /**
  * Manipulations of polynomials with where `coff[i]` multiplies x^i
  * @internal
@@ -1027,6 +1027,29 @@ export class PowerPolynomial {
  * @internal
  */
 export class TrigPolynomial {
+    // tolerance for small angle decision.
+    static { this._smallAngle = 1.0e-11; }
+    // see core\geometry\internaldocs\unitCircleEllipseIntersection.md for derivation of these coefficients.
+    /** Standard Basis coefficients for the numerator of the y-coordinate y(t) = S(t)/W(t) in the rational semicircle parameterization. */
+    static { this.S = Float64Array.from([0.0, 2.0, -2.0]); }
+    /** Standard Basis coefficients for the numerator of the x-coordinate x(t) = C(t)/W(t) in the rational semicircle parameterization. */
+    static { this.C = Float64Array.from([1.0, -2.0]); }
+    /** Standard Basis coefficients for the denominator of x(t) and y(t) in the rational semicircle parameterization. */
+    static { this.W = Float64Array.from([1.0, -2.0, 2.0]); }
+    /** Standard Basis coefficients for C(t) * W(t). */
+    static { this.CW = Float64Array.from([1.0, -4.0, 6.0, -4.0]); }
+    /** Standard Basis coefficients for S(t) * W(t). */
+    static { this.SW = Float64Array.from([0.0, 2.0, -6.0, 8.0, -4.0]); }
+    /** Standard Basis coefficients for S(t) * C(t). */
+    static { this.SC = Float64Array.from([0.0, 2.0, -6.0, 4.0]); }
+    /** Standard Basis coefficients for S(t) * S(t). */
+    static { this.SS = Float64Array.from([0.0, 0.0, 4.0, -8.0, 4.0]); }
+    /** Standard Basis coefficients for C(t) * C(t). */
+    static { this.CC = Float64Array.from([1.0, -4.0, 4.0]); }
+    /** Standard Basis coefficients for W(t) * W(t). */
+    static { this.WW = Float64Array.from([1.0, -4.0, 8.0, -8.0, 4.0]); }
+    /** Standard Basis coefficients for C(t) * C(t) - S(t) * S(t). */
+    static { this.CCMinusSS = Float64Array.from([1.0, -4.0, 0.0, 8.0, -4.0]); } // eslint-disable-line @typescript-eslint/naming-convention
     /**
      * Find the roots of a univariate polynomial created from substituting the rational parameterization of the unit
      * circle into a trigonometric polynomial. Roots are returned as radian angles.
@@ -1057,7 +1080,7 @@ export class TrigPolynomial {
             degree--;
         const roots = new GrowableFloat64Array();
         if (degree === -1) {
-            // do nothing
+            // p(t) is identically zero (degenerate); do nothing.
         }
         else {
             if (degree === 0) {
@@ -1096,6 +1119,7 @@ export class TrigPolynomial {
         }
         return radians.length > 0;
     }
+    static { this._coefficientRelTol = 1.0e-12; }
     /**
      * Compute intersections of the unit circle `x^2 + y^2 = 1` with the general quadric (conic)
      * `axx x^2 + axy xy + ayy y^2 + ax x + ay y + a = 0`.
@@ -1207,30 +1231,6 @@ export class TrigPolynomial {
         return status;
     }
 }
-// tolerance for small angle decision.
-TrigPolynomial._smallAngle = 1.0e-11;
-// see core\geometry\internaldocs\unitCircleEllipseIntersection.md for derivation of these coefficients.
-/** Standard Basis coefficients for the numerator of the y-coordinate y(t) = S(t)/W(t) in the rational semicircle parameterization. */
-TrigPolynomial.S = Float64Array.from([0.0, 2.0, -2.0]);
-/** Standard Basis coefficients for the numerator of the x-coordinate x(t) = C(t)/W(t) in the rational semicircle parameterization. */
-TrigPolynomial.C = Float64Array.from([1.0, -2.0]);
-/** Standard Basis coefficients for the denominator of x(t) and y(t) in the rational semicircle parameterization. */
-TrigPolynomial.W = Float64Array.from([1.0, -2.0, 2.0]);
-/** Standard Basis coefficients for C(t) * W(t). */
-TrigPolynomial.CW = Float64Array.from([1.0, -4.0, 6.0, -4.0]);
-/** Standard Basis coefficients for S(t) * W(t). */
-TrigPolynomial.SW = Float64Array.from([0.0, 2.0, -6.0, 8.0, -4.0]);
-/** Standard Basis coefficients for S(t) * C(t). */
-TrigPolynomial.SC = Float64Array.from([0.0, 2.0, -6.0, 4.0]);
-/** Standard Basis coefficients for S(t) * S(t). */
-TrigPolynomial.SS = Float64Array.from([0.0, 0.0, 4.0, -8.0, 4.0]);
-/** Standard Basis coefficients for C(t) * C(t). */
-TrigPolynomial.CC = Float64Array.from([1.0, -4.0, 4.0]);
-/** Standard Basis coefficients for W(t) * W(t). */
-TrigPolynomial.WW = Float64Array.from([1.0, -4.0, 8.0, -8.0, 4.0]);
-/** Standard Basis coefficients for C(t) * C(t) - S(t) * S(t). */
-TrigPolynomial.CCMinusSS = Float64Array.from([1.0, -4.0, 0.0, 8.0, -4.0]); // eslint-disable-line @typescript-eslint/naming-convention
-TrigPolynomial._coefficientRelTol = 1.0e-12;
 /**
  * * bilinear expression
  * * `f(u,v) = a + b * u * c * v + d * u * v`