@faicad/faijs 0.13.0

package/dist/vendored/brepjs/gear/gearFns.js

@@ -0,0 +1,479 @@
+/**
+ * Public gear API — external + internal spur gears, planetary assemblies.
+ *
+ * @example Single external gear
+ * ```typescript
+ * await initOC();
+ * const result = makeExternalGear({ teeth: 24, moduleSize: 2, thickness: 8, bore: 6 });
+ * if (isOk(result)) writeSTEP(result.value.solid, 'gear.step');
+ * ```
+ *
+ * @example Planetary assembly
+ * ```typescript
+ * const planetary = makePlanetaryGear({ thickness: 10 });
+ * if (isOk(planetary)) {
+ *   const { sun, planets, ring, contactRatio } = planetary.value;
+ * }
+ * ```
+ */
+var __addDisposableResource = (this && this.__addDisposableResource) || function (env, value, async) {
+    if (value !== null && value !== void 0) {
+        if (typeof value !== "object" && typeof value !== "function") throw new TypeError("Object expected.");
+        var dispose, inner;
+        if (async) {
+            if (!Symbol.asyncDispose) throw new TypeError("Symbol.asyncDispose is not defined.");
+            dispose = value[Symbol.asyncDispose];
+        }
+        if (dispose === void 0) {
+            if (!Symbol.dispose) throw new TypeError("Symbol.dispose is not defined.");
+            dispose = value[Symbol.dispose];
+            if (async) inner = dispose;
+        }
+        if (typeof dispose !== "function") throw new TypeError("Object not disposable.");
+        if (inner) dispose = function() { try { inner.call(this); } catch (e) { return Promise.reject(e); } };
+        env.stack.push({ value: value, dispose: dispose, async: async });
+    }
+    else if (async) {
+        env.stack.push({ async: true });
+    }
+    return value;
+};
+var __disposeResources = (this && this.__disposeResources) || (function (SuppressedError) {
+    return function (env) {
+        function fail(e) {
+            env.error = env.hasError ? new SuppressedError(e, env.error, "An error was suppressed during disposal.") : e;
+            env.hasError = true;
+        }
+        var r, s = 0;
+        function next() {
+            while (r = env.stack.pop()) {
+                try {
+                    if (!r.async && s === 1) return s = 0, env.stack.push(r), Promise.resolve().then(next);
+                    if (r.dispose) {
+                        var result = r.dispose.call(r.value);
+                        if (r.async) return s |= 2, Promise.resolve(result).then(next, function(e) { fail(e); return next(); });
+                    }
+                    else s |= 1;
+                }
+                catch (e) {
+                    fail(e);
+                }
+            }
+            if (s === 1) return env.hasError ? Promise.reject(env.error) : Promise.resolve();
+            if (env.hasError) throw env.error;
+        }
+        return next();
+    };
+})(typeof SuppressedError === "function" ? SuppressedError : function (error, suppressed, message) {
+    var e = new Error(message);
+    return e.name = "SuppressedError", e.error = error, e.suppressed = suppressed, e;
+});
+import { ok, err, isErr } from '../core/result.js';
+import { validationError } from '../core/errors.js';
+import { makeCircle, assembleWire } from '../topology/curveBuilders.js';
+import { makeFace } from '../topology/surfaceBuilders.js';
+import { DisposalScope } from '../core/disposal.js';
+import { rotate, translate } from '../topology/transformFns.js';
+import { cut } from '../topology/booleanFns.js';
+import { extrude } from '../operations/extrudeFns.js';
+import { DEFAULT_CLEARANCE, DEFAULT_PRESSURE_ANGLE_DEG, backlashHalf, evenToothPhaseOffset, externalExternalContactRatio, externalInternalContactRatio, filletStressConcentrationFactor, gearGeometry, lewisRootStressCorrected, planetSelfRotationAngle, ringTeeth, solvePlanetRingWorkingPressureAngle, solveSunPlanetWorkingPressureAngle, undercutDeficit, validatePlanetary, workingCenterDistance, } from './gearMath.js';
+import { makeExternalGearProfileWire, makeInternalGearProfileWire } from './gearProfile.js';
+export function makeExternalGear(params) {
+    const { teeth, moduleSize, thickness, pressureAngleDeg = DEFAULT_PRESSURE_ANGLE_DEG, shift = 0, clearance = DEFAULT_CLEARANCE, flankThinning = 0, bore = 0, samples, } = params;
+    if (thickness <= 0)
+        return err(validationError('GEAR_THICKNESS_NONPOSITIVE', 'thickness must be > 0'));
+    if (samples !== undefined && (!Number.isInteger(samples) || samples < 1))
+        return err(validationError('GEAR_SAMPLES_INVALID', `samples must be a positive integer; got ${String(samples)}`));
+    const alpha = (pressureAngleDeg * Math.PI) / 180;
+    const geom = gearGeometry(teeth, moduleSize, alpha, shift, clearance, flankThinning, false);
+    if (bore > 0 && bore >= 2 * geom.rRoot)
+        return err(validationError('GEAR_BORE_TOO_LARGE', `bore diameter ${bore.toFixed(2)} ≥ root diameter ${(2 * geom.rRoot).toFixed(2)} — would erase the gear teeth`));
+    const wireResult = makeExternalGearProfileWire({
+        teeth,
+        moduleSize,
+        pressureAngle: alpha,
+        shift,
+        clearance,
+        flankThinning,
+        ...(samples !== undefined ? { samples } : {}),
+    });
+    if (isErr(wireResult))
+        return wireResult;
+    const diagnostics = externalGearDiagnostics(teeth, alpha, shift);
+    return finalizeExternalSolid(wireResult.value, thickness, bore, geom, diagnostics);
+}
+export function makeInternalGear(params) {
+    const { teeth, moduleSize, thickness, pressureAngleDeg = DEFAULT_PRESSURE_ANGLE_DEG, shift = 0, clearance = DEFAULT_CLEARANCE, flankThinning = 0, ringWallThickness = 2 * params.moduleSize, samples, } = params;
+    if (thickness <= 0)
+        return err(validationError('GEAR_THICKNESS_NONPOSITIVE', 'thickness must be > 0'));
+    if (ringWallThickness <= 0)
+        return err(validationError('GEAR_WALL_NONPOSITIVE', 'ringWallThickness must be > 0'));
+    if (samples !== undefined && (!Number.isInteger(samples) || samples < 1))
+        return err(validationError('GEAR_SAMPLES_INVALID', `samples must be a positive integer; got ${String(samples)}`));
+    const alpha = (pressureAngleDeg * Math.PI) / 180;
+    const innerWireResult = makeInternalGearProfileWire({
+        teeth,
+        moduleSize,
+        pressureAngle: alpha,
+        shift,
+        clearance,
+        flankThinning,
+        ...(samples !== undefined ? { samples } : {}),
+    });
+    if (isErr(innerWireResult))
+        return innerWireResult;
+    const geom = gearGeometry(teeth, moduleSize, alpha, shift, clearance, flankThinning, true);
+    const outerRadius = geom.rPitch + ringWallThickness;
+    const outerWireResult = makeOuterCircleWire(outerRadius);
+    if (isErr(outerWireResult))
+        return outerWireResult;
+    // Internal gears use a different undercut criterion (involute-cutter geometry,
+    // shift sign inverted), so the rack-cut formula that drives external-gear
+    // diagnostics doesn't apply here. Leave diagnostics empty until a ring-specific
+    // check exists.
+    return finalizeInternalSolid(outerWireResult.value, innerWireResult.value, thickness, geom, []);
+}
+/** Diagnostics for an external (rack-cut) gear in isolation. Not valid for internal gears. */
+function externalGearDiagnostics(teeth, alpha, shift) {
+    const deficit = undercutDeficit(teeth, alpha, shift);
+    if (deficit <= 0)
+        return [];
+    return [
+        {
+            code: 'UNDERCUT_RISK',
+            severity: 'warning',
+            message: `gear is undercut: increase shift by ${deficit.toFixed(3)} to avoid (z=${teeth})`,
+            context: { deficit, teeth },
+        },
+    ];
+}
+export function makePlanetaryGear(params) {
+    const resolved = resolvePlanetaryParams(params);
+    if (isErr(resolved))
+        return resolved;
+    const cfg = resolved.value;
+    const stages = buildPlanetaryStages(cfg);
+    if (isErr(stages))
+        return stages;
+    const { sun, planet, ring } = stages.value;
+    // Compute metrics before disposing `planet.solid` — the function only reads
+    // scalar fields off GearResult so the order is technically irrelevant today,
+    // but it leaves a footgun for anyone adding a solid-touching read later.
+    const metrics = computeMeshMetrics(cfg, sun, planet, ring);
+    const planets = placePlanets(planet.solid, cfg);
+    planet.solid.delete();
+    const ringPhased = applyRingPhase(ring.solid, cfg.zr);
+    const diagnostics = collectDiagnostics(cfg, metrics);
+    return ok({
+        sun: sun.solid,
+        planets,
+        ring: ringPhased,
+        ringTeeth: cfg.zr,
+        workingPressureAngle: cfg.alphaW_sp,
+        centerDistance: cfg.centerDistance,
+        contactRatio: { sunPlanet: metrics.crSunPlanet, planetRing: metrics.crPlanetRing },
+        undercutDeficit: { sun: metrics.undercutSun, planet: metrics.undercutPlanet },
+        ...(metrics.lewisStress ? { lewisStress: metrics.lewisStress } : {}),
+        ...(metrics.stressConcentrationFactor
+            ? { stressConcentrationFactor: metrics.stressConcentrationFactor }
+            : {}),
+        diagnostics,
+    });
+}
+function buildPlanetaryStages(cfg) {
+    const samplesField = cfg.samples !== undefined ? { samples: cfg.samples } : {};
+    const common = {
+        moduleSize: cfg.moduleSize,
+        thickness: cfg.thickness,
+        pressureAngleDeg: cfg.pressureAngleDeg,
+        clearance: cfg.clearance,
+        flankThinning: cfg.bHalf,
+        ...samplesField,
+    };
+    const sun = makeExternalGear({
+        ...common,
+        teeth: cfg.sunTeeth,
+        shift: cfg.sunShift,
+        bore: cfg.sunBore,
+    });
+    if (isErr(sun))
+        return sun;
+    const planet = makeExternalGear({
+        ...common,
+        teeth: cfg.planetTeeth,
+        shift: cfg.planetShift,
+        bore: cfg.planetBore,
+    });
+    if (isErr(planet)) {
+        sun.value.solid.delete();
+        return planet;
+    }
+    const ring = makeInternalGear({
+        ...common,
+        teeth: cfg.zr,
+        shift: cfg.ringShift,
+        ringWallThickness: cfg.ringWallThickness,
+    });
+    if (isErr(ring)) {
+        sun.value.solid.delete();
+        planet.value.solid.delete();
+        return ring;
+    }
+    return ok({ sun: sun.value, planet: planet.value, ring: ring.value });
+}
+function resolvePlanetaryParams(params) {
+    const moduleSize = params.moduleSize ?? 3;
+    const sunTeeth = params.sunTeeth ?? 15;
+    const planetTeeth = params.planetTeeth ?? 12;
+    const numPlanets = params.numPlanets ?? 3;
+    const pressureAngleDeg = params.pressureAngleDeg ?? DEFAULT_PRESSURE_ANGLE_DEG;
+    const alpha = (pressureAngleDeg * Math.PI) / 180;
+    const sunShift = params.sunShift ?? 0;
+    const planetShift = params.planetShift ?? 0;
+    const ringShift = params.ringShift ?? 0;
+    if (params.thickness <= 0)
+        return err(validationError('GEAR_THICKNESS_NONPOSITIVE', 'thickness must be > 0'));
+    const validation = validatePlanetary(sunTeeth, planetTeeth, numPlanets, planetShift);
+    if (isErr(validation))
+        return validation;
+    const zr = ringTeeth(sunTeeth, planetTeeth);
+    const sp = solveSunPlanetWorkingPressureAngle(alpha, sunShift, planetShift, sunTeeth, planetTeeth);
+    if (isErr(sp))
+        return sp;
+    const pr = solvePlanetRingWorkingPressureAngle(alpha, planetShift, ringShift, planetTeeth, zr);
+    if (isErr(pr))
+        return pr;
+    const centerDistance = workingCenterDistance(sunTeeth, planetTeeth, moduleSize, alpha, sp.value);
+    return ok({
+        moduleSize,
+        sunTeeth,
+        planetTeeth,
+        numPlanets,
+        pressureAngleDeg,
+        alpha,
+        alphaW_sp: sp.value,
+        alphaW_pr: pr.value,
+        clearance: params.clearance ?? DEFAULT_CLEARANCE,
+        bHalf: backlashHalf(params.backlash ?? 0),
+        sunShift,
+        planetShift,
+        ringShift,
+        ringWallThickness: params.ringWallThickness ?? 2 * moduleSize,
+        thickness: params.thickness,
+        sunBore: params.sunBore ?? 0,
+        planetBore: params.planetBore ?? 0,
+        ...(params.appliedTorque !== undefined ? { appliedTorque: params.appliedTorque } : {}),
+        ...(params.samples !== undefined ? { samples: params.samples } : {}),
+        zr,
+        centerDistance,
+    });
+}
+function makeOuterCircleWire(radius) {
+    const env_1 = { stack: [], error: void 0, hasError: false };
+    try {
+        const scope = __addDisposableResource(env_1, new DisposalScope(), false);
+        const circleEdge = scope.register(makeCircle(radius, [0, 0, 0], [0, 0, 1]));
+        const wireResult = assembleWire([circleEdge]);
+        if (isErr(wireResult))
+            return wireResult;
+        // A single-edge circle wire is closed and planar by construction; brand once.
+        return ok(wireResult.value);
+    }
+    catch (e_1) {
+        env_1.error = e_1;
+        env_1.hasError = true;
+    }
+    finally {
+        __disposeResources(env_1);
+    }
+}
+function finalizeExternalSolid(wire, thickness, bore, geom, diagnostics) {
+    const env_2 = { stack: [], error: void 0, hasError: false };
+    try {
+        // The profile wire and its face are consumed by makeFace/extrude (which share
+        // their refcounted TShape); dispose them once the solid is built. The solid is
+        // NOT registered here — it is the returned value on the no-bore path.
+        const profileScope = __addDisposableResource(env_2, new DisposalScope(), false);
+        profileScope.register(wire);
+        const faceResult = makeFace(wire);
+        if (isErr(faceResult))
+            return faceResult;
+        profileScope.register(faceResult.value);
+        const solidResult = extrude(faceResult.value, [0, 0, thickness]);
+        if (isErr(solidResult))
+            return solidResult;
+        if (bore <= 0)
+            return ok(buildGearResult(solidResult.value, geom, diagnostics));
+        // Bore overshoots both ends (z = -0.5 to thickness + 0.5) so no faces are coplanar.
+        const boreScope = __addDisposableResource(env_2, new DisposalScope(), false);
+        boreScope.register(solidResult.value); // un-bored gear, consumed by cut below
+        const boreFaceResult = makeBoreFace(bore / 2);
+        if (isErr(boreFaceResult))
+            return boreFaceResult;
+        boreScope.register(boreFaceResult.value);
+        const boreRaw = extrude(boreFaceResult.value, [0, 0, thickness + 1]);
+        if (isErr(boreRaw))
+            return boreRaw;
+        boreScope.register(boreRaw.value);
+        const boreSolid = boreScope.register(translate(boreRaw.value, [0, 0, -0.5]));
+        const cutResult = cut(solidResult.value, boreSolid);
+        if (isErr(cutResult))
+            return cutResult;
+        return ok(buildGearResult(cutResult.value, geom, diagnostics));
+    }
+    catch (e_2) {
+        env_2.error = e_2;
+        env_2.hasError = true;
+    }
+    finally {
+        __disposeResources(env_2);
+    }
+}
+function makeBoreFace(radius) {
+    const wireResult = makeOuterCircleWire(radius);
+    if (isErr(wireResult))
+        return wireResult;
+    return makeFace(wireResult.value);
+}
+function finalizeInternalSolid(outerWire, innerToothedWire, thickness, geom, diagnostics) {
+    const faceResult = makeFace(outerWire, [innerToothedWire]);
+    if (isErr(faceResult))
+        return faceResult;
+    const solidResult = extrude(faceResult.value, [0, 0, thickness]);
+    if (isErr(solidResult))
+        return solidResult;
+    return ok(buildGearResult(solidResult.value, geom, diagnostics));
+}
+function buildGearResult(solid, geom, diagnostics = []) {
+    return {
+        solid,
+        diagnostics,
+        pitchDiameter: 2 * geom.rPitch,
+        baseDiameter: 2 * geom.rb,
+        tipDiameter: 2 * geom.rTip,
+        rootDiameter: 2 * geom.rRoot,
+    };
+}
+function placePlanets(planetProto, cfg) {
+    const planets = [];
+    try {
+        for (let i = 0; i < cfg.numPlanets; i++) {
+            const orbital = (i * 2 * Math.PI) / cfg.numPlanets;
+            const selfRot = planetSelfRotationAngle(orbital, cfg.sunTeeth, cfg.planetTeeth);
+            const rotated = rotate(planetProto, (selfRot * 180) / Math.PI);
+            try {
+                const offset = [
+                    cfg.centerDistance * Math.cos(orbital),
+                    cfg.centerDistance * Math.sin(orbital),
+                    0,
+                ];
+                planets.push(translate(rotated, offset));
+            }
+            finally {
+                rotated.delete();
+            }
+        }
+        return planets;
+    }
+    catch (e) {
+        // Mid-loop failure: clean up the planets that did succeed before rethrowing.
+        for (const p of planets)
+            p.delete();
+        throw e;
+    }
+}
+function applyRingPhase(ring, zr) {
+    const phaseRad = evenToothPhaseOffset(zr);
+    if (phaseRad === 0)
+        return ring;
+    const rotated = rotate(ring, (phaseRad * 180) / Math.PI);
+    ring.delete();
+    return rotated;
+}
+function computeMeshMetrics(cfg, sun, planet, ring) {
+    const crSunPlanet = externalExternalContactRatio(sun.tipDiameter / 2, sun.baseDiameter / 2, planet.tipDiameter / 2, planet.baseDiameter / 2, cfg.centerDistance, cfg.moduleSize, cfg.alpha, cfg.alphaW_sp);
+    const crPlanetRing = externalInternalContactRatio(planet.tipDiameter / 2, planet.baseDiameter / 2, ring.tipDiameter / 2, ring.baseDiameter / 2, cfg.centerDistance, cfg.moduleSize, cfg.alpha, cfg.alphaW_pr);
+    const undercutSun = undercutDeficit(cfg.sunTeeth, cfg.alpha, cfg.sunShift);
+    const undercutPlanet = undercutDeficit(cfg.planetTeeth, cfg.alpha, cfg.planetShift);
+    const metrics = { crSunPlanet, crPlanetRing, undercutSun, undercutPlanet };
+    if (cfg.appliedTorque !== undefined) {
+        // Force balance: W_t is shared at each mesh, so each gear's own-pitch torque
+        // is T_eff = T_sun · z / z_sun (preserves W_t through Lewis's 2T/(z·m) term).
+        const tSun = cfg.appliedTorque;
+        metrics.lewisStress = {
+            sun: lewisRootStressCorrected(tSun, cfg.moduleSize, cfg.thickness, cfg.sunTeeth, cfg.alpha),
+            planet: lewisRootStressCorrected((tSun * cfg.planetTeeth) / cfg.sunTeeth, cfg.moduleSize, cfg.thickness, cfg.planetTeeth, cfg.alpha),
+            ring: lewisRootStressCorrected((tSun * cfg.zr) / cfg.sunTeeth, cfg.moduleSize, cfg.thickness, cfg.zr, cfg.alpha, true),
+        };
+        metrics.stressConcentrationFactor = {
+            sun: filletStressConcentrationFactor(cfg.sunTeeth, cfg.alpha),
+            planet: filletStressConcentrationFactor(cfg.planetTeeth, cfg.alpha),
+            ring: filletStressConcentrationFactor(cfg.zr, cfg.alpha, true),
+        };
+    }
+    return metrics;
+}
+function collectDiagnostics(cfg, metrics) {
+    const diagnostics = [];
+    if (metrics.crSunPlanet < 1.2) {
+        diagnostics.push({
+            code: 'CONTACT_RATIO_LOW_SUN_PLANET',
+            severity: 'warning',
+            message: `sun-planet contact ratio ${metrics.crSunPlanet.toFixed(2)} is below 1.2 — may run unevenly`,
+            context: { value: metrics.crSunPlanet },
+        });
+    }
+    if (metrics.crPlanetRing < 1.2) {
+        diagnostics.push({
+            code: 'CONTACT_RATIO_LOW_PLANET_RING',
+            severity: 'warning',
+            message: `planet-ring contact ratio ${metrics.crPlanetRing.toFixed(2)} is below 1.2 — may run unevenly`,
+            context: { value: metrics.crPlanetRing },
+        });
+    }
+    if (metrics.undercutSun > 0) {
+        diagnostics.push({
+            code: 'UNDERCUT_RISK_SUN',
+            severity: 'warning',
+            message: `sun gear is undercut: increase sunShift by ${metrics.undercutSun.toFixed(3)} to avoid`,
+            context: { deficit: metrics.undercutSun, sunTeeth: cfg.sunTeeth },
+        });
+    }
+    if (metrics.undercutPlanet > 0) {
+        diagnostics.push({
+            code: 'UNDERCUT_RISK_PLANET',
+            severity: 'warning',
+            message: `planet gear is undercut: increase planetShift by ${metrics.undercutPlanet.toFixed(3)} to avoid`,
+            context: { deficit: metrics.undercutPlanet, planetTeeth: cfg.planetTeeth },
+        });
+    }
+    // Only emit when shifts are large enough to matter. The diagnostic exists to
+    // warn that Lewis Y(z) uses unshifted geometry — at ~0.05 of profile shift
+    // the error is under 2.5%, well within engineering noise for a strength
+    // estimate, so don't pester users on canonical setups with small shifts.
+    const totalShiftMagnitude = Math.abs(cfg.sunShift) + Math.abs(cfg.planetShift) + Math.abs(cfg.ringShift);
+    if (cfg.appliedTorque !== undefined && totalShiftMagnitude > 0.05) {
+        diagnostics.push({
+            code: 'LEWIS_Y_SHIFT_UNCORRECTED',
+            severity: 'info',
+            message: 'Lewis stress uses unshifted Y(z) and K_f approximations; expect ±5% per 0.1 of profile shift',
+        });
+    }
+    // Kinematic compatibility: x_ring should equal x_sun + 2·x_planet for both meshes to share αw.
+    const kinematicError = cfg.ringShift - (cfg.sunShift + 2 * cfg.planetShift);
+    if (Math.abs(kinematicError) > 1e-6) {
+        diagnostics.push({
+            code: 'PLANETARY_SHIFT_KINEMATIC_MISMATCH',
+            severity: 'warning',
+            message: `ringShift should equal sunShift + 2·planetShift for both meshes to share working PA; off by ${kinematicError.toFixed(3)}`,
+            context: {
+                kinematicError,
+                sunShift: cfg.sunShift,
+                planetShift: cfg.planetShift,
+                ringShift: cfg.ringShift,
+            },
+        });
+    }
+    return diagnostics;
+}
+//# sourceMappingURL=gearFns.js.map

