@theerapat-s28/ts-structural-eng-tools 1.5.0 → 2.0.1

package/dist/esm/utils/matrix.js

@@ -0,0 +1,145 @@
+/** Returns an m×n zero matrix. */
+export const zeros = (m, n) => Array.from({ length: m }, () => new Array(n).fill(0));
+/** Returns the n×n identity matrix. */
+export const identity = (n) => {
+    const I = zeros(n, n);
+    for (let i = 0; i < n; i++)
+        I[i][i] = 1;
+    return I;
+};
+/** Returns the transpose of a matrix. */
+export const transpose = (A) => {
+    const m = A.length;
+    const n = A[0].length;
+    const T = zeros(n, m);
+    for (let i = 0; i < m; i++) {
+        for (let j = 0; j < n; j++)
+            T[j][i] = A[i][j];
+    }
+    return T;
+};
+/** Returns the product A·B of an m×n and an n×p matrix. */
+export const matMul = (A, B) => {
+    const m = A.length;
+    const n = B.length;
+    const p = B[0].length;
+    const C = zeros(m, p);
+    for (let i = 0; i < m; i++) {
+        for (let k = 0; k < n; k++) {
+            const aik = A[i][k];
+            if (aik === 0)
+                continue;
+            for (let j = 0; j < p; j++)
+                C[i][j] += aik * B[k][j];
+        }
+    }
+    return C;
+};
+/** Returns the product A·x of an m×n matrix and an n-vector. */
+export const matVec = (A, x) => A.map((row) => row.reduce((sum, a, j) => sum + a * x[j], 0));
+/**
+ * Cholesky factorisation A = L·Lᵀ of a symmetric positive-definite matrix.
+ *
+ * @param A - Symmetric matrix (only the lower triangle is read)
+ * @param tolerance - A pivot at or below `tolerance` times its own diagonal term is treated as zero
+ * @returns The lower-triangular factor L, or `null` when A is not positive definite (singular or a mechanism)
+ */
+export const choleskyFactor = (A, tolerance = 1e-12) => {
+    const n = A.length;
+    const L = zeros(n, n);
+    for (let i = 0; i < n; i++) {
+        for (let j = 0; j <= i; j++) {
+            let sum = A[i][j];
+            for (let k = 0; k < j; k++)
+                sum -= L[i][k] * L[j][k];
+            if (i === j) {
+                if (!(sum > tolerance * A[i][i]))
+                    return null;
+                L[i][i] = Math.sqrt(sum);
+            }
+            else {
+                L[i][j] = sum / L[j][j];
+            }
+        }
+    }
+    return L;
+};
+/**
+ * Solves A·x = b from the Cholesky factor L of A (forward then back substitution), so one
+ * factorisation serves any number of right-hand sides.
+ *
+ * @param L - Lower-triangular factor from `choleskyFactor`
+ * @param b - Right-hand side
+ */
+export const choleskySolveFactored = (L, b) => {
+    const n = L.length;
+    const y = new Array(n).fill(0);
+    for (let i = 0; i < n; i++) {
+        let sum = b[i];
+        for (let k = 0; k < i; k++)
+            sum -= L[i][k] * y[k];
+        y[i] = sum / L[i][i];
+    }
+    const x = new Array(n).fill(0);
+    for (let i = n - 1; i >= 0; i--) {
+        let sum = y[i];
+        for (let k = i + 1; k < n; k++)
+            sum -= L[k][i] * x[k];
+        x[i] = sum / L[i][i];
+    }
+    return x;
+};
+/**
+ * Solves A·x = b for a symmetric positive-definite A by Cholesky factorisation.
+ *
+ * @param A - Symmetric matrix (only the lower triangle is read)
+ * @param b - Right-hand side
+ * @param tolerance - A pivot at or below `tolerance` times its own diagonal term is treated as zero
+ * @returns The solution, or `null` when A is not positive definite (singular or a mechanism)
+ */
+export const choleskySolve = (A, b, tolerance = 1e-12) => {
+    const L = choleskyFactor(A, tolerance);
+    return L ? choleskySolveFactored(L, b) : null;
+};
+/**
+ * Solves A·X = B for a general square A by Gaussian elimination with partial pivoting.
+ *
+ * @param A - n×n matrix
+ * @param B - n×p right-hand-side matrix
+ * @param tolerance - A pivot at or below `tolerance` times the largest entry of A is treated as zero
+ * @returns The n×p solution, or `null` when A is singular
+ */
+export const solveLinearSystem = (A, B, tolerance = 1e-12) => {
+    const n = A.length;
+    const p = B[0].length;
+    const M = A.map((row, i) => [...row, ...B[i]]);
+    const scale = Math.max(...A.map((row) => Math.max(...row.map(Math.abs))));
+    for (let col = 0; col < n; col++) {
+        let pivot = col;
+        for (let r = col + 1; r < n; r++) {
+            if (Math.abs(M[r][col]) > Math.abs(M[pivot][col]))
+                pivot = r;
+        }
+        if (!(Math.abs(M[pivot][col]) > tolerance * scale))
+            return null;
+        [M[col], M[pivot]] = [M[pivot], M[col]];
+        for (let r = col + 1; r < n; r++) {
+            const factor = M[r][col] / M[col][col];
+            if (factor === 0)
+                continue;
+            for (let c = col; c < n + p; c++)
+                M[r][c] -= factor * M[col][c];
+        }
+    }
+    const X = zeros(n, p);
+    for (let c = 0; c < p; c++) {
+        for (let i = n - 1; i >= 0; i--) {
+            let sum = M[i][n + c];
+            for (let k = i + 1; k < n; k++)
+                sum -= M[i][k] * X[k][c];
+            X[i][c] = sum / M[i][i];
+        }
+    }
+    return X;
+};
+//# sourceMappingURL=matrix.js.map

