@graphty/graph-samples 0.0.1 → 0.1.2

package/dist/src/generators/growth.js

@@ -0,0 +1,445 @@
+/**
+ * Sequential growth models: the random recursive tree, forest fire, duplication-divergence,
+ * Newman-Watts, Bianconi-Barabasi fitness, the random Apollonian network and Wilson's uniform
+ * spanning tree of a grid (a perfect maze). Each draws from one stream (seed, domain, 0): the
+ * processes are sequential, so there is no row-parallel form.
+ */
+import { detLog } from "../random/log.js";
+import { RandomStream, resolveSeed } from "../random/stream.js";
+import { checkEdgeCount, checkInt, checkProbability, EdgeBuffer, toGraph } from "./util.js";
+import { applyWeights } from "./weights.js";
+/**
+ * The random recursive tree (R. T. Smythe and H. M. Mahmoud, "A survey of recursive trees",
+ * Theory Probab. Math. Statist. 51, 1-27, 1995): node t = 1 .. n - 1 attaches to a uniform earlier
+ * node, drawn as nextBelow(t) from the stream (seed, "recursive-tree", 0). Edge t - 1 is
+ * (parent, t), in t order. O(n). Depth about ln n, root degree about ln n.
+ * @param options - n and seed
+ * @returns the undirected tree
+ */
+export function randomRecursiveTreeGraph(options) {
+    const { n } = options;
+    checkInt("n", n, 1);
+    const stream = new RandomStream(resolveSeed(options.seed), "recursive-tree", 0);
+    const src = new Uint32Array(n - 1);
+    const dst = new Uint32Array(n - 1);
+    for (let t = 1; t < n; t++) {
+        src[t - 1] = stream.nextBelow(t);
+        dst[t - 1] = t;
+    }
+    return applyWeights({ directed: false, nodeCount: n, src, dst }, options);
+}
+/**
+ * Pick `count` of `candidates` uniformly without replacement (a partial Fisher-Yates shuffle in
+ * place): the picks are candidates[0 .. count).
+ * @param stream - the stream
+ * @param candidates - the candidates, reordered in place
+ * @param count - how many to pick, <= candidates.length
+ */
+function pickPrefix(stream, candidates, count) {
+    for (let i = 0; i < count; i++) {
+        const j = i + stream.nextBelow(candidates.length - i);
+        const t = candidates[i];
+        candidates[i] = candidates[j];
+        candidates[j] = t;
+    }
+}
+/**
+ * The forest fire model (J. Leskovec, J. Kleinberg and C. Faloutsos, "Graphs over time:
+ * densification laws, shrinking diameters and possible explanations", KDD 2005,
+ * doi:10.1145/1081870.1081893). Directed; heavy-tailed degrees, densification, community locality.
+ *
+ * Node 0 starts alone. Each new node v = 1 .. n - 1 draws an ambassador w = nextBelow(v) and burns
+ * breadth-first from it: for each burning node x in burn order, draw x_out = floor(ln(1 - u) / ln p)
+ * (geometric, P(x_out >= k) = p^k, mean p / (1 - p)) and then x_in the same way with r p; then pick
+ * x_out of x's out-neighbours and then x_in of its in-neighbours uniformly among those not yet
+ * burned (adjacency order, partial Fisher-Yates; all of them when there are fewer), and burn them.
+ * Burning stops at `maxBurn` nodes. v then links to every burned node: arcs (v, t) in burn order.
+ * Edges are listed by v. All draws come from the stream (seed, "forest-fire", 0).
+ *
+ * Cost: O(sum over burning nodes of their degree), since each burning node's lists are scanned;
+ * the edge count explodes for p above about 0.37, which `maxBurn` bounds at n maxBurn arcs.
+ * @param options - n, forward, backward, maxBurn and seed
+ * @returns the directed graph
+ */
+export function forestFireGraph(options) {
+    const { n, forward: p, backward: r } = options;
+    const maxBurn = options.maxBurn ?? 1000;
+    checkInt("n", n, 1);
+    if (!(p >= 0 && p < 1)) {
+        throw new RangeError(`forward must be a probability in [0, 1), got ${String(p)}`);
+    }
+    if (!(r >= 0 && r * p < 1)) {
+        throw new RangeError(`backward must be >= 0 with backward * forward < 1, got ${String(r)}`);
+    }
+    checkInt("maxBurn", maxBurn, 1);
+    const stream = new RandomStream(resolveSeed(options.seed), "forest-fire", 0);
+    const pb = r * p;
+    const logForward = detLog(p);
+    const logBackward = detLog(pb);
+    // ponytail: number[][] adjacency, fine to a few million arcs; a growable CSR beyond
+    const out = Array.from({ length: n }, () => []);
+    const inn = Array.from({ length: n }, () => []);
+    const stamp = new Uint32Array(n);
+    const edges = new EdgeBuffer(4 * n);
+    const burned = [];
+    const burn = (v, list, count) => {
+        const candidates = list.filter((t) => stamp[t] !== v);
+        const take = Math.min(count, candidates.length, maxBurn - burned.length);
+        pickPrefix(stream, candidates, take);
+        for (let i = 0; i < take; i++) {
+            stamp[candidates[i]] = v;
+            burned.push(candidates[i]);
+        }
+    };
+    for (let v = 1; v < n; v++) {
+        burned.length = 0;
+        const w = stream.nextBelow(v);
+        stamp[w] = v;
+        burned.push(w);
+        for (let head = 0; head < burned.length && burned.length < maxBurn; head++) {
+            const x = burned[head];
+            const xOut = p === 0 ? 0 : stream.nextSkip(logForward);
+            const xIn = pb === 0 ? 0 : stream.nextSkip(logBackward);
+            burn(v, out[x], xOut);
+            burn(v, inn[x], xIn);
+        }
+        for (const t of burned) {
+            edges.push(v, t);
+            out[v].push(t);
+            inn[t].push(v);
+        }
+    }
+    return applyWeights(toGraph(n, edges, true), options);
+}
+/**
+ * The duplication-divergence model (I. Ispolatov, P. L. Krapivsky and A. Yuryev,
+ * "Duplication-divergence model of protein interaction network", Phys. Rev. E 71, 061911, 2005,
+ * doi:10.1103/PhysRevE.71.061911), as networkx's `duplication_divergence_graph`: protein-network
+ * like, with many nodes sharing neighbours.
+ *
+ * Start with the edge (0, 1). While there are fewer than n nodes, draw a node x = nextBelow(i) for
+ * the next index i, then one float per neighbour of x (in the order they were joined to x): below
+ * `retention`, the edge (i, neighbour). If i got at least one edge it is kept, otherwise the attempt
+ * is discarded. All draws come from the stream (seed, "duplication-divergence", 0). Edges are
+ * listed by new node. Simple and connected.
+ *
+ * Cost: O(sum of the copied nodes' degrees). Each attempt succeeds with probability at least
+ * `retention`; after 1000 n attempts in total it throws a RangeError rather than loop on.
+ * @param options - n, retention and seed
+ * @returns the undirected graph
+ */
+export function duplicationDivergenceGraph(options) {
+    const { n, retention } = options;
+    checkInt("n", n, 2);
+    if (!(retention > 0 && retention <= 1)) {
+        throw new RangeError(`retention must be a probability in (0, 1], got ${String(retention)}`);
+    }
+    const stream = new RandomStream(resolveSeed(options.seed), "duplication-divergence", 0);
+    // ponytail: number[][] adjacency, fine to a few million edges; a growable CSR beyond
+    const adj = Array.from({ length: n }, () => []);
+    const edges = new EdgeBuffer(4 * n);
+    edges.push(0, 1);
+    adj[0].push(1);
+    adj[1].push(0);
+    const maxAttempts = 1000 * n;
+    let attempts = 0;
+    for (let i = 2; i < n;) {
+        if (++attempts > maxAttempts) {
+            throw new RangeError(`duplicationDivergenceGraph gave up after ${maxAttempts} attempts at ${i} of ${n} nodes: retention ${retention} is too small`);
+        }
+        const copied = adj[stream.nextBelow(i)];
+        for (const t of copied) {
+            if (stream.nextFloat() < retention) {
+                edges.push(i, t);
+                adj[i].push(t);
+            }
+        }
+        for (const t of adj[i]) {
+            adj[t].push(i);
+        }
+        if (adj[i].length > 0) {
+            i++;
+        }
+    }
+    return applyWeights(toGraph(n, edges, false), options);
+}
+/**
+ * The Newman-Watts small world (M. E. J. Newman and D. J. Watts, "Renormalization group analysis
+ * of the small-world network model", Phys. Lett. A 263, 341-346, 1999,
+ * doi:10.1016/S0375-9601(99)00757-4): Watts-Strogatz with shortcuts added instead of rewired, so
+ * the graph stays connected.
+ *
+ * The ring lattice of `wattsStrogatzGraph` comes first: edge slot u k/2 + (j - 1) is
+ * (u, (u + j) mod n). Then for each slot in order, with source u: draw a float; if it is below p and
+ * u is not joined to every other node, redraw w = nextBelow(n) until w is neither u nor a
+ * neighbour of u, and append the shortcut (u, w). All draws come from the stream (seed,
+ * "newman-watts", 0). O(n k) expected.
+ * @param options - n, k, p and seed
+ * @returns the undirected graph
+ */
+export function newmanWattsGraph(options) {
+    const { n, k, p } = options;
+    checkInt("k", k, 2);
+    if (k % 2 !== 0) {
+        throw new RangeError(`k must be even, got ${k}`);
+    }
+    checkInt("n", n, k + 1);
+    checkProbability("p", p);
+    const half = k / 2;
+    const ring = n * half;
+    checkEdgeCount(2 * ring);
+    const edges = new EdgeBuffer(ring + Math.ceil(ring * p + 6 * Math.sqrt(ring) + 16));
+    const degree = new Uint32Array(n).fill(k);
+    // ponytail: a Set of u * n + v keys, like wattsStrogatzGraph; a hashed typed table beyond millions
+    const present = new Set();
+    const key = (a, b) => (a < b ? a * n + b : b * n + a);
+    for (let u = 0; u < n; u++) {
+        for (let j = 1; j <= half; j++) {
+            const v = (u + j) % n;
+            edges.push(u, v);
+            present.add(key(u, v));
+        }
+    }
+    const stream = new RandomStream(resolveSeed(options.seed), "newman-watts", 0);
+    for (let slot = 0; slot < ring; slot++) {
+        const u = edges.src[slot];
+        if (stream.nextFloat() >= p || degree[u] >= n - 1) {
+            continue;
+        }
+        let w = stream.nextBelow(n);
+        while (w === u || present.has(key(u, w))) {
+            w = stream.nextBelow(n);
+        }
+        edges.push(u, w);
+        present.add(key(u, w));
+        degree[u]++;
+        degree[w]++;
+    }
+    return applyWeights(toGraph(n, edges, false), options);
+}
+/**
+ * The Bianconi-Barabasi fitness model (G. Bianconi and A.-L. Barabasi, "Competition and
+ * multiscaling in evolving networks", Europhys. Lett. 54, 436-442, 2001,
+ * doi:10.1209/epl/i2001-00260-6): preferential attachment weighted by fitness, so a fit latecomer
+ * can overtake older hubs ("fit get richer").
+ *
+ * Start from the clique on nodes 0 .. m, edges (a, b) for a < b in lexicographic order. Each new
+ * node s = m + 1 .. n - 1 picks m distinct existing nodes, one at a time, with probability
+ * proportional to fitness_i * degree_i among those not yet picked: x = nextFloat() * total, then a
+ * descent of a Fenwick tree over the Float64 weights. Edges (s, t) follow in pick order. Draws
+ * come from the stream (seed, "bianconi-barabasi", 0). O(n m log n). Emits the Float64 node column
+ * `fitness`.
+ * @param options - n, m, fitness and seed
+ * @returns the undirected graph
+ */
+export function bianconiBarabasiGraph(options) {
+    const { n, m } = options;
+    checkInt("m", m, 1);
+    checkInt("n", n, m + 1);
+    const seed = resolveSeed(options.seed);
+    checkEdgeCount((m * (m + 1)) / 2 + (n - m - 1) * m);
+    const fitness = new Float64Array(n);
+    if (options.fitness === undefined) {
+        const draw = new RandomStream(seed, "fitness", 0);
+        for (let i = 0; i < n; i++) {
+            draw.reset(i);
+            fitness[i] = 1 - draw.nextFloat();
+        }
+    }
+    else {
+        if (options.fitness.length !== n) {
+            throw new RangeError(`fitness must have n = ${n} values, got ${options.fitness.length}`);
+        }
+        for (let i = 0; i < n; i++) {
+            const f = options.fitness[i];
+            if (!(f > 0 && Number.isFinite(f))) {
+                throw new RangeError(`fitness[${i}] must be finite and positive, got ${String(f)}`);
+            }
+            fitness[i] = f;
+        }
+    }
+    const edges = new EdgeBuffer((m * (m + 1)) / 2 + (n - m - 1) * m);
+    const degree = new Uint32Array(n);
+    for (let a = 0; a <= m; a++) {
+        for (let b = a + 1; b <= m; b++) {
+            edges.push(a, b);
+        }
+        degree[a] = m;
+    }
+    // Fenwick tree over weight[i] = fitness[i] * degree[i] (1-based internally)
+    const weight = new Float64Array(n);
+    const tree = new Float64Array(n + 1);
+    const add = (i, delta) => {
+        weight[i] += delta;
+        for (let j = i + 1; j <= n; j += j & -j) {
+            tree[j] += delta;
+        }
+    };
+    let top = 1;
+    while (top * 2 <= n) {
+        top *= 2;
+    }
+    const find = (x, limit) => {
+        let pos = 0;
+        let rest = x;
+        for (let step = top; step >= 1; step = Math.floor(step / 2)) {
+            if (pos + step <= n && tree[pos + step] <= rest) {
+                pos += step;
+                rest -= tree[pos];
+            }
+        }
+        // rounding guard: land on a node with weight, below the newcomer
+        pos = Math.min(pos, limit - 1);
+        while (weight[pos] === 0) {
+            pos--;
+        }
+        return pos;
+    };
+    let total = 0;
+    for (let i = 0; i <= m; i++) {
+        add(i, fitness[i] * m);
+        total += fitness[i] * m;
+    }
+    const stream = new RandomStream(seed, "bianconi-barabasi", 0);
+    const picked = new Uint32Array(m);
+    for (let s = m + 1; s < n; s++) {
+        let remaining = total;
+        for (let i = 0; i < m; i++) {
+            const t = find(stream.nextFloat() * remaining, s);
+            picked[i] = t;
+            remaining -= weight[t];
+            add(t, -weight[t]);
+        }
+        for (let i = 0; i < m; i++) {
+            const t = picked[i];
+            edges.push(s, t);
+            degree[t]++;
+            add(t, fitness[t] * degree[t]);
+        }
+        degree[s] = m;
+        add(s, fitness[s] * m);
+        total = 0;
+        // recompute the total from the tree's roots so it never drifts from the tree
+        for (let j = n; j > 0; j -= j & -j) {
+            total += tree[j];
+        }
+    }
+    return applyWeights(toGraph(n, edges, false, { fitness }), options);
+}
+/**
+ * The random Apollonian network (T. Zhou, G. Yan and B.-H. Wang, "Maximal planar networks with
+ * large clustering coefficient and power-law degree distribution", Phys. Rev. E 71, 046141, 2005,
+ * doi:10.1103/PhysRevE.71.046141; the deterministic one is J. S. Andrade et al., Phys. Rev. Lett. 94,
+ * 018702, 2005): maximal planar, power-law degrees, high clustering.
+ *
+ * Start from the triangle 0, 1, 2 (edges (0, 1), (0, 2), (1, 2)) whose inner face is face 0; the
+ * outer face is never split. Each new node v = 3 .. n - 1 picks face f = nextBelow(faceCount) =
+ * (a, b, c), adds the edges (v, a), (v, b), (v, c), and splits f: face f becomes (a, b, v), and
+ * (b, c, v) and (a, c, v) are appended. Draws come from the stream (seed, "apollonian", 0).
+ * 3n - 6 edges. O(n).
+ * @param options - n and seed
+ * @returns the undirected graph
+ */
+export function randomApollonianGraph(options) {
+    const { n } = options;
+    checkInt("n", n, 3);
+    const m = 3 * n - 6;
+    checkEdgeCount(m);
+    const stream = new RandomStream(resolveSeed(options.seed), "apollonian", 0);
+    const faces = 2 * n - 5;
+    const fa = new Uint32Array(faces);
+    const fb = new Uint32Array(faces);
+    const fc = new Uint32Array(faces);
+    fa[0] = 0;
+    fb[0] = 1;
+    fc[0] = 2;
+    let faceCount = 1;
+    const src = new Uint32Array(m);
+    const dst = new Uint32Array(m);
+    src.set([0, 0, 1]);
+    dst.set([1, 2, 2]);
+    let e = 3;
+    for (let v = 3; v < n; v++) {
+        const f = stream.nextBelow(faceCount);
+        const a = fa[f];
+        const b = fb[f];
+        const c = fc[f];
+        src.fill(v, e, e + 3);
+        dst[e++] = a;
+        dst[e++] = b;
+        dst[e++] = c;
+        fc[f] = v;
+        fa[faceCount] = b;
+        fb[faceCount] = c;
+        fc[faceCount++] = v;
+        fa[faceCount] = a;
+        fb[faceCount] = c;
+        fc[faceCount++] = v;
+    }
+    return applyWeights({ directed: false, nodeCount: n, src, dst }, options);
+}
+/**
+ * A uniform spanning tree of the rows x cols grid, a perfect maze, by Wilson's algorithm
+ * (D. B. Wilson, "Generating random spanning trees more quickly than the cover time", STOC 1996,
+ * doi:10.1145/237814.237880): loop-erased random walks.
+ *
+ * Node (r, c) is index r cols + c. The root is node 0. For each start node in index order that is
+ * not yet in the tree, walk until the tree is hit, each step to a uniform grid neighbour
+ * (candidates in the order left, right, up, down; nextBelow(count)), remembering only the last exit
+ * from each node (which erases the loops); then add the path: edges (u, next(u)) from the start to
+ * the tree, in that order. Draws come from the stream (seed, "wilson", 0). rows cols - 1 edges.
+ * Float64 node columns `x` = column and `y` = row. Expected cost is the mean hitting time of the
+ * root, about O(n log n) on a grid.
+ * @param options - rows, cols and seed
+ * @returns the undirected tree
+ */
+export function wilsonMazeGraph(options) {
+    const { rows, cols } = options;
+    checkInt("rows", rows, 1);
+    checkInt("cols", cols, 1);
+    const n = rows * cols;
+    checkInt("rows * cols", n, 1);
+    const stream = new RandomStream(resolveSeed(options.seed), "wilson", 0);
+    const inTree = new Uint8Array(n);
+    const next = new Uint32Array(n);
+    const src = new Uint32Array(n - 1);
+    const dst = new Uint32Array(n - 1);
+    const options4 = new Uint32Array(4);
+    inTree[0] = 1;
+    let e = 0;
+    for (let start = 1; start < n; start++) {
+        let u = start;
+        while (inTree[u] === 0) {
+            const r = Math.floor(u / cols);
+            const c = u - r * cols;
+            let count = 0;
+            if (c > 0) {
+                options4[count++] = u - 1;
+            }
+            if (c + 1 < cols) {
+                options4[count++] = u + 1;
+            }
+            if (r > 0) {
+                options4[count++] = u - cols;
+            }
+            if (r + 1 < rows) {
+                options4[count++] = u + cols;
+            }
+            next[u] = options4[stream.nextBelow(count)];
+            u = next[u];
+        }
+        for (u = start; inTree[u] === 0; u = next[u]) {
+            inTree[u] = 1;
+            src[e] = u;
+            dst[e++] = next[u];
+        }
+    }
+    const x = new Float64Array(n);
+    const y = new Float64Array(n);
+    for (let i = 0; i < n; i++) {
+        y[i] = Math.floor(i / cols);
+        x[i] = i - y[i] * cols;
+    }
+    return applyWeights({ directed: false, nodeCount: n, src, dst, nodeColumns: { x, y } }, options);
+}
+//# sourceMappingURL=growth.js.map

