@graphty/layout 1.2.3 → 1.2.4

package/dist/src/algorithms/index.d.ts

@@ -0,0 +1,5 @@
+/**
+ * Re-export all algorithm utilities
+ */
+export * from './planarity';
+export * from './optimization';

package/dist/src/algorithms/index.js

@@ -0,0 +1,6 @@
+/**
+ * Re-export all algorithm utilities
+ */
+export * from './planarity';
+export * from './optimization';
+//# sourceMappingURL=index.js.map

package/dist/src/algorithms/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/algorithms/index.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,cAAc,aAAa,CAAC;AAC5B,cAAc,gBAAgB,CAAC"}

package/dist/src/algorithms/optimization/index.d.ts

@@ -0,0 +1,7 @@
+/**
+ * Re-export all optimization algorithms
+ */
+export type { DistanceMap } from './types';
+export { _computeShortestPathDistances, _kamadaKawaiSolve, _kamadaKawaiCostfn } from './kamada-kawai-solver';
+export { _lbfgsDirection } from './lbfgs';
+export { _backtrackingLineSearch } from './line-search';

package/dist/src/algorithms/optimization/index.js

@@ -0,0 +1,7 @@
+/**
+ * Re-export all optimization algorithms
+ */
+export { _computeShortestPathDistances, _kamadaKawaiSolve, _kamadaKawaiCostfn } from './kamada-kawai-solver';
+export { _lbfgsDirection } from './lbfgs';
+export { _backtrackingLineSearch } from './line-search';
+//# sourceMappingURL=index.js.map

package/dist/src/algorithms/optimization/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../../src/algorithms/optimization/index.ts"],"names":[],"mappings":"AAAA;;GAEG;AAGH,OAAO,EACL,6BAA6B,EAC7B,iBAAiB,EACjB,kBAAkB,EACnB,MAAM,uBAAuB,CAAC;AAC/B,OAAO,EAAE,eAAe,EAAE,MAAM,SAAS,CAAC;AAC1C,OAAO,EAAE,uBAAuB,EAAE,MAAM,eAAe,CAAC"}

package/dist/src/algorithms/optimization/kamada-kawai-solver.d.ts

@@ -0,0 +1,32 @@
+/**
+ * Kamada-Kawai layout algorithm optimization functions
+ */
+import { Graph } from '../../types';
+import { DistanceMap } from './types';
+/**
+ * Compute all-pairs shortest path distances for the graph
+ *
+ * @param G - NetworkX graph
+ * @param weight - Edge attribute for weight
+ * @returns Dictionary of dictionaries of shortest path distances
+ */
+export declare function _computeShortestPathDistances(G: Graph, weight: string): DistanceMap;
+/**
+ * Solve the Kamada-Kawai layout optimization problem
+ *
+ * @param distMatrix - Matrix of desired distances between nodes
+ * @param positions - Initial node positions
+ * @param dim - Dimension of layout
+ * @returns Optimized node positions
+ */
+export declare function _kamadaKawaiSolve(distMatrix: number[][], positions: number[][], dim: number): number[][];
+/**
+ * Cost function and gradient for Kamada-Kawai layout algorithm
+ *
+ * @param posVec - Flattened position array
+ * @param invDist - Inverse distance matrix
+ * @param meanWeight - Weight for centering positions
+ * @param dim - Dimension of layout
+ * @returns Array with [cost, gradient]
+ */
+export declare function _kamadaKawaiCostfn(posVec: number[], invDist: number[][], meanWeight: number, dim: number): [number, number[]];

package/dist/src/algorithms/optimization/kamada-kawai-solver.js