package/dist/esm/utils/matrix.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"matrix.js","sourceRoot":"","sources":["../../../src/utils/matrix.ts"],"names":[],"mappings":"AAEA,kCAAkC;AAClC,MAAM,CAAC,MAAM,KAAK,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE,CACpD,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;AAEhE,uCAAuC;AACvC,MAAM,CAAC,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAU,EAAE;IAC5C,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACxC,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF,yCAAyC;AACzC,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,CAAS,EAAU,EAAE;IAC7C,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACtB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF,2DAA2D;AAC3D,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE;IACrD,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACtB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpB,IAAI,GAAG,KAAK,CAAC;gBAAE,SAAS;YACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvD,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF,gEAAgE;AAChE,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,CAAS,EAAE,CAAW,EAAY,EAAE,CACzD,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAE/D;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,CAAS,EAAE,SAAS,GAAG,KAAK,EAAiB,EAAE;IAC5E,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrD,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;gBACZ,IAAI,CAAC,CAAC,GAAG,GAAG,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAAE,OAAO,IAAI,CAAC;gBAC9C,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YAC3B,CAAC;iBAAM,CAAC;gBACN,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1B,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,qBAAqB,GAAG,CAAC,CAAS,EAAE,CAAW,EAAY,EAAE;IACxE,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACvC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtD,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAE,CAAW,EAAE,SAAS,GAAG,KAAK,EAAmB,EAAE;IAC1F,MAAM,CAAC,GAAG,cAAc,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACvC,OAAO,CAAC,CAAC,CAAC,CAAC,qBAAqB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;AAChD,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,SAAS,GAAG,KAAK,EAAiB,EAAE;IAC1F,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACtB,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1E,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,IAAI,KAAK,GAAG,GAAG,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;gBAAE,KAAK,GAAG,CAAC,CAAC;QAC/D,CAAC;QACD,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,GAAG,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC;QAChE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACxC,KAAK,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACvC,IAAI,MAAM,KAAK,CAAC;gBAAE,SAAS;YAC3B,KAAK,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAClE,CAAC;IACH,CAAC;IACD,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAChC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACtB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzD,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC"}

package/dist/esm/utils/merge-warning.js

@@ -0,0 +1,11 @@
+/**
+ * Merges two warning arrays into a single combined array.
+ *
+ * @param existingWarnings - Warnings already accumulated
+ * @param newWarnings - Additional warnings to append
+ * @returns Combined warning array
+ */
+export const mergeWarnings = (existingWarnings, newWarnings) => {
+    return [...existingWarnings, ...newWarnings];
+};
+//# sourceMappingURL=merge-warning.js.map

package/dist/esm/utils/merge-warning.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"merge-warning.js","sourceRoot":"","sources":["../../../src/utils/merge-warning.ts"],"names":[],"mappings":"AAEA;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,gBAA0B,EAAE,WAAqB,EAAY,EAAE;IAC3F,OAAO,CAAC,GAAG,gBAAgB,EAAE,GAAG,WAAW,CAAC,CAAC;AAC/C,CAAC,CAAC"}

package/dist/esm/utils/skyline.js

