@theerapat-s28/ts-structural-eng-tools 3.3.0 → 3.4.0

package/dist/fea/analysis/post/wall-forces.js

@@ -0,0 +1,305 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.presentSpandrelForces = exports.presentPierForces = exports.prepareWalls = exports.WALL_COMPONENTS = void 0;
+const fea_constant_1 = require("../../../core/constants/fea.constant");
+const member_diagram_1 = require("../../frame/post/member-diagram");
+const global_assembly_1 = require("../../shared/solver/global-assembly");
+const math_1 = require("../../../utils/math");
+const wall_geometry_1 = require("./wall-geometry");
+/** The resultants, in the order of the raw values. */
+exports.WALL_COMPONENTS = ["P", "V2", "V3", "T", "M2", "M3"];
+const FORCE_DECIMALS = 4;
+const LENGTH_DECIMALS = 3;
+const round = (value, decimals = FORCE_DECIMALS) => (0, math_1.roundToDecimalPlaces)(value, decimals) + 0; // + 0: no −0
+const roundPoint = (p) => [round(p[0], 3), round(p[1], 3), round(p[2], 3)];
+const range = (values) => ({ min: Math.min(...values), max: Math.max(...values) });
+/**
+ * One cut: the shell elements and member pieces on its side that reach it, the members that
+ * cross it, its section (a leg per area, a point per member) and the section's centroid.
+ */
+const prepareCut = (system, shape, e1, s, face, tol) => {
+    const negative = (0, wall_geometry_1.fromNegativeSide)(face);
+    const sign = negative ? 1 : -1;
+    const at = (n) => {
+        const node = system.nodes[n];
+        return [node.x, node.y, node.z];
+    };
+    const shells = [];
+    const legs = [];
+    for (const { areaId, shells: list } of shape.areas) {
+        const onCut = [];
+        let normal;
+        let thickness = 0;
+        for (const index of list) {
+            const shell = system.shells[index];
+            const points = shell.nodeIndices.map(at);
+            const r = range(points.map((p) => (0, wall_geometry_1.dot)(p, e1)));
+            if (r.max - r.min <= tol)
+                continue;
+            if (Math.abs((negative ? r.max : r.min) - s) > tol)
+                continue;
+            const corners = points.flatMap((p, i) => (Math.abs((0, wall_geometry_1.dot)(p, e1) - s) <= tol ? [i] : []));
+            shells.push({ shell: index, corners, sign });
+            onCut.push(...corners.map((i) => points[i]));
+            normal ?? (normal = shell.element.axes.lambda[2]);
+            thickness = shell.section.thickness;
+        }
+        if (!normal)
+            continue;
+        const d = (0, wall_geometry_1.cross)(e1, normal);
+        const dn = Math.hypot(d[0], d[1], d[2]);
+        if (dn === 0)
+            continue;
+        // legs run up (spandrel), else toward +X (+Y for a leg parallel to Y)
+        let along = [d[0] / dn, d[1] / dn, d[2] / dn];
+        const flip = Math.abs(along[2]) > fea_constant_1.VERTICAL_MEMBER_TOLERANCE
+            ? along[2] < 0
+            : Math.abs(along[0]) > fea_constant_1.VERTICAL_MEMBER_TOLERANCE
+                ? along[0] < 0
+                : along[1] < 0;
+        if (flip)
+            along = [-along[0], -along[1], -along[2]];
+        let start = onCut[0];
+        let end = onCut[0];
+        for (const p of onCut) {
+            if ((0, wall_geometry_1.dot)(p, along) < (0, wall_geometry_1.dot)(start, along))
+                start = p;
+            if ((0, wall_geometry_1.dot)(p, along) > (0, wall_geometry_1.dot)(end, along))
+                end = p;
+        }
+        const length = (0, wall_geometry_1.dot)((0, wall_geometry_1.sub)(end, start), along);
+        if (length > tol)
+            legs.push({ areaId, start, end, length, thickness });
+    }
+    const ends = [];
+    const crossings = [];
+    const members = [];
+    for (const { memberId, pieces } of shape.members) {
+        for (const piece of pieces) {
+            const member = system.members[piece];
+            const p0 = at(member.startIndex);
+            const p1 = at(member.endIndex);
+            const s0 = (0, wall_geometry_1.dot)(p0, e1);
+            const s1 = (0, wall_geometry_1.dot)(p1, e1);
+            if (Math.abs(s1 - s0) <= tol)
+                continue; // lies in the cut plane
+            const r = range([s0, s1]);
+            if (r.min < s - tol && r.max > s + tol) {
+                const t = (s - s0) / (s1 - s0);
+                const point = [
+                    p0[0] + t * (p1[0] - p0[0]),
+                    p0[1] + t * (p1[1] - p0[1]),
+                    p0[2] + t * (p1[2] - p0[2]),
+                ];
+                crossings.push({ piece, x: t * member.length, point, sign: s1 > s0 ? 1 : -1 });
+                members.push({ memberId, point, A: member.section.A });
+                break;
+            }
+            if (Math.abs((negative ? r.max : r.min) - s) <= tol) {
+                // the piece's end on the cut: the one with the larger s from below, the smaller from above
+                const end = negative === s1 > s0 ? 1 : 0;
+                const point = end === 1 ? p1 : p0;
+                ends.push({ piece, end, point, sign });
+                members.push({ memberId, point, A: member.section.A });
+                break;
+            }
+        }
+    }
+    // centroid: legs by length × thickness at their midpoints, members by A
+    let area = 0;
+    const sum = [0, 0, 0];
+    const addWeighted = (p, w) => {
+        area += w;
+        for (let k = 0; k < 3; k++)
+            sum[k] += w * p[k];
+    };
+    for (const leg of legs) {
+        const mid = [0, 1, 2].map((k) => (leg.start[k] + leg.end[k]) / 2);
+        addWeighted(mid, leg.length * leg.thickness);
+    }
+    for (const m of members)
+        addWeighted(m.point, m.A);
+    let centroid;
+    if (area > 0) {
+        centroid = sum.map((v) => v / area);
+    }
+    else {
+        // only corners touch the cut: the mean of the points on it
+        const points = [
+            ...shells.flatMap((t) => t.corners.map((i) => at(system.shells[t.shell].nodeIndices[i]))),
+            ...ends.map((t) => t.point),
+            ...crossings.map((t) => t.point),
+        ];
+        centroid = [0, 1, 2].map((k) => points.reduce((acc, p) => acc + p[k], 0) / Math.max(points.length, 1));
+    }
+    return {
+        face,
+        centroid,
+        section: {
+            legs: legs.map((leg) => ({
+                areaId: leg.areaId,
+                start: roundPoint(leg.start),
+                end: roundPoint(leg.end),
+                length: round(leg.length, LENGTH_DECIMALS),
+                thickness: leg.thickness,
+            })),
+            members: members.map((m) => ({ memberId: m.memberId, point: roundPoint(m.point), A: m.A })),
+            length: round(legs.reduce((acc, leg) => acc + leg.length, 0), LENGTH_DECIMALS),
+            area: round(area, LENGTH_DECIMALS),
+            centroid: roundPoint(centroid),
+        },
+        shells,
+        ends,
+        crossings,
+    };
+};
+/**
+ * Prepares the piers and spandrels of a model that passed `validateModel`: their axes, cuts
+ * (see `wallCutPositions`) and, for each cut, the shell elements and member pieces whose nodal
+ * forces make up its resultant.
+ *
+ * @param walls - The model's piers and spandrels
+ */
+const prepareWalls = (walls, system) => {
+    const at = (n) => {
+        const node = system.nodes[n];
+        return [node.x, node.y, node.z];
+    };
+    const prepare = (kind, wall) => {
+        const areas = (wall.areaIds ?? []).map((areaId) => ({
+            areaId,
+            shells: system.shellsByArea.get(areaId) ?? [],
+        }));
+        const members = (wall.memberIds ?? []).map((memberId) => ({
+            memberId,
+            pieces: system.memberGroups[system.memberIndex.get(memberId)].pieces,
+        }));
+        const shape = {
+            areas: areas.map(({ areaId, shells }) => ({
+                areaId,
+                elements: shells.map((s) => system.shells[s].nodeIndices.map(at)),
+            })),
+            members: members.map(({ memberId, pieces }) => ({
+                memberId,
+                pieces: pieces.map((m) => [
+                    at(system.members[m].startIndex),
+                    at(system.members[m].endIndex),
+                ]),
+            })),
+        };
+        const axes = (0, wall_geometry_1.wallAxes)(kind, shape, {
+            angle: wall.angle,
+            reverse: wall.reverse,
+        });
+        const e1 = axes[0];
+        const extent = (0, wall_geometry_1.wallExtent)(shape, e1);
+        const tol = (0, wall_geometry_1.wallTolerance)(extent);
+        const offset = kind === "pier" ? 0 : extent.min;
+        const extra = (wall.cuts ?? []).map((c) => c + offset);
+        const cuts = (0, wall_geometry_1.wallCutPositions)(kind, shape, e1, extra).map(({ s, face }) => ({
+            position: s - offset,
+            ...prepareCut(system, { areas, members }, e1, s, face, tol),
+        }));
+        return { id: wall.id, axes, length: extent.max - extent.min, cuts };
+    };
+    return {
+        piers: (walls.piers ?? []).map((p) => prepare("pier", p)),
+        spandrels: (walls.spandrels ?? []).map((s) => prepare("spandrel", s)),
+    };
+};
+exports.prepareWalls = prepareWalls;
+/**
+ * The unrounded global force (N) and moment about the cut's centroid (N·mm) on the +1 face of
+ * the negative part of one cut: the nodal forces `k·u − f` of the shell elements on the cut's
+ * side at their corners on it, the end forces of the member pieces ending on it from that side
+ * (both negated when the side is the positive one), and the internal forces of the members
+ * crossing it. Linear in the state.
+ */
+const cutRawForces = (system, state, cut) => {
+    const F = [0, 0, 0];
+    const M = [0, 0, 0];
+    const c = cut.centroid;
+    const addAt = (p, f, m, sign) => {
+        const r = (0, wall_geometry_1.sub)(p, c);
+        const rf = (0, wall_geometry_1.cross)(r, f);
+        for (let k = 0; k < 3; k++) {
+            F[k] += sign * f[k];
+            M[k] += sign * (m[k] + rf[k]);
+        }
+    };
+    for (const { shell: index, corners, sign } of cut.shells) {
+        const shell = system.shells[index];
+        const dofs = (0, global_assembly_1.elementDofs)(shell.nodeIndices);
+        const u = dofs.map((d) => state.U[d]);
+        const p = state.shellPressures[index];
+        for (const i of corners) {
+            const f = [0, 0, 0, 0, 0, 0];
+            for (let r = 0; r < 6; r++) {
+                const row = shell.element.kGlobal[6 * i + r];
+                let v = 0;
+                for (let j = 0; j < 24; j++)
+                    v += row[j] * u[j];
+                f[r] = v;
+            }
+            for (let k = 0; k < 3; k++)
+                f[k] -= p[k] * 1e-3 * shell.cornerAreas[i];
+            const node = system.nodes[shell.nodeIndices[i]];
+            addAt([node.x, node.y, node.z], [f[0], f[1], f[2]], [f[3], f[4], f[5]], sign);
+        }
+    }
+    const toGlobal = (lambda, local) => [
+        lambda[0][0] * local[0] + lambda[1][0] * local[1] + lambda[2][0] * local[2],
+        lambda[0][1] * local[0] + lambda[1][1] * local[1] + lambda[2][1] * local[2],
+        lambda[0][2] * local[0] + lambda[1][2] * local[1] + lambda[2][2] * local[2],
+    ];
+    for (const { piece, end, point, sign } of cut.ends) {
+        const { lambda } = system.members[piece];
+        const f = state.memberForces[piece].slice(6 * end, 6 * end + 6);
+        addAt(point, toGlobal(lambda, f.slice(0, 3)), toGlobal(lambda, f.slice(3, 6)), sign);
+    }
+    for (const { piece, x, point, sign } of cut.crossings) {
+        const { lambda } = system.members[piece];
+        const f = (0, member_diagram_1.memberInternalForces)(state.memberForces[piece].slice(0, 6), state.memberLoads[piece], x);
+        addAt(point, toGlobal(lambda, f.slice(0, 3)), toGlobal(lambda, f.slice(3, 6)), sign);
+    }
+    return { F, M };
+};
+/** The presented resultant of one cut in the wall's axes (kN, kN·m). */
+const cutResultants = (system, state, axes, cut) => {
+    const { F, M } = cutRawForces(system, state, cut);
+    return {
+        P: round((0, wall_geometry_1.dot)(F, axes[0]) * 1e-3),
+        V2: round((0, wall_geometry_1.dot)(F, axes[1]) * 1e-3),
+        V3: round((0, wall_geometry_1.dot)(F, axes[2]) * 1e-3),
+        T: round((0, wall_geometry_1.dot)(M, axes[0]) * 1e-6),
+        M2: round((0, wall_geometry_1.dot)(M, axes[1]) * 1e-6),
+        M3: round((0, wall_geometry_1.dot)(M, axes[2]) * 1e-6),
+    };
+};
+const presentAxes = (axes) => axes.map((e) => e.map((v) => round(v, 9)));
+/** The presented forces of one pier in a raw state. */
+const presentPierForces = (system, state, pier) => ({
+    pierId: pier.id,
+    axes: presentAxes(pier.axes),
+    cuts: pier.cuts.map((cut) => ({
+        z: round(cut.position, LENGTH_DECIMALS),
+        face: cut.face,
+        ...cutResultants(system, state, pier.axes, cut),
+        section: cut.section,
+    })),
+});
+exports.presentPierForces = presentPierForces;
+/** The presented forces of one spandrel in a raw state. */
+const presentSpandrelForces = (system, state, spandrel) => ({
+    spandrelId: spandrel.id,
+    axes: presentAxes(spandrel.axes),
+    length: round(spandrel.length, LENGTH_DECIMALS),
+    cuts: spandrel.cuts.map((cut) => ({
+        x: round(cut.position, LENGTH_DECIMALS),
+        face: cut.face,
+        ...cutResultants(system, state, spandrel.axes, cut),
+        section: cut.section,
+    })),
+});
+exports.presentSpandrelForces = presentSpandrelForces;
+//# sourceMappingURL=wall-forces.js.map