package/dist/src/generators/growth.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"growth.js","sourceRoot":"","sources":["../../../src/generators/growth.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,kBAAkB,CAAC;AAC1C,OAAO,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,qBAAqB,CAAC;AAEhE,OAAO,EAAE,cAAc,EAAE,QAAQ,EAAE,gBAAgB,EAAE,UAAU,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAC5F,OAAO,EAAE,YAAY,EAAsB,MAAM,cAAc,CAAC;AAUhE;;;;;;;GAOG;AACH,MAAM,UAAU,wBAAwB,CAAC,OAAmC;IACxE,MAAM,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC;IACtB,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,gBAAgB,EAAE,CAAC,CAAC,CAAC;IAChF,MAAM,GAAG,GAAG,IAAI,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnC,MAAM,GAAG,GAAG,IAAI,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACzB,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACjC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACnB,CAAC;IACD,OAAO,YAAY,CAAC,EAAE,QAAQ,EAAE,KAAK,EAAE,SAAS,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,OAAO,CAAC,CAAC;AAC9E,CAAC;AAgBD;;;;;;GAMG;AACH,SAAS,UAAU,CAAC,MAAoB,EAAE,UAAoB,EAAE,KAAa;IACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACtD,MAAM,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACxB,UAAU,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QAC9B,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACtB,CAAC;AACL,CAAC;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,MAAM,UAAU,eAAe,CAAC,OAA0B;IACtD,MAAM,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC;IAC/C,MAAM,OAAO,GAAG,OAAO,CAAC,OAAO,IAAI,IAAI,CAAC;IACxC,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QACrB,MAAM,IAAI,UAAU,CAAC,gDAAgD,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACtF,CAAC;IACD,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QACzB,MAAM,IAAI,UAAU,CAAC,0DAA0D,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAChG,CAAC;IACD,QAAQ,CAAC,SAAS,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;IAChC,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,aAAa,EAAE,CAAC,CAAC,CAAC;IAC7E,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,UAAU,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,WAAW,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC;IAC/B,oFAAoF;IACpF,MAAM,GAAG,GAAe,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC5D,MAAM,GAAG,GAAe,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC5D,MAAM,KAAK,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpC,MAAM,MAAM,GAAa,EAAE,CAAC;IAC5B,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,IAAc,EAAE,KAAa,EAAQ,EAAE;QAC5D,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;QACtD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,UAAU,CAAC,MAAM,EAAE,OAAO,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;QACzE,UAAU,CAAC,MAAM,EAAE,UAAU,EAAE,IAAI,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,KAAK,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACzB,MAAM,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/B,CAAC;IACL,CAAC,CAAC;IACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACzB,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC;QAClB,MAAM,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,MAAM,CAAC,MAAM,IAAI,MAAM,CAAC,MAAM,GAAG,OAAO,EAAE,IAAI,EAAE,EAAE,CAAC;YACzE,MAAM,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;YACvB,MAAM,IAAI,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,CAAC;YACvD,MAAM,GAAG,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,CAAC;YACxD,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;YACtB,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACzB,CAAC;QACD,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;YACrB,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACjB,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACf,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,CAAC;IACL,CAAC;IACD,OAAO,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,EAAE,OAAO,CAAC,CAAC;AAC1D,CAAC;AAYD;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,UAAU,0BAA0B,CAAC,OAAqC;IAC5E,MAAM,EAAE,CAAC,EAAE,SAAS,EAAE,GAAG,OAAO,CAAC;IACjC,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,IAAI,CAAC,CAAC,SAAS,GAAG,CAAC,IAAI,SAAS,IAAI,CAAC,CAAC,EAAE,CAAC;QACrC,MAAM,IAAI,UAAU,CAAC,kDAAkD,MAAM,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC;IAChG,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,wBAAwB,EAAE,CAAC,CAAC,CAAC;IACxF,qFAAqF;IACrF,MAAM,GAAG,GAAe,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC5D,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpC,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACjB,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACf,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACf,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,CAAC;IAC7B,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAI,CAAC;QACtB,IAAI,EAAE,QAAQ,GAAG,WAAW,EAAE,CAAC;YAC3B,MAAM,IAAI,UAAU,CAChB,4CAA4C,WAAW,gBAAgB,CAAC,OAAO,CAAC,qBAAqB,SAAS,eAAe,CAChI,CAAC;QACN,CAAC;QACD,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;YACrB,IAAI,MAAM,CAAC,SAAS,EAAE,GAAG,SAAS,EAAE,CAAC;gBACjC,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACjB,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACnB,CAAC;QACL,CAAC;QACD,KAAK,MAAM,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC;YACrB,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,CAAC;QACD,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACpB,CAAC,EAAE,CAAC;QACR,CAAC;IACL,CAAC;IACD,OAAO,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,EAAE,OAAO,CAAC,CAAC;AAC3D,CAAC;AAcD;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,gBAAgB,CAAC,OAA2B;IACxD,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC;IAC5B,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;QACd,MAAM,IAAI,UAAU,CAAC,uBAAuB,CAAC,EAAE,CAAC,CAAC;IACrD,CAAC;IACD,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACxB,gBAAgB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACzB,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC;IACtB,cAAc,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IACzB,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACpF,MAAM,MAAM,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,mGAAmG;IACnG,MAAM,OAAO,GAAG,IAAI,GAAG,EAAU,CAAC;IAClC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAC9E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACtB,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACjB,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC;IACL,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,cAAc,EAAE,CAAC,CAAC,CAAC;IAC9E,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC;QACrC,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC1B,IAAI,MAAM,CAAC,SAAS,EAAE,IAAI,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YAChD,SAAS;QACb,CAAC;QACD,IAAI,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAC5B,OAAO,CAAC,KAAK,CAAC,IAAI,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;YACvC,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAC5B,CAAC;QACD,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACjB,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;QACZ,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;IAChB,CAAC;IACD,OAAO,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,EAAE,OAAO,CAAC,CAAC;AAC3D,CAAC;AAiBD;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,qBAAqB,CAAC,OAAgC;IAClE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC;IACzB,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACxB,MAAM,IAAI,GAAG,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;IACvC,cAAc,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IACpC,IAAI,OAAO,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;QAChC,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,EAAE,SAAS,EAAE,CAAC,CAAC,CAAC;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACzB,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACd,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;QACtC,CAAC;IACL,CAAC;SAAM,CAAC;QACJ,IAAI,OAAO,CAAC,OAAO,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC/B,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,gBAAgB,OAAO,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC;QAC7F,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACzB,MAAM,CAAC,GAAG,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;YAC7B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACjC,MAAM,IAAI,UAAU,CAAC,WAAW,CAAC,sCAAsC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACxF,CAAC;YACD,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACnB,CAAC;IACL,CAAC;IACD,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAClE,MAAM,MAAM,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACrB,CAAC;QACD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAClB,CAAC;IACD,4EAA4E;IAC5E,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,KAAa,EAAQ,EAAE;QAC3C,MAAM,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACtC,IAAI,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;QACrB,CAAC;IACL,CAAC,CAAC;IACF,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,OAAO,GAAG,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;QAClB,GAAG,IAAI,CAAC,CAAC;IACb,CAAC;IACD,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,KAAa,EAAU,EAAE;QAC9C,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,KAAK,IAAI,IAAI,GAAG,GAAG,EAAE,IAAI,IAAI,CAAC,EAAE,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC;YAC1D,IAAI,GAAG,GAAG,IAAI,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,IAAI,EAAE,CAAC;gBAC9C,GAAG,IAAI,IAAI,CAAC;gBACZ,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC;YACtB,CAAC;QACL,CAAC;QACD,iEAAiE;QACjE,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC;QAC/B,OAAO,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;YACvB,GAAG,EAAE,CAAC;QACV,CAAC;QACD,OAAO,GAAG,CAAC;IACf,CAAC,CAAC;IACF,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1B,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACvB,KAAK,IAAI,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC5B,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,EAAE,mBAAmB,EAAE,CAAC,CAAC,CAAC;IAC9D,MAAM,MAAM,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,IAAI,SAAS,GAAG,KAAK,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACzB,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,GAAG,SAAS,EAAE,CAAC,CAAC,CAAC;YAClD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACd,SAAS,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC;YACvB,GAAG,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACvB,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACzB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACpB,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACjB,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;YACZ,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACnC,CAAC;QACD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACd,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACvB,KAAK,GAAG,CAAC,CAAC;QACV,6EAA6E;QAC7E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACjC,KAAK,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC;QACrB,CAAC;IACL,CAAC;IACD,OAAO,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;AACxE,CAAC;AAUD;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,qBAAqB,CAAC,OAAgC;IAClE,MAAM,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC;IACtB,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,cAAc,CAAC,CAAC,CAAC,CAAC;IAClB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,YAAY,EAAE,CAAC,CAAC,CAAC;IAC5E,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,IAAI,WAAW,CAAC,KAAK,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,IAAI,WAAW,CAAC,KAAK,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,IAAI,WAAW,CAAC,KAAK,CAAC,CAAC;IAClC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACV,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACV,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,MAAM,GAAG,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IAC/B,MAAM,GAAG,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IAC/B,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACzB,MAAM,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,SAAS,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAChB,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACtB,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;QACb,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;QACb,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;QACb,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACV,EAAE,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;QAClB,EAAE,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;QAClB,EAAE,CAAC,SAAS,EAAE,CAAC,GAAG,CAAC,CAAC;QACpB,EAAE,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;QAClB,EAAE,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;QAClB,EAAE,CAAC,SAAS,EAAE,CAAC,GAAG,CAAC,CAAC;IACxB,CAAC;IACD,OAAO,YAAY,CAAC,EAAE,QAAQ,EAAE,KAAK,EAAE,SAAS,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,OAAO,CAAC,CAAC;AAC9E,CAAC;AAYD;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,eAAe,CAAC,OAA0B;IACtD,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,OAAO,CAAC;IAC/B,QAAQ,CAAC,MAAM,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;IAC1B,QAAQ,CAAC,MAAM,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC;IACtB,QAAQ,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC9B,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;IACxE,MAAM,MAAM,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IAChC,MAAM,GAAG,GAAG,IAAI,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnC,MAAM,GAAG,GAAG,IAAI,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnC,MAAM,QAAQ,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACd,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;QACrC,IAAI,CAAC,GAAG,KAAK,CAAC;QACd,OAAO,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;YACrB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;YAC/B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;YACvB,IAAI,KAAK,GAAG,CAAC,CAAC;YACd,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACR,QAAQ,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC9B,CAAC;YACD,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,EAAE,CAAC;gBACf,QAAQ,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC9B,CAAC;YACD,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACR,QAAQ,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;YACjC,CAAC;YACD,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,EAAE,CAAC;gBACf,QAAQ,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;YACjC,CAAC;YACD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC;YAC5C,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QAChB,CAAC;QACD,KAAK,CAAC,GAAG,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;YAC3C,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACd,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACX,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACvB,CAAC;IACL,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACzB,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QAC5B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;IAC3B,CAAC;IACD,OAAO,YAAY,CAAC,EAAE,QAAQ,EAAE,KAAK,EAAE,SAAS,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,OAAO,CAAC,CAAC;AACrG,CAAC"}