@@ -0,0 +1,151 @@
+/** Creates a zero skyline matrix with the given first row of each column (`first[j] ≤ j`). */
+export const createSkyline = (first) => {
+    const n = first.length;
+    const firstRows = Int32Array.from(first);
+    const columns = Array.from({ length: n }, (_, j) => new Float64Array(j - firstRows[j] + 1));
+    return { n, first: firstRows, columns };
+};
+/** Adds `value` to A[i][j] (and so to A[j][i]); the entry must lie inside the profile. */
+export const skylineAdd = (S, i, j, value) => {
+    const row = Math.min(i, j);
+    const col = Math.max(i, j);
+    const offset = row - S.first[col];
+    if (offset < 0)
+        throw new Error(`Skyline entry (${row}, ${col}) is outside the profile.`);
+    S.columns[col][offset] += value;
+};
+/** Reads A[i][j] (0 outside the profile). */
+export const skylineGet = (S, i, j) => {
+    const row = Math.min(i, j);
+    const col = Math.max(i, j);
+    const offset = row - S.first[col];
+    return offset < 0 ? 0 : S.columns[col][offset];
+};
+/** Number of stored entries, a measure of memory and of factorisation work. */
+export const skylineSize = (S) => S.columns.reduce((sum, column) => sum + column.length, 0);
+/**
+ * Factorises a symmetric positive-definite skyline matrix in place as A = Uᵀ·U (U upper
+ * triangular, same profile).
+ *
+ * @param S - Matrix, overwritten by U
+ * @param tolerance - A pivot at or below `tolerance` times its own diagonal term is treated as zero
+ * @returns `{ ok: true }`, or `{ ok: false, failedPivot }` with the index of the first pivot that
+ *   failed (A is singular or indefinite there — for a structure, a mechanism involving that DOF)
+ */
+export const skylineFactor = (S, tolerance = 1e-12) => {
+    const { n, first, columns } = S;
+    for (let j = 0; j < n; j++) {
+        const colJ = columns[j];
+        const fj = first[j];
+        for (let i = fj; i < j; i++) {
+            const colI = columns[i];
+            const fi = first[i];
+            const k0 = Math.max(fi, fj);
+            let sum = colJ[i - fj];
+            for (let k = k0; k < i; k++)
+                sum -= colI[k - fi] * colJ[k - fj];
+            colJ[i - fj] = sum / colI[i - fi];
+        }
+        const diagonal = colJ[j - fj];
+        let sum = diagonal;
+        for (let k = fj; k < j; k++)
+            sum -= colJ[k - fj] * colJ[k - fj];
+        if (!(sum > tolerance * diagonal))
+            return { ok: false, failedPivot: j };
+        colJ[j - fj] = Math.sqrt(sum);
+    }
+    return { ok: true };
+};
+/**
+ * Solves A·x = b from the factor produced by `skylineFactor` (forward then back substitution).
+ * One factorisation serves any number of right-hand sides.
+ */
+export const skylineSolve = (U, b) => {
+    const { n, first, columns } = U;
+    const y = Array.from(b);
+    // Uᵀ·y = b
+    for (let j = 0; j < n; j++) {
+        const col = columns[j];
+        const fj = first[j];
+        let sum = y[j];
+        for (let k = fj; k < j; k++)
+            sum -= col[k - fj] * y[k];
+        y[j] = sum / col[j - fj];
+    }
+    // U·x = y
+    for (let j = n - 1; j >= 0; j--) {
+        const col = columns[j];
+        const fj = first[j];
+        y[j] /= col[j - fj];
+        const xj = y[j];
+        for (let k = fj; k < j; k++)
+            y[k] -= col[k - fj] * xj;
+    }
+    return y;
+};
+/**
+ * Reverse Cuthill–McKee ordering of a graph, which keeps connected vertices close together and
+ * so narrows the band of the matrix built on it. Each connected component starts from a
+ * pseudo-peripheral vertex.
+ *
+ * @param adjacency - Neighbours of each vertex (undirected)
+ * @returns The vertices in their new order (`order[k]` = original index of the k-th vertex)
+ */
+export const reverseCuthillMcKee = (adjacency) => {
+    const n = adjacency.length;
+    const degree = adjacency.map((neighbours) => neighbours.length);
+    const visited = new Uint8Array(n);
+    const order = [];
+    /** Breadth-first levels from `start`, within unvisited vertices. */
+    const levels = (start) => {
+        const level = new Int32Array(n).fill(-1);
+        level[start] = 0;
+        const queue = [start];
+        let last = start;
+        for (let q = 0; q < queue.length; q++) {
+            const v = queue[q];
+            for (const w of adjacency[v]) {
+                if (!visited[w] && level[w] < 0) {
+                    level[w] = level[v] + 1;
+                    queue.push(w);
+                    if (level[w] > level[last] || (level[w] === level[last] && degree[w] < degree[last])) {
+                        last = w;
+                    }
+                }
+            }
+        }
+        return { depth: level[last], last };
+    };
+    const byDegree = Array.from({ length: n }, (_, i) => i).sort((a, b) => degree[a] - degree[b]);
+    for (const seed of byDegree) {
+        if (visited[seed])
+            continue;
+        // pseudo-peripheral start: move to the far end of the BFS until the depth stops growing
+        let start = seed;
+        let { depth, last } = levels(start);
+        for (let iteration = 0; iteration < 8; iteration++) {
+            const next = levels(last);
+            if (next.depth <= depth)
+                break;
+            start = last;
+            depth = next.depth;
+            last = next.last;
+        }
+        const component = [start];
+        visited[start] = 1;
+        for (let q = 0; q < component.length; q++) {
+            const neighbours = adjacency[component[q]]
+                .filter((w) => !visited[w])
+                .sort((a, b) => degree[a] - degree[b]);
+            for (const w of neighbours) {
+                if (!visited[w]) {
+                    visited[w] = 1;
+                    component.push(w);
+                }
+            }
+        }
+        order.push(...component);
+    }
+    return order.reverse();
+};
+//# sourceMappingURL=skyline.js.map

