@pirireis/webglobeplugins 0.9.6 → 0.9.8

package/{Math/frustum/from-globeinfo.ts → dist/Math/frustum/from-globeinfo.js}

@@ -1,63 +1,48 @@
-// import { Plane } from "../plane";
-// import { Vector3D } from "../vector3d";
-// import { WORLD_RADIUS_3D } from "../constants";
-// import { Quaternion } from "../quaternion";
-
-// import { FrustumPlanes } from "./types";
-
-// const _quaternion = new Quaternion();
-
-// export function getFrustum(globe: any, fieldOfView: number = 50 * Math.PI / 180,
-//     out: FrustumPlanes): void {
-
-//     const lookinfo = globe.api_GetCurrentLookInfo();
-//     lookinfo.CenterLong *= Math.PI / 180; // convert degrees to radians
-//     lookinfo.CenterLat *= Math.PI / 180; // convert degrees to radians
-//     lookinfo.NorthAng *= Math.PI / 180; // convert degrees to radians
-
-//     const lookAtPosition = Vector3D.fromLongLatRadians(lookinfo.CenterLong, lookinfo.CenterLat);
-//     const cameraPosition = new Vector3D(globe.Fp.x, globe.Fp.y, globe.Fp.z).divideByScaler(WORLD_RADIUS_3D);
-//     const lookFromCamera = lookAtPosition.clone().subtract(cameraPosition).normalize();
-
-
-//     const cameraZ = new Vector3D(globe.Fu.x, globe.Fu.y, globe.Fu.z);
-//     const cameraY = cameraZ.clone().cross(lookFromCamera).normalize();
-
-//     const bottomNormal = cameraZ.clone().applyQuaternion(
-//         _quaternion.setFromAxisAngle(cameraY, -fieldOfView)
-//     ).normalize();
-
-//     const topNormal = cameraZ.clone().applyQuaternion(
-//         _quaternion.setFromAxisAngle(cameraY, fieldOfView)
-//     ).negate().normalize();
-//     const leftNormal = cameraY.clone().applyQuaternion(
-//         _quaternion.setFromAxisAngle(cameraZ, fieldOfView)
-//     ).negate().normalize();
-
-//     const rightNormal = cameraY.clone().applyQuaternion(
-//         _quaternion.setFromAxisAngle(cameraZ, -fieldOfView)
-//     ).normalize();
-
-
-//     out.top.set(topNormal, topNormal.dot(cameraPosition));
-//     out.bottom.set(bottomNormal, bottomNormal.dot(cameraPosition));
-//     out.left.set(leftNormal, leftNormal.dot(cameraPosition));
-//     out.right.set(rightNormal, rightNormal.dot(cameraPosition));
-//     out.far.set(cameraPosition.normalize(), 0.01);
-//     out.near = null; // Near plane is not calculated in this function, set to null
-//     vectorInfo(globe, [out]);
-
-// }
-
-// let interval = 0;
-
-// const vectorInfo = (globe, args: null | Array<any>) => {
-//     interval++;
-//     if (interval % 100 === 0) {
-//         console.log("Fp:", globe.Fp.x / WORLD_RADIUS_3D, globe.Fp.y / WORLD_RADIUS_3D, globe.Fp.z / WORLD_RADIUS_3D);
-//         console.log("Fu:", globe.Fu.x, globe.Fu.y, globe.Fu.z);
-
-//         if (args) console.log(...args)
-//         interval = 0;
-//     }
-// }
+"use strict";
+// import { Plane } from "../plane";
+// import { Vector3D } from "../vector3d";
+// import { WORLD_RADIUS_3D } from "../constants";
+// import { Quaternion } from "../quaternion";
+// import { FrustumPlanes } from "./types";
+// const _quaternion = new Quaternion();
+// export function getFrustum(globe: any, fieldOfView: number = 50 * Math.PI / 180,
+//     out: FrustumPlanes): void {
+//     const lookinfo = globe.api_GetCurrentLookInfo();
+//     lookinfo.CenterLong *= Math.PI / 180; // convert degrees to radians
+//     lookinfo.CenterLat *= Math.PI / 180; // convert degrees to radians
+//     lookinfo.NorthAng *= Math.PI / 180; // convert degrees to radians
+//     const lookAtPosition = Vector3D.fromLongLatRadians(lookinfo.CenterLong, lookinfo.CenterLat);
+//     const cameraPosition = new Vector3D(globe.Fp.x, globe.Fp.y, globe.Fp.z).divideByScaler(WORLD_RADIUS_3D);
+//     const lookFromCamera = lookAtPosition.clone().subtract(cameraPosition).normalize();
+//     const cameraZ = new Vector3D(globe.Fu.x, globe.Fu.y, globe.Fu.z);
+//     const cameraY = cameraZ.clone().cross(lookFromCamera).normalize();
+//     const bottomNormal = cameraZ.clone().applyQuaternion(
+//         _quaternion.setFromAxisAngle(cameraY, -fieldOfView)
+//     ).normalize();
+//     const topNormal = cameraZ.clone().applyQuaternion(
+//         _quaternion.setFromAxisAngle(cameraY, fieldOfView)
+//     ).negate().normalize();
+//     const leftNormal = cameraY.clone().applyQuaternion(
+//         _quaternion.setFromAxisAngle(cameraZ, fieldOfView)
+//     ).negate().normalize();
+//     const rightNormal = cameraY.clone().applyQuaternion(
+//         _quaternion.setFromAxisAngle(cameraZ, -fieldOfView)
+//     ).normalize();
+//     out.top.set(topNormal, topNormal.dot(cameraPosition));
+//     out.bottom.set(bottomNormal, bottomNormal.dot(cameraPosition));
+//     out.left.set(leftNormal, leftNormal.dot(cameraPosition));
+//     out.right.set(rightNormal, rightNormal.dot(cameraPosition));
+//     out.far.set(cameraPosition.normalize(), 0.01);
+//     out.near = null; // Near plane is not calculated in this function, set to null
+//     vectorInfo(globe, [out]);
+// }
+// let interval = 0;
+// const vectorInfo = (globe, args: null | Array<any>) => {
+//     interval++;
+//     if (interval % 100 === 0) {
+//         console.log("Fp:", globe.Fp.x / WORLD_RADIUS_3D, globe.Fp.y / WORLD_RADIUS_3D, globe.Fp.z / WORLD_RADIUS_3D);
+//         console.log("Fu:", globe.Fu.x, globe.Fu.y, globe.Fu.z);
+//         if (args) console.log(...args)
+//         interval = 0;
+//     }
+// }