package/dist/fea/analysis/post/wall-forces.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"wall-forces.js","sourceRoot":"","sources":["../../../../src/fea/analysis/post/wall-forces.ts"],"names":[],"mappings":";;;AAAA,mEAA6E;AAC7E,uEAA0E;AAE1E,4EAAqE;AACrE,0CAAuD;AAYvD,mDAcyB;AAEzB,sDAAsD;AACzC,QAAA,eAAe,GAAkC,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;AAEjG,MAAM,cAAc,GAAG,CAAC,CAAC;AACzB,MAAM,eAAe,GAAG,CAAC,CAAC;AA6C1B,MAAM,KAAK,GAAG,CAAC,KAAa,EAAE,QAAQ,GAAG,cAAc,EAAE,EAAE,CACzD,IAAA,2BAAoB,EAAC,KAAK,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,aAAa;AAC1D,MAAM,UAAU,GAAG,CAAC,CAAO,EAAQ,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAEvF,MAAM,KAAK,GAAG,CAAC,MAAgB,EAAE,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC;AAE7F;;;GAGG;AACH,MAAM,UAAU,GAAG,CACjB,MAAsB,EACtB,KAGC,EACD,EAAQ,EACR,CAAS,EACT,IAAiB,EACjB,GAAW,EACoB,EAAE;IACjC,MAAM,QAAQ,GAAG,IAAA,gCAAgB,EAAC,IAAI,CAAC,CAAC;IACxC,MAAM,IAAI,GAAW,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,MAAM,EAAE,GAAG,CAAC,CAAS,EAAQ,EAAE;QAC7B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC7B,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC;IAEF,MAAM,MAAM,GAAgB,EAAE,CAAC;IAC/B,MAAM,IAAI,GAA2B,EAAE,CAAC;IACxC,KAAK,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;QACnD,MAAM,KAAK,GAAW,EAAE,CAAC;QACzB,IAAI,MAAwB,CAAC;QAC7B,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,KAAK,MAAM,KAAK,IAAI,IAAI,EAAE,CAAC;YACzB,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;YACnC,MAAM,MAAM,GAAG,KAAK,CAAC,WAAW,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;YACzC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG;gBAAE,SAAS;YACnC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG;gBAAE,SAAS;YAC7D,MAAM,OAAO,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACvF,MAAM,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC;YAC7C,KAAK,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7C,MAAM,KAAN,MAAM,GAAK,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAS,EAAC;YAChD,SAAS,GAAG,KAAK,CAAC,OAAO,CAAC,SAAS,CAAC;QACtC,CAAC;QACD,IAAI,CAAC,MAAM;YAAE,SAAS;QACtB,MAAM,CAAC,GAAG,IAAA,qBAAK,EAAC,EAAE,EAAE,MAAM,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,IAAI,EAAE,KAAK,CAAC;YAAE,SAAS;QACvB,sEAAsE;QACtE,IAAI,KAAK,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACpD,MAAM,IAAI,GACR,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,wCAAyB;YAC5C,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC;YACd,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,wCAAyB;gBAC9C,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC;gBACd,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACrB,IAAI,IAAI;YAAE,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACpD,IAAI,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACrB,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACnB,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC;YACtB,IAAI,IAAA,mBAAG,EAAC,CAAC,EAAE,KAAK,CAAC,GAAG,IAAA,mBAAG,EAAC,KAAK,EAAE,KAAK,CAAC;gBAAE,KAAK,GAAG,CAAC,CAAC;YACjD,IAAI,IAAA,mBAAG,EAAC,CAAC,EAAE,KAAK,CAAC,GAAG,IAAA,mBAAG,EAAC,GAAG,EAAE,KAAK,CAAC;gBAAE,GAAG,GAAG,CAAC,CAAC;QAC/C,CAAC;QACD,MAAM,MAAM,GAAG,IAAA,mBAAG,EAAC,IAAA,mBAAG,EAAC,GAAG,EAAE,KAAK,CAAC,EAAE,KAAK,CAAC,CAAC;QAC3C,IAAI,MAAM,GAAG,GAAG;YAAE,IAAI,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE,SAAS,EAAE,CAAC,CAAC;IACzE,CAAC;IAED,MAAM,IAAI,GAAc,EAAE,CAAC;IAC3B,MAAM,SAAS,GAAgB,EAAE,CAAC;IAClC,MAAM,OAAO,GAA8B,EAAE,CAAC;IAC9C,KAAK,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,IAAI,KAAK,CAAC,OAAO,EAAE,CAAC;QACjD,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;YAC3B,MAAM,MAAM,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;YACjC,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;YAC/B,MAAM,EAAE,GAAG,IAAA,mBAAG,EAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,IAAA,mBAAG,EAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YACvB,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,IAAI,GAAG;gBAAE,SAAS,CAAC,wBAAwB;YAChE,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YAC1B,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC;gBACvC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;gBAC/B,MAAM,KAAK,GAAS;oBAClB,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBAC3B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBAC3B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;iBAC5B,CAAC;gBACF,SAAS,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBAC/E,OAAO,CAAC,IAAI,CAAC,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC;gBACvD,MAAM;YACR,CAAC;YACD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;gBACpD,2FAA2F;gBAC3F,MAAM,GAAG,GAAU,QAAQ,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChD,MAAM,KAAK,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;gBAClC,IAAI,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,CAAC,CAAC;gBACvC,OAAO,CAAC,IAAI,CAAC,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC;gBACvD,MAAM;YACR,CAAC;QACH,CAAC;IACH,CAAC;IAED,wEAAwE;IACxE,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,MAAM,GAAG,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5B,MAAM,WAAW,GAAG,CAAC,CAAO,EAAE,CAAS,EAAE,EAAE;QACzC,IAAI,IAAI,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,CAAC,CAAC;IACF,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE,CAAC;QACvB,MAAM,GAAG,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAS,CAAC;QAChF,WAAW,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC;IAC/C,CAAC;IACD,KAAK,MAAM,CAAC,IAAI,OAAO;QAAE,WAAW,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,IAAI,QAAc,CAAC;IACnB,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;QACb,QAAQ,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAS,CAAC;IAC9C,CAAC;SAAM,CAAC;QACN,2DAA2D;QAC3D,MAAM,MAAM,GAAG;YACb,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzF,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;YAC3B,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;SACjC,CAAC;QACF,QAAQ,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CACtB,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CACrE,CAAC;IACZ,CAAC;IAED,OAAO;QACL,IAAI;QACJ,QAAQ;QACR,OAAO,EAAE;YACP,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;gBACvB,MAAM,EAAE,GAAG,CAAC,MAAM;gBAClB,KAAK,EAAE,UAAU,CAAC,GAAG,CAAC,KAAK,CAAC;gBAC5B,GAAG,EAAE,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC;gBACxB,MAAM,EAAE,KAAK,CAAC,GAAG,CAAC,MAAM,EAAE,eAAe,CAAC;gBAC1C,SAAS,EAAE,GAAG,CAAC,SAAS;aACzB,CAAC,CAAC;YACH,OAAO,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,KAAK,EAAE,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC3F,MAAM,EAAE,KAAK,CACX,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,EAC9C,eAAe,CAChB;YACD,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,eAAe,CAAC;YAClC,QAAQ,EAAE,UAAU,CAAC,QAAQ,CAAC;SAC/B;QACD,MAAM;QACN,IAAI;QACJ,SAAS;KACV,CAAC;AACJ,CAAC,CAAC;AAEF;;;;;;GAMG;AACI,MAAM,YAAY,GAAG,CAC1B,KAAyD,EACzD,MAAsB,EACgC,EAAE;IACxD,MAAM,EAAE,GAAG,CAAC,CAAS,EAAQ,EAAE;QAC7B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC7B,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC;IACF,MAAM,OAAO,GAAG,CAAC,IAAc,EAAE,IAA6B,EAAgB,EAAE;QAC9E,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,OAAO,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;YAClD,MAAM;YACN,MAAM,EAAE,MAAM,CAAC,YAAY,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,EAAE;SAC9C,CAAC,CAAC,CAAC;QACJ,MAAM,OAAO,GAAG,CAAC,IAAI,CAAC,SAAS,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC;YACxD,QAAQ;YACR,MAAM,EAAE,MAAM,CAAC,YAAY,CAAC,MAAM,CAAC,WAAW,CAAC,GAAG,CAAC,QAAQ,CAAW,CAAC,CAAC,MAAM;SAC/E,CAAC,CAAC,CAAC;QACJ,MAAM,KAAK,GAAc;YACvB,KAAK,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,CAAC;gBACxC,MAAM;gBACN,QAAQ,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;aAClE,CAAC,CAAC;YACH,OAAO,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,CAAC;gBAC9C,QAAQ;gBACR,MAAM,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAgB,EAAE,CAAC;oBACtC,EAAE,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC;oBAChC,EAAE,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;iBAC/B,CAAC;aACH,CAAC,CAAC;SACJ,CAAC;QACF,MAAM,IAAI,GAAG,IAAA,wBAAQ,EAAC,IAAI,EAAE,KAAK,EAAE;YACjC,KAAK,EAAG,IAAiB,CAAC,KAAK;YAC/B,OAAO,EAAG,IAAqB,CAAC,OAAO;SACxC,CAAC,CAAC;QACH,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,MAAM,GAAG,IAAA,0BAAU,EAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACrC,MAAM,GAAG,GAAG,IAAA,6BAAa,EAAC,MAAM,CAAC,CAAC;QAClC,MAAM,MAAM,GAAG,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;QAChD,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;QACvD,MAAM,IAAI,GAAG,IAAA,gCAAgB,EAAC,IAAI,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,CAAC;YAC1E,QAAQ,EAAE,CAAC,GAAG,MAAM;YACpB,GAAG,UAAU,CAAC,MAAM,EAAE,EAAE,KAAK,EAAE,OAAO,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,CAAC;SAC5D,CAAC,CAAC,CAAC;QACJ,OAAO,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,IAAI,EAAE,CAAC;IACtE,CAAC,CAAC;IACF,OAAO;QACL,KAAK,EAAE,CAAC,KAAK,CAAC,KAAK,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;QACzD,SAAS,EAAE,CAAC,KAAK,CAAC,SAAS,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;KACtE,CAAC;AACJ,CAAC,CAAC;AAjDW,QAAA,YAAY,gBAiDvB;AAEF;;;;;;GAMG;AACH,MAAM,YAAY,GAAG,CAAC,MAAsB,EAAE,KAAe,EAAE,GAAgB,EAAE,EAAE;IACjF,MAAM,CAAC,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC;IACvB,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,CAAO,EAAE,IAAY,EAAE,EAAE;QACxD,MAAM,CAAC,GAAG,IAAA,mBAAG,EAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACpB,MAAM,EAAE,GAAG,IAAA,qBAAK,EAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACpB,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAChC,CAAC;IACH,CAAC,CAAC;IAEF,KAAK,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,IAAI,GAAG,CAAC,MAAM,EAAE,CAAC;QACzD,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,IAAA,6BAAW,EAAC,KAAK,CAAC,WAAW,CAAC,CAAC;QAC5C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,KAAK,CAAC,cAAc,CAAC,KAAK,CAAC,CAAC;QACtC,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,GAAG,GAAG,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC7C,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;oBAAE,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChD,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACX,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;YACvE,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC;YAChD,KAAK,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;QAChF,CAAC;IACH,CAAC;IAED,MAAM,QAAQ,GAAG,CAAC,MAAkB,EAAE,KAAe,EAAQ,EAAE,CAAC;QAC9D,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QAC3E,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QAC3E,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;KAC5E,CAAC;IACF,KAAK,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC;QACnD,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,KAAK,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QAChE,KAAK,CAAC,KAAK,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACvF,CAAC;IACD,KAAK,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,GAAG,CAAC,SAAS,EAAE,CAAC;QACtD,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,IAAA,qCAAoB,EAC5B,KAAK,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACrC,KAAK,CAAC,WAAW,CAAC,KAAK,CAAC,EACxB,CAAC,CACF,CAAC;QACF,KAAK,CAAC,KAAK,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACvF,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAClB,CAAC,CAAC;AAEF,wEAAwE;AACxE,MAAM,aAAa,GAAG,CACpB,MAAsB,EACtB,KAAe,EACf,IAAc,EACd,GAAgB,EACA,EAAE;IAClB,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,YAAY,CAAC,MAAM,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;IAClD,OAAO;QACL,CAAC,EAAE,KAAK,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAChC,EAAE,EAAE,KAAK,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACjC,EAAE,EAAE,KAAK,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACjC,CAAC,EAAE,KAAK,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAChC,EAAE,EAAE,KAAK,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACjC,EAAE,EAAE,KAAK,CAAC,IAAA,mBAAG,EAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;KAClC,CAAC;AACJ,CAAC,CAAC;AAEF,MAAM,WAAW,GAAG,CAAC,IAAc,EAAsB,EAAE,CACzD,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAS,CAAuB,CAAC;AAE3E,uDAAuD;AAChD,MAAM,iBAAiB,GAAG,CAC/B,MAAsB,EACtB,KAAe,EACf,IAAkB,EACN,EAAE,CAAC,CAAC;IAChB,MAAM,EAAE,IAAI,CAAC,EAAE;IACf,IAAI,EAAE,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC;IAC5B,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAC5B,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,QAAQ,EAAE,eAAe,CAAC;QACvC,IAAI,EAAE,GAAG,CAAC,IAA0C;QACpD,GAAG,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,IAAI,EAAE,GAAG,CAAC;QAC/C,OAAO,EAAE,GAAG,CAAC,OAAO;KACrB,CAAC,CAAC;CACJ,CAAC,CAAC;AAbU,QAAA,iBAAiB,qBAa3B;AAEH,2DAA2D;AACpD,MAAM,qBAAqB,GAAG,CACnC,MAAsB,EACtB,KAAe,EACf,QAAsB,EACN,EAAE,CAAC,CAAC;IACpB,UAAU,EAAE,QAAQ,CAAC,EAAE;IACvB,IAAI,EAAE,WAAW,CAAC,QAAQ,CAAC,IAAI,CAAC;IAChC,MAAM,EAAE,KAAK,CAAC,QAAQ,CAAC,MAAM,EAAE,eAAe,CAAC;IAC/C,IAAI,EAAE,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAChC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,QAAQ,EAAE,eAAe,CAAC;QACvC,IAAI,EAAE,GAAG,CAAC,IAA8C;QACxD,GAAG,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,QAAQ,CAAC,IAAI,EAAE,GAAG,CAAC;QACnD,OAAO,EAAE,GAAG,CAAC,OAAO;KACrB,CAAC,CAAC;CACJ,CAAC,CAAC;AAdU,QAAA,qBAAqB,yBAc/B"}