@@ -0,0 +1,176 @@
+/**
+ * Kamada-Kawai layout algorithm optimization functions
+ */
+import { _lbfgsDirection } from './lbfgs';
+import { _backtrackingLineSearch } from './line-search';
+import { getNodesFromGraph, getEdgesFromGraph } from '../../utils/graph';
+/**
+ * Compute all-pairs shortest path distances for the graph
+ *
+ * @param G - NetworkX graph
+ * @param weight - Edge attribute for weight
+ * @returns Dictionary of dictionaries of shortest path distances
+ */
+export function _computeShortestPathDistances(G, weight) {
+    const distances = {};
+    const nodes = getNodesFromGraph(G);
+    const edges = getEdgesFromGraph(G);
+    // Initialize distances with direct edges
+    for (const node of nodes) {
+        distances[node] = {};
+        distances[node][node] = 0;
+        for (const other of nodes) {
+            if (node !== other) {
+                distances[node][other] = Infinity;
+            }
+        }
+    }
+    // Add direct edges
+    for (const [source, target] of edges) {
+        // In a real implementation, we would get the weight from the graph
+        // For now, assume weight = 1 or use weight attribute if available
+        let edgeWeight = 1;
+        if (G.getEdgeData) {
+            edgeWeight = G.getEdgeData(source, target, weight) || 1;
+        }
+        distances[source][target] = edgeWeight;
+        distances[target][source] = edgeWeight; // Assuming undirected graph
+    }
+    // Floyd-Warshall algorithm for all-pairs shortest paths
+    for (const k of nodes) {
+        for (const i of nodes) {
+            for (const j of nodes) {
+                if (distances[i][k] + distances[k][j] < distances[i][j]) {
+                    distances[i][j] = distances[i][k] + distances[k][j];
+                }
+            }
+        }
+    }
+    return distances;
+}
+/**
+ * Solve the Kamada-Kawai layout optimization problem
+ *
+ * @param distMatrix - Matrix of desired distances between nodes
+ * @param positions - Initial node positions
+ * @param dim - Dimension of layout
+ * @returns Optimized node positions
+ */
+export function _kamadaKawaiSolve(distMatrix, positions, dim) {
+    // Implementation of L-BFGS optimization for Kamada-Kawai
+    const nNodes = positions.length;
+    const meanWeight = 1e-3;
+    // Convert distances to inverse distances (with protection against division by zero)
+    const invDistMatrix = distMatrix.map(row => row.map(d => d === 0 ? 0 : 1 / (d + 1e-3)));
+    // Flatten positions for optimization
+    let posVec = positions.flat();
+    // Optimization parameters
+    const maxIter = 500;
+    const gtol = 1e-5;
+    const m = 10; // L-BFGS memory size
+    // Implement a simplified L-BFGS-B algorithm
+    let alpha = 1.0;
+    const oldValues = [];
+    const oldGrads = [];
+    for (let iter = 0; iter < maxIter; iter++) {
+        // Calculate cost and gradient
+        const [cost, grad] = _kamadaKawaiCostfn(posVec, invDistMatrix, meanWeight, dim);
+        // Compute search direction using L-BFGS approximation
+        const direction = _lbfgsDirection(grad, oldValues, oldGrads, m);
+        // Simple line search for step size
+        alpha = _backtrackingLineSearch(posVec, direction, cost, grad, (x) => _kamadaKawaiCostfn(x, invDistMatrix, meanWeight, dim)[0], alpha);
+        // Save current position and gradient for next iteration
+        const oldPos = [...posVec];
+        // Update position
+        for (let i = 0; i < posVec.length; i++) {
+            posVec[i] += alpha * direction[i];
+        }
+        // Calculate new gradient
+        const [, newGrad] = _kamadaKawaiCostfn(posVec, invDistMatrix, meanWeight, dim);
+        // Update L-BFGS memory
+        oldValues.push(posVec.map((val, i) => val - oldPos[i]));
+        oldGrads.push(newGrad.map((val, i) => val - grad[i]));
+        // Keep only m most recent updates
+        if (oldValues.length > m) {
+            oldValues.shift();
+            oldGrads.shift();
+        }
+        // Check convergence
+        const gradNorm = Math.sqrt(newGrad.reduce((sum, g) => sum + g * g, 0));
+        if (gradNorm < gtol) {
+            break;
+        }
+    }
+    // Reshape result back into positions array
+    const result = [];
+    for (let i = 0; i < nNodes; i++) {
+        result.push(posVec.slice(i * dim, (i + 1) * dim));
+    }
+    return result;
+}
+/**
+ * Cost function and gradient for Kamada-Kawai layout algorithm
+ *
+ * @param posVec - Flattened position array
+ * @param invDist - Inverse distance matrix
+ * @param meanWeight - Weight for centering positions
+ * @param dim - Dimension of layout
+ * @returns Array with [cost, gradient]
+ */
+export function _kamadaKawaiCostfn(posVec, invDist, meanWeight, dim) {
+    const nNodes = invDist.length;
+    const positions = [];
+    // Reshape flat vector into positions array
+    for (let i = 0; i < nNodes; i++) {
+        positions.push(posVec.slice(i * dim, (i + 1) * dim));
+    }
+    // Calculate cost
+    let cost = 0;
+    // Add mean position penalty term
+    const sumPos = Array(dim).fill(0);
+    for (let i = 0; i < nNodes; i++) {
+        for (let d = 0; d < dim; d++) {
+            sumPos[d] += positions[i][d];
+        }
+    }
+    cost += 0.5 * meanWeight * sumPos.reduce((sum, val) => sum + val * val, 0);
+    // Add distance penalty terms
+    for (let i = 0; i < nNodes; i++) {
+        for (let j = i + 1; j < nNodes; j++) {
+            // Calculate actual distance
+            const diff = positions[i].map((val, d) => val - positions[j][d]);
+            const distance = Math.sqrt(diff.reduce((sum, d) => sum + d * d, 0));
+            // Add penalty for difference between actual and ideal distance
+            const idealInvDist = invDist[i][j];
+            const offset = distance * idealInvDist - 1.0;
+            cost += 0.5 * offset * offset;
+        }
+    }
+    // Calculate gradient
+    const grad = new Array(posVec.length).fill(0);
+    // Add gradient of mean position penalty
+    for (let i = 0; i < nNodes; i++) {
+        for (let d = 0; d < dim; d++) {
+            grad[i * dim + d] += meanWeight * sumPos[d];
+        }
+    }
+    // Add gradient of distance penalties
+    for (let i = 0; i < nNodes; i++) {
+        for (let j = i + 1; j < nNodes; j++) {
+            // Calculate actual distance and direction
+            const diff = positions[i].map((val, d) => val - positions[j][d]);
+            const distance = Math.sqrt(diff.reduce((sum, d) => sum + d * d, 0)) || 1e-10;
+            const direction = diff.map(d => d / distance);
+            // Calculate contribution to gradient
+            const idealInvDist = invDist[i][j];
+            const offset = distance * idealInvDist - 1.0;
+            for (let d = 0; d < dim; d++) {
+                const force = idealInvDist * offset * direction[d];
+                grad[i * dim + d] += force;
+                grad[j * dim + d] -= force;
+            }
+        }
+    }
+    return [cost, grad];
+}
+//# sourceMappingURL=kamada-kawai-solver.js.map