package/dist/esm/utils/skyline.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"skyline.js","sourceRoot":"","sources":["../../../src/utils/skyline.ts"],"names":[],"mappings":"AAaA,8FAA8F;AAC9F,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,KAAwB,EAAiB,EAAE;IACvE,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IACvB,MAAM,SAAS,GAAG,UAAU,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACzC,MAAM,OAAO,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,YAAY,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC5F,OAAO,EAAE,CAAC,EAAE,KAAK,EAAE,SAAS,EAAE,OAAO,EAAE,CAAC;AAC1C,CAAC,CAAC;AAEF,0FAA0F;AAC1F,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,CAAgB,EAAE,CAAS,EAAE,CAAS,EAAE,KAAa,EAAQ,EAAE;IACxF,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IAClC,IAAI,MAAM,GAAG,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,kBAAkB,GAAG,KAAK,GAAG,2BAA2B,CAAC,CAAC;IAC1F,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,IAAI,KAAK,CAAC;AAClC,CAAC,CAAC;AAEF,6CAA6C;AAC7C,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,CAAgB,EAAE,CAAS,EAAE,CAAS,EAAU,EAAE;IAC3E,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IAClC,OAAO,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;AACjD,CAAC,CAAC;AAEF,+EAA+E;AAC/E,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAgB,EAAU,EAAE,CACtD,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,MAAM,EAAE,EAAE,CAAC,GAAG,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;AAE5D;;;;;;;;GAQG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAC3B,CAAgB,EAChB,SAAS,GAAG,KAAK,EACkC,EAAE;IACrD,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,CAAC,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACpB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAC5B,IAAI,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,IAAI,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;YAChE,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACpC,CAAC;QACD,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QAC9B,IAAI,GAAG,GAAG,QAAQ,CAAC;QACnB,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QAChE,IAAI,CAAC,CAAC,GAAG,GAAG,SAAS,GAAG,QAAQ,CAAC;YAAE,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,WAAW,EAAE,CAAC,EAAE,CAAC;QACxE,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAChC,CAAC;IACD,OAAO,EAAE,EAAE,EAAE,IAAI,EAAE,CAAC;AACtB,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,CAAgB,EAAE,CAAoB,EAAY,EAAE;IAC/E,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,CAAC,CAAC;IAChC,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,WAAW;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpB,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACvD,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IAC3B,CAAC;IACD,UAAU;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpB,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACpB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;IACxD,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAG,CAAC,SAAqB,EAAY,EAAE;IACrE,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,CAAC;IAC3B,MAAM,MAAM,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,UAAU,EAAE,EAAE,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;IAChE,MAAM,OAAO,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,oEAAoE;IACpE,MAAM,MAAM,GAAG,CAAC,KAAa,EAAE,EAAE;QAC/B,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,KAAK,GAAG,CAAC,KAAK,CAAC,CAAC;QACtB,IAAI,IAAI,GAAG,KAAK,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,MAAM,CAAC,IAAI,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC;gBAC7B,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;oBAChC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBACxB,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;oBACd,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,IAAI,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;wBACrF,IAAI,GAAG,CAAC,CAAC;oBACX,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;QACD,OAAO,EAAE,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,CAAC;IACtC,CAAC,CAAC;IAEF,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9F,KAAK,MAAM,IAAI,IAAI,QAAQ,EAAE,CAAC;QAC5B,IAAI,OAAO,CAAC,IAAI,CAAC;YAAE,SAAS;QAC5B,wFAAwF;QACxF,IAAI,KAAK,GAAG,IAAI,CAAC;QACjB,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;QACpC,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE,CAAC;YACnD,MAAM,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;YAC1B,IAAI,IAAI,CAAC,KAAK,IAAI,KAAK;gBAAE,MAAM;YAC/B,KAAK,GAAG,IAAI,CAAC;YACb,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC;YACnB,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC;QACnB,CAAC;QAED,MAAM,SAAS,GAAa,CAAC,KAAK,CAAC,CAAC;QACpC,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,UAAU,GAAG,SAAS,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;iBACvC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;iBAC1B,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;YACzC,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE,CAAC;gBAC3B,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC;oBAChB,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBACf,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBACpB,CAAC;YACH,CAAC;QACH,CAAC;QACD,KAAK,CAAC,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC;IAC3B,CAAC;IACD,OAAO,KAAK,CAAC,OAAO,EAAE,CAAC;AACzB,CAAC,CAAC"}