package/dist/Math/frustum/types.js

@@ -0,0 +1,2 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });

package/{Math/globe-util/horizon-plane.ts → dist/Math/globe-util/horizon-plane.js}

@@ -1,137 +1,112 @@
-// TODO: REMOVE
-
-// // import { Vector3D, Quaternion, Plane, Ray, Matrix4 } from 'your-3d-math-library'; // <-- ADAPT THIS IMPORT
-
-// import { Vector3D } from '../vector3d';
-// import { Plane } from '../plane';
-// import { Line } from '../line';
-// import { Matrix4 } from '../matrix4';
-// import { Quaternion } from '../quaternion';
-
-// // --- Assume these static helper variables exist within your class ---
-// const _lookFromOrigin = new Vector3D();
-// const _cameraPosition = new Vector3D();
-// const _finalForward = new Vector3D();
-// const _approxUp = new Vector3D();
-// const _finalRight = new Vector3D();
-// const _finalUp = new Vector3D();
-// const _tempVec = new Vector3D(); // General purpose helper
-// const _initialBottomRayDir = new Vector3D();
-// const _qOrient = new Quaternion();
-// const _rollQuaternion = new Quaternion();
-// const _rotationMatrix = new Matrix4();
-// const _fieldOfViewBottomRay = new Line(new Vector3D(), new Vector3D());
-// const _origin3d = new Vector3D(0, 0, 0); // Globe Center
-
-// // --- Assume these constants exist ---
-// const WORLD_RADIUS_3D = 6371000; // Example: Earth radius
-// const FIELD_OF_VIEW = (25 * Math.PI) / 180; // Example: 25 degrees half-FOV in radians
-
-// // // --- Assume Plane class has these methods ---
-// // interface Plane {
-// //     setFromNormalAndPoint(normal: Vector3D, point: Vector3D): this;
-// //     set(normal: Vector3D, constant: number): this;
-// //     // ... other plane methods
-// // }
-
-// // // --- Assume Ray class has these methods ---
-// // interface Ray {
-// //     origin: Vector3D;
-// //     direction: Vector3D;
-// //     intersectionSphere(center: Vector3D, radius: number): Vector3D[] | null; // Or similar
-// //     // ... other ray methods
-// // }
-
-// // // --- Assume Vector3D has setFromLongLat ---
-// // interface Vector3D {
-// //     setFromLongLat(lon: number, lat: number): this;
-// //     // ... other vector methods
-// // }
-// // // --- END OF ASSUMPTIONS ---
-
-
-// /**
-//  * Calculates a culling plane representing the bottom edge of the camera's field of view.
-//  * @param globe - The globe object, expected to have api_GetCurrentLookInfo and Fp properties.
-//  * @param out - The Plane object to store the result.
-//  * @param fieldOfView - The vertical half-angle of the field of view in radians.
-//  * @returns The calculated Plane.
-//  */
-// function calculateHorizonPlane(globe: any, out: Plane, fieldOfView: number = FIELD_OF_VIEW): Plane {
-
-//     const cameraLookInfo = globe.api_GetCurrentLookInfo();
-//     cameraLookInfo.CenterLong *= Math.PI / 180;
-//     cameraLookInfo.CenterLat *= Math.PI / 180;
-//     // cameraLookInfo.Tilt *= Math.PI / 180; // Tilt is now implicitly handled
-//     cameraLookInfo.NorthAng *= Math.PI / 180;
-
-//     // Set _lookFromOrigin based on camera target (Lon/Lat)
-//     (_lookFromOrigin as any).setFromLongLat(cameraLookInfo.CenterLong, cameraLookInfo.CenterLat); // Cast to any if setFromLongLat is custom
-
-//     // 1. Calculate _cameraPosition and _finalForward (_lookFromCamera)
-//     _cameraPosition.set(globe.Fp.x, globe.Fp.y, globe.Fp.z); // Using your X negation
-//     _cameraPosition.divideScalar(WORLD_RADIUS_3D); // Scale to unit sphere
-
-//     Vector3D.subVectors(_lookFromOrigin, _cameraPosition, _finalForward).normalize();
-
-//     // 2. Calculate Q_orient (Final Camera Orientation Quaternion)
-
-//     // Use the vector from origin to camera as an initial 'up' reference
-//     _approxUp.copy(_cameraPosition).normalize();
-
-//     // Calculate Right vector using cross product
-//     Vector3D.crossVectors(_finalForward, _approxUp, _finalRight).normalize();
-
-//     // Handle edge case: If looking straight along (or opposite to) _approxUp
-//     if (_finalRight.lengthSq() < 0.0001) {
-//         // If vectors are parallel, _approxUp is not good.
-//         // Choose a different 'up', e.g., if _approxUp is near Z, use Y.
-//         if (Math.abs(_approxUp.z) > 0.9) {
-//             _approxUp.set(0, 1, 0);
-//         } else {
-//             _approxUp.set(0, 0, 1);
-//         }
-//         Vector3D.crossVectors(_finalForward, _approxUp, _finalRight).normalize();
-//     }
-
-//     // Calculate the 'true' Up vector, orthogonal to Right and Forward
-//     Vector3D.crossVectors(_finalRight, _finalForward, _finalUp).normalize();
-
-//     // Apply NorthAng (Roll) around the Forward axis
-//     _rollQuaternion.setFromAxisAngle(_finalForward, cameraLookInfo.NorthAng);
-//     _finalUp.applyQuaternion(_rollQuaternion);
-//     _finalRight.applyQuaternion(_rollQuaternion); // Must rotate Right as well!
-
-//     // Create Q_orient from the final basis vectors (Right, Up, -Forward)
-//     // This assumes a Matrix4.makeBasis(x, y, z) exists and builds a rotation matrix.
-//     // It also assumes local camera looks down -Z, hence negating _finalForward for the Z basis.
-//     _tempVec.copy(_finalForward).negate();
-//     _rotationMatrix.makeBasis(_finalRight, _finalUp, _tempVec); // <-- ADAPT: Ensure your library has makeBasis or equivalent
-//     _qOrient.setFromRotationMatrix(_rotationMatrix); // <-- ADAPT: Ensure this function exists
-
-//     // 3. Calculate Initial Bottom Ray (Local: -Z forward, Y up, X right)
-//     // Rotates (0, 0, -1) around (1, 0, 0) by fieldOfView
-//     _initialBottomRayDir.set(0, -Math.sin(fieldOfView), -Math.cos(fieldOfView));
-
-//     // 4. Calculate Final Bottom Ray (World)
-//     _fieldOfViewBottomRay.direction.copy(_initialBottomRayDir).applyQuaternion(_qOrient).normalize();
-//     _fieldOfViewBottomRay.origin.copy(_cameraPosition);
-
-//     // 5. Intersection and Plane
-//     const intersection = (_fieldOfViewBottomRay as any).intersectionSphere(_origin3d, 1); // Cast to any if method isn't standard
-
-//     if (intersection && intersection.length > 0) {
-//         // If ray intersects, use the closest point for the plane
-//         out.setFromNormalAndPoint(_lookFromOrigin, intersection[0]);
-//     } else {
-//         // If ray does not intersect (FOV covers globe or misses)
-//         // WARNING: This sets a plane through the globe's center.
-//         // This might not be the desired behavior for culling.
-//         // Consider what should happen if the whole globe *is* visible.
-//         out.set(_lookFromOrigin, 0.0);
-//     }
-
-//     return out;
-// }
-
-// export { calculateHorizonPlane };
+"use strict";
+// TODO: REMOVE
+// // import { Vector3D, Quaternion, Plane, Ray, Matrix4 } from 'your-3d-math-library'; // <-- ADAPT THIS IMPORT
+// import { Vector3D } from '../vector3d';
+// import { Plane } from '../plane';
+// import { Line } from '../line';
+// import { Matrix4 } from '../matrix4';
+// import { Quaternion } from '../quaternion';
+// // --- Assume these static helper variables exist within your class ---
+// const _lookFromOrigin = new Vector3D();
+// const _cameraPosition = new Vector3D();
+// const _finalForward = new Vector3D();
+// const _approxUp = new Vector3D();
+// const _finalRight = new Vector3D();
+// const _finalUp = new Vector3D();
+// const _tempVec = new Vector3D(); // General purpose helper
+// const _initialBottomRayDir = new Vector3D();
+// const _qOrient = new Quaternion();
+// const _rollQuaternion = new Quaternion();
+// const _rotationMatrix = new Matrix4();
+// const _fieldOfViewBottomRay = new Line(new Vector3D(), new Vector3D());
+// const _origin3d = new Vector3D(0, 0, 0); // Globe Center
+// // --- Assume these constants exist ---
+// const WORLD_RADIUS_3D = 6371000; // Example: Earth radius
+// const FIELD_OF_VIEW = (25 * Math.PI) / 180; // Example: 25 degrees half-FOV in radians
+// // // --- Assume Plane class has these methods ---
+// // interface Plane {
+// //     setFromNormalAndPoint(normal: Vector3D, point: Vector3D): this;
+// //     set(normal: Vector3D, constant: number): this;
+// //     // ... other plane methods
+// // }
+// // // --- Assume Ray class has these methods ---
+// // interface Ray {
+// //     origin: Vector3D;
+// //     direction: Vector3D;
+// //     intersectionSphere(center: Vector3D, radius: number): Vector3D[] | null; // Or similar
+// //     // ... other ray methods
+// // }
+// // // --- Assume Vector3D has setFromLongLat ---
+// // interface Vector3D {
+// //     setFromLongLat(lon: number, lat: number): this;
+// //     // ... other vector methods
+// // }
+// // // --- END OF ASSUMPTIONS ---
+// /**
+//  * Calculates a culling plane representing the bottom edge of the camera's field of view.
+//  * @param globe - The globe object, expected to have api_GetCurrentLookInfo and Fp properties.
+//  * @param out - The Plane object to store the result.
+//  * @param fieldOfView - The vertical half-angle of the field of view in radians.
+//  * @returns The calculated Plane.
+//  */
+// function calculateHorizonPlane(globe: any, out: Plane, fieldOfView: number = FIELD_OF_VIEW): Plane {
+//     const cameraLookInfo = globe.api_GetCurrentLookInfo();
+//     cameraLookInfo.CenterLong *= Math.PI / 180;
+//     cameraLookInfo.CenterLat *= Math.PI / 180;
+//     // cameraLookInfo.Tilt *= Math.PI / 180; // Tilt is now implicitly handled
+//     cameraLookInfo.NorthAng *= Math.PI / 180;
+//     // Set _lookFromOrigin based on camera target (Lon/Lat)
+//     (_lookFromOrigin as any).setFromLongLat(cameraLookInfo.CenterLong, cameraLookInfo.CenterLat); // Cast to any if setFromLongLat is custom
+//     // 1. Calculate _cameraPosition and _finalForward (_lookFromCamera)
+//     _cameraPosition.set(globe.Fp.x, globe.Fp.y, globe.Fp.z); // Using your X negation
+//     _cameraPosition.divideScalar(WORLD_RADIUS_3D); // Scale to unit sphere
+//     Vector3D.subVectors(_lookFromOrigin, _cameraPosition, _finalForward).normalize();
+//     // 2. Calculate Q_orient (Final Camera Orientation Quaternion)
+//     // Use the vector from origin to camera as an initial 'up' reference
+//     _approxUp.copy(_cameraPosition).normalize();
+//     // Calculate Right vector using cross product
+//     Vector3D.crossVectors(_finalForward, _approxUp, _finalRight).normalize();
+//     // Handle edge case: If looking straight along (or opposite to) _approxUp
+//     if (_finalRight.lengthSq() < 0.0001) {
+//         // If vectors are parallel, _approxUp is not good.
+//         // Choose a different 'up', e.g., if _approxUp is near Z, use Y.
+//         if (Math.abs(_approxUp.z) > 0.9) {
+//             _approxUp.set(0, 1, 0);
+//         } else {
+//             _approxUp.set(0, 0, 1);
+//         }
+//         Vector3D.crossVectors(_finalForward, _approxUp, _finalRight).normalize();
+//     }
+//     // Calculate the 'true' Up vector, orthogonal to Right and Forward
+//     Vector3D.crossVectors(_finalRight, _finalForward, _finalUp).normalize();
+//     // Apply NorthAng (Roll) around the Forward axis
+//     _rollQuaternion.setFromAxisAngle(_finalForward, cameraLookInfo.NorthAng);
+//     _finalUp.applyQuaternion(_rollQuaternion);
+//     _finalRight.applyQuaternion(_rollQuaternion); // Must rotate Right as well!
+//     // Create Q_orient from the final basis vectors (Right, Up, -Forward)
+//     // This assumes a Matrix4.makeBasis(x, y, z) exists and builds a rotation matrix.
+//     // It also assumes local camera looks down -Z, hence negating _finalForward for the Z basis.
+//     _tempVec.copy(_finalForward).negate();
+//     _rotationMatrix.makeBasis(_finalRight, _finalUp, _tempVec); // <-- ADAPT: Ensure your library has makeBasis or equivalent
+//     _qOrient.setFromRotationMatrix(_rotationMatrix); // <-- ADAPT: Ensure this function exists
+//     // 3. Calculate Initial Bottom Ray (Local: -Z forward, Y up, X right)
+//     // Rotates (0, 0, -1) around (1, 0, 0) by fieldOfView
+//     _initialBottomRayDir.set(0, -Math.sin(fieldOfView), -Math.cos(fieldOfView));
+//     // 4. Calculate Final Bottom Ray (World)
+//     _fieldOfViewBottomRay.direction.copy(_initialBottomRayDir).applyQuaternion(_qOrient).normalize();
+//     _fieldOfViewBottomRay.origin.copy(_cameraPosition);
+//     // 5. Intersection and Plane
+//     const intersection = (_fieldOfViewBottomRay as any).intersectionSphere(_origin3d, 1); // Cast to any if method isn't standard
+//     if (intersection && intersection.length > 0) {
+//         // If ray intersects, use the closest point for the plane
+//         out.setFromNormalAndPoint(_lookFromOrigin, intersection[0]);
+//     } else {
+//         // If ray does not intersect (FOV covers globe or misses)
+//         // WARNING: This sets a plane through the globe's center.
+//         // This might not be the desired behavior for culling.
+//         // Consider what should happen if the whole globe *is* visible.
+//         out.set(_lookFromOrigin, 0.0);
+//     }
+//     return out;
+// }
+// export { calculateHorizonPlane };