package/dist/src/algorithms/optimization/kamada-kawai-solver.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"kamada-kawai-solver.js","sourceRoot":"","sources":["../../../../src/algorithms/optimization/kamada-kawai-solver.ts"],"names":[],"mappings":"AAAA;;GAEG;AAIH,OAAO,EAAE,eAAe,EAAE,MAAM,SAAS,CAAC;AAC1C,OAAO,EAAE,uBAAuB,EAAE,MAAM,eAAe,CAAC;AACxD,OAAO,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,MAAM,mBAAmB,CAAC;AAEzE;;;;;;GAMG;AACH,MAAM,UAAU,6BAA6B,CAC3C,CAAQ,EACR,MAAc;IAEd,MAAM,SAAS,GAAgB,EAAE,CAAC;IAClC,MAAM,KAAK,GAAG,iBAAiB,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,iBAAiB,CAAC,CAAC,CAAC,CAAC;IAEnC,yCAAyC;IACzC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,SAAS,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC;QACrB,SAAS,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAE1B,KAAK,MAAM,KAAK,IAAI,KAAK,EAAE,CAAC;YAC1B,IAAI,IAAI,KAAK,KAAK,EAAE,CAAC;gBACnB,SAAS,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,GAAG,QAAQ,CAAC;YACpC,CAAC;QACH,CAAC;IACH,CAAC;IAED,mBAAmB;IACnB,KAAK,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC;QACrC,mEAAmE;QACnE,kEAAkE;QAClE,IAAI,UAAU,GAAG,CAAC,CAAC;QACnB,IAAI,CAAC,CAAC,WAAW,EAAE,CAAC;YAClB,UAAU,GAAG,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC;QAC1D,CAAC;QAED,SAAS,CAAC,MAAM,CAAC,CAAC,MAAM,CAAC,GAAG,UAAU,CAAC;QACvC,SAAS,CAAC,MAAM,CAAC,CAAC,MAAM,CAAC,GAAG,UAAU,CAAC,CAAE,4BAA4B;IACvE,CAAC;IAED,wDAAwD;IACxD,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC;QACtB,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC;YACtB,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC;gBACtB,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;oBACxD,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACtD,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,iBAAiB,CAC/B,UAAsB,EACtB,SAAqB,EACrB,GAAW;IAEX,yDAAyD;IACzD,MAAM,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC;IAChC,MAAM,UAAU,GAAG,IAAI,CAAC;IAExB,oFAAoF;IACpF,MAAM,aAAa,GAAG,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CACzC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAC3C,CAAC;IAEF,qCAAqC;IACrC,IAAI,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;IAE9B,0BAA0B;IAC1B,MAAM,OAAO,GAAG,GAAG,CAAC;IACpB,MAAM,IAAI,GAAG,IAAI,CAAC;IAClB,MAAM,CAAC,GAAG,EAAE,CAAC,CAAE,qBAAqB;IAEpC,4CAA4C;IAC5C,IAAI,KAAK,GAAG,GAAG,CAAC;IAChB,MAAM,SAAS,GAAe,EAAE,CAAC;IACjC,MAAM,QAAQ,GAAe,EAAE,CAAC;IAEhC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,OAAO,EAAE,IAAI,EAAE,EAAE,CAAC;QAC1C,8BAA8B;QAC9B,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,kBAAkB,CAAC,MAAM,EAAE,aAAa,EAAE,UAAU,EAAE,GAAG,CAAC,CAAC;QAEhF,sDAAsD;QACtD,MAAM,SAAS,GAAG,eAAe,CAAC,IAAI,EAAE,SAAS,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;QAEhE,mCAAmC;QACnC,KAAK,GAAG,uBAAuB,CAC7B,MAAM,EAAE,SAAS,EAAE,IAAI,EAAE,IAAI,EAC7B,CAAC,CAAW,EAAE,EAAE,CAAC,kBAAkB,CAAC,CAAC,EAAE,aAAa,EAAE,UAAU,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EACzE,KAAK,CACN,CAAC;QAEF,wDAAwD;QACxD,MAAM,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;QAE3B,kBAAkB;QAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,MAAM,CAAC,CAAC,CAAC,IAAI,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QACpC,CAAC;QAED,yBAAyB;QACzB,MAAM,CAAC,EAAE,OAAO,CAAC,GAAG,kBAAkB,CAAC,MAAM,EAAE,aAAa,EAAE,UAAU,EAAE,GAAG,CAAC,CAAC;QAE/E,uBAAuB;QACvB,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxD,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAEtD,kCAAkC;QAClC,IAAI,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACzB,SAAS,CAAC,KAAK,EAAE,CAAC;YAClB,QAAQ,CAAC,KAAK,EAAE,CAAC;QACnB,CAAC;QAED,oBAAoB;QACpB,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACvE,IAAI,QAAQ,GAAG,IAAI,EAAE,CAAC;YACpB,MAAM;QACR,CAAC;IACH,CAAC;IAED,2CAA2C;IAC3C,MAAM,MAAM,GAAe,EAAE,CAAC;IAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACpD,CAAC;IAED,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,kBAAkB,CAChC,MAAgB,EAChB,OAAmB,EACnB,UAAkB,EAClB,GAAW;IAEX,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;IAC9B,MAAM,SAAS,GAAe,EAAE,CAAC;IAEjC,2CAA2C;IAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACvD,CAAC;IAED,iBAAiB;IACjB,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,iCAAiC;IACjC,MAAM,MAAM,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,MAAM,CAAC,CAAC,CAAC,IAAI,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/B,CAAC;IACH,CAAC;IACD,IAAI,IAAI,GAAG,GAAG,UAAU,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;IAE3E,6BAA6B;IAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,4BAA4B;YAC5B,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACjE,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAEpE,+DAA+D;YAC/D,MAAM,YAAY,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnC,MAAM,MAAM,GAAG,QAAQ,GAAG,YAAY,GAAG,GAAG,CAAC;YAC7C,IAAI,IAAI,GAAG,GAAG,MAAM,GAAG,MAAM,CAAC;QAChC,CAAC;IACH,CAAC;IAED,qBAAqB;IACrB,MAAM,IAAI,GAAG,IAAI,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAE9C,wCAAwC;IACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,IAAI,UAAU,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAC9C,CAAC;IACH,CAAC;IAED,qCAAqC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,0CAA0C;YAC1C,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACjE,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;YAC7E,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC;YAE9C,qCAAqC;YACrC,MAAM,YAAY,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnC,MAAM,MAAM,GAAG,QAAQ,GAAG,YAAY,GAAG,GAAG,CAAC;YAE7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,MAAM,KAAK,GAAG,YAAY,GAAG,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBACnD,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,IAAI,KAAK,CAAC;gBAC3B,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,IAAI,KAAK,CAAC;YAC7B,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;AACtB,CAAC"}