package/dist/fea/frame/elements/beam-3d-stiffness.d.ts

@@ -0,0 +1,29 @@
+import { type Matrix } from "../../../utils/matrix";
+/**
+ * Local stiffness matrix of a 3D prismatic beam element, exactly as in the reference
+ * "Complete Beam Stiffness Matrix".
+ *
+ * DOF order for node i then node j: [U1, U2, U3, R1, R2, R3]. Local 1 runs along the member.
+ * The reference matrix measures the bending rotations R2 and R3 opposite to the right-hand
+ * rule; use `beam3dLocalStiffnessRightHanded` for the right-handed form.
+ *
+ * The reference uses one `EI`; here the U2–R3 terms use `E·I33` and the U3–R2 terms use
+ * `E·I22`, which reduces to the reference when `I22 = I33`.
+ *
+ * @param E - Young's modulus in MPa
+ * @param G - Shear modulus in MPa
+ * @param A - Area in mm²
+ * @param I22 - Second moment about local 2 in mm⁴
+ * @param I33 - Second moment about local 3 in mm⁴
+ * @param J - Torsion constant in mm⁴
+ * @param L - Member length in mm
+ * @returns 12×12 symmetric matrix in N and mm
+ */
+export declare const beam3dLocalStiffness: (E: number, G: number, A: number, I22: number, I33: number, J: number, L: number) => Matrix;
+/**
+ * Converts a local matrix between the reference rotation convention and the right-hand rule
+ * (`D·k·D`, with `D` flipping the sign of R2 and R3 at both ends). The conversion is its own inverse.
+ */
+export declare const flipBendingRotationSigns: (k: Matrix) => Matrix;
+/** Same element as `beam3dLocalStiffness` with all rotations following the right-hand rule. */
+export declare const beam3dLocalStiffnessRightHanded: (E: number, G: number, A: number, I22: number, I33: number, J: number, L: number) => Matrix;

package/dist/fea/frame/elements/beam-3d-stiffness.js

@@ -0,0 +1,90 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.beam3dLocalStiffnessRightHanded = exports.flipBendingRotationSigns = exports.beam3dLocalStiffness = void 0;
+const matrix_1 = require("../../../utils/matrix");
+/**
+ * Local DOF indices (0-based) whose reference-matrix rotation sign is opposite to the
+ * right-hand rule: R2 and R3 at the start and end nodes.
+ */
+const FLIPPED_DOFS = [4, 5, 10, 11];
+/**
+ * Local stiffness matrix of a 3D prismatic beam element, exactly as in the reference
+ * "Complete Beam Stiffness Matrix".
+ *
+ * DOF order for node i then node j: [U1, U2, U3, R1, R2, R3]. Local 1 runs along the member.
+ * The reference matrix measures the bending rotations R2 and R3 opposite to the right-hand
+ * rule; use `beam3dLocalStiffnessRightHanded` for the right-handed form.
+ *
+ * The reference uses one `EI`; here the U2–R3 terms use `E·I33` and the U3–R2 terms use
+ * `E·I22`, which reduces to the reference when `I22 = I33`.
+ *
+ * @param E - Young's modulus in MPa
+ * @param G - Shear modulus in MPa
+ * @param A - Area in mm²
+ * @param I22 - Second moment about local 2 in mm⁴
+ * @param I33 - Second moment about local 3 in mm⁴
+ * @param J - Torsion constant in mm⁴
+ * @param L - Member length in mm
+ * @returns 12×12 symmetric matrix in N and mm
+ */
+const beam3dLocalStiffness = (E, G, A, I22, I33, J, L) => {
+    const k = (0, matrix_1.zeros)(12, 12);
+    // 1-based row/column, upper triangle as written in the reference; mirrored for symmetry.
+    const set = (row, col, value) => {
+        k[row - 1][col - 1] = value;
+        k[col - 1][row - 1] = value;
+    };
+    const ax = (E * A) / L;
+    const tor = (G * J) / L;
+    // U2–R3 bending (I33)
+    const a3 = (12 * E * I33) / L ** 3;
+    const b3 = (6 * E * I33) / L ** 2;
+    const c3 = (4 * E * I33) / L;
+    const d3 = (2 * E * I33) / L;
+    // U3–R2 bending (I22)
+    const a2 = (12 * E * I22) / L ** 3;
+    const b2 = (6 * E * I22) / L ** 2;
+    const c2 = (4 * E * I22) / L;
+    const d2 = (2 * E * I22) / L;
+    set(1, 1, ax);
+    set(1, 7, -ax);
+    set(7, 7, ax);
+    set(2, 2, a3);
+    set(2, 6, -b3);
+    set(2, 8, -a3);
+    set(2, 12, -b3);
+    set(6, 6, c3);
+    set(6, 8, b3);
+    set(6, 12, d3);
+    set(8, 8, a3);
+    set(8, 12, b3);
+    set(12, 12, c3);
+    set(3, 3, a2);
+    set(3, 5, b2);
+    set(3, 9, -a2);
+    set(3, 11, b2);
+    set(5, 5, c2);
+    set(5, 9, -b2);
+    set(5, 11, d2);
+    set(9, 9, a2);
+    set(9, 11, -b2);
+    set(11, 11, c2);
+    set(4, 4, tor);
+    set(4, 10, -tor);
+    set(10, 10, tor);
+    return k;
+};
+exports.beam3dLocalStiffness = beam3dLocalStiffness;
+/**
+ * Converts a local matrix between the reference rotation convention and the right-hand rule
+ * (`D·k·D`, with `D` flipping the sign of R2 and R3 at both ends). The conversion is its own inverse.
+ */
+const flipBendingRotationSigns = (k) => {
+    const flipped = new Set(FLIPPED_DOFS);
+    return k.map((row, i) => row.map((value, j) => (flipped.has(i) !== flipped.has(j) ? -value : value)));
+};
+exports.flipBendingRotationSigns = flipBendingRotationSigns;
+/** Same element as `beam3dLocalStiffness` with all rotations following the right-hand rule. */
+const beam3dLocalStiffnessRightHanded = (E, G, A, I22, I33, J, L) => (0, exports.flipBendingRotationSigns)((0, exports.beam3dLocalStiffness)(E, G, A, I22, I33, J, L));
+exports.beam3dLocalStiffnessRightHanded = beam3dLocalStiffnessRightHanded;
+//# sourceMappingURL=beam-3d-stiffness.js.map

