@itwin/core-geometry 4.0.0-dev.22 → 4.0.0-dev.24

package/lib/cjs/geometry3d/BarycentricTriangle.js

@@ -7,11 +7,100 @@
  * @module CartesianGeometry
  */
 Object.defineProperty(exports, "__esModule", { value: true });
-exports.BarycentricTriangle = void 0;
+exports.BarycentricTriangle = exports.TriangleLocationDetail = void 0;
+const core_bentley_1 = require("@itwin/core-bentley");
 const Geometry_1 = require("../Geometry");
+const Matrix3d_1 = require("./Matrix3d");
 const Point3dVector3d_1 = require("./Point3dVector3d");
+const Ray3d_1 = require("./Ray3d");
 /**
- * 3 points defining a triangle to be evaluated with Barycentric coordinates.
+ * Carries data about a location in the plane of a triangle.
+ * * Each instance carries both world and barycentric coordinates for the point, and provides query services on the latter.
+ * * No tolerance is used when querying barycentric coordinates (e.g., `isInsideOrOn`, `classify`). Use `snapLocalToEdge`
+ * to adjust the barycentric coordinates to a triangle edge if they lie within a parametric tolerance.
+ *
+ * Properties of the barycentric coordinates (b0, b1, b2) of a point p in the plane of a triangle T with vertices p0, p1, p2:
+ * * 1 = b0 + b1 + b2
+ * * p = b0 * p0 + b1 * p1 + b2 * p2
+ * * If T is spanned by the vectors U=p1-p0 and V=p2-p0, then the vector P=p-p0 can be written P = b1 * U + b2 * V.
+ * * The coordinates are all nonnegative if and only if p is inside or on T.
+ * * Exactly one coordinate is zero if and only if p lies on an (infinitely extended) edge of T.
+ * * Exactly two coordinates are zero if and only if p coincides with a vertex of T.
+ * @public
+ */
+class TriangleLocationDetail {
+    constructor() {
+        this.world = new Point3dVector3d_1.Point3d();
+        this.local = new Point3dVector3d_1.Point3d();
+        this.a = 0.0;
+        this.closestEdgeIndex = 0;
+        this.closestEdgeParam = 0.0;
+    }
+    /** Invalidate this detail. */
+    invalidate() {
+        this.world.setZero();
+        this.local.setZero();
+        this.a = 0.0;
+        this.closestEdgeIndex = 0;
+        this.closestEdgeParam = 0.0;
+    }
+    /** Create an invalid detail.
+     * @param result optional pre-allocated object to fill and return
+     */
+    static create(result) {
+        if (undefined === result)
+            result = new TriangleLocationDetail();
+        else
+            result.invalidate();
+        return result;
+    }
+    /** Set the instance contents from the other detail.
+     * @param other detail to clone
+     */
+    copyContentsFrom(other) {
+        this.world.setFrom(other.world);
+        this.local.setFrom(other.local);
+        this.a = other.a;
+        this.closestEdgeIndex = other.closestEdgeIndex;
+        this.closestEdgeParam = other.closestEdgeParam;
+    }
+    /** Whether this detail is invalid. */
+    get isValid() {
+        return !this.local.isZero;
+    }
+    /** Queries the barycentric coordinates to determine whether this instance specifies a location inside or on the triangle.
+     * @see classify
+     */
+    get isInsideOrOn() {
+        return this.isValid && this.local.x >= 0.0 && this.local.y >= 0.0 && this.local.z >= 0.0;
+    }
+    /** Queries this detail to classify the location of this instance with respect to the triangle.
+     * @returns location code
+     * @see isInsideOrOn
+     */
+    get classify() {
+        if (!this.isValid)
+            return Geometry_1.PolygonLocation.Unknown;
+        if (this.isInsideOrOn) {
+            let numZero = 0;
+            if (Math.abs(this.local.x) === 0.0)
+                ++numZero;
+            if (Math.abs(this.local.y) === 0.0)
+                ++numZero;
+            if (Math.abs(this.local.z) === 0.0)
+                ++numZero;
+            if (2 === numZero)
+                return Geometry_1.PolygonLocation.OnPolygonVertex;
+            if (1 === numZero)
+                return Geometry_1.PolygonLocation.OnPolygonEdgeInterior;
+            return Geometry_1.PolygonLocation.InsidePolygonProjectsToEdgeInterior;
+        }
+        return (this.closestEdgeParam === 0.0) ? Geometry_1.PolygonLocation.OutsidePolygonProjectsToVertex : Geometry_1.PolygonLocation.OutsidePolygonProjectsToEdgeInterior;
+    }
+}
+exports.TriangleLocationDetail = TriangleLocationDetail;
+/**
+ * 3 points defining a triangle to be evaluated with barycentric coordinates.
  * @public
  */
 class BarycentricTriangle {
@@ -23,6 +112,10 @@ class BarycentricTriangle {
         this.points.push(point0);
         this.points.push(point1);
         this.points.push(point2);
+        this.edgeLength2 = [];
+        this.edgeLength2.push(point1.distanceSquared(point2));
+        this.edgeLength2.push(point0.distanceSquared(point2));
+        this.edgeLength2.push(point0.distanceSquared(point1));
     }
     /**
      * Return a `BarycentricTriangle` with coordinates given by enumerated x,y,z of the 3 points.
@@ -49,31 +142,331 @@ class BarycentricTriangle {
     }
     /** Return area divided by sum of squared lengths. */
     get aspectRatio() {
-        return Geometry_1.Geometry.safeDivideFraction(0.5 * this.points[0].crossProductToPointsMagnitude(this.points[1], this.points[2]), (this.points[0].distanceSquared(this.points[1]) + this.points[1].distanceSquared(this.points[2]) + this.points[2].distanceSquared(this.points[0])), 0);
+        return Geometry_1.Geometry.safeDivideFraction(this.area, this.edgeLengthSquared(0) + this.edgeLengthSquared(1) + this.edgeLengthSquared(2), 0);
     }
     /** Return the area of the triangle. */
     get area() {
         return 0.5 * this.points[0].crossProductToPointsMagnitude(this.points[1], this.points[2]);
     }
+    /** Return the perimeter of the triangle */
+    get perimeter() {
+        return this.edgeLength(0) + this.edgeLength(1) + this.edgeLength(2);
+    }
     /** Sum the points with given scales.
-     * * In normal use, the scales will add to 1 and the result point is in the plane of the triangle.
-     * * If scales do not add to 1, the point is in the triangle scaled (by the scale sum) from the origin.
+     * * If the scales sum to 1, they are barycentric coordinates, and hence the result point is in the plane of the triangle.
+     * * If the scales do not sum to 1, the point is in the triangle scaled (by the scale sum) from the origin.
+     * @param b0 scale to apply to vertex 0
+     * @param b1 scale to apply to vertex 1
+     * @param b2 scale to apply to vertex 2
+     * @param result optional pre-allocated point to fill and return
+     * @return linear combination of the vertices of this triangle
+     * @see pointToFraction
+     */
+    fractionToPoint(b0, b1, b2, result) {
+        return Point3dVector3d_1.Point3d.createAdd3Scaled(this.points[0], b0, this.points[1], b1, this.points[2], b2, result);
+    }
+    /** Compute the projection of the given point onto the plane of this triangle.
+     * @param point point p to project
+     * @param result optional pre-allocated object to fill and return
+     * @returns details d of the projection point P = `d.point`:
+     * * `d.isValid` returns true if and only if `this.normal()` is defined.
+     * * `d.classify` can be used to determine where P lies with respect to the triangle.
+     * * `d.a` is the signed projection distance: P = p + a * `this.normal()`.
+     * @see fractionToPoint
+     */
+    pointToFraction(point, result) {
+        const normal = BarycentricTriangle._workVector0 = this.normal(BarycentricTriangle._workVector0);
+        if (undefined === normal)
+            return TriangleLocationDetail.create(result);
+        const ray = BarycentricTriangle._workRay = Ray3d_1.Ray3d.create(point, normal, BarycentricTriangle._workRay);
+        return this.intersectRay3d(ray, result); // is free to use workVector0
+    }
+    /** Convert from opposite-vertex to start-vertex edge indexing. */
+    static edgeOppositeVertexIndexToStartVertexIndex(edgeIndex) {
+        return Geometry_1.Geometry.cyclic3dAxis(edgeIndex + 1);
+    }
+    /** Convert from start-vertex to opposite-vertex edge indexing. */
+    static edgeStartVertexIndexToOppositeVertexIndex(startVertexIndex) {
+        return Geometry_1.Geometry.cyclic3dAxis(startVertexIndex - 1);
+    }
+    /** Examine a point's barycentric coordinates to determine if it lies inside the triangle but not on an edge/vertex.
+     * * No parametric tolerance is used.
+     * * It is assumed b0 + b1 + b2 = 1.
+     * @returns whether the point with barycentric coordinates is strictly inside the triangle.
+     */
+    static isInsideTriangle(b0, b1, b2) {
+        return b0 > 0 && b1 > 0 && b2 > 0;
+    }
+    /** Examine a point's barycentric coordinates to determine if it lies inside the triangle or on an edge/vertex.
+     * * No parametric tolerance is used.
+     * * It is assumed b0 + b1 + b2 = 1.
+     * @returns whether the point with barycentric coordinates is inside or on the triangle.
+     */
+    static isInsideOrOnTriangle(b0, b1, b2) {
+        return b0 >= 0 && b1 >= 0 && b2 >= 0;
+    }
+    /** Examine a point's barycentric coordinates to determine if it lies "outside" an edge of the triangle.
+     * * No parametric tolerance is used.
+     * * It is assumed b0 + b1 + b2 = 1.
+     * @returns index of vertex/edge i for which b_i < 0 and b_j >= 0 and b_k >= 0, or -1
      */
-    fractionToPoint(a0, a1, a2, result) {
-        return Point3dVector3d_1.Point3d.createAdd3Scaled(this.points[0], a0, this.points[1], a1, this.points[2], a2, result);
+    static isInRegionBeyondEdge(b0, b1, b2) {
+        if (b0 < 0 && b1 >= 0 && b2 >= 0)
+            return 0;
+        if (b0 >= 0 && b1 < 0 && b2 >= 0)
+            return 1;
+        if (b0 >= 0 && b1 >= 0 && b2 < 0)
+            return 2;
+        return -1;
+    }
+    /** Examine a point's barycentric coordinates to determine if it lies "outside" a vertex of the triangle.
+     * * No parametric tolerance is used.
+     * * It is assumed b0 + b1 + b2 = 1.
+     * @returns index of vertex i for which and b_j < 0 and b_k < 0, or -1
+     */
+    static isInRegionBeyondVertex(b0, b1, b2) {
+        if (b1 < 0 && b2 < 0)
+            return 0;
+        if (b0 < 0 && b2 < 0)
+            return 1;
+        if (b0 < 0 && b1 < 0)
+            return 2;
+        return -1;
+    }
+    /** Examine a point's barycentric coordinates to determine if it lies on a vertex of the triangle.
+     * * No parametric tolerance is used.
+     * * It is assumed b0 + b1 + b2 = 1.
+     * @returns index of vertex i for which b_i = 1 and b_j = b_k = 0, or -1
+     */
+    static isOnVertex(b0, b1, b2) {
+        if (b0 === 1 && b1 === 0 && b2 === 0)
+            return 0;
+        if (b0 === 0 && b1 === 1 && b2 === 0)
+            return 1;
+        if (b0 === 0 && b1 === 0 && b2 === 1)
+            return 2;
+        return -1;
+    }
+    /** Examine a point's barycentric coordinates to determine if it lies on a bounded edge of the triangle.
+     * * No parametric tolerance is used.
+     * * It is assumed b0 + b1 + b2 = 1.
+     * @returns index of vertex/edge i for which b_i = 0 and b_j > 0 and b_k > 0, or -1
+     */
+    static isOnBoundedEdge(b0, b1, b2) {
+        if (b0 === 0 && b1 > 0 && b2 > 0)
+            return 0;
+        if (b0 > 0 && b1 === 0 && b2 > 0)
+            return 1;
+        if (b0 > 0 && b1 > 0 && b2 === 0)
+            return 2;
+        return -1;
+    }
+    /** @returns edge/vertex index (0,1,2) for which the function has a minimum value */
+    static indexOfMinimum(fn) {
+        let i = 0;
+        let min = fn(0);
+        const val = fn(1);
+        if (min > val) {
+            i = 1;
+            min = val;
+        }
+        if (min > fn(2))
+            i = 2;
+        return i;
+    }
+    /** Return the index of the closest triangle vertex to the point given by its barycentric coordinates. */
+    closestVertexIndex(b0, b1, b2) {
+        return BarycentricTriangle.indexOfMinimum((i) => {
+            const a = BarycentricTriangle._workPoint = Point3dVector3d_1.Point3d.createZero(BarycentricTriangle._workPoint);
+            a.setAt(i, 1.0);
+            return this.distanceSquared(a.x, a.y, a.z, b0, b1, b2);
+        });
+    }
+    /** Compute the projection of barycentric point p onto the (unbounded) edge e_k(v_i,v_j) of the triangle T(v_i,v_j,v_k).
+     * @param k vertex v_k is opposite the edge e_k
+     * @param b barycentric coordinates of point to project
+     * @returns parameter s along e_k, such that:
+     * * the projection point is q = v_i + s * (v_j - v_i)
+     * * the barycentric coords of the projection are q_ijk = (1 - s, s, 0)
+     */
+    computeProjectionToEdge(k, b) {
+        // Let U=v_j-v_i and V=v_k-v_i. Then
+        // 0 = (p-q).(v_j-v_i) = ((p-v_i)-(q-v_i)).(v_j-v_i) = (b[j]U + b[k]V - sU).U = b[j]U.U + b[k]U.V - sU.U
+        // Thus s = b[j] + b[k]U.V/U.U
+        k = Geometry_1.Geometry.cyclic3dAxis(k);
+        const i = Geometry_1.Geometry.cyclic3dAxis(k + 1);
+        const j = Geometry_1.Geometry.cyclic3dAxis(i + 1);
+        return b[j] + b[k] * this.dotProductOfEdgeVectorsAtVertex(i) / this.edgeLengthSquared(k);
+    }
+    /** Compute the projection of a barycentric point p to the triangle T(v_0,v_1,v_2).
+     * @param b0 barycentric coordinate of p corresponding to v_0
+     * @param b1 barycentric coordinate of p corresponding to v_1
+     * @param b2 barycentric coordinate of p corresponding to v_2
+     * @returns closest edge start vertex index i and projection parameter u such that the projection q = v_i + u * (v_j - v_i).
+     */
+    closestPoint(b0, b1, b2) {
+        const b = [b0, b1, b2];
+        let edgeIndex = -1; // opposite-vertex index
+        let edgeParam = 0.0;
+        if (BarycentricTriangle.isInsideTriangle(b0, b1, b2)) { // projects to any edge
+            edgeIndex = BarycentricTriangle.indexOfMinimum((i) => {
+                // We want smallest projection distance d_i of p to e_i.
+                // Since b[i]=d_i|e_i|/2A we can compare quantities b[i]/|e_i|.
+                return b[i] * b[i] / this.edgeLengthSquared(i); // avoid sqrt
+            });
+            edgeParam = this.computeProjectionToEdge(edgeIndex, b);
+        }
+        else if ((edgeIndex = BarycentricTriangle.isInRegionBeyondVertex(b0, b1, b2)) >= 0) { // projects to other edges, or any vertex
+            edgeIndex = Geometry_1.Geometry.cyclic3dAxis(edgeIndex + 1);
+            edgeParam = this.computeProjectionToEdge(edgeIndex, b);
+            if (edgeParam < 0 || edgeParam > 1) {
+                edgeIndex = Geometry_1.Geometry.cyclic3dAxis(edgeIndex + 1);
+                edgeParam = this.computeProjectionToEdge(edgeIndex, b);
+                if (edgeParam < 0 || edgeParam > 1) {
+                    edgeParam = 0.0;
+                    edgeIndex = BarycentricTriangle.edgeStartVertexIndexToOppositeVertexIndex(this.closestVertexIndex(b0, b1, b2));
+                }
+            }
+        }
+        else if ((edgeIndex = BarycentricTriangle.isInRegionBeyondEdge(b0, b1, b2)) >= 0) { // projects to the edge or its vertices
+            edgeParam = this.computeProjectionToEdge(edgeIndex, b);
+            if (edgeParam < 0) {
+                edgeParam = 0.0; // start of this edge
+            }
+            else if (edgeParam > 1) {
+                edgeParam = 0.0;
+                edgeIndex = Geometry_1.Geometry.cyclic3dAxis(edgeIndex + 1); // end of this edge = start of next edge
+            }
+        }
+        else if ((edgeIndex = BarycentricTriangle.isOnBoundedEdge(b0, b1, b2)) >= 0) {
+            edgeParam = 1 - b[BarycentricTriangle.edgeOppositeVertexIndexToStartVertexIndex(edgeIndex)];
+        }
+        else if ((edgeIndex = BarycentricTriangle.isOnVertex(b0, b1, b2)) >= 0) {
+            edgeParam = 0.0;
+            edgeIndex = BarycentricTriangle.edgeStartVertexIndexToOppositeVertexIndex(edgeIndex);
+        }
+        // invalid edgeIndex shouldn't happen, but propagate it anyway
+        (0, core_bentley_1.assert)(edgeIndex === 0 || edgeIndex === 1 || edgeIndex === 2);
+        return {
+            closestEdgeIndex: (edgeIndex < 0) ? -1 : BarycentricTriangle.edgeOppositeVertexIndexToStartVertexIndex(edgeIndex),
+            closestEdgeParam: edgeParam,
+        };
+    }
+    /** Compute the intersection of a line (parameterized as a ray) with the plane of this triangle.
+     * @param ray infinite line to intersect, as a ray
+     * @param result optional pre-allocated object to fill and return
+     * @returns details d of the line-plane intersection `d.point`:
+     * * `d.isValid` returns true if and only if the line intersects the plane.
+     * * `d.classify` can be used to determine where the intersection lies with respect to the triangle.
+     * * `d.a` is the intersection parameter. If `d.a` >= 0, the ray intersects the plane of the triangle.
+     * @see pointToFraction
+    */
+    intersectRay3d(ray, result) {
+        result = TriangleLocationDetail.create(result);
+        // Let r0 = ray.origin, d = ray.direction. Write intersection point p two ways for unknown scalars s,b0,b1,b2:
+        //    r0 + s*d = p = b0*p0 + b1*p1 + b2*p2
+        // Subtract p0 from both ends, let u=p1-p0, v=p2-p0, c=r0-p0, and enforce b0+b1+b2=1:
+        //    b1*u + b2*v - s*d = c
+        // This is a linear system Mx=c where M has columns u,v,d and solution x=(b1,b2,-s).
+        const r0 = ray.origin;
+        const d = ray.direction;
+        const u = BarycentricTriangle._workVector0 = Point3dVector3d_1.Vector3d.createStartEnd(this.points[0], this.points[1], BarycentricTriangle._workVector0);
+        const v = BarycentricTriangle._workVector1 = Point3dVector3d_1.Vector3d.createStartEnd(this.points[0], this.points[2], BarycentricTriangle._workVector1);
+        const M = BarycentricTriangle._workMatrix = Matrix3d_1.Matrix3d.createColumns(u, v, d, BarycentricTriangle._workMatrix);
+        const c = Point3dVector3d_1.Vector3d.createStartEnd(this.points[0], r0, BarycentricTriangle._workVector0); // reuse workVector0
+        const solution = BarycentricTriangle._workVector1; // reuse workVector1
+        if (undefined === M.multiplyInverse(c, solution))
+            return result; // invalid
+        result.a = -solution.z;
+        ray.fractionToPoint(result.a, result.world);
+        result.local.set(1.0 - solution.x - solution.y, solution.x, solution.y);
+        const proj = this.closestPoint(result.local.x, result.local.y, result.local.z);
+        result.closestEdgeIndex = proj.closestEdgeIndex;
+        result.closestEdgeParam = proj.closestEdgeParam;
+        return result;
+    }
+    /** Compute the intersection of a line (parameterized as a line segment) with the plane of this triangle.
+     * @param point0 start point of segment on line to intersect
+     * @param point1 end point of segment on line to intersect
+     * @param result optional pre-allocated object to fill and return
+     * @returns details d of the line-plane intersection `d.point`:
+     * * `d.isValid` returns true if and only if the line intersects the plane.
+     * * `d.classify` can be used to determine where the intersection lies with respect to the triangle.
+     * * `d.a` is the intersection parameter. If `d.a` is in [0,1], the segment intersects the plane of the triangle.
+     * @see intersectRay3d
+    */
+    intersectSegment(point0, point1, result) {
+        BarycentricTriangle._workRay = Ray3d_1.Ray3d.createStartEnd(point0, point1, BarycentricTriangle._workRay);
+        return this.intersectRay3d(BarycentricTriangle._workRay, result);
+    }
+    /** Adjust the location to an edge of the triangle if within either given tolerance.
+     * @param location details of a point in the plane of the triangle, `location.local` and `location.world` possibly updated to lie on a triangle edge
+     * @param distanceTolerance absolute distance tolerance (or zero to ignore)
+     * @param parameterTolerance barycentric coordinate fractional tolerance (or zero to ignore)
+     * @return whether the location was adjusted
+     */
+    snapLocationToEdge(location, distanceTolerance = Geometry_1.Geometry.smallMetricDistance, parameterTolerance = Geometry_1.Geometry.smallFraction) {
+        if (!location.isValid)
+            return false;
+        // first try parametric tol to zero barycentric coordinate (no vertices or world distances used!)
+        if (parameterTolerance > 0.0) {
+            let numSnapped = 0;
+            let newSum = 0.0;
+            for (let i = 0; i < 3; i++) {
+                const barycentricDist = Math.abs(location.local.at(i));
+                if (barycentricDist > 0.0 && barycentricDist < parameterTolerance) {
+                    location.local.setAt(i, 0.0);
+                    numSnapped++;
+                }
+                newSum += location.local.at(i);
+            }
+            if (numSnapped > 0 && newSum > 0.0) {
+                location.local.scaleInPlace(1.0 / newSum);
+                if (1 === numSnapped) {
+                    location.closestEdgeIndex = BarycentricTriangle.edgeOppositeVertexIndexToStartVertexIndex(BarycentricTriangle.isOnBoundedEdge(location.local.x, location.local.y, location.local.z));
+                    location.closestEdgeParam = 1.0 - location.local.at(location.closestEdgeIndex);
+                }
+                else { // 2 snapped, at vertex
+                    location.closestEdgeIndex = BarycentricTriangle.isOnVertex(location.local.x, location.local.y, location.local.z);
+                    location.closestEdgeParam = 0.0;
+                }
+                this.fractionToPoint(location.local.x, location.local.y, location.local.z, location.world);
+                return true;
+            }
+        }
+        // failing that, try distance tol to closest edge projection
+        if (distanceTolerance > 0.0) {
+            const i = location.closestEdgeIndex;
+            const j = (i + 1) % 3;
+            const k = (j + 1) % 3;
+            const edgeProjection = BarycentricTriangle._workPoint = this.points[i].interpolate(location.closestEdgeParam, this.points[j], BarycentricTriangle._workPoint);
+            const dist = location.world.distance(edgeProjection);
+            if (dist > 0.0 && dist < distanceTolerance) {
+                location.local.setAt(i, 1.0 - location.closestEdgeParam);
+                location.local.setAt(j, location.closestEdgeParam);
+                location.local.setAt(k, 0.0);
+                location.world.setFrom(edgeProjection);
+                return true;
+            }
+        }
+        return false;
     }
     /** Copy all values from `other`
      */
     setFrom(other) {
-        this.points[0].setFromPoint3d(other.points[0]);
-        this.points[1].setFromPoint3d(other.points[1]);
-        this.points[2].setFromPoint3d(other.points[2]);
+        for (let i = 0; i < 3; ++i) {
+            this.points[i].setFromPoint3d(other.points[i]);
+            this.edgeLength2[i] = other.edgeLength2[i];
+        }
     }
-    /** copy contents of (not pointers to) the given points. */
+    /** Copy contents of (not pointers to) the given points. A vertex is zeroed if its corresponding input point is undefined. */
     set(point0, point1, point2) {
         this.points[0].setFromPoint3d(point0);
         this.points[1].setFromPoint3d(point1);
         this.points[2].setFromPoint3d(point2);
+        this.edgeLength2[0] = this.points[1].distanceSquared(this.points[2]);
+        this.edgeLength2[1] = this.points[0].distanceSquared(this.points[2]);
+        this.edgeLength2[2] = this.points[0].distanceSquared(this.points[1]);
     }
     /**
      * * For `this` and `other` BarycentricTriangles, compute cross products of vectors from point0 to point1 and from point0 to point2.
@@ -90,12 +483,67 @@ class BarycentricTriangle {
         // Do the scale as true division (rather than multiply by precomputed 1/3).  This might protect one bit of result.
         return Point3dVector3d_1.Point3d.create((this.points[0].x + this.points[1].x + this.points[2].x) / 3.0, (this.points[0].y + this.points[1].y + this.points[2].y) / 3.0, (this.points[0].z + this.points[1].z + this.points[2].z) / 3.0, result);
     }
+    /** Return the incenter of the triangle. */
+    incenter(result) {
+        const a = this.edgeLength(0);
+        const b = this.edgeLength(1);
+        const c = this.edgeLength(2);
+        const scale = Geometry_1.Geometry.safeDivideFraction(1.0, a + b + c, 0.0);
+        return this.fractionToPoint(scale * a, scale * b, scale * c, result);
+    }
+    /** Return the circumcenter of the triangle. */
+    circumcenter(result) {
+        const a2 = this.edgeLengthSquared(0);
+        const b2 = this.edgeLengthSquared(1);
+        const c2 = this.edgeLengthSquared(2);
+        const x = a2 * (b2 + c2 - a2);
+        const y = b2 * (c2 + a2 - b2);
+        const z = c2 * (a2 + b2 - c2);
+        const scale = Geometry_1.Geometry.safeDivideFraction(1.0, x + y + z, 0.0);
+        return this.fractionToPoint(scale * x, scale * y, scale * z, result);
+    }
+    /** Return the unit normal of the triangle.
+     * @param result optional pre-allocated vector to fill and return.
+     * @returns unit normal, or undefined if cross product length is too small.
+     */
+    normal(result) {
+        const cross = this.points[0].crossProductToPoints(this.points[1], this.points[2], result);
+        if (cross.tryNormalizeInPlace())
+            return cross;
+        return undefined;
+    }
     /** test for point-by-point `isAlmostEqual` relationship. */
     isAlmostEqual(other) {
         return this.points[0].isAlmostEqual(other.points[0])
             && this.points[1].isAlmostEqual(other.points[1])
             && this.points[2].isAlmostEqual(other.points[2]);
     }
+    /** Compute length of the triangle edge opposite the vertex with the given index.
+     * @see edgeStartVertexIndexToOppositeVertexIndex
+     */
+    edgeLength(oppositeVertexIndex) {
+        return Math.sqrt(this.edgeLengthSquared(oppositeVertexIndex));
+    }
+    /** Compute squared length of the triangle edge opposite the vertex with the given index.
+     * @see edgeStartVertexIndexToOppositeVertexIndex
+    */
+    edgeLengthSquared(oppositeVertexIndex) {
+        return this.edgeLength2[Geometry_1.Geometry.cyclic3dAxis(oppositeVertexIndex)];
+    }
+    /** Compute dot product of the edge vectors based at the vertex with the given index. */
+    dotProductOfEdgeVectorsAtVertex(baseVertexIndex) {
+        const i = Geometry_1.Geometry.cyclic3dAxis(baseVertexIndex);
+        const j = Geometry_1.Geometry.cyclic3dAxis(i + 1);
+        const k = Geometry_1.Geometry.cyclic3dAxis(j + 1);
+        return Geometry_1.Geometry.dotProductXYZXYZ(this.points[j].x - this.points[i].x, this.points[j].y - this.points[i].y, this.points[j].z - this.points[i].z, this.points[k].x - this.points[i].x, this.points[k].y - this.points[i].y, this.points[k].z - this.points[i].z);
+    }
+    /** Compute the squared distance between two points given by their barycentric coordinates.
+     * * It is assumed that a0 + a1 + a2 = b0 + b1 + b2 = 1.
+     */
+    distanceSquared(a0, a1, a2, b0, b1, b2) {
+        // the barycentric displacement vector distance formula
+        return -this.edgeLengthSquared(0) * (b1 - a1) * (b2 - a2) - this.edgeLengthSquared(1) * (b2 - a2) * (b0 - a0) - this.edgeLengthSquared(2) * (b0 - a0) * (b1 - a1);
+    }
 }
 exports.BarycentricTriangle = BarycentricTriangle;
 //# sourceMappingURL=BarycentricTriangle.js.map