package/dist/src/algorithms/optimization/lbfgs.d.ts

@@ -0,0 +1,13 @@
+/**
+ * L-BFGS optimization utilities
+ */
+/**
+ * Compute the search direction using L-BFGS approximation
+ *
+ * @param grad - Current gradient
+ * @param sList - List of position differences (s_k)
+ * @param yList - List of gradient differences (y_k)
+ * @param m - Memory size
+ * @returns Direction vector
+ */
+export declare function _lbfgsDirection(grad: number[], sList: number[][], yList: number[][], m: number): number[];

package/dist/src/algorithms/optimization/lbfgs.js

@@ -0,0 +1,56 @@
+/**
+ * L-BFGS optimization utilities
+ */
+/**
+ * Compute the search direction using L-BFGS approximation
+ *
+ * @param grad - Current gradient
+ * @param sList - List of position differences (s_k)
+ * @param yList - List of gradient differences (y_k)
+ * @param m - Memory size
+ * @returns Direction vector
+ */
+export function _lbfgsDirection(grad, sList, yList, m) {
+    if (sList.length === 0) {
+        // First iteration - use negative gradient
+        return grad.map(g => -g);
+    }
+    const q = grad.slice();
+    const alpha = Array(sList.length).fill(0);
+    const rho = [];
+    // Compute rho values
+    for (let i = 0; i < sList.length; i++) {
+        const s = sList[i];
+        const y = yList[i];
+        rho.push(1 / y.reduce((sum, val, j) => sum + val * s[j], 0));
+    }
+    // Forward pass
+    for (let i = sList.length - 1; i >= 0; i--) {
+        const s = sList[i];
+        alpha[i] = rho[i] * s.reduce((sum, val, j) => sum + val * q[j], 0);
+        for (let j = 0; j < q.length; j++) {
+            q[j] -= alpha[i] * yList[i][j];
+        }
+    }
+    // Scale initial Hessian approximation
+    let gamma = 1;
+    if (sList.length > 0 && yList.length > 0) {
+        const y = yList[yList.length - 1];
+        const s = sList[sList.length - 1];
+        gamma = s.reduce((sum, val, i) => sum + val * y[i], 0) /
+            y.reduce((sum, val) => sum + val * val, 0);
+    }
+    // Initialize direction with scaled negative gradient
+    const direction = q.map(val => -gamma * val);
+    // Backward pass
+    for (let i = 0; i < sList.length; i++) {
+        const s = sList[i];
+        const y = yList[i];
+        const beta = rho[i] * y.reduce((sum, val, j) => sum + val * direction[j], 0);
+        for (let j = 0; j < direction.length; j++) {
+            direction[j] += s[j] * (alpha[i] - beta);
+        }
+    }
+    return direction;
+}
+//# sourceMappingURL=lbfgs.js.map