package/dist/vendored/brepjs/gear/gearFns.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gearFns.js","sourceRoot":"","sources":["../../../../src/vendored/brepjs/gear/gearFns.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;GAiBG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAGH,OAAO,EAAe,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,MAAM,mBAAmB,CAAC;AAChE,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AAEpD,OAAO,EAAE,UAAU,EAAE,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACxE,OAAO,EAAE,QAAQ,EAAE,MAAM,gCAAgC,CAAC;AAC1D,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AACpD,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,6BAA6B,CAAC;AAChE,OAAO,EAAE,GAAG,EAAE,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAE,OAAO,EAAE,MAAM,6BAA6B,CAAC;AACtD,OAAO,EAEL,iBAAiB,EACjB,0BAA0B,EAC1B,YAAY,EACZ,oBAAoB,EACpB,4BAA4B,EAC5B,4BAA4B,EAC5B,+BAA+B,EAC/B,YAAY,EACZ,wBAAwB,EACxB,uBAAuB,EACvB,SAAS,EACT,mCAAmC,EACnC,kCAAkC,EAClC,eAAe,EACf,iBAAiB,EACjB,qBAAqB,GACtB,MAAM,eAAe,CAAC;AACvB,OAAO,EAAE,2BAA2B,EAAE,2BAA2B,EAAE,MAAM,kBAAkB,CAAC;AAoH5F,MAAM,UAAU,gBAAgB,CAAC,MAA0B;IACzD,MAAM,EACJ,KAAK,EACL,UAAU,EACV,SAAS,EACT,gBAAgB,GAAG,0BAA0B,EAC7C,KAAK,GAAG,CAAC,EACT,SAAS,GAAG,iBAAiB,EAC7B,aAAa,GAAG,CAAC,EACjB,IAAI,GAAG,CAAC,EACR,OAAO,GACR,GAAG,MAAM,CAAC;IACX,IAAI,SAAS,IAAI,CAAC;QAChB,OAAO,GAAG,CAAC,eAAe,CAAC,4BAA4B,EAAE,uBAAuB,CAAC,CAAC,CAAC;IACrF,IAAI,OAAO,KAAK,SAAS,IAAI,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,IAAI,OAAO,GAAG,CAAC,CAAC;QACtE,OAAO,GAAG,CACR,eAAe,CACb,sBAAsB,EACtB,2CAA2C,MAAM,CAAC,OAAO,CAAC,EAAE,CAC7D,CACF,CAAC;IAEJ,MAAM,KAAK,GAAG,CAAC,gBAAgB,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACjD,MAAM,IAAI,GAAG,YAAY,CAAC,KAAK,EAAE,UAAU,EAAE,KAAK,EAAE,KAAK,EAAE,SAAS,EAAE,aAAa,EAAE,KAAK,CAAC,CAAC;IAC5F,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK;QACpC,OAAO,GAAG,CACR,eAAe,CACb,qBAAqB,EACrB,iBAAiB,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,oBAAoB,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,+BAA+B,CAC/G,CACF,CAAC;IAEJ,MAAM,UAAU,GAAG,2BAA2B,CAAC;QAC7C,KAAK;QACL,UAAU;QACV,aAAa,EAAE,KAAK;QACpB,KAAK;QACL,SAAS;QACT,aAAa;QACb,GAAG,CAAC,OAAO,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;KAC9C,CAAC,CAAC;IACH,IAAI,KAAK,CAAC,UAAU,CAAC;QAAE,OAAO,UAAU,CAAC;IACzC,MAAM,WAAW,GAAG,uBAAuB,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;IACjE,OAAO,qBAAqB,CAAC,UAAU,CAAC,KAAK,EAAE,SAAS,EAAE,IAAI,EAAE,IAAI,EAAE,WAAW,CAAC,CAAC;AACrF,CAAC;AAED,MAAM,UAAU,gBAAgB,CAAC,MAA0B;IACzD,MAAM,EACJ,KAAK,EACL,UAAU,EACV,SAAS,EACT,gBAAgB,GAAG,0BAA0B,EAC7C,KAAK,GAAG,CAAC,EACT,SAAS,GAAG,iBAAiB,EAC7B,aAAa,GAAG,CAAC,EACjB,iBAAiB,GAAG,CAAC,GAAG,MAAM,CAAC,UAAU,EACzC,OAAO,GACR,GAAG,MAAM,CAAC;IACX,IAAI,SAAS,IAAI,CAAC;QAChB,OAAO,GAAG,CAAC,eAAe,CAAC,4BAA4B,EAAE,uBAAuB,CAAC,CAAC,CAAC;IACrF,IAAI,iBAAiB,IAAI,CAAC;QACxB,OAAO,GAAG,CAAC,eAAe,CAAC,uBAAuB,EAAE,+BAA+B,CAAC,CAAC,CAAC;IACxF,IAAI,OAAO,KAAK,SAAS,IAAI,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,IAAI,OAAO,GAAG,CAAC,CAAC;QACtE,OAAO,GAAG,CACR,eAAe,CACb,sBAAsB,EACtB,2CAA2C,MAAM,CAAC,OAAO,CAAC,EAAE,CAC7D,CACF,CAAC;IAEJ,MAAM,KAAK,GAAG,CAAC,gBAAgB,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACjD,MAAM,eAAe,GAAG,2BAA2B,CAAC;QAClD,KAAK;QACL,UAAU;QACV,aAAa,EAAE,KAAK;QACpB,KAAK;QACL,SAAS;QACT,aAAa;QACb,GAAG,CAAC,OAAO,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;KAC9C,CAAC,CAAC;IACH,IAAI,KAAK,CAAC,eAAe,CAAC;QAAE,OAAO,eAAe,CAAC;IAEnD,MAAM,IAAI,GAAG,YAAY,CAAC,KAAK,EAAE,UAAU,EAAE,KAAK,EAAE,KAAK,EAAE,SAAS,EAAE,aAAa,EAAE,IAAI,CAAC,CAAC;IAC3F,MAAM,WAAW,GAAG,IAAI,CAAC,MAAM,GAAG,iBAAiB,CAAC;IACpD,MAAM,eAAe,GAAG,mBAAmB,CAAC,WAAW,CAAC,CAAC;IACzD,IAAI,KAAK,CAAC,eAAe,CAAC;QAAE,OAAO,eAAe,CAAC;IAEnD,+EAA+E;IAC/E,0EAA0E;IAC1E,gFAAgF;IAChF,gBAAgB;IAChB,OAAO,qBAAqB,CAAC,eAAe,CAAC,KAAK,EAAE,eAAe,CAAC,KAAK,EAAE,SAAS,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC;AAClG,CAAC;AAED,8FAA8F;AAC9F,SAAS,uBAAuB,CAAC,KAAa,EAAE,KAAa,EAAE,KAAa;IAC1E,MAAM,OAAO,GAAG,eAAe,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;IACrD,IAAI,OAAO,IAAI,CAAC;QAAE,OAAO,EAAE,CAAC;IAC5B,OAAO;QACL;YACE,IAAI,EAAE,eAAe;YACrB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,uCAAuC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,gBAAgB,KAAK,GAAG;YAC1F,OAAO,EAAE,EAAE,OAAO,EAAE,KAAK,EAAE;SAC5B;KACF,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,iBAAiB,CAAC,MAA2B;IAC3D,MAAM,QAAQ,GAAG,sBAAsB,CAAC,MAAM,CAAC,CAAC;IAChD,IAAI,KAAK,CAAC,QAAQ,CAAC;QAAE,OAAO,QAAQ,CAAC;IACrC,MAAM,GAAG,GAAG,QAAQ,CAAC,KAAK,CAAC;IAE3B,MAAM,MAAM,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC;IACzC,IAAI,KAAK,CAAC,MAAM,CAAC;QAAE,OAAO,MAAM,CAAC;IACjC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC;IAE3C,4EAA4E;IAC5E,6EAA6E;IAC7E,yEAAyE;IACzE,MAAM,OAAO,GAAG,kBAAkB,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,CAAC,CAAC;IAC3D,MAAM,OAAO,GAAG,YAAY,CAAC,MAAM,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;IAChD,MAAM,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;IACtB,MAAM,UAAU,GAAG,cAAc,CAAC,IAAI,CAAC,KAAK,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC;IACtD,MAAM,WAAW,GAAG,kBAAkB,CAAC,GAAG,EAAE,OAAO,CAAC,CAAC;IAErD,OAAO,EAAE,CAAC;QACR,GAAG,EAAE,GAAG,CAAC,KAAK;QACd,OAAO;QACP,IAAI,EAAE,UAAU;QAChB,SAAS,EAAE,GAAG,CAAC,EAAE;QACjB,oBAAoB,EAAE,GAAG,CAAC,SAAS;QACnC,cAAc,EAAE,GAAG,CAAC,cAAc;QAClC,YAAY,EAAE,EAAE,SAAS,EAAE,OAAO,CAAC,WAAW,EAAE,UAAU,EAAE,OAAO,CAAC,YAAY,EAAE;QAClF,eAAe,EAAE,EAAE,GAAG,EAAE,OAAO,CAAC,WAAW,EAAE,MAAM,EAAE,OAAO,CAAC,cAAc,EAAE;QAC7E,GAAG,CAAC,OAAO,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,EAAE,OAAO,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACpE,GAAG,CAAC,OAAO,CAAC,yBAAyB;YACnC,CAAC,CAAC,EAAE,yBAAyB,EAAE,OAAO,CAAC,yBAAyB,EAAE;YAClE,CAAC,CAAC,EAAE,CAAC;QACP,WAAW;KACZ,CAAC,CAAC;AACL,CAAC;AAED,SAAS,oBAAoB,CAC3B,GAAsB;IAEtB,MAAM,YAAY,GAAG,GAAG,CAAC,OAAO,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,GAAG,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;IAC/E,MAAM,MAAM,GAAG;QACb,UAAU,EAAE,GAAG,CAAC,UAAU;QAC1B,SAAS,EAAE,GAAG,CAAC,SAAS;QACxB,gBAAgB,EAAE,GAAG,CAAC,gBAAgB;QACtC,SAAS,EAAE,GAAG,CAAC,SAAS;QACxB,aAAa,EAAE,GAAG,CAAC,KAAK;QACxB,GAAG,YAAY;KAChB,CAAC;IACF,MAAM,GAAG,GAAG,gBAAgB,CAAC;QAC3B,GAAG,MAAM;QACT,KAAK,EAAE,GAAG,CAAC,QAAQ;QACnB,KAAK,EAAE,GAAG,CAAC,QAAQ;QACnB,IAAI,EAAE,GAAG,CAAC,OAAO;KAClB,CAAC,CAAC;IACH,IAAI,KAAK,CAAC,GAAG,CAAC;QAAE,OAAO,GAAG,CAAC;IAC3B,MAAM,MAAM,GAAG,gBAAgB,CAAC;QAC9B,GAAG,MAAM;QACT,KAAK,EAAE,GAAG,CAAC,WAAW;QACtB,KAAK,EAAE,GAAG,CAAC,WAAW;QACtB,IAAI,EAAE,GAAG,CAAC,UAAU;KACrB,CAAC,CAAC;IACH,IAAI,KAAK,CAAC,MAAM,CAAC,EAAE,CAAC;QAClB,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QACzB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,MAAM,IAAI,GAAG,gBAAgB,CAAC;QAC5B,GAAG,MAAM;QACT,KAAK,EAAE,GAAG,CAAC,EAAE;QACb,KAAK,EAAE,GAAG,CAAC,SAAS;QACpB,iBAAiB,EAAE,GAAG,CAAC,iBAAiB;KACzC,CAAC,CAAC;IACH,IAAI,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC;QAChB,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QACzB,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,KAAK,EAAE,CAAC,CAAC;AACxE,CAAC;AA4BD,SAAS,sBAAsB,CAAC,MAA2B;IACzD,MAAM,UAAU,GAAG,MAAM,CAAC,UAAU,IAAI,CAAC,CAAC;IAC1C,MAAM,QAAQ,GAAG,MAAM,CAAC,QAAQ,IAAI,EAAE,CAAC;IACvC,MAAM,WAAW,GAAG,MAAM,CAAC,WAAW,IAAI,EAAE,CAAC;IAC7C,MAAM,UAAU,GAAG,MAAM,CAAC,UAAU,IAAI,CAAC,CAAC;IAC1C,MAAM,gBAAgB,GAAG,MAAM,CAAC,gBAAgB,IAAI,0BAA0B,CAAC;IAC/E,MAAM,KAAK,GAAG,CAAC,gBAAgB,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACjD,MAAM,QAAQ,GAAG,MAAM,CAAC,QAAQ,IAAI,CAAC,CAAC;IACtC,MAAM,WAAW,GAAG,MAAM,CAAC,WAAW,IAAI,CAAC,CAAC;IAC5C,MAAM,SAAS,GAAG,MAAM,CAAC,SAAS,IAAI,CAAC,CAAC;IAExC,IAAI,MAAM,CAAC,SAAS,IAAI,CAAC;QACvB,OAAO,GAAG,CAAC,eAAe,CAAC,4BAA4B,EAAE,uBAAuB,CAAC,CAAC,CAAC;IAErF,MAAM,UAAU,GAAG,iBAAiB,CAAC,QAAQ,EAAE,WAAW,EAAE,UAAU,EAAE,WAAW,CAAC,CAAC;IACrF,IAAI,KAAK,CAAC,UAAU,CAAC;QAAE,OAAO,UAAU,CAAC;IAEzC,MAAM,EAAE,GAAG,SAAS,CAAC,QAAQ,EAAE,WAAW,CAAC,CAAC;IAE5C,MAAM,EAAE,GAAG,kCAAkC,CAC3C,KAAK,EACL,QAAQ,EACR,WAAW,EACX,QAAQ,EACR,WAAW,CACZ,CAAC;IACF,IAAI,KAAK,CAAC,EAAE,CAAC;QAAE,OAAO,EAAE,CAAC;IACzB,MAAM,EAAE,GAAG,mCAAmC,CAAC,KAAK,EAAE,WAAW,EAAE,SAAS,EAAE,WAAW,EAAE,EAAE,CAAC,CAAC;IAC/F,IAAI,KAAK,CAAC,EAAE,CAAC;QAAE,OAAO,EAAE,CAAC;IACzB,MAAM,cAAc,GAAG,qBAAqB,CAAC,QAAQ,EAAE,WAAW,EAAE,UAAU,EAAE,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC;IAEjG,OAAO,EAAE,CAAC;QACR,UAAU;QACV,QAAQ;QACR,WAAW;QACX,UAAU;QACV,gBAAgB;QAChB,KAAK;QACL,SAAS,EAAE,EAAE,CAAC,KAAK;QACnB,SAAS,EAAE,EAAE,CAAC,KAAK;QACnB,SAAS,EAAE,MAAM,CAAC,SAAS,IAAI,iBAAiB;QAChD,KAAK,EAAE,YAAY,CAAC,MAAM,CAAC,QAAQ,IAAI,CAAC,CAAC;QACzC,QAAQ;QACR,WAAW;QACX,SAAS;QACT,iBAAiB,EAAE,MAAM,CAAC,iBAAiB,IAAI,CAAC,GAAG,UAAU;QAC7D,SAAS,EAAE,MAAM,CAAC,SAAS;QAC3B,OAAO,EAAE,MAAM,CAAC,OAAO,IAAI,CAAC;QAC5B,UAAU,EAAE,MAAM,CAAC,UAAU,IAAI,CAAC;QAClC,GAAG,CAAC,MAAM,CAAC,aAAa,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,aAAa,EAAE,MAAM,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACtF,GAAG,CAAC,MAAM,CAAC,OAAO,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,MAAM,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACpE,EAAE;QACF,cAAc;KACf,CAAC,CAAC;AACL,CAAC;AAED,SAAS,mBAAmB,CAAC,MAAc;;;QACzC,MAAM,KAAK,kCAAG,IAAI,aAAa,EAAE,QAAA,CAAC;QAClC,MAAM,UAAU,GAAG,KAAK,CAAC,QAAQ,CAAC,UAAU,CAAC,MAAM,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5E,MAAM,UAAU,GAAG,YAAY,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC;QAC9C,IAAI,KAAK,CAAC,UAAU,CAAC;YAAE,OAAO,UAAU,CAAC;QACzC,8EAA8E;QAC9E,OAAO,EAAE,CAAC,UAAU,CAAC,KAAgC,CAAC,CAAC;;;;;;;;;CACxD;AAED,SAAS,qBAAqB,CAC5B,IAA6B,EAC7B,SAAiB,EACjB,IAAY,EACZ,IAAqC,EACrC,WAA6B;;;QAE7B,8EAA8E;QAC9E,+EAA+E;QAC/E,sEAAsE;QACtE,MAAM,YAAY,kCAAG,IAAI,aAAa,EAAE,QAAA,CAAC;QACzC,YAAY,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QAC5B,MAAM,UAAU,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;QAClC,IAAI,KAAK,CAAC,UAAU,CAAC;YAAE,OAAO,UAAU,CAAC;QACzC,YAAY,CAAC,QAAQ,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;QACxC,MAAM,WAAW,GAAG,OAAO,CAAC,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC;QACjE,IAAI,KAAK,CAAC,WAAW,CAAC;YAAE,OAAO,WAAW,CAAC;QAC3C,IAAI,IAAI,IAAI,CAAC;YAAE,OAAO,EAAE,CAAC,eAAe,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,EAAE,WAAW,CAAC,CAAC,CAAC;QAEhF,oFAAoF;QACpF,MAAM,SAAS,kCAAG,IAAI,aAAa,EAAE,QAAA,CAAC;QACtC,SAAS,CAAC,QAAQ,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC,uCAAuC;QAC9E,MAAM,cAAc,GAAG,YAAY,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QAC9C,IAAI,KAAK,CAAC,cAAc,CAAC;YAAE,OAAO,cAAc,CAAC;QACjD,SAAS,CAAC,QAAQ,CAAC,cAAc,CAAC,KAAK,CAAC,CAAC;QACzC,MAAM,OAAO,GAAG,OAAO,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC;QACrE,IAAI,KAAK,CAAC,OAAO,CAAC;YAAE,OAAO,OAAO,CAAC;QACnC,SAAS,CAAC,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;QAClC,MAAM,SAAS,GAAG,SAAS,CAAC,QAAQ,CAAC,SAAS,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC7E,MAAM,SAAS,GAAG,GAAG,CAAC,WAAW,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;QACpD,IAAI,KAAK,CAAC,SAAS,CAAC;YAAE,OAAO,SAAS,CAAC;QACvC,OAAO,EAAE,CAAC,eAAe,CAAC,SAAS,CAAC,KAAK,EAAE,IAAI,EAAE,WAAW,CAAC,CAAC,CAAC;;;;;;;;;CAChE;AAED,SAAS,YAAY,CAAC,MAAc;IAClC,MAAM,UAAU,GAAG,mBAAmB,CAAC,MAAM,CAAC,CAAC;IAC/C,IAAI,KAAK,CAAC,UAAU,CAAC;QAAE,OAAO,UAAU,CAAC;IACzC,OAAO,QAAQ,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;AACpC,CAAC;AAED,SAAS,qBAAqB,CAC5B,SAAkC,EAClC,gBAAyC,EACzC,SAAiB,EACjB,IAAqC,EACrC,WAA6B;IAE7B,MAAM,UAAU,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC,gBAAgB,CAAC,CAAC,CAAC;IAC3D,IAAI,KAAK,CAAC,UAAU,CAAC;QAAE,OAAO,UAAU,CAAC;IACzC,MAAM,WAAW,GAAG,OAAO,CAAC,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC;IACjE,IAAI,KAAK,CAAC,WAAW,CAAC;QAAE,OAAO,WAAW,CAAC;IAC3C,OAAO,EAAE,CAAC,eAAe,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,EAAE,WAAW,CAAC,CAAC,CAAC;AACnE,CAAC;AAED,SAAS,eAAe,CACtB,KAAiB,EACjB,IAAqC,EACrC,cAAgC,EAAE;IAElC,OAAO;QACL,KAAK;QACL,WAAW;QACX,aAAa,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM;QAC9B,YAAY,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE;QACzB,WAAW,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI;QAC1B,YAAY,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK;KAC7B,CAAC;AACJ,CAAC;AAED,SAAS,YAAY,CAAC,WAAuB,EAAE,GAAsB;IACnE,MAAM,OAAO,GAAiB,EAAE,CAAC;IACjC,IAAI,CAAC;QACH,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;YACxC,MAAM,OAAO,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,UAAU,CAAC;YACnD,MAAM,OAAO,GAAG,uBAAuB,CAAC,OAAO,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,WAAW,CAAC,CAAC;YAChF,MAAM,OAAO,GAAG,MAAM,CAAC,WAAW,EAAE,CAAC,OAAO,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;YAC/D,IAAI,CAAC;gBACH,MAAM,MAAM,GAAS;oBACnB,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC;oBACtC,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC;oBACtC,CAAC;iBACF,CAAC;gBACF,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC,CAAC;YAC3C,CAAC;oBAAS,CAAC;gBACT,OAAO,CAAC,MAAM,EAAE,CAAC;YACnB,CAAC;QACH,CAAC;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;IAAC,OAAO,CAAC,EAAE,CAAC;QACX,6EAA6E;QAC7E,KAAK,MAAM,CAAC,IAAI,OAAO;YAAE,CAAC,CAAC,MAAM,EAAE,CAAC;QACpC,MAAM,CAAC,CAAC;IACV,CAAC;AACH,CAAC;AAED,SAAS,cAAc,CAAC,IAAgB,EAAE,EAAU;IAClD,MAAM,QAAQ,GAAG,oBAAoB,CAAC,EAAE,CAAC,CAAC;IAC1C,IAAI,QAAQ,KAAK,CAAC;QAAE,OAAO,IAAI,CAAC;IAChC,MAAM,OAAO,GAAG,MAAM,CAAC,IAAI,EAAE,CAAC,QAAQ,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IACzD,IAAI,CAAC,MAAM,EAAE,CAAC;IACd,OAAO,OAAO,CAAC;AACjB,CAAC;AAWD,SAAS,kBAAkB,CACzB,GAAsB,EACtB,GAAe,EACf,MAAkB,EAClB,IAAgB;IAEhB,MAAM,WAAW,GAAG,4BAA4B,CAC9C,GAAG,CAAC,WAAW,GAAG,CAAC,EACnB,GAAG,CAAC,YAAY,GAAG,CAAC,EACpB,MAAM,CAAC,WAAW,GAAG,CAAC,EACtB,MAAM,CAAC,YAAY,GAAG,CAAC,EACvB,GAAG,CAAC,cAAc,EAClB,GAAG,CAAC,UAAU,EACd,GAAG,CAAC,KAAK,EACT,GAAG,CAAC,SAAS,CACd,CAAC;IACF,MAAM,YAAY,GAAG,4BAA4B,CAC/C,MAAM,CAAC,WAAW,GAAG,CAAC,EACtB,MAAM,CAAC,YAAY,GAAG,CAAC,EACvB,IAAI,CAAC,WAAW,GAAG,CAAC,EACpB,IAAI,CAAC,YAAY,GAAG,CAAC,EACrB,GAAG,CAAC,cAAc,EAClB,GAAG,CAAC,UAAU,EACd,GAAG,CAAC,KAAK,EACT,GAAG,CAAC,SAAS,CACd,CAAC;IACF,MAAM,WAAW,GAAG,eAAe,CAAC,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC;IAC3E,MAAM,cAAc,GAAG,eAAe,CAAC,GAAG,CAAC,WAAW,EAAE,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,WAAW,CAAC,CAAC;IAEpF,MAAM,OAAO,GAAgB,EAAE,WAAW,EAAE,YAAY,EAAE,WAAW,EAAE,cAAc,EAAE,CAAC;IACxF,IAAI,GAAG,CAAC,aAAa,KAAK,SAAS,EAAE,CAAC;QACpC,6EAA6E;QAC7E,8EAA8E;QAC9E,MAAM,IAAI,GAAG,GAAG,CAAC,aAAa,CAAC;QAC/B,OAAO,CAAC,WAAW,GAAG;YACpB,GAAG,EAAE,wBAAwB,CAAC,IAAI,EAAE,GAAG,CAAC,UAAU,EAAE,GAAG,CAAC,SAAS,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,KAAK,CAAC;YAC3F,MAAM,EAAE,wBAAwB,CAC9B,CAAC,IAAI,GAAG,GAAG,CAAC,WAAW,CAAC,GAAG,GAAG,CAAC,QAAQ,EACvC,GAAG,CAAC,UAAU,EACd,GAAG,CAAC,SAAS,EACb,GAAG,CAAC,WAAW,EACf,GAAG,CAAC,KAAK,CACV;YACD,IAAI,EAAE,wBAAwB,CAC5B,CAAC,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,QAAQ,EAC9B,GAAG,CAAC,UAAU,EACd,GAAG,CAAC,SAAS,EACb,GAAG,CAAC,EAAE,EACN,GAAG,CAAC,KAAK,EACT,IAAI,CACL;SACF,CAAC;QACF,OAAO,CAAC,yBAAyB,GAAG;YAClC,GAAG,EAAE,+BAA+B,CAAC,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,KAAK,CAAC;YAC7D,MAAM,EAAE,+BAA+B,CAAC,GAAG,CAAC,WAAW,EAAE,GAAG,CAAC,KAAK,CAAC;YACnE,IAAI,EAAE,+BAA+B,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC;SAC/D,CAAC;IACJ,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AAED,SAAS,kBAAkB,CAAC,GAAsB,EAAE,OAAoB;IACtE,MAAM,WAAW,GAAqB,EAAE,CAAC;IACzC,IAAI,OAAO,CAAC,WAAW,GAAG,GAAG,EAAE,CAAC;QAC9B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,8BAA8B;YACpC,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,4BAA4B,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,kCAAkC;YACrG,OAAO,EAAE,EAAE,KAAK,EAAE,OAAO,CAAC,WAAW,EAAE;SACxC,CAAC,CAAC;IACL,CAAC;IACD,IAAI,OAAO,CAAC,YAAY,GAAG,GAAG,EAAE,CAAC;QAC/B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,+BAA+B;YACrC,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,6BAA6B,OAAO,CAAC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,kCAAkC;YACvG,OAAO,EAAE,EAAE,KAAK,EAAE,OAAO,CAAC,YAAY,EAAE;SACzC,CAAC,CAAC;IACL,CAAC;IACD,IAAI,OAAO,CAAC,WAAW,GAAG,CAAC,EAAE,CAAC;QAC5B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,mBAAmB;YACzB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,8CAA8C,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,WAAW;YAChG,OAAO,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,WAAW,EAAE,QAAQ,EAAE,GAAG,CAAC,QAAQ,EAAE;SAClE,CAAC,CAAC;IACL,CAAC;IACD,IAAI,OAAO,CAAC,cAAc,GAAG,CAAC,EAAE,CAAC;QAC/B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,sBAAsB;YAC5B,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,oDAAoD,OAAO,CAAC,cAAc,CAAC,OAAO,CAAC,CAAC,CAAC,WAAW;YACzG,OAAO,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,cAAc,EAAE,WAAW,EAAE,GAAG,CAAC,WAAW,EAAE;SAC3E,CAAC,CAAC;IACL,CAAC;IACD,6EAA6E;IAC7E,2EAA2E;IAC3E,wEAAwE;IACxE,yEAAyE;IACzE,MAAM,mBAAmB,GACvB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,WAAW,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IAC/E,IAAI,GAAG,CAAC,aAAa,KAAK,SAAS,IAAI,mBAAmB,GAAG,IAAI,EAAE,CAAC;QAClE,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,2BAA2B;YACjC,QAAQ,EAAE,MAAM;YAChB,OAAO,EACL,8FAA8F;SACjG,CAAC,CAAC;IACL,CAAC;IAED,+FAA+F;IAC/F,MAAM,cAAc,GAAG,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,GAAG,CAAC,WAAW,CAAC,CAAC;IAC5E,IAAI,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,GAAG,IAAI,EAAE,CAAC;QACpC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,oCAAoC;YAC1C,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,+FAA+F,cAAc,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE;YACnI,OAAO,EAAE;gBACP,cAAc;gBACd,QAAQ,EAAE,GAAG,CAAC,QAAQ;gBACtB,WAAW,EAAE,GAAG,CAAC,WAAW;gBAC5B,SAAS,EAAE,GAAG,CAAC,SAAS;aACzB;SACF,CAAC,CAAC;IACL,CAAC;IACD,OAAO,WAAW,CAAC;AACrB,CAAC","sourcesContent":["/**\n * Public gear API — external + internal spur gears, planetary assemblies.\n *\n * @example Single external gear\n * ```typescript\n * await initOC();\n * const result = makeExternalGear({ teeth: 24, moduleSize: 2, thickness: 8, bore: 6 });\n * if (isOk(result)) writeSTEP(result.value.solid, 'gear.step');\n * ```\n *\n * @example Planetary assembly\n * ```typescript\n * const planetary = makePlanetaryGear({ thickness: 10 });\n * if (isOk(planetary)) {\n *   const { sun, planets, ring, contactRatio } = planetary.value;\n * }\n * ```\n */\n\nimport type { Vec3 } from '../core/types.js';\nimport { type Result, ok, err, isErr } from '../core/result.js';\nimport { validationError } from '../core/errors.js';\nimport type { ClosedWire, PlanarWire, ValidSolid } from '../core/shapeTypes.js';\nimport { makeCircle, assembleWire } from '../topology/curveBuilders.js';\nimport { makeFace } from '../topology/surfaceBuilders.js';\nimport { DisposalScope } from '../core/disposal.js';\nimport { rotate, translate } from '../topology/transformFns.js';\nimport { cut } from '../topology/booleanFns.js';\nimport { extrude } from '../operations/extrudeFns.js';\nimport {\n  type GearDiagnostic,\n  DEFAULT_CLEARANCE,\n  DEFAULT_PRESSURE_ANGLE_DEG,\n  backlashHalf,\n  evenToothPhaseOffset,\n  externalExternalContactRatio,\n  externalInternalContactRatio,\n  filletStressConcentrationFactor,\n  gearGeometry,\n  lewisRootStressCorrected,\n  planetSelfRotationAngle,\n  ringTeeth,\n  solvePlanetRingWorkingPressureAngle,\n  solveSunPlanetWorkingPressureAngle,\n  undercutDeficit,\n  validatePlanetary,\n  workingCenterDistance,\n} from './gearMath.js';\nimport { makeExternalGearProfileWire, makeInternalGearProfileWire } from './gearProfile.js';\n\nexport interface ExternalGearParams {\n  teeth: number;\n  moduleSize: number;\n  thickness: number;\n  pressureAngleDeg?: number;\n  shift?: number;\n  clearance?: number;\n  /**\n   * Per-gear flank thinning (mm); 0 = theoretical involute. For an equal-pair\n   * mesh, total mesh backlash is 2× this value.\n   */\n  flankThinning?: number;\n  /** Diameter (mm) of central bore; 0 or omitted = no bore. */\n  bore?: number;\n  /**\n   * Sample count per involute flank. Lower = faster mesh, coarser surface.\n   * Defaults to `adaptiveSampleCount(moduleSize)` (≈ 16 at module 2). Try\n   * 8 for previews, 24+ for 3D printing.\n   */\n  samples?: number;\n}\n\nexport interface InternalGearParams {\n  teeth: number;\n  moduleSize: number;\n  thickness: number;\n  pressureAngleDeg?: number;\n  shift?: number;\n  clearance?: number;\n  flankThinning?: number;\n  /** Wall thickness from pitch radius outward; defaults to 2·moduleSize. */\n  ringWallThickness?: number;\n  /** See {@link ExternalGearParams.samples}. */\n  samples?: number;\n}\n\nexport interface PlanetaryGearParams {\n  /** Extrusion thickness (mm). */\n  thickness: number;\n  moduleSize?: number;\n  sunTeeth?: number;\n  planetTeeth?: number;\n  numPlanets?: number;\n  pressureAngleDeg?: number;\n  clearance?: number;\n  /** Total mesh backlash (mm); split as b/2 per gear. */\n  backlash?: number;\n  sunShift?: number;\n  planetShift?: number;\n  ringShift?: number;\n  ringWallThickness?: number;\n  /** Diameter (mm) of central bore in the sun gear; 0 or omitted = no bore. */\n  sunBore?: number;\n  /** Diameter (mm) of central bore in each planet gear; 0 or omitted = no bore. */\n  planetBore?: number;\n  /**\n   * Applied torque on the SUN (input) shaft, in N·m. Planet and ring stresses\n   * are derived via force balance (shared tangential force at the mesh).\n   * When supplied, `lewisStress` and `stressConcentrationFactor` are computed.\n   */\n  appliedTorque?: number;\n  /** See {@link ExternalGearParams.samples}. Applied to sun, planet, and ring. */\n  samples?: number;\n}\n\nexport interface GearResult {\n  solid: ValidSolid;\n  pitchDiameter: number;\n  baseDiameter: number;\n  tipDiameter: number;\n  rootDiameter: number;\n  /**\n   * Diagnostics specific to this gear in isolation (e.g. undercut risk for\n   * low-tooth-count gears with no profile shift). Empty when the gear is\n   * geometrically clean. Planetary assemblies emit additional cross-mesh\n   * diagnostics on `PlanetaryGearAssembly.diagnostics`.\n   */\n  diagnostics: GearDiagnostic[];\n}\n\nexport interface PlanetaryGearAssembly {\n  sun: ValidSolid;\n  planets: ValidSolid[];\n  ring: ValidSolid;\n  ringTeeth: number;\n  /** Working pressure angle (radians) under the supplied profile shifts. */\n  workingPressureAngle: number;\n  /** True (working) center distance between sun and planet axes (mm). */\n  centerDistance: number;\n  /** Transverse contact ratios per mesh; ≥ 1.2 is industry-acceptable. */\n  contactRatio: { sunPlanet: number; planetRing: number };\n  /**\n   * Additional profile shift each gear needs to clear the undercut threshold.\n   * Positive = undercut risk; zero = clear.\n   */\n  undercutDeficit: { sun: number; planet: number };\n  /**\n   * Root bending stress (MPa) at each gear, with the fillet stress\n   * concentration factor applied: `σ = σ_Lewis · K_f`. Only present when\n   * `appliedTorque` was supplied. See {@link stressConcentrationFactor} for\n   * the K_f values used.\n   */\n  lewisStress?: { sun: number; planet: number; ring: number };\n  /**\n   * Dolan-Broghamer stress concentration factor `K_f` at each gear's root\n   * fillet, the multiplier already folded into {@link lewisStress}. Exposed\n   * so callers can recover the raw Lewis stress (σ_Lewis = lewisStress / K_f)\n   * or compare against ISO 6336-3 `Y_F·Y_S` from an external check. The ring\n   * value uses the Niemann internal-gear reduction (concave fillet).\n   */\n  stressConcentrationFactor?: { sun: number; planet: number; ring: number };\n  diagnostics: GearDiagnostic[];\n}\n\nexport function makeExternalGear(params: ExternalGearParams): Result<GearResult> {\n  const {\n    teeth,\n    moduleSize,\n    thickness,\n    pressureAngleDeg = DEFAULT_PRESSURE_ANGLE_DEG,\n    shift = 0,\n    clearance = DEFAULT_CLEARANCE,\n    flankThinning = 0,\n    bore = 0,\n    samples,\n  } = params;\n  if (thickness <= 0)\n    return err(validationError('GEAR_THICKNESS_NONPOSITIVE', 'thickness must be > 0'));\n  if (samples !== undefined && (!Number.isInteger(samples) || samples < 1))\n    return err(\n      validationError(\n        'GEAR_SAMPLES_INVALID',\n        `samples must be a positive integer; got ${String(samples)}`\n      )\n    );\n\n  const alpha = (pressureAngleDeg * Math.PI) / 180;\n  const geom = gearGeometry(teeth, moduleSize, alpha, shift, clearance, flankThinning, false);\n  if (bore > 0 && bore >= 2 * geom.rRoot)\n    return err(\n      validationError(\n        'GEAR_BORE_TOO_LARGE',\n        `bore diameter ${bore.toFixed(2)} ≥ root diameter ${(2 * geom.rRoot).toFixed(2)} — would erase the gear teeth`\n      )\n    );\n\n  const wireResult = makeExternalGearProfileWire({\n    teeth,\n    moduleSize,\n    pressureAngle: alpha,\n    shift,\n    clearance,\n    flankThinning,\n    ...(samples !== undefined ? { samples } : {}),\n  });\n  if (isErr(wireResult)) return wireResult;\n  const diagnostics = externalGearDiagnostics(teeth, alpha, shift);\n  return finalizeExternalSolid(wireResult.value, thickness, bore, geom, diagnostics);\n}\n\nexport function makeInternalGear(params: InternalGearParams): Result<GearResult> {\n  const {\n    teeth,\n    moduleSize,\n    thickness,\n    pressureAngleDeg = DEFAULT_PRESSURE_ANGLE_DEG,\n    shift = 0,\n    clearance = DEFAULT_CLEARANCE,\n    flankThinning = 0,\n    ringWallThickness = 2 * params.moduleSize,\n    samples,\n  } = params;\n  if (thickness <= 0)\n    return err(validationError('GEAR_THICKNESS_NONPOSITIVE', 'thickness must be > 0'));\n  if (ringWallThickness <= 0)\n    return err(validationError('GEAR_WALL_NONPOSITIVE', 'ringWallThickness must be > 0'));\n  if (samples !== undefined && (!Number.isInteger(samples) || samples < 1))\n    return err(\n      validationError(\n        'GEAR_SAMPLES_INVALID',\n        `samples must be a positive integer; got ${String(samples)}`\n      )\n    );\n\n  const alpha = (pressureAngleDeg * Math.PI) / 180;\n  const innerWireResult = makeInternalGearProfileWire({\n    teeth,\n    moduleSize,\n    pressureAngle: alpha,\n    shift,\n    clearance,\n    flankThinning,\n    ...(samples !== undefined ? { samples } : {}),\n  });\n  if (isErr(innerWireResult)) return innerWireResult;\n\n  const geom = gearGeometry(teeth, moduleSize, alpha, shift, clearance, flankThinning, true);\n  const outerRadius = geom.rPitch + ringWallThickness;\n  const outerWireResult = makeOuterCircleWire(outerRadius);\n  if (isErr(outerWireResult)) return outerWireResult;\n\n  // Internal gears use a different undercut criterion (involute-cutter geometry,\n  // shift sign inverted), so the rack-cut formula that drives external-gear\n  // diagnostics doesn't apply here. Leave diagnostics empty until a ring-specific\n  // check exists.\n  return finalizeInternalSolid(outerWireResult.value, innerWireResult.value, thickness, geom, []);\n}\n\n/** Diagnostics for an external (rack-cut) gear in isolation. Not valid for internal gears. */\nfunction externalGearDiagnostics(teeth: number, alpha: number, shift: number): GearDiagnostic[] {\n  const deficit = undercutDeficit(teeth, alpha, shift);\n  if (deficit <= 0) return [];\n  return [\n    {\n      code: 'UNDERCUT_RISK',\n      severity: 'warning',\n      message: `gear is undercut: increase shift by ${deficit.toFixed(3)} to avoid (z=${teeth})`,\n      context: { deficit, teeth },\n    },\n  ];\n}\n\nexport function makePlanetaryGear(params: PlanetaryGearParams): Result<PlanetaryGearAssembly> {\n  const resolved = resolvePlanetaryParams(params);\n  if (isErr(resolved)) return resolved;\n  const cfg = resolved.value;\n\n  const stages = buildPlanetaryStages(cfg);\n  if (isErr(stages)) return stages;\n  const { sun, planet, ring } = stages.value;\n\n  // Compute metrics before disposing `planet.solid` — the function only reads\n  // scalar fields off GearResult so the order is technically irrelevant today,\n  // but it leaves a footgun for anyone adding a solid-touching read later.\n  const metrics = computeMeshMetrics(cfg, sun, planet, ring);\n  const planets = placePlanets(planet.solid, cfg);\n  planet.solid.delete();\n  const ringPhased = applyRingPhase(ring.solid, cfg.zr);\n  const diagnostics = collectDiagnostics(cfg, metrics);\n\n  return ok({\n    sun: sun.solid,\n    planets,\n    ring: ringPhased,\n    ringTeeth: cfg.zr,\n    workingPressureAngle: cfg.alphaW_sp,\n    centerDistance: cfg.centerDistance,\n    contactRatio: { sunPlanet: metrics.crSunPlanet, planetRing: metrics.crPlanetRing },\n    undercutDeficit: { sun: metrics.undercutSun, planet: metrics.undercutPlanet },\n    ...(metrics.lewisStress ? { lewisStress: metrics.lewisStress } : {}),\n    ...(metrics.stressConcentrationFactor\n      ? { stressConcentrationFactor: metrics.stressConcentrationFactor }\n      : {}),\n    diagnostics,\n  });\n}\n\nfunction buildPlanetaryStages(\n  cfg: ResolvedPlanetary\n): Result<{ sun: GearResult; planet: GearResult; ring: GearResult }> {\n  const samplesField = cfg.samples !== undefined ? { samples: cfg.samples } : {};\n  const common = {\n    moduleSize: cfg.moduleSize,\n    thickness: cfg.thickness,\n    pressureAngleDeg: cfg.pressureAngleDeg,\n    clearance: cfg.clearance,\n    flankThinning: cfg.bHalf,\n    ...samplesField,\n  };\n  const sun = makeExternalGear({\n    ...common,\n    teeth: cfg.sunTeeth,\n    shift: cfg.sunShift,\n    bore: cfg.sunBore,\n  });\n  if (isErr(sun)) return sun;\n  const planet = makeExternalGear({\n    ...common,\n    teeth: cfg.planetTeeth,\n    shift: cfg.planetShift,\n    bore: cfg.planetBore,\n  });\n  if (isErr(planet)) {\n    sun.value.solid.delete();\n    return planet;\n  }\n  const ring = makeInternalGear({\n    ...common,\n    teeth: cfg.zr,\n    shift: cfg.ringShift,\n    ringWallThickness: cfg.ringWallThickness,\n  });\n  if (isErr(ring)) {\n    sun.value.solid.delete();\n    planet.value.solid.delete();\n    return ring;\n  }\n  return ok({ sun: sun.value, planet: planet.value, ring: ring.value });\n}\n\ninterface ResolvedPlanetary {\n  moduleSize: number;\n  sunTeeth: number;\n  planetTeeth: number;\n  numPlanets: number;\n  pressureAngleDeg: number;\n  alpha: number;\n  /** Working PA for the sun-planet mesh. */\n  alphaW_sp: number;\n  /** Working PA for the planet-ring mesh; equals alphaW_sp iff xr = xs + 2·xp. */\n  alphaW_pr: number;\n  clearance: number;\n  bHalf: number;\n  sunShift: number;\n  planetShift: number;\n  ringShift: number;\n  ringWallThickness: number;\n  thickness: number;\n  sunBore: number;\n  planetBore: number;\n  appliedTorque?: number;\n  samples?: number;\n  zr: number;\n  centerDistance: number;\n}\n\nfunction resolvePlanetaryParams(params: PlanetaryGearParams): Result<ResolvedPlanetary> {\n  const moduleSize = params.moduleSize ?? 3;\n  const sunTeeth = params.sunTeeth ?? 15;\n  const planetTeeth = params.planetTeeth ?? 12;\n  const numPlanets = params.numPlanets ?? 3;\n  const pressureAngleDeg = params.pressureAngleDeg ?? DEFAULT_PRESSURE_ANGLE_DEG;\n  const alpha = (pressureAngleDeg * Math.PI) / 180;\n  const sunShift = params.sunShift ?? 0;\n  const planetShift = params.planetShift ?? 0;\n  const ringShift = params.ringShift ?? 0;\n\n  if (params.thickness <= 0)\n    return err(validationError('GEAR_THICKNESS_NONPOSITIVE', 'thickness must be > 0'));\n\n  const validation = validatePlanetary(sunTeeth, planetTeeth, numPlanets, planetShift);\n  if (isErr(validation)) return validation;\n\n  const zr = ringTeeth(sunTeeth, planetTeeth);\n\n  const sp = solveSunPlanetWorkingPressureAngle(\n    alpha,\n    sunShift,\n    planetShift,\n    sunTeeth,\n    planetTeeth\n  );\n  if (isErr(sp)) return sp;\n  const pr = solvePlanetRingWorkingPressureAngle(alpha, planetShift, ringShift, planetTeeth, zr);\n  if (isErr(pr)) return pr;\n  const centerDistance = workingCenterDistance(sunTeeth, planetTeeth, moduleSize, alpha, sp.value);\n\n  return ok({\n    moduleSize,\n    sunTeeth,\n    planetTeeth,\n    numPlanets,\n    pressureAngleDeg,\n    alpha,\n    alphaW_sp: sp.value,\n    alphaW_pr: pr.value,\n    clearance: params.clearance ?? DEFAULT_CLEARANCE,\n    bHalf: backlashHalf(params.backlash ?? 0),\n    sunShift,\n    planetShift,\n    ringShift,\n    ringWallThickness: params.ringWallThickness ?? 2 * moduleSize,\n    thickness: params.thickness,\n    sunBore: params.sunBore ?? 0,\n    planetBore: params.planetBore ?? 0,\n    ...(params.appliedTorque !== undefined ? { appliedTorque: params.appliedTorque } : {}),\n    ...(params.samples !== undefined ? { samples: params.samples } : {}),\n    zr,\n    centerDistance,\n  });\n}\n\nfunction makeOuterCircleWire(radius: number): Result<ClosedWire & PlanarWire> {\n  using scope = new DisposalScope();\n  const circleEdge = scope.register(makeCircle(radius, [0, 0, 0], [0, 0, 1]));\n  const wireResult = assembleWire([circleEdge]);\n  if (isErr(wireResult)) return wireResult;\n  // A single-edge circle wire is closed and planar by construction; brand once.\n  return ok(wireResult.value as ClosedWire & PlanarWire);\n}\n\nfunction finalizeExternalSolid(\n  wire: ClosedWire & PlanarWire,\n  thickness: number,\n  bore: number,\n  geom: ReturnType<typeof gearGeometry>,\n  diagnostics: GearDiagnostic[]\n): Result<GearResult> {\n  // The profile wire and its face are consumed by makeFace/extrude (which share\n  // their refcounted TShape); dispose them once the solid is built. The solid is\n  // NOT registered here — it is the returned value on the no-bore path.\n  using profileScope = new DisposalScope();\n  profileScope.register(wire);\n  const faceResult = makeFace(wire);\n  if (isErr(faceResult)) return faceResult;\n  profileScope.register(faceResult.value);\n  const solidResult = extrude(faceResult.value, [0, 0, thickness]);\n  if (isErr(solidResult)) return solidResult;\n  if (bore <= 0) return ok(buildGearResult(solidResult.value, geom, diagnostics));\n\n  // Bore overshoots both ends (z = -0.5 to thickness + 0.5) so no faces are coplanar.\n  using boreScope = new DisposalScope();\n  boreScope.register(solidResult.value); // un-bored gear, consumed by cut below\n  const boreFaceResult = makeBoreFace(bore / 2);\n  if (isErr(boreFaceResult)) return boreFaceResult;\n  boreScope.register(boreFaceResult.value);\n  const boreRaw = extrude(boreFaceResult.value, [0, 0, thickness + 1]);\n  if (isErr(boreRaw)) return boreRaw;\n  boreScope.register(boreRaw.value);\n  const boreSolid = boreScope.register(translate(boreRaw.value, [0, 0, -0.5]));\n  const cutResult = cut(solidResult.value, boreSolid);\n  if (isErr(cutResult)) return cutResult;\n  return ok(buildGearResult(cutResult.value, geom, diagnostics));\n}\n\nfunction makeBoreFace(radius: number): Result<Parameters<typeof extrude>[0]> {\n  const wireResult = makeOuterCircleWire(radius);\n  if (isErr(wireResult)) return wireResult;\n  return makeFace(wireResult.value);\n}\n\nfunction finalizeInternalSolid(\n  outerWire: ClosedWire & PlanarWire,\n  innerToothedWire: ClosedWire & PlanarWire,\n  thickness: number,\n  geom: ReturnType<typeof gearGeometry>,\n  diagnostics: GearDiagnostic[]\n): Result<GearResult> {\n  const faceResult = makeFace(outerWire, [innerToothedWire]);\n  if (isErr(faceResult)) return faceResult;\n  const solidResult = extrude(faceResult.value, [0, 0, thickness]);\n  if (isErr(solidResult)) return solidResult;\n  return ok(buildGearResult(solidResult.value, geom, diagnostics));\n}\n\nfunction buildGearResult(\n  solid: ValidSolid,\n  geom: ReturnType<typeof gearGeometry>,\n  diagnostics: GearDiagnostic[] = []\n): GearResult {\n  return {\n    solid,\n    diagnostics,\n    pitchDiameter: 2 * geom.rPitch,\n    baseDiameter: 2 * geom.rb,\n    tipDiameter: 2 * geom.rTip,\n    rootDiameter: 2 * geom.rRoot,\n  };\n}\n\nfunction placePlanets(planetProto: ValidSolid, cfg: ResolvedPlanetary): ValidSolid[] {\n  const planets: ValidSolid[] = [];\n  try {\n    for (let i = 0; i < cfg.numPlanets; i++) {\n      const orbital = (i * 2 * Math.PI) / cfg.numPlanets;\n      const selfRot = planetSelfRotationAngle(orbital, cfg.sunTeeth, cfg.planetTeeth);\n      const rotated = rotate(planetProto, (selfRot * 180) / Math.PI);\n      try {\n        const offset: Vec3 = [\n          cfg.centerDistance * Math.cos(orbital),\n          cfg.centerDistance * Math.sin(orbital),\n          0,\n        ];\n        planets.push(translate(rotated, offset));\n      } finally {\n        rotated.delete();\n      }\n    }\n    return planets;\n  } catch (e) {\n    // Mid-loop failure: clean up the planets that did succeed before rethrowing.\n    for (const p of planets) p.delete();\n    throw e;\n  }\n}\n\nfunction applyRingPhase(ring: ValidSolid, zr: number): ValidSolid {\n  const phaseRad = evenToothPhaseOffset(zr);\n  if (phaseRad === 0) return ring;\n  const rotated = rotate(ring, (phaseRad * 180) / Math.PI);\n  ring.delete();\n  return rotated;\n}\n\ninterface MeshMetrics {\n  crSunPlanet: number;\n  crPlanetRing: number;\n  undercutSun: number;\n  undercutPlanet: number;\n  lewisStress?: { sun: number; planet: number; ring: number };\n  stressConcentrationFactor?: { sun: number; planet: number; ring: number };\n}\n\nfunction computeMeshMetrics(\n  cfg: ResolvedPlanetary,\n  sun: GearResult,\n  planet: GearResult,\n  ring: GearResult\n): MeshMetrics {\n  const crSunPlanet = externalExternalContactRatio(\n    sun.tipDiameter / 2,\n    sun.baseDiameter / 2,\n    planet.tipDiameter / 2,\n    planet.baseDiameter / 2,\n    cfg.centerDistance,\n    cfg.moduleSize,\n    cfg.alpha,\n    cfg.alphaW_sp\n  );\n  const crPlanetRing = externalInternalContactRatio(\n    planet.tipDiameter / 2,\n    planet.baseDiameter / 2,\n    ring.tipDiameter / 2,\n    ring.baseDiameter / 2,\n    cfg.centerDistance,\n    cfg.moduleSize,\n    cfg.alpha,\n    cfg.alphaW_pr\n  );\n  const undercutSun = undercutDeficit(cfg.sunTeeth, cfg.alpha, cfg.sunShift);\n  const undercutPlanet = undercutDeficit(cfg.planetTeeth, cfg.alpha, cfg.planetShift);\n\n  const metrics: MeshMetrics = { crSunPlanet, crPlanetRing, undercutSun, undercutPlanet };\n  if (cfg.appliedTorque !== undefined) {\n    // Force balance: W_t is shared at each mesh, so each gear's own-pitch torque\n    // is T_eff = T_sun · z / z_sun (preserves W_t through Lewis's 2T/(z·m) term).\n    const tSun = cfg.appliedTorque;\n    metrics.lewisStress = {\n      sun: lewisRootStressCorrected(tSun, cfg.moduleSize, cfg.thickness, cfg.sunTeeth, cfg.alpha),\n      planet: lewisRootStressCorrected(\n        (tSun * cfg.planetTeeth) / cfg.sunTeeth,\n        cfg.moduleSize,\n        cfg.thickness,\n        cfg.planetTeeth,\n        cfg.alpha\n      ),\n      ring: lewisRootStressCorrected(\n        (tSun * cfg.zr) / cfg.sunTeeth,\n        cfg.moduleSize,\n        cfg.thickness,\n        cfg.zr,\n        cfg.alpha,\n        true\n      ),\n    };\n    metrics.stressConcentrationFactor = {\n      sun: filletStressConcentrationFactor(cfg.sunTeeth, cfg.alpha),\n      planet: filletStressConcentrationFactor(cfg.planetTeeth, cfg.alpha),\n      ring: filletStressConcentrationFactor(cfg.zr, cfg.alpha, true),\n    };\n  }\n  return metrics;\n}\n\nfunction collectDiagnostics(cfg: ResolvedPlanetary, metrics: MeshMetrics): GearDiagnostic[] {\n  const diagnostics: GearDiagnostic[] = [];\n  if (metrics.crSunPlanet < 1.2) {\n    diagnostics.push({\n      code: 'CONTACT_RATIO_LOW_SUN_PLANET',\n      severity: 'warning',\n      message: `sun-planet contact ratio ${metrics.crSunPlanet.toFixed(2)} is below 1.2 — may run unevenly`,\n      context: { value: metrics.crSunPlanet },\n    });\n  }\n  if (metrics.crPlanetRing < 1.2) {\n    diagnostics.push({\n      code: 'CONTACT_RATIO_LOW_PLANET_RING',\n      severity: 'warning',\n      message: `planet-ring contact ratio ${metrics.crPlanetRing.toFixed(2)} is below 1.2 — may run unevenly`,\n      context: { value: metrics.crPlanetRing },\n    });\n  }\n  if (metrics.undercutSun > 0) {\n    diagnostics.push({\n      code: 'UNDERCUT_RISK_SUN',\n      severity: 'warning',\n      message: `sun gear is undercut: increase sunShift by ${metrics.undercutSun.toFixed(3)} to avoid`,\n      context: { deficit: metrics.undercutSun, sunTeeth: cfg.sunTeeth },\n    });\n  }\n  if (metrics.undercutPlanet > 0) {\n    diagnostics.push({\n      code: 'UNDERCUT_RISK_PLANET',\n      severity: 'warning',\n      message: `planet gear is undercut: increase planetShift by ${metrics.undercutPlanet.toFixed(3)} to avoid`,\n      context: { deficit: metrics.undercutPlanet, planetTeeth: cfg.planetTeeth },\n    });\n  }\n  // Only emit when shifts are large enough to matter. The diagnostic exists to\n  // warn that Lewis Y(z) uses unshifted geometry — at ~0.05 of profile shift\n  // the error is under 2.5%, well within engineering noise for a strength\n  // estimate, so don't pester users on canonical setups with small shifts.\n  const totalShiftMagnitude =\n    Math.abs(cfg.sunShift) + Math.abs(cfg.planetShift) + Math.abs(cfg.ringShift);\n  if (cfg.appliedTorque !== undefined && totalShiftMagnitude > 0.05) {\n    diagnostics.push({\n      code: 'LEWIS_Y_SHIFT_UNCORRECTED',\n      severity: 'info',\n      message:\n        'Lewis stress uses unshifted Y(z) and K_f approximations; expect ±5% per 0.1 of profile shift',\n    });\n  }\n\n  // Kinematic compatibility: x_ring should equal x_sun + 2·x_planet for both meshes to share αw.\n  const kinematicError = cfg.ringShift - (cfg.sunShift + 2 * cfg.planetShift);\n  if (Math.abs(kinematicError) > 1e-6) {\n    diagnostics.push({\n      code: 'PLANETARY_SHIFT_KINEMATIC_MISMATCH',\n      severity: 'warning',\n      message: `ringShift should equal sunShift + 2·planetShift for both meshes to share working PA; off by ${kinematicError.toFixed(3)}`,\n      context: {\n        kinematicError,\n        sunShift: cfg.sunShift,\n        planetShift: cfg.planetShift,\n        ringShift: cfg.ringShift,\n      },\n    });\n  }\n  return diagnostics;\n}\n"]}

package/dist/vendored/brepjs/gear/gearMath.d.ts

@@ -0,0 +1,89 @@
+import type { Vec3 } from '../core/types.js';
+import { type Result } from '../core/result.js';
+export declare const DEFAULT_PRESSURE_ANGLE_DEG = 20;
+export declare const DEFAULT_CLEARANCE = 0.25;
+export type GearDiagnosticSeverity = 'warning' | 'info';
+export type GearDiagnosticCode = 'CONTACT_RATIO_LOW_SUN_PLANET' | 'CONTACT_RATIO_LOW_PLANET_RING' | 'UNDERCUT_RISK' | 'UNDERCUT_RISK_SUN' | 'UNDERCUT_RISK_PLANET' | 'LEWIS_Y_SHIFT_UNCORRECTED' | 'PLANETARY_SHIFT_KINEMATIC_MISMATCH';
+export interface GearDiagnostic {
+    code: GearDiagnosticCode;
+    severity: GearDiagnosticSeverity;
+    message: string;
+    context?: Record<string, number | string>;
+}
+export declare const inv: (alpha: number) => number;
+export declare function involutePoint(rb: number, alpha: number, theta0: number, sign: 1 | -1): Vec3;
+export declare function cosineSpaceFlankSamples(rb: number, alphaMax: number, theta0: number, count: number, sign: 1 | -1): Vec3[];
+export declare function adaptiveSampleCount(moduleSize: number): number;
+export interface GearGeometry {
+    rPitch: number;
+    rb: number;
+    /** Outer radius for external, inner for internal. */
+    rTip: number;
+    /** Inner radius for external, outer for internal. */
+    rRoot: number;
+    alphaPitch: number;
+    /** Half-tooth angular width at the pitch circle, with shift and flank thinning applied. */
+    halfToothAngle: number;
+    /** Involute parameter α at the radius the flank reaches (rTip ext, rRoot int). */
+    alphaTip: number;
+    toothPitch: number;
+    isInternal: boolean;
+}
+export declare function gearGeometry(z: number, moduleSize: number, alpha: number, shift: number, clearance: number, backlashHalf: number, isInternal: boolean): GearGeometry;
+/**
+ * Solve `inv(αw) = inv(α) + 2·(summedShift)·tan α / totalTeeth` by bisection.
+ *
+ * For external-external mesh: summedShift = x1 + x2, totalTeeth = z1 + z2.
+ * For external-internal mesh: summedShift = xInternal − xExternal, totalTeeth = zInternal − zExternal.
+ *
+ * Returns Err for shifts that push αw outside (epsilon, π/2).
+ */
+export declare function solveWorkingPressureAngle(alpha: number, summedShift: number, totalTeeth: number): Result<number>;
+/** Convenience: working PA for an external-external mesh (sun-planet). */
+export declare function solveSunPlanetWorkingPressureAngle(alpha: number, xSun: number, xPlanet: number, zSun: number, zPlanet: number): Result<number>;
+/** Convenience: working PA for an external-internal mesh (planet-ring). */
+export declare function solvePlanetRingWorkingPressureAngle(alpha: number, xPlanet: number, xRing: number, zPlanet: number, zRing: number): Result<number>;
+export declare function workingCenterDistance(zs: number, zp: number, moduleSize: number, alpha: number, alphaW: number): number;
+export declare function validatePlanetary(zs: number, zp: number, n: number, planetShift: number): Result<void>;
+export declare function externalExternalContactRatio(ra1: number, rb1: number, ra2: number, rb2: number, centerDistance: number, moduleSize: number, alpha: number, alphaW: number): number;
+export declare function externalInternalContactRatio(ra_p: number, rb_p: number, ra_r: number, rb_r: number, centerDistance: number, moduleSize: number, alpha: number, alphaW: number): number;
+export declare function undercutMinimumShift(z: number, alpha: number): number;
+export declare function undercutDeficit(z: number, alpha: number, shift: number): number;
+export declare function lewisYFactor(z: number): number;
+export declare function lewisRootStress(appliedTorqueNm: number, moduleSize: number, faceWidth: number, z: number): number;
+/**
+ * Dolan-Broghamer fit is for *external* spur teeth. Set `isInternal=true` on
+ * ring gears to apply Niemann's 0.85× reduction — the concave root fillet has
+ * a gentler geometric riser (matches AGMA 908-B89 trend).
+ */
+export declare function filletStressConcentrationFactor(z: number, alpha: number, isInternal?: boolean): number;
+export declare function lewisRootStressCorrected(appliedTorqueNm: number, moduleSize: number, faceWidth: number, z: number, alpha: number, isInternal?: boolean): number;
+export declare function backlashHalf(totalBacklash: number): number;
+export declare function ringTeeth(zs: number, zp: number): number;
+export declare function evenToothPhaseOffset(z: number): number;
+export declare function planetSelfRotationAngle(orbitalAngle: number, zs: number, zp: number): number;
+/** Pure-kinematics inputs for {@link planetPlacements}; subset of `PlanetaryGearParams`. */
+export interface PlanetPlacementParams {
+    moduleSize?: number;
+    sunTeeth?: number;
+    planetTeeth?: number;
+    numPlanets?: number;
+    pressureAngleDeg?: number;
+    sunShift?: number;
+    planetShift?: number;
+}
+export interface PlanetPlacement {
+    /** Self-rotation around the planet's own Z axis, in degrees. */
+    rotationDeg: number;
+    /** Center of the planet in the assembly frame: `[x, y, 0]`. */
+    position: Vec3;
+}
+/**
+ * Compute the rigid placement of each planet without building any solid.
+ *
+ * Pairs with one materialized planet from `makeExternalGear` to support
+ * GPU-instanced rendering: mesh once, draw N times under these transforms.
+ * Defaults match {@link PlanetaryGearParams}.
+ */
+export declare function planetPlacements(params?: PlanetPlacementParams): Result<PlanetPlacement[]>;
+//# sourceMappingURL=gearMath.d.ts.map