package/dist/fea/frame/elements/beam-3d-stiffness.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"beam-3d-stiffness.js","sourceRoot":"","sources":["../../../../src/fea/frame/elements/beam-3d-stiffness.ts"],"names":[],"mappings":";;;AAAA,8CAAuD;AAEvD;;;GAGG;AACH,MAAM,YAAY,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;AAEpC;;;;;;;;;;;;;;;;;;;GAmBG;AACI,MAAM,oBAAoB,GAAG,CAClC,CAAS,EACT,CAAS,EACT,CAAS,EACT,GAAW,EACX,GAAW,EACX,CAAS,EACT,CAAS,EACD,EAAE;IACV,MAAM,CAAC,GAAG,IAAA,cAAK,EAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACxB,yFAAyF;IACzF,MAAM,GAAG,GAAG,CAAC,GAAW,EAAE,GAAW,EAAE,KAAa,EAAE,EAAE;QACtD,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QAC5B,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;IAC9B,CAAC,CAAC;IAEF,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IAExB,sBAAsB;IACtB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,sBAAsB;IACtB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAE7B,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IAEd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAEhB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACd,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAEhB,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACjB,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IAEjB,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AA7DW,QAAA,oBAAoB,wBA6D/B;AAEF;;;GAGG;AACI,MAAM,wBAAwB,GAAG,CAAC,CAAS,EAAU,EAAE;IAC5D,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,YAAY,CAAC,CAAC;IACtC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CACtB,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAC5E,CAAC;AACJ,CAAC,CAAC;AALW,QAAA,wBAAwB,4BAKnC;AAEF,+FAA+F;AACxF,MAAM,+BAA+B,GAAG,CAC7C,CAAS,EACT,CAAS,EACT,CAAS,EACT,GAAW,EACX,GAAW,EACX,CAAS,EACT,CAAS,EACD,EAAE,CAAC,IAAA,gCAAwB,EAAC,IAAA,4BAAoB,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AARxE,QAAA,+BAA+B,mCAQyC"}

package/dist/fea/frame/elements/fixed-end-forces.d.ts

@@ -0,0 +1,22 @@
+/** A member load resolved into local axes, in N and mm. */
+export type LocalMemberLoad = {
+    type: "udl";
+    axis: 1 | 2 | 3;
+    w: number;
+    start: number;
+    end: number;
+} | {
+    type: "point";
+    axis: 1 | 2 | 3;
+    P: number;
+    a: number;
+};
+/**
+ * Fixed-end forces of a member under local loads: the end forces (right-hand local axes,
+ * N and N·mm) that hold both ends fully fixed. Equivalent nodal loads are the negative.
+ *
+ * @param loads - Loads in local axes, N and mm
+ * @param L - Member length in mm
+ * @throws {FrameAnalysisError} 308 if a load lies outside the member
+ */
+export declare const fixedEndForces: (loads: LocalMemberLoad[], L: number) => number[];