package/dist/fea/analysis/post/wall-geometry.d.ts

@@ -0,0 +1,73 @@
+import type { FeaId } from "../../shared/types/model.type";
+import type { PierCutFace, SpandrelCutFace } from "../types/results.type";
+/** A point or vector in global axes: X, Y, Z in mm. */
+export type Vec3 = [number, number, number];
+/** Pier or spandrel axes 1, 2, 3 (rows) in global components. */
+export type WallAxes = [Vec3, Vec3, Vec3];
+export type WallKind = "pier" | "spandrel";
+export type WallCutFace = PierCutFace | SpandrelCutFace;
+/** The geometry of a pier or spandrel: its areas' shell elements and its members' pieces. */
+export interface WallShape {
+    areas: {
+        areaId: FeaId;
+        elements: Vec3[][];
+    }[];
+    members: {
+        memberId: FeaId;
+        pieces: [Vec3, Vec3][];
+    }[];
+}
+/** One cut of a pier or spandrel: its coordinate along axis 1 and the side it is taken from. */
+export interface WallCutPosition {
+    s: number;
+    face: WallCutFace;
+    /** Index into the user's `cuts` for an extra cut. */
+    cut?: number;
+}
+export declare const sub: (p: Vec3, q: Vec3) => Vec3;
+export declare const dot: (p: Vec3, q: Vec3) => number;
+export declare const cross: (p: Vec3, q: Vec3) => Vec3;
+/** Unit normal of a polygon by Newell's method (zero for a degenerate polygon). */
+export declare const polygonNormal: (points: readonly Vec3[]) => Vec3;
+/** The horizontal direction in a vertical plane of normal `n`, toward +X (+Y if along Y). */
+export declare const horizontalInPlane: (n: Vec3) => Vec3;
+/**
+ * Axes of a pier or spandrel (see `WallPier`, `WallSpandrel`).
+ *
+ * - Pier: 1 = +Z; 2 = global X rotated by `angle` about +Z, else the horizontal direction of
+ *   the area with the largest horizontal extent (the first wins a tie), else global X; 3 = 1 × 2.
+ * - Spandrel: 1 = horizontal in the plane of its first area, toward +X (+Y for a plane parallel
+ *   to Y), reversed by `reverse`; 2 = +Z; 3 = 1 × 2.
+ */
+export declare const wallAxes: (kind: WallKind, shape: WallShape, options: {
+    angle?: number;
+    reverse?: boolean;
+}) => WallAxes;
+/** Smallest and largest coordinate along `e1` of every element corner and member end. */
+export declare const wallExtent: (shape: WallShape, e1: Vec3) => {
+    min: number;
+    max: number;
+};
+/** The distance within which a point counts as on a cut (see `WALL_CUT_TOLERANCE`). */
+export declare const wallTolerance: (extent: {
+    min: number;
+    max: number;
+}) => number;
+/**
+ * The cuts of a pier or spandrel, in order along axis 1.
+ *
+ * - Pier: at each elevation where one of its areas or members starts or ends, a `"top"` cut
+ *   when something of the pier ends there from below and a `"bottom"` cut when something starts
+ *   there going up (the top first); then each extra cut not at such an elevation, as `"cut"`.
+ * - Spandrel: `"left"` at its smallest coordinate, `"right"` at its largest, and each extra cut
+ *   strictly between them as `"cut"`.
+ *
+ * @param extra - The extra cuts as coordinates along axis 1 (mm)
+ */
+export declare const wallCutPositions: (kind: WallKind, shape: WallShape, e1: Vec3, extra?: readonly number[]) => WallCutPosition[];
+/** Whether a face takes the elements on the negative side of its cut (below, or to the left). */
+export declare const fromNegativeSide: (face: WallCutFace) => face is "right" | "top" | "cut";
+/** Areas with a shell element on both sides of the cut at s (no mesh line there). */
+export declare const areasAcrossCut: (shape: WallShape, e1: Vec3, s: number, tol: number) => FeaId[];
+/** Whether some element or member of the shape reaches the cut at s. */
+export declare const reachesCut: (shape: WallShape, e1: Vec3, s: number, tol: number) => boolean;