package/dist/Math/index.js

@@ -0,0 +1 @@
+"use strict";

package/dist/Math/juction/arc-plane.js

@@ -0,0 +1,83 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.arcSlice = arcSlice;
+const constants_1 = require("../constants");
+const vec3_1 = require("../vec3");
+const line_1 = require("../line");
+const plane_1 = require("../plane");
+const arc_1 = require("../arc");
+const plane_plane_1 = require("./plane-plane");
+const line_sphere_1 = require("./line-sphere");
+const _intersectionLine = /*@__PURE__*/ line_1.line.create();
+const _originPlane = /*@__PURE__*/ plane_1.plane.create();
+_originPlane.distance = 0;
+const _originSphere = /*@__PURE__*/ { center: vec3_1.vec3.create(0, 0, 0), radius: 1 };
+const _intersectionPoints = /*@__PURE__*/ [vec3_1.vec3.create(), vec3_1.vec3.create()];
+// TODO: out must be [Arc, Arc] there is a case where split into three arcs, visible by points or arcs but middle is not visible
+/**
+ *
+ * @param out \
+ * @param inArc
+ * @param juctionPlane
+ * @returns number of arcs segments in the junction divided by the junction plane.
+ */
+function arcSlice(out, inArc, juctionPlane) {
+    // arc coverPlane and junctionPlane intersection exist in the range of unit sphere
+    const coverRadiusAngle = plane_1.plane.getUnitSphereRadiusAngle(inArc.coverPlane);
+    const visibleRadiusAngle = plane_1.plane.getUnitSphereRadiusAngle(juctionPlane);
+    // TODO: RESEARCH efficient approach --
+    const angleBetween_Cover_Visible = Math.acos(vec3_1.vec3.dot(inArc.coverPlane.normal, juctionPlane.normal));
+    if (coverRadiusAngle + visibleRadiusAngle < angleBetween_Cover_Visible - constants_1.EPSILON) { // case A: out of range
+        return 0; // No intersection
+    }
+    // ------------------------------------
+    // the case when the arc is completely covered by the juction plane
+    if (visibleRadiusAngle + constants_1.EPSILON >= angleBetween_Cover_Visible + coverRadiusAngle) { // case B: fully visible
+        arc_1.arc.copy(out[0], inArc);
+        return 1;
+    }
+    // plane-plane intersection line should be calculated for the rest of the calculations
+    vec3_1.vec3.copy(_originPlane.normal, inArc.normal);
+    const isPlaneJunctions = (0, plane_plane_1.planePlaneJuction)(_intersectionLine, _originPlane, juctionPlane);
+    if (!isPlaneJunctions) { // case C: planes are parallel. 
+        // case A covers opposite directions
+        // case B Covers they face each other
+        throw new Error("Unexpected case: planes are parallel, case A and B should Cover it");
+    }
+    // --- read until here ---
+    // calculate the intersection points
+    const isSphereIntersection = (0, line_sphere_1.lineSphereIntersection)(_intersectionPoints, _intersectionLine, _originSphere);
+    if (!isSphereIntersection) {
+        // other edge caes should be covered by now
+        return 0; // No intersection
+    }
+    const i0IsCovered = plane_1.plane.distanceToPoint(inArc.coverPlane, _intersectionPoints[0]) > -constants_1.EPSILON;
+    const i1IsCovered = plane_1.plane.distanceToPoint(inArc.coverPlane, _intersectionPoints[1]) > -constants_1.EPSILON;
+    const p0IsVisible = plane_1.plane.distanceToPoint(juctionPlane, inArc.p0) > -constants_1.EPSILON;
+    const p1IsVisible = plane_1.plane.distanceToPoint(juctionPlane, inArc.p1) > -constants_1.EPSILON;
+    if (!p0IsVisible && !p1IsVisible && !i0IsCovered && !i1IsCovered) {
+        return 0; // No intersection
+    }
+    if (i0IsCovered && i1IsCovered && p0IsVisible && p1IsVisible) {
+        // calculate which points are closer
+        const p0i0DistanceSquared = vec3_1.vec3.distanceSquared(inArc.p0, _intersectionPoints[0]);
+        const p0i1DistanceSquared = vec3_1.vec3.distanceSquared(inArc.p0, _intersectionPoints[1]);
+        const case0 = p0i0DistanceSquared < p0i1DistanceSquared;
+        arc_1.arc.set(out[0], inArc.p0, case0 ? _intersectionPoints[0] : _intersectionPoints[1]);
+        arc_1.arc.set(out[1], inArc.p1, !case0 ? _intersectionPoints[0] : _intersectionPoints[1]);
+        return 2;
+    }
+    if (i0IsCovered && i1IsCovered) {
+        arc_1.arc.set(out[0], _intersectionPoints[0], _intersectionPoints[1]);
+        return 1;
+    }
+    if (p0IsVisible && p1IsVisible) {
+        arc_1.arc.copy(out[0], inArc);
+        return 1;
+    }
+    if ((p0IsVisible || p1IsVisible) !== (i0IsCovered || i1IsCovered)) {
+        throw new Error("Unexpected case: one covered and one visible point must be present");
+    }
+    arc_1.arc.set(out[0], p0IsVisible ? inArc.p0 : inArc.p1, i0IsCovered ? _intersectionPoints[0] : _intersectionPoints[1]);
+    return 1;
+}