package/dist/fea/frame/elements/fixed-end-forces.js

@@ -0,0 +1,79 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.fixedEndForces = void 0;
+const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.error");
+// 4-point Gauss–Legendre rule, exact for the cubic integrands used here
+const GAUSS_POINTS = [
+    -0.8611363115940526, -0.3399810435848563, 0.3399810435848563, 0.8611363115940526,
+];
+const GAUSS_WEIGHTS = [
+    0.3478548451374538, 0.6521451548625461, 0.6521451548625461, 0.3478548451374538,
+];
+/**
+ * Equivalent nodal loads of a unit transverse/axial load at distance `x`, from the
+ * Euler–Bernoulli shape functions: the exact fixed-end solution for a prismatic member.
+ * Returned per unit of load, right-hand local DOF order [i: U1..R3, j: U1..R3].
+ */
+const unitEquivalentLoad = (axis, x, L) => {
+    const xi = x / L;
+    const out = new Array(12).fill(0);
+    if (axis === 1) {
+        out[0] = 1 - xi;
+        out[6] = xi;
+        return out;
+    }
+    const n1 = 1 - 3 * xi ** 2 + 2 * xi ** 3;
+    const n2 = L * (xi - 2 * xi ** 2 + xi ** 3);
+    const n3 = 3 * xi ** 2 - 2 * xi ** 3;
+    const n4 = L * (-(xi ** 2) + xi ** 3);
+    if (axis === 2) {
+        // load along local 2: force U2, moment about 3 (R3 = +dv/dx)
+        out[1] = n1;
+        out[5] = n2;
+        out[7] = n3;
+        out[11] = n4;
+    }
+    else {
+        // load along local 3: force U3, moment about 2 (R2 = −dw/dx)
+        out[2] = n1;
+        out[4] = -n2;
+        out[8] = n3;
+        out[10] = -n4;
+    }
+    return out;
+};
+/**
+ * Fixed-end forces of a member under local loads: the end forces (right-hand local axes,
+ * N and N·mm) that hold both ends fully fixed. Equivalent nodal loads are the negative.
+ *
+ * @param loads - Loads in local axes, N and mm
+ * @param L - Member length in mm
+ * @throws {FrameAnalysisError} 308 if a load lies outside the member
+ */
+const fixedEndForces = (loads, L) => {
+    const equivalent = new Array(12).fill(0);
+    const add = (unit, factor) => {
+        for (let i = 0; i < 12; i++)
+            equivalent[i] += unit[i] * factor;
+    };
+    for (const load of loads) {
+        if (load.type === "point") {
+            if (!(load.a >= 0 && load.a <= L))
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.LOAD_OUTSIDE_MEMBER);
+            add(unitEquivalentLoad(load.axis, load.a, L), load.P);
+        }
+        else {
+            if (!(load.start >= 0 && load.end <= L && load.start < load.end)) {
+                throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.LOAD_OUTSIDE_MEMBER);
+            }
+            const half = (load.end - load.start) / 2;
+            const mid = (load.end + load.start) / 2;
+            GAUSS_POINTS.forEach((p, index) => {
+                add(unitEquivalentLoad(load.axis, mid + half * p, L), load.w * half * GAUSS_WEIGHTS[index]);
+            });
+        }
+    }
+    return equivalent.map((v) => -v);
+};
+exports.fixedEndForces = fixedEndForces;
+//# sourceMappingURL=fixed-end-forces.js.map

package/dist/fea/frame/elements/fixed-end-forces.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"fixed-end-forces.js","sourceRoot":"","sources":["../../../../src/fea/frame/elements/fixed-end-forces.ts"],"names":[],"mappings":";;;AAAA,gFAAwF;AAOxF,wEAAwE;AACxE,MAAM,YAAY,GAAG;IACnB,CAAC,kBAAkB,EAAE,CAAC,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CACjF,CAAC;AACF,MAAM,aAAa,GAAG;IACpB,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CAC/E,CAAC;AAEF;;;;GAIG;AACH,MAAM,kBAAkB,GAAG,CAAC,IAAe,EAAE,CAAS,EAAE,CAAS,EAAY,EAAE;IAC7E,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QAChB,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,OAAO,GAAG,CAAC;IACb,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IAC5C,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACrC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC;IACtC,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;IACf,CAAC;SAAM,CAAC;QACN,6DAA6D;QAC7D,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QACb,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC;IAChB,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAEF;;;;;;;GAOG;AACI,MAAM,cAAc,GAAG,CAAC,KAAwB,EAAE,CAAS,EAAY,EAAE;IAC9E,MAAM,UAAU,GAAG,IAAI,KAAK,CAAS,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,GAAG,GAAG,CAAC,IAAc,EAAE,MAAc,EAAE,EAAE;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,UAAU,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;IACjE,CAAC,CAAC;IAEF,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;YAC1B,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;gBAC/B,MAAM,IAAI,yCAAkB,CAAC,kCAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,GAAG,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACjE,MAAM,IAAI,yCAAkB,CAAC,kCAAW,CAAC,mBAAmB,CAAC,CAAC;YAChE,CAAC;YACD,MAAM,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACzC,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YACxC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE;gBAChC,GAAG,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,GAAG,GAAG,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;YAC9F,CAAC,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,OAAO,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AACnC,CAAC,CAAC;AAvBW,QAAA,cAAc,kBAuBzB"}