package/dist/src/algorithms/optimization/lbfgs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"lbfgs.js","sourceRoot":"","sources":["../../../../src/algorithms/optimization/lbfgs.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH;;;;;;;;GAQG;AACH,MAAM,UAAU,eAAe,CAC7B,IAAc,EACd,KAAiB,EACjB,KAAiB,EACjB,CAAS;IAET,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACvB,0CAA0C;QAC1C,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3B,CAAC;IAED,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;IACvB,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,GAAG,GAAa,EAAE,CAAC;IAEzB,qBAAqB;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACnB,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC/D,CAAC;IAED,eAAe;IACf,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3C,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACnB,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACnE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,CAAC;IACH,CAAC;IAED,sCAAsC;IACtC,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACzC,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAClC,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAClC,KAAK,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACpD,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;IAC/C,CAAC;IAED,qDAAqD;IACrD,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC;IAE7C,gBAAgB;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QAC3C,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/dist/src/algorithms/optimization/line-search.d.ts

@@ -0,0 +1,15 @@
+/**
+ * Line search optimization utilities
+ */
+/**
+ * Backtracking line search to find step size
+ *
+ * @param x - Current position
+ * @param direction - Search direction
+ * @param f - Function value at current position
+ * @param grad - Gradient at current position
+ * @param func - Function to evaluate cost
+ * @param alpha0 - Initial step size
+ * @returns Optimal step size
+ */
+export declare function _backtrackingLineSearch(x: number[], direction: number[], f: number, grad: number[], func: (x: number[]) => number, alpha0: number): number;

package/dist/src/algorithms/optimization/line-search.js

@@ -0,0 +1,37 @@
+/**
+ * Line search optimization utilities
+ */
+/**
+ * Backtracking line search to find step size
+ *
+ * @param x - Current position
+ * @param direction - Search direction
+ * @param f - Function value at current position
+ * @param grad - Gradient at current position
+ * @param func - Function to evaluate cost
+ * @param alpha0 - Initial step size
+ * @returns Optimal step size
+ */
+export function _backtrackingLineSearch(x, direction, f, grad, func, alpha0) {
+    const c1 = 1e-4;
+    const c2 = 0.9;
+    const initialSlope = grad.reduce((sum, g, i) => sum + g * direction[i], 0);
+    if (initialSlope >= 0) {
+        return 1e-8; // Direction is not a descent direction
+    }
+    let alpha = alpha0;
+    const maxIter = 20;
+    for (let i = 0; i < maxIter; i++) {
+        // Try step
+        const newX = x.map((val, i) => val + alpha * direction[i]);
+        const newF = func(newX);
+        // Check sufficient decrease condition (Armijo condition)
+        if (newF <= f + c1 * alpha * initialSlope) {
+            return alpha;
+        }
+        // Reduce step size
+        alpha *= c2;
+    }
+    return alpha; // Return last alpha even if not optimal
+}
+//# sourceMappingURL=line-search.js.map

package/dist/src/algorithms/optimization/line-search.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"line-search.js","sourceRoot":"","sources":["../../../../src/algorithms/optimization/line-search.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH;;;;;;;;;;GAUG;AACH,MAAM,UAAU,uBAAuB,CACrC,CAAW,EACX,SAAmB,EACnB,CAAS,EACT,IAAc,EACd,IAA6B,EAC7B,MAAc;IAEd,MAAM,EAAE,GAAG,IAAI,CAAC;IAChB,MAAM,EAAE,GAAG,GAAG,CAAC;IACf,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAE3E,IAAI,YAAY,IAAI,CAAC,EAAE,CAAC;QACtB,OAAO,IAAI,CAAC,CAAE,uCAAuC;IACvD,CAAC;IAED,IAAI,KAAK,GAAG,MAAM,CAAC;IACnB,MAAM,OAAO,GAAG,EAAE,CAAC;IAEnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,WAAW;QACX,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3D,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;QAExB,yDAAyD;QACzD,IAAI,IAAI,IAAI,CAAC,GAAG,EAAE,GAAG,KAAK,GAAG,YAAY,EAAE,CAAC;YAC1C,OAAO,KAAK,CAAC;QACf,CAAC;QAED,mBAAmB;QACnB,KAAK,IAAI,EAAE,CAAC;IACd,CAAC;IAED,OAAO,KAAK,CAAC,CAAE,wCAAwC;AACzD,CAAC"}

package/dist/src/algorithms/optimization/types.d.ts

@@ -0,0 +1,5 @@
+/**
+ * Type definitions for optimization algorithms
+ */
+import { Node } from '../../types';
+export type DistanceMap = Record<Node, Record<Node, number>>;

package/dist/src/algorithms/optimization/types.js

@@ -0,0 +1,5 @@
+/**
+ * Type definitions for optimization algorithms
+ */
+export {};
+//# sourceMappingURL=types.js.map

package/dist/src/algorithms/optimization/types.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"types.js","sourceRoot":"","sources":["../../../../src/algorithms/optimization/types.ts"],"names":[],"mappings":"AAAA;;GAEG"}