package/dist/Math/juction/line-sphere.js

@@ -0,0 +1,25 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.lineSphereIntersection = lineSphereIntersection;
+const vec3_1 = require("../vec3");
+const line_1 = require("../line");
+const _0vector = /*@__PURE__*/ vec3_1.vec3.create(0, 0, 0);
+function lineSphereIntersection(out, inLine, inSphere) {
+    vec3_1.vec3.subtract(_0vector, inLine.origin, inSphere.center);
+    const distanceSquared = vec3_1.vec3.lengthSquared(_0vector);
+    const radiusSquared = inSphere.radius * inSphere.radius;
+    const dot = vec3_1.vec3.dot(_0vector, inLine.direction);
+    const dotSquared = dot * dot;
+    const _a = dotSquared + radiusSquared - distanceSquared;
+    if (_a < 0) {
+        return false; // no intersection
+    }
+    else {
+        const a = Math.sqrt(_a);
+        const t1 = dot - a;
+        const t2 = dot + a;
+        line_1.line.at(out[0], inLine, t1);
+        line_1.line.at(out[1], inLine, t2);
+        return true;
+    }
+}

package/{Math/juction/plane-plane.ts → dist/Math/juction/plane-plane.js}

@@ -1,66 +1,58 @@
-import { Plane, Line, Vec3 } from "../types";
-import { EPSILON } from "../constants";
-import { plane } from "../plane";
-import { line } from "../line";
-import { vec3 } from "../vec3";
-
-const _normal1 = vec3.create();
-const _normal2 = vec3.create();
-
-export function planePlaneJuction(out: Line, plane1: Plane, plane2: Plane): boolean {
-    vec3.copy(_normal1, plane1.normal);
-    vec3.copy(_normal2, plane2.normal);
-    const distance1 = plane1.distance;
-    const distance2 = plane2.distance;
-
-    const { origin, direction } = out;
-
-    const dot = vec3.dot(_normal1, _normal2);
-
-    if (Math.abs(dot) > 1 - EPSILON) {
-        return false; // Planes are parallel, no intersection
-    }
-
-    vec3.cross(direction, _normal1, _normal2);
-
-    const magnitudeSquired = vec3.lengthSquared(out.direction);
-    if (magnitudeSquired < EPSILON) {
-        return false; // No valid intersection line
-    }
-
-    const magnitude = Math.sqrt(magnitudeSquired);
-    direction[0] /= magnitude;
-    direction[1] /= magnitude;
-    direction[2] /= magnitude;
-
-
-
-    let determinant: number;
-    // Calculate the intersection point
-    // set z = 0
-    determinant = _normal1[0] * _normal2[1] - _normal1[1] * _normal2[0];
-    if (Math.abs(determinant) > EPSILON) {
-        origin[0] = (distance1 * _normal2[1] - distance2 * _normal1[1]) / determinant;
-        origin[1] = (distance2 * _normal1[0] - distance1 * _normal2[0]) / determinant;
-        origin[2] = 0;
-    } else {
-        // set y = 0
-        determinant = _normal1[0] * _normal2[2] - _normal1[2] * _normal2[0];
-        if (Math.abs(determinant) > EPSILON) {
-            origin[0] = (distance1 * _normal2[2] - distance2 * _normal1[2]) / determinant;
-            origin[1] = 0;
-            origin[2] = (distance2 * _normal1[0] - distance1 * _normal2[0]) / determinant;
-        } else {
-            // set x = 0
-            determinant = _normal1[1] * _normal2[2] - _normal1[2] * _normal2[1];
-            if (Math.abs(determinant) > EPSILON) {
-                origin[0] = 0;
-                origin[1] = (distance1 * _normal2[2] - distance2 * _normal1[2]) / determinant;
-                origin[2] = (distance2 * _normal1[1] - distance1 * _normal2[1]) / determinant;
-            } else {
-                return false; // No valid intersection point
-            }
-        }
-    }
-    return true;
-}
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.planePlaneJuction = planePlaneJuction;
+const constants_1 = require("../constants");
+const vec3_1 = require("../vec3");
+const _normal1 = vec3_1.vec3.create();
+const _normal2 = vec3_1.vec3.create();
+function planePlaneJuction(out, plane1, plane2) {
+    vec3_1.vec3.copy(_normal1, plane1.normal);
+    vec3_1.vec3.copy(_normal2, plane2.normal);
+    const distance1 = plane1.distance;
+    const distance2 = plane2.distance;
+    const { origin, direction } = out;
+    const dot = vec3_1.vec3.dot(_normal1, _normal2);
+    if (Math.abs(dot) > 1 - constants_1.EPSILON) {
+        return false; // Planes are parallel, no intersection
+    }
+    vec3_1.vec3.cross(direction, _normal1, _normal2);
+    const magnitudeSquired = vec3_1.vec3.lengthSquared(out.direction);
+    if (magnitudeSquired < constants_1.EPSILON) {
+        return false; // No valid intersection line
+    }
+    const magnitude = Math.sqrt(magnitudeSquired);
+    direction[0] /= magnitude;
+    direction[1] /= magnitude;
+    direction[2] /= magnitude;
+    let determinant;
+    // Calculate the intersection point
+    // set z = 0
+    determinant = _normal1[0] * _normal2[1] - _normal1[1] * _normal2[0];
+    if (Math.abs(determinant) > constants_1.EPSILON) {
+        origin[0] = (distance1 * _normal2[1] - distance2 * _normal1[1]) / determinant;
+        origin[1] = (distance2 * _normal1[0] - distance1 * _normal2[0]) / determinant;
+        origin[2] = 0;
+    }
+    else {
+        // set y = 0
+        determinant = _normal1[0] * _normal2[2] - _normal1[2] * _normal2[0];
+        if (Math.abs(determinant) > constants_1.EPSILON) {
+            origin[0] = (distance1 * _normal2[2] - distance2 * _normal1[2]) / determinant;
+            origin[1] = 0;
+            origin[2] = (distance2 * _normal1[0] - distance1 * _normal2[0]) / determinant;
+        }
+        else {
+            // set x = 0
+            determinant = _normal1[1] * _normal2[2] - _normal1[2] * _normal2[1];
+            if (Math.abs(determinant) > constants_1.EPSILON) {
+                origin[0] = 0;
+                origin[1] = (distance1 * _normal2[2] - distance2 * _normal1[2]) / determinant;
+                origin[2] = (distance2 * _normal1[1] - distance1 * _normal2[1]) / determinant;
+            }
+            else {
+                return false; // No valid intersection point
+            }
+        }
+    }
+    return true;
+}