package/dist/src/generators/hyperbolic.d.ts

@@ -0,0 +1,84 @@
+/**
+ * Random hyperbolic graphs (D. Krioukov, F. Papadopoulos, M. Kitsak, A. Vahdat and M. Boguna,
+ * "Hyperbolic geometry of complex networks", Phys. Rev. E 82, 036106, 2010,
+ * doi:10.1103/PhysRevE.82.036106): power-law degrees, high clustering and a native 2D embedding.
+ */
+import { type SampleGraph } from "../types.js";
+import { EdgeBuffer } from "./util.js";
+import { type WeightOptions } from "./weights.js";
+/** Options of {@link hyperbolicGraph}. */
+export interface HyperbolicOptions extends WeightOptions {
+    /** The node count, in [0, 20000]. */
+    n: number;
+    /** The target mean degree, a finite number > 0 (approximate: see {@link hyperbolicGraph}). */
+    averageDegree: number;
+    /** The power-law exponent gamma of the degree distribution, a finite number > 2. */
+    exponent: number;
+    /** The temperature T in [0, 1); 0 (default) is the threshold model, higher T lowers clustering. */
+    temperature?: number | undefined;
+    /** The seed, an integer in [0, 2^53); default 0. */
+    seed?: number | undefined;
+}
+/** The checked options, the disk radius and the points. */
+interface HyperbolicPlan {
+    readonly n: number;
+    readonly seed: number;
+    readonly temperature: number;
+    /** The disk radius R. */
+    readonly R: number;
+    readonly coshR: number;
+    readonly cosT: Float64Array<ArrayBuffer>;
+    readonly sinT: Float64Array<ArrayBuffer>;
+    readonly coshr: Float64Array<ArrayBuffer>;
+    readonly sinhr: Float64Array<ArrayBuffer>;
+    readonly x: Float64Array<ArrayBuffer>;
+    readonly y: Float64Array<ArrayBuffer>;
+    readonly radius: Float64Array<ArrayBuffer>;
+}
+/**
+ * Check the options, compute the disk radius and draw the points.
+ * @param options - the options
+ * @returns the plan
+ */
+export declare function planHyperbolic(options: HyperbolicOptions): HyperbolicPlan;
+/**
+ * Rows [start, end) of the hyperbolic graph: row u holds the pairs (u, v), v = u + 1 .. n - 1 in
+ * ascending v, with cosh d = cosh r_u cosh r_v - sinh r_u sinh r_v cos(theta_u - theta_v) (the
+ * cosine of the angle difference expanded from the stored cosines and sines). At T = 0 the pair is
+ * an edge iff cosh d < cosh R (no draws). At T > 0 row u draws one nextFloat() t per pair from
+ * stream (seed, "hyperbolic", u) and the pair is an edge iff t < 1 / (1 + exp((d - R) / (2 T))),
+ * decided as acosh(cosh d) < R + 2 T ln((1 - t) / t) with no exp. Any split of [0, n) into
+ * consecutive ranges gives the whole graph.
+ * @param plan - the checked options and points
+ * @param start - the first row
+ * @param end - one past the last row
+ * @param out - receives the edges
+ */
+export declare function hyperbolicRows(plan: HyperbolicPlan, start: number, end: number, out: EdgeBuffer): void;
+/**
+ * A random hyperbolic graph (Krioukov et al. 2010): n points in a hyperbolic disk of radius R,
+ * angles uniform and radii with density proportional to sinh(alpha r), alpha = (gamma - 1) / 2,
+ * so the degrees follow a power law with exponent gamma. At temperature 0 two nodes are joined iff
+ * their hyperbolic distance is below R; at T > 0 with probability 1 / (1 + exp((d - R) / (2 T))).
+ *
+ * R comes from the paper's approximation of the mean degree, R = 2 ln(2 xi^2 n / (pi k)) with
+ * xi = alpha / (alpha - 1/2), times (pi T / sin(pi T)) inside the logarithm for T > 0. It is an
+ * approximation, so the measured mean degree only lands near `averageDegree`: measured at n = 3000
+ * and 10000, within about 10% for exponent >= 2.5 and T <= 0.7, but about 0.6x at exponent 2.2 and
+ * about 0.5x at T = 0.9 (the approximation degrades as gamma approaches 2 and T approaches 1).
+ *
+ * Node i draws its angle theta = 2 pi u and then its radius r = acosh(1 + (cosh(alpha R) - 1) u') /
+ * alpha (inverse CDF) from stream (seed, "hyperbolic-points", i); the pairs are decided by
+ * {@link hyperbolicRows}. Edges are the pairs (u, v), u < v, ascending u then v. Node columns (all
+ * Float64Array): `x`, `y` (the native polar embedding r cos theta, r sin theta) and `radius` (r).
+ *
+ * Cost: O(n^2) pairs (and one draw per pair at T > 0), so n is capped at 20000.
+ * ponytail: naive pair loop, n <= 20000; the band algorithm of M. von Looz, H. Meyerhenke and
+ * R. Prutkin, "Generating random hyperbolic graphs in subquadratic time", ISAAC 2015,
+ * doi:10.1007/978-3-662-48971-0_40, is the upgrade path.
+ * @param options - n, averageDegree, exponent, temperature and seed
+ * @returns the undirected graph
+ */
+export declare function hyperbolicGraph(options: HyperbolicOptions): SampleGraph;
+export {};
+//# sourceMappingURL=hyperbolic.d.ts.map