package/dist/fea/analysis/post/wall-geometry.js

@@ -0,0 +1,188 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.reachesCut = exports.areasAcrossCut = exports.fromNegativeSide = exports.wallCutPositions = exports.wallTolerance = exports.wallExtent = exports.wallAxes = exports.horizontalInPlane = exports.polygonNormal = exports.cross = exports.dot = exports.sub = void 0;
+const fea_constant_1 = require("../../../core/constants/fea.constant");
+const sub = (p, q) => [p[0] - q[0], p[1] - q[1], p[2] - q[2]];
+exports.sub = sub;
+const dot = (p, q) => p[0] * q[0] + p[1] * q[1] + p[2] * q[2];
+exports.dot = dot;
+const cross = (p, q) => [
+    p[1] * q[2] - p[2] * q[1],
+    p[2] * q[0] - p[0] * q[2],
+    p[0] * q[1] - p[1] * q[0],
+];
+exports.cross = cross;
+const norm = (p) => Math.hypot(p[0], p[1], p[2]);
+const unit = (p) => {
+    const n = norm(p);
+    return n > 0 ? [p[0] / n, p[1] / n, p[2] / n] : [0, 0, 0];
+};
+/** Unit normal of a polygon by Newell's method (zero for a degenerate polygon). */
+const polygonNormal = (points) => {
+    const n = [0, 0, 0];
+    points.forEach((p, i) => {
+        const q = points[(i + 1) % points.length];
+        n[0] += (p[1] - q[1]) * (p[2] + q[2]);
+        n[1] += (p[2] - q[2]) * (p[0] + q[0]);
+        n[2] += (p[0] - q[0]) * (p[1] + q[1]);
+    });
+    return unit(n);
+};
+exports.polygonNormal = polygonNormal;
+/** The horizontal direction in a vertical plane of normal `n`, toward +X (+Y if along Y). */
+const horizontalInPlane = (n) => {
+    const d = unit([-n[1], n[0], 0]);
+    const flip = Math.abs(d[0]) > fea_constant_1.VERTICAL_MEMBER_TOLERANCE ? d[0] < 0 : d[1] < 0;
+    return flip ? [-d[0], -d[1], -d[2]] : d;
+};
+exports.horizontalInPlane = horizontalInPlane;
+/**
+ * Axes of a pier or spandrel (see `WallPier`, `WallSpandrel`).
+ *
+ * - Pier: 1 = +Z; 2 = global X rotated by `angle` about +Z, else the horizontal direction of
+ *   the area with the largest horizontal extent (the first wins a tie), else global X; 3 = 1 × 2.
+ * - Spandrel: 1 = horizontal in the plane of its first area, toward +X (+Y for a plane parallel
+ *   to Y), reversed by `reverse`; 2 = +Z; 3 = 1 × 2.
+ */
+const wallAxes = (kind, shape, options) => {
+    const up = [0, 0, 1];
+    if (kind === "pier") {
+        let e2 = [1, 0, 0];
+        if (options.angle !== undefined) {
+            e2 = [Math.cos(options.angle), Math.sin(options.angle), 0];
+        }
+        else {
+            let best = 0;
+            for (const area of shape.areas) {
+                if (area.elements.length === 0)
+                    continue;
+                const d = (0, exports.horizontalInPlane)((0, exports.polygonNormal)(area.elements[0]));
+                const s = area.elements.flat().map((p) => (0, exports.dot)(p, d));
+                const extent = Math.max(...s) - Math.min(...s);
+                if (extent > best) {
+                    best = extent;
+                    e2 = d;
+                }
+            }
+        }
+        return [up, e2, (0, exports.cross)(up, e2)];
+    }
+    const first = shape.areas.find((a) => a.elements.length > 0);
+    let e1 = first ? (0, exports.horizontalInPlane)((0, exports.polygonNormal)(first.elements[0])) : [1, 0, 0];
+    if (options.reverse)
+        e1 = [-e1[0], -e1[1], -e1[2]];
+    return [e1, up, (0, exports.cross)(e1, up)];
+};
+exports.wallAxes = wallAxes;
+/** Smallest and largest coordinate along `e1` of every element corner and member end. */
+const wallExtent = (shape, e1) => {
+    const s = [
+        ...shape.areas.flatMap((a) => a.elements.flat()),
+        ...shape.members.flatMap((m) => m.pieces.flat()),
+    ].map((p) => (0, exports.dot)(p, e1));
+    return { min: Math.min(...s), max: Math.max(...s) };
+};
+exports.wallExtent = wallExtent;
+/** The distance within which a point counts as on a cut (see `WALL_CUT_TOLERANCE`). */
+const wallTolerance = (extent) => fea_constant_1.WALL_CUT_TOLERANCE * Math.max(extent.max - extent.min, 1);
+exports.wallTolerance = wallTolerance;
+const range = (points, e1) => {
+    const s = points.map((p) => (0, exports.dot)(p, e1));
+    return { min: Math.min(...s), max: Math.max(...s) };
+};
+/** Whether anything of the shape lies on the negative / positive side touching the cut at s. */
+const touching = (shape, e1, s, tol) => {
+    let negative = false;
+    let positive = false;
+    const check = (points) => {
+        const r = range(points, e1);
+        if (r.max - r.min <= tol)
+            return; // lies in the cut plane
+        if (Math.abs(r.max - s) <= tol)
+            negative = true;
+        if (Math.abs(r.min - s) <= tol)
+            positive = true;
+    };
+    shape.areas.forEach((a) => a.elements.forEach(check));
+    shape.members.forEach((m) => m.pieces.forEach(check));
+    return { negative, positive };
+};
+/**
+ * The cuts of a pier or spandrel, in order along axis 1.
+ *
+ * - Pier: at each elevation where one of its areas or members starts or ends, a `"top"` cut
+ *   when something of the pier ends there from below and a `"bottom"` cut when something starts
+ *   there going up (the top first); then each extra cut not at such an elevation, as `"cut"`.
+ * - Spandrel: `"left"` at its smallest coordinate, `"right"` at its largest, and each extra cut
+ *   strictly between them as `"cut"`.
+ *
+ * @param extra - The extra cuts as coordinates along axis 1 (mm)
+ */
+const wallCutPositions = (kind, shape, e1, extra = []) => {
+    const extent = (0, exports.wallExtent)(shape, e1);
+    const tol = (0, exports.wallTolerance)(extent);
+    const auto = [];
+    const addAuto = (s) => {
+        if (!auto.some((a) => Math.abs(a - s) <= tol))
+            auto.push(s);
+    };
+    if (kind === "pier") {
+        for (const area of shape.areas) {
+            if (area.elements.length === 0)
+                continue;
+            const r = range(area.elements.flat(), e1);
+            addAuto(r.min);
+            addAuto(r.max);
+        }
+        for (const member of shape.members) {
+            const r = range(member.pieces.flat(), e1);
+            addAuto(r.min);
+            addAuto(r.max);
+        }
+    }
+    else {
+        addAuto(extent.min);
+        addAuto(extent.max);
+    }
+    auto.sort((a, b) => a - b);
+    const cuts = [];
+    for (const s of auto) {
+        if (kind === "spandrel") {
+            cuts.push({ s, face: s === extent.min ? "left" : "right" });
+            continue;
+        }
+        const side = touching(shape, e1, s, tol);
+        if (side.negative)
+            cuts.push({ s, face: "top" });
+        if (side.positive)
+            cuts.push({ s, face: "bottom" });
+    }
+    extra.forEach((s, cut) => {
+        if (auto.some((a) => Math.abs(a - s) <= tol))
+            return;
+        if (cuts.some((c) => c.face === "cut" && Math.abs(c.s - s) <= tol))
+            return;
+        cuts.push({ s, face: "cut", cut });
+    });
+    // stable: a "top" stays before the "bottom" at the same level
+    return cuts.sort((a, b) => a.s - b.s);
+};
+exports.wallCutPositions = wallCutPositions;
+/** Whether a face takes the elements on the negative side of its cut (below, or to the left). */
+const fromNegativeSide = (face) => face === "top" || face === "right" || face === "cut";
+exports.fromNegativeSide = fromNegativeSide;
+/** Areas with a shell element on both sides of the cut at s (no mesh line there). */
+const areasAcrossCut = (shape, e1, s, tol) => shape.areas
+    .filter((a) => a.elements.some((el) => {
+    const r = range(el, e1);
+    return r.min < s - tol && r.max > s + tol;
+}))
+    .map((a) => a.areaId);
+exports.areasAcrossCut = areasAcrossCut;
+/** Whether some element or member of the shape reaches the cut at s. */
+const reachesCut = (shape, e1, s, tol) => [...shape.areas.flatMap((a) => a.elements), ...shape.members.flatMap((m) => m.pieces)].some((points) => {
+    const r = range(points, e1);
+    return r.min <= s + tol && r.max >= s - tol;
+});
+exports.reachesCut = reachesCut;
+//# sourceMappingURL=wall-geometry.js.map