package/dist/fea/frame/elements/static-condensation.d.ts

@@ -0,0 +1,19 @@
+import { type Matrix } from "../../../utils/matrix";
+export interface CondensedMember {
+    k: Matrix;
+    /** Condenses a 12-term fixed-end force vector; the result is zero at released DOFs. */
+    condenseFixedEndForces: (fixed: number[]) => number[];
+}
+/**
+ * Statically condenses released local DOFs out of a member, so those DOFs carry no force.
+ *
+ * With released set r and retained set c:
+ * `k* = k_cc − k_cr·k_rr⁻¹·k_rc`, `f* = f_c − k_cr·k_rr⁻¹·f_r`.
+ * The stiffness is condensed once; the returned function condenses the fixed-end forces of
+ * each load case.
+ *
+ * @param k - 12×12 local stiffness (right-hand)
+ * @param released - Local DOF indices (0–11) to release
+ * @throws {FrameAnalysisError} 307 if the released set makes the member a mechanism
+ */
+export declare const condenseMember: (k: Matrix, released: number[]) => CondensedMember;

package/dist/fea/frame/elements/static-condensation.js

@@ -0,0 +1,59 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.condenseMember = void 0;
+const frame_analysis_error_1 = require("../../../core/errors/frame-analysis.error");
+const matrix_1 = require("../../../utils/matrix");
+/**
+ * Statically condenses released local DOFs out of a member, so those DOFs carry no force.
+ *
+ * With released set r and retained set c:
+ * `k* = k_cc − k_cr·k_rr⁻¹·k_rc`, `f* = f_c − k_cr·k_rr⁻¹·f_r`.
+ * The stiffness is condensed once; the returned function condenses the fixed-end forces of
+ * each load case.
+ *
+ * @param k - 12×12 local stiffness (right-hand)
+ * @param released - Local DOF indices (0–11) to release
+ * @throws {FrameAnalysisError} 307 if the released set makes the member a mechanism
+ */
+const condenseMember = (k, released) => {
+    if (released.length === 0) {
+        return { k, condenseFixedEndForces: (fixed) => [...fixed] };
+    }
+    const retained = Array.from({ length: 12 }, (_, i) => i).filter((i) => !released.includes(i));
+    const sub = (rows, cols) => rows.map((r) => cols.map((c) => k[r][c]));
+    const krr = sub(released, released);
+    const krc = sub(released, retained);
+    const kcr = sub(retained, released);
+    const kcc = sub(retained, retained);
+    // [k_rr⁻¹·k_rc | k_rr⁻¹]
+    const rhs = released.map((_, i) => [...krc[i], ...released.map((__, j) => (i === j ? 1 : 0))]);
+    const solved = (0, matrix_1.solveLinearSystem)(krr, rhs);
+    if (!solved)
+        throw new frame_analysis_error_1.FrameAnalysisError(frame_analysis_error_1.FrameErrors.INVALID_RELEASE);
+    const nc = retained.length;
+    const kStar = k.map((row) => row.map(() => 0));
+    retained.forEach((r, ri) => {
+        retained.forEach((c, ci) => {
+            let value = kcc[ri][ci];
+            released.forEach((_, j) => {
+                value -= kcr[ri][j] * solved[j][ci];
+            });
+            kStar[r][c] = value;
+        });
+    });
+    // k_cr·k_rr⁻¹, applied to f_r for every load case
+    const kcrKrrInv = retained.map((_, ri) => released.map((__, m) => released.reduce((s, ___, j) => s + kcr[ri][j] * solved[j][nc + m], 0)));
+    return {
+        k: kStar,
+        condenseFixedEndForces: (fixed) => {
+            const fStar = new Array(12).fill(0);
+            const fr = released.map((r) => fixed[r]);
+            retained.forEach((r, ri) => {
+                fStar[r] = fixed[r] - kcrKrrInv[ri].reduce((s, v, m) => s + v * fr[m], 0);
+            });
+            return fStar;
+        },
+    };
+};
+exports.condenseMember = condenseMember;
+//# sourceMappingURL=static-condensation.js.map