package/dist/src/algorithms/planarity/check.d.ts

@@ -0,0 +1,17 @@
+/**
+ * Main planarity checking function
+ */
+import { Graph, Node, Edge, Embedding } from '../../types';
+/**
+ * Check if graph is planar using a simplified version of Boyer-Myrvold algorithm.
+ * Returns planarity and embedding information.
+ *
+ * @param G - Graph
+ * @param nodes - List of nodes
+ * @param edges - List of edges
+ * @returns Object containing isPlanar flag and embedding
+ */
+export declare function checkPlanarity(G: Graph, nodes: Node[], edges: Edge[]): {
+    isPlanar: boolean;
+    embedding: Embedding | null;
+};

package/dist/src/algorithms/planarity/check.js

@@ -0,0 +1,30 @@
+/**
+ * Main planarity checking function
+ */
+import { isK5, isK33 } from './special-graphs';
+import { lrPlanarityTest } from './lr-test';
+import { createTriangulationEmbedding } from './embedding';
+/**
+ * Check if graph is planar using a simplified version of Boyer-Myrvold algorithm.
+ * Returns planarity and embedding information.
+ *
+ * @param G - Graph
+ * @param nodes - List of nodes
+ * @param edges - List of edges
+ * @returns Object containing isPlanar flag and embedding
+ */
+export function checkPlanarity(G, nodes, edges) {
+    // For small graphs (n <= 4), all are planar
+    if (nodes.length <= 4) {
+        return { isPlanar: true, embedding: createTriangulationEmbedding(nodes, edges) };
+    }
+    // For K5 (complete graph with 5 nodes) and K3,3 (complete bipartite with 3,3 nodes)
+    // these are not planar by Kuratowski's theorem
+    if (isK5(nodes, edges) || isK33(nodes, edges)) {
+        return { isPlanar: false, embedding: null };
+    }
+    // For other graphs, use LR algorithm (Left-Right Planarity Test)
+    const result = lrPlanarityTest(nodes, edges);
+    return result;
+}
+//# sourceMappingURL=check.js.map