package/dist/fea/analysis/post/wall-geometry.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"wall-geometry.js","sourceRoot":"","sources":["../../../../src/fea/analysis/post/wall-geometry.ts"],"names":[],"mappings":";;;AAAA,mEAAiG;AA6B1F,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAA1E,QAAA,GAAG,OAAuE;AAChF,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAApE,QAAA,GAAG,OAAiE;AAC1E,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC;IAC/C,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;CAC1B,CAAC;AAJW,QAAA,KAAK,SAIhB;AACF,MAAM,IAAI,GAAG,CAAC,CAAO,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACvD,MAAM,IAAI,GAAG,CAAC,CAAO,EAAQ,EAAE;IAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IAClB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAC5D,CAAC,CAAC;AAEF,mFAAmF;AAC5E,MAAM,aAAa,GAAG,CAAC,MAAuB,EAAQ,EAAE;IAC7D,MAAM,CAAC,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACtB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;QAC1C,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxC,CAAC,CAAC,CAAC;IACH,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC;AACjB,CAAC,CAAC;AATW,QAAA,aAAa,iBASxB;AAEF,6FAA6F;AACtF,MAAM,iBAAiB,GAAG,CAAC,CAAO,EAAQ,EAAE;IACjD,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,wCAAyB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC9E,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1C,CAAC,CAAC;AAJW,QAAA,iBAAiB,qBAI5B;AAEF;;;;;;;GAOG;AACI,MAAM,QAAQ,GAAG,CACtB,IAAc,EACd,KAAgB,EAChB,OAA8C,EACpC,EAAE;IACZ,MAAM,EAAE,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;QACpB,IAAI,EAAE,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACzB,IAAI,OAAO,CAAC,KAAK,KAAK,SAAS,EAAE,CAAC;YAChC,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7D,CAAC;aAAM,CAAC;YACN,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,KAAK,MAAM,IAAI,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;gBAC/B,IAAI,IAAI,CAAC,QAAQ,CAAC,MAAM,KAAK,CAAC;oBAAE,SAAS;gBACzC,MAAM,CAAC,GAAG,IAAA,yBAAiB,EAAC,IAAA,qBAAa,EAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC7D,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAA,WAAG,EAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBACrD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC/C,IAAI,MAAM,GAAG,IAAI,EAAE,CAAC;oBAClB,IAAI,GAAG,MAAM,CAAC;oBACd,EAAE,GAAG,CAAC,CAAC;gBACT,CAAC;YACH,CAAC;QACH,CAAC;QACD,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,IAAA,aAAK,EAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IACjC,CAAC;IACD,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC7D,IAAI,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,IAAA,yBAAiB,EAAC,IAAA,qBAAa,EAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAU,CAAC;IAC3F,IAAI,OAAO,CAAC,OAAO;QAAE,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,IAAA,aAAK,EAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;AACjC,CAAC,CAAC;AA7BW,QAAA,QAAQ,YA6BnB;AAEF,yFAAyF;AAClF,MAAM,UAAU,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAgC,EAAE;IACrF,MAAM,CAAC,GAAG;QACR,GAAG,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,EAAE,CAAC;QAChD,GAAG,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC;KACjD,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAA,WAAG,EAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACzB,OAAO,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;AACtD,CAAC,CAAC;AANW,QAAA,UAAU,cAMrB;AAEF,uFAAuF;AAChF,MAAM,aAAa,GAAG,CAAC,MAAoC,EAAE,EAAE,CACpE,iCAAkB,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;AAD/C,QAAA,aAAa,iBACkC;AAE5D,MAAM,KAAK,GAAG,CAAC,MAAuB,EAAE,EAAQ,EAAE,EAAE;IAClD,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAA,WAAG,EAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACxC,OAAO,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;AACtD,CAAC,CAAC;AAEF,gGAAgG;AAChG,MAAM,QAAQ,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAE,CAAS,EAAE,GAAW,EAAE,EAAE;IACtE,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,MAAM,KAAK,GAAG,CAAC,MAAuB,EAAE,EAAE;QACxC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;QAC5B,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG;YAAE,OAAO,CAAC,wBAAwB;QAC1D,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,GAAG;YAAE,QAAQ,GAAG,IAAI,CAAC;QAChD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,GAAG;YAAE,QAAQ,GAAG,IAAI,CAAC;IAClD,CAAC,CAAC;IACF,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;IACtD,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;IACtD,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,CAAC;AAChC,CAAC,CAAC;AAEF;;;;;;;;;;GAUG;AACI,MAAM,gBAAgB,GAAG,CAC9B,IAAc,EACd,KAAgB,EAChB,EAAQ,EACR,QAA2B,EAAE,EACV,EAAE;IACrB,MAAM,MAAM,GAAG,IAAA,kBAAU,EAAC,KAAK,EAAE,EAAE,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,IAAA,qBAAa,EAAC,MAAM,CAAC,CAAC;IAClC,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,MAAM,OAAO,GAAG,CAAC,CAAS,EAAE,EAAE;QAC5B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9D,CAAC,CAAC;IACF,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;QACpB,KAAK,MAAM,IAAI,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;YAC/B,IAAI,IAAI,CAAC,QAAQ,CAAC,MAAM,KAAK,CAAC;gBAAE,SAAS;YACzC,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,CAAC;YAC1C,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACf,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,CAAC;QACD,KAAK,MAAM,MAAM,IAAI,KAAK,CAAC,OAAO,EAAE,CAAC;YACnC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,CAAC;YAC1C,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACf,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,CAAC;IACH,CAAC;SAAM,CAAC;QACN,OAAO,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;QACpB,OAAO,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACtB,CAAC;IACD,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAE3B,MAAM,IAAI,GAAsB,EAAE,CAAC;IACnC,KAAK,MAAM,CAAC,IAAI,IAAI,EAAE,CAAC;QACrB,IAAI,IAAI,KAAK,UAAU,EAAE,CAAC;YACxB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,KAAK,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;YAC5D,SAAS;QACX,CAAC;QACD,MAAM,IAAI,GAAG,QAAQ,CAAC,KAAK,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QACzC,IAAI,IAAI,CAAC,QAAQ;YAAE,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,CAAC;QACjD,IAAI,IAAI,CAAC,QAAQ;YAAE,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,CAAC,CAAC;IACtD,CAAC;IACD,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE;QACvB,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAAE,OAAO;QACrD,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAAE,OAAO;QAC3E,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC;IACrC,CAAC,CAAC,CAAC;IACH,8DAA8D;IAC9D,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AACxC,CAAC,CAAC;AA/CW,QAAA,gBAAgB,oBA+C3B;AAEF,iGAAiG;AAC1F,MAAM,gBAAgB,GAAG,CAAC,IAAiB,EAAE,EAAE,CACpD,IAAI,KAAK,KAAK,IAAI,IAAI,KAAK,OAAO,IAAI,IAAI,KAAK,KAAK,CAAC;AAD1C,QAAA,gBAAgB,oBAC0B;AAEvD,qFAAqF;AAC9E,MAAM,cAAc,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAE,CAAS,EAAE,GAAW,EAAW,EAAE,CAC5F,KAAK,CAAC,KAAK;KACR,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CACZ,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE;IACrB,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;AAC5C,CAAC,CAAC,CACH;KACA,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;AARb,QAAA,cAAc,kBAQD;AAE1B,wEAAwE;AACjE,MAAM,UAAU,GAAG,CAAC,KAAgB,EAAE,EAAQ,EAAE,CAAS,EAAE,GAAW,EAAW,EAAE,CACxF,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,EAAE,GAAG,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CACzF,CAAC,MAAM,EAAE,EAAE;IACT,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IAC5B,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC;AAC9C,CAAC,CACF,CAAC;AANS,QAAA,UAAU,cAMnB"}

package/dist/fea/analysis/solver/analysis-system.js

@@ -1,7 +1,7 @@
 "use strict";
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.prepareAnalysisSystem = void 0;
-const frame_constant_1 = require("../../../core/constants/frame.constant");
+const fea_constant_1 = require("../../../core/constants/fea.constant");
 const fea_error_1 = require("../../../core/errors/fea.error");
 const member_assembly_1 = require("../../frame/solver/member-assembly");
 const system_member_1 = require("../../frame/solver/system-member");
@@ -144,7 +144,7 @@ const prepareAnalysisSystem = (model) => {
                 }
             }
         }
-        const result = (0, skyline_1.skylineFactor)(K, frame_constant_1.SINGULARITY_TOLERANCE);
+        const result = (0, skyline_1.skylineFactor)(K, fea_constant_1.SINGULARITY_TOLERANCE);
         if (!result.ok) {
             const dof = free[result.failedPivot];
             (0, issues_1.throwIssues)([