package/dist/src/generators/hyperbolic.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"hyperbolic.d.ts","sourceRoot":"","sources":["../../../src/generators/hyperbolic.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAMH,OAAO,EAAE,KAAK,WAAW,EAAE,MAAM,aAAa,CAAC;AAC/C,OAAO,EAAY,UAAU,EAAW,MAAM,WAAW,CAAC;AAC1D,OAAO,EAAgB,KAAK,aAAa,EAAE,MAAM,cAAc,CAAC;AAKhE,0CAA0C;AAC1C,MAAM,WAAW,iBAAkB,SAAQ,aAAa;IACpD,qCAAqC;IACrC,CAAC,EAAE,MAAM,CAAC;IACV,8FAA8F;IAC9F,aAAa,EAAE,MAAM,CAAC;IACtB,oFAAoF;IACpF,QAAQ,EAAE,MAAM,CAAC;IACjB,mGAAmG;IACnG,WAAW,CAAC,EAAE,MAAM,GAAG,SAAS,CAAC;IACjC,oDAAoD;IACpD,IAAI,CAAC,EAAE,MAAM,GAAG,SAAS,CAAC;CAC7B;AAED,2DAA2D;AAC3D,UAAU,cAAc;IACpB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,yBAAyB;IACzB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,IAAI,EAAE,YAAY,CAAC,WAAW,CAAC,CAAC;IACzC,QAAQ,CAAC,IAAI,EAAE,YAAY,CAAC,WAAW,CAAC,CAAC;IACzC,QAAQ,CAAC,KAAK,EAAE,YAAY,CAAC,WAAW,CAAC,CAAC;IAC1C,QAAQ,CAAC,KAAK,EAAE,YAAY,CAAC,WAAW,CAAC,CAAC;IAC1C,QAAQ,CAAC,CAAC,EAAE,YAAY,CAAC,WAAW,CAAC,CAAC;IACtC,QAAQ,CAAC,CAAC,EAAE,YAAY,CAAC,WAAW,CAAC,CAAC;IACtC,QAAQ,CAAC,MAAM,EAAE,YAAY,CAAC,WAAW,CAAC,CAAC;CAC9C;AAoBD;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,OAAO,EAAE,iBAAiB,GAAG,cAAc,CAwDzE;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,cAAc,CAAC,IAAI,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,UAAU,GAAG,IAAI,CAmBtG;AAED;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,wBAAgB,eAAe,CAAC,OAAO,EAAE,iBAAiB,GAAG,WAAW,CAKvE"}