package/dist/src/algorithms/planarity/check.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"check.js","sourceRoot":"","sources":["../../../../src/algorithms/planarity/check.ts"],"names":[],"mappings":"AAAA;;GAEG;AAGH,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,MAAM,WAAW,CAAC;AAC5C,OAAO,EAAE,4BAA4B,EAAE,MAAM,aAAa,CAAC;AAE3D;;;;;;;;GAQG;AACH,MAAM,UAAU,cAAc,CAC5B,CAAQ,EACR,KAAa,EACb,KAAa;IAEb,4CAA4C;IAC5C,IAAI,KAAK,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACtB,OAAO,EAAE,QAAQ,EAAE,IAAI,EAAE,SAAS,EAAE,4BAA4B,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,CAAC;IACnF,CAAC;IAED,oFAAoF;IACpF,+CAA+C;IAC/C,IAAI,IAAI,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,CAAC;QAC9C,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,SAAS,EAAE,IAAI,EAAE,CAAC;IAC9C,CAAC;IAED,iEAAiE;IACjE,MAAM,MAAM,GAAG,eAAe,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IAC7C,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/dist/src/algorithms/planarity/embedding.d.ts

@@ -0,0 +1,29 @@
+/**
+ * Embedding functions for planar graphs
+ */
+import { Node, Edge, Embedding, PositionMap } from '../../types';
+/**
+ * Create a triangulation-based embedding for a planar graph
+ *
+ * @param nodes - List of nodes
+ * @param edges - List of edges
+ * @returns Embedding object
+ */
+export declare function createTriangulationEmbedding(nodes: Node[], edges: Edge[]): Embedding;
+/**
+ * Find a simple cycle in the graph (for outer face)
+ *
+ * @param nodes - List of nodes
+ * @param edges - List of edges
+ * @param adjMap - Adjacency map
+ * @returns Cycle as array of nodes, or null if none found
+ */
+export declare function findCycle(nodes: Node[], edges: Edge[], adjMap: Record<Node, Set<Node>>): Node[] | null;
+/**
+ * Convert a combinatorial embedding to node positions
+ *
+ * @param embedding - The embedding object
+ * @param nodes - List of nodes
+ * @returns Dictionary mapping nodes to positions
+ */
+export declare function combinatorialEmbeddingToPos(embedding: Embedding, nodes: Node[]): PositionMap;