react-cheminfo 0.9.0 → 0.11.0

package/lib/scatter/ui/useScatterSelection.js

@@ -0,0 +1,72 @@
+import { useCallback, useMemo, useState } from 'react';
+import { mergeScatterSelection, scatterSelectedIndices, scatterSelectionMask, } from "../core/scatterSelection.js";
+/**
+ * The set of points a scatter has picked out, controlled or not.
+ *
+ * A caller may own the selection and pass it back through `selected`, or leave
+ * it here and only listen. Both work through one hook, and the controlled
+ * value wins on every render: a parent that answers a change by passing
+ * something else — a filter that keeps only the rows still on screen — sees
+ * what it asked for and not what the gesture asked for. The hook keeps its own
+ * copy up to date all the same, so a caller that stops passing `selected`
+ * carries on from where the reader left off rather than from nothing.
+ *
+ * A preview is deliberately not a change. What a half-drawn lasso would catch
+ * is merged in for the render and never reported, which is what lets a
+ * controlled parent stay out of a drag it would otherwise be re-rendered by
+ * sixty times a second.
+ * @param options - See {@link ScatterSelectionOptions}.
+ * @returns The selection. See {@link ScatterSelectionApi}.
+ */
+export function useScatterSelection(options) {
+    const { count, selected, defaultSelected, onSelectionChange } = options;
+    const [own, setOwn] = useState(defaultSelected ?? EMPTY);
+    const [draft, setDraft] = useState(null);
+    const committed = selected ?? own;
+    const mask = useMemo(() => {
+        const base = scatterSelectionMask(committed, count);
+        if (draft === null)
+            return base;
+        return mergeScatterSelection(base, draft.hits, draft.mode, base);
+    }, [committed, count, draft]);
+    const indices = useMemo(() => scatterSelectedIndices(mask), [mask]);
+    const commit = useCallback((rows, mode, source) => {
+        setDraft(null);
+        setOwn(rows);
+        onSelectionChange?.({ indices: rows, mode, source });
+    }, [onSelectionChange]);
+    const selectMask = useCallback((hits, mode = 'replace', source = 'lasso') => {
+        const base = scatterSelectionMask(committed, count);
+        const merged = mergeScatterSelection(base, hits, mode, base);
+        commit(scatterSelectedIndices(merged), mode, source);
+    }, [commit, committed, count]);
+    const select = useCallback((rows, mode = 'replace', source = 'point') => {
+        selectMask(scatterSelectionMask(rows, count), mode, source);
+    }, [count, selectMask]);
+    const preview = useCallback((hits, mode = 'replace') => {
+        setDraft(hits === null ? null : { hits, mode });
+    }, []);
+    const clear = useCallback(() => {
+        commit([], 'replace', 'clear');
+    }, [commit]);
+    const selectAll = useCallback(() => {
+        const rows = new Array(Math.max(0, count));
+        for (let index = 0; index < rows.length; index++)
+            rows[index] = index;
+        commit(rows, 'replace', 'all');
+    }, [commit, count]);
+    return {
+        mask,
+        indices,
+        selectedCount: indices.length,
+        previewing: draft !== null,
+        select,
+        selectMask,
+        preview,
+        clear,
+        selectAll,
+    };
+}
+/** Shared, so an uncontrolled scatter with nothing picked allocates nothing. */
+const EMPTY = [];
+//# sourceMappingURL=useScatterSelection.js.map

package/lib/scatter/ui/useScatterSelection.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"useScatterSelection.js","sourceRoot":"","sources":["../../../src/scatter/ui/useScatterSelection.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,WAAW,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,OAAO,CAAC;AAGvD,OAAO,EACL,qBAAqB,EACrB,sBAAsB,EACtB,oBAAoB,GACrB,MAAM,6BAA6B,CAAC;AAsErC;;;;;;;;;;;;;;;;;GAiBG;AACH,MAAM,UAAU,mBAAmB,CACjC,OAAgC;IAEhC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,eAAe,EAAE,iBAAiB,EAAE,GAAG,OAAO,CAAC;IAExE,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,QAAQ,CAAoB,eAAe,IAAI,KAAK,CAAC,CAAC;IAC5E,MAAM,CAAC,KAAK,EAAE,QAAQ,CAAC,GAAG,QAAQ,CAAwB,IAAI,CAAC,CAAC;IAChE,MAAM,SAAS,GAAG,QAAQ,IAAI,GAAG,CAAC;IAElC,MAAM,IAAI,GAAG,OAAO,CAAC,GAAG,EAAE;QACxB,MAAM,IAAI,GAAG,oBAAoB,CAAC,SAAS,EAAE,KAAK,CAAC,CAAC;QACpD,IAAI,KAAK,KAAK,IAAI;YAAE,OAAO,IAAI,CAAC;QAChC,OAAO,qBAAqB,CAAC,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IACnE,CAAC,EAAE,CAAC,SAAS,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC,CAAC;IAE9B,MAAM,OAAO,GAAG,OAAO,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC;IAEpE,MAAM,MAAM,GAAG,WAAW,CACxB,CACE,IAAc,EACd,IAA0B,EAC1B,MAAiC,EACjC,EAAE;QACF,QAAQ,CAAC,IAAI,CAAC,CAAC;QACf,MAAM,CAAC,IAAI,CAAC,CAAC;QACb,iBAAiB,EAAE,CAAC,EAAE,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,CAAC;IACvD,CAAC,EACD,CAAC,iBAAiB,CAAC,CACpB,CAAC;IAEF,MAAM,UAAU,GAAG,WAAW,CAC5B,CACE,IAAgB,EAChB,OAA6B,SAAS,EACtC,SAAoC,OAAO,EAC3C,EAAE;QACF,MAAM,IAAI,GAAG,oBAAoB,CAAC,SAAS,EAAE,KAAK,CAAC,CAAC;QACpD,MAAM,MAAM,GAAG,qBAAqB,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC7D,MAAM,CAAC,sBAAsB,CAAC,MAAM,CAAC,EAAE,IAAI,EAAE,MAAM,CAAC,CAAC;IACvD,CAAC,EACD,CAAC,MAAM,EAAE,SAAS,EAAE,KAAK,CAAC,CAC3B,CAAC;IAEF,MAAM,MAAM,GAAG,WAAW,CACxB,CACE,IAAsB,EACtB,OAA6B,SAAS,EACtC,SAAoC,OAAO,EAC3C,EAAE;QACF,UAAU,CAAC,oBAAoB,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,IAAI,EAAE,MAAM,CAAC,CAAC;IAC9D,CAAC,EACD,CAAC,KAAK,EAAE,UAAU,CAAC,CACpB,CAAC;IAEF,MAAM,OAAO,GAAG,WAAW,CACzB,CAAC,IAAuB,EAAE,OAA6B,SAAS,EAAE,EAAE;QAClE,QAAQ,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,CAAC;IAClD,CAAC,EACD,EAAE,CACH,CAAC;IAEF,MAAM,KAAK,GAAG,WAAW,CAAC,GAAG,EAAE;QAC7B,MAAM,CAAC,EAAE,EAAE,SAAS,EAAE,OAAO,CAAC,CAAC;IACjC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC;IAEb,MAAM,SAAS,GAAG,WAAW,CAAC,GAAG,EAAE;QACjC,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;QACnD,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE;YAAE,IAAI,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC;QACtE,MAAM,CAAC,IAAI,EAAE,SAAS,EAAE,KAAK,CAAC,CAAC;IACjC,CAAC,EAAE,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC;IAEpB,OAAO;QACL,IAAI;QACJ,OAAO;QACP,aAAa,EAAE,OAAO,CAAC,MAAM;QAC7B,UAAU,EAAE,KAAK,KAAK,IAAI;QAC1B,MAAM;QACN,UAAU;QACV,OAAO;QACP,KAAK;QACL,SAAS;KACV,CAAC;AACJ,CAAC;AAED,gFAAgF;AAChF,MAAM,KAAK,GAAsB,EAAE,CAAC"}

package/lib/scatter3d/core/cloudGesture.d.ts

@@ -0,0 +1,17 @@
+/**
+ * What a drag over a cloud does.
+ *
+ * Two gestures want the same drag. Turning the box is the one a reader reaches
+ * for first and must stay a plain drag; drawing a lasso is the one the flat map
+ * has already taught them. Neither can be given a modifier without taking one
+ * that already means "add to the selection" or "take out of it", so the choice
+ * is a control the reader can see rather than a shortcut somebody has to tell
+ * them about — and inside a lasso the modifiers keep exactly the meaning they
+ * have on the map.
+ */
+export type CloudGesture = 
+/** Turn the box. A tap still picks the sample under it. */
+'turn'
+/** Draw a lasso, exactly as on the map. */
+ | 'select';
+//# sourceMappingURL=cloudGesture.d.ts.map

package/lib/scatter3d/core/cloudGesture.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cloudGesture.d.ts","sourceRoot":"","sources":["../../../src/scatter3d/core/cloudGesture.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AACH,MAAM,MAAM,YAAY;AACtB,2DAA2D;AACzD,MAAM;AACR,2CAA2C;GACzC,QAAQ,CAAC"}

package/lib/scatter3d/core/cloudGesture.js

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=cloudGesture.js.map

package/lib/scatter3d/core/cloudGesture.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cloudGesture.js","sourceRoot":"","sources":["../../../src/scatter3d/core/cloudGesture.ts"],"names":[],"mappings":""}

package/lib/scatter3d/core/confidenceEllipsoid.d.ts

@@ -0,0 +1,63 @@
+/**
+ * The shell that summarises where a group of samples sits in the cloud.
+ *
+ * It is the outline of `scatter/core` with one dimension more, and it is
+ * described the same way — a centre and the axes of the group's own spread,
+ * in the data's units — so that turning it into something an SVG can draw
+ * stays a separate step. Here that step matters more than it does on a map:
+ * a shell has no `<ellipse>` to fall back on and has to be tessellated,
+ * projected and painted back to front.
+ */
+import type { EllipseSize } from '../../scatter/core/confidenceEllipse.ts';
+import type { Vector3 } from './orbitCamera.ts';
+/** A group's shell, in the cloud's own units. */
+export interface ConfidenceEllipsoid {
+    /** Centre of the shell, which is the group's average. */
+    center: Vector3;
+    /**
+     * Three orthogonal semi-axes, longest first, each vector as long as its own
+     * radius. A direction the group has no spread in is the zero vector, which
+     * is what flattens the shell rather than giving it a false thickness.
+     */
+    axes: readonly [Vector3, Vector3, Vector3];
+    /** How many finite points went into it. */
+    count: number;
+}
+/** How {@link confidenceEllipsoid} measures a group. */
+export interface ConfidenceEllipsoidOptions {
+    /**
+     * How large to draw the shell. A share is read in three dimensions, so the
+     * same share gives a wider shell here than the same call gives on the map.
+     * @default { kind: 'coverage', probability: 0.95 }
+     */
+    size?: EllipseSize;
+    /**
+     * How many finite points a group needs before it is given a shell at all.
+     * Below four a group has no volume to measure — three points lie on a plane
+     * however they are arranged — so fewer than four is read as four.
+     * @default 4
+     */
+    minimumPoints?: number;
+}
+/**
+ * The shell that summarises a group of samples.
+ *
+ * Points that are not finite are skipped rather than poisoning the average,
+ * and the spread is measured with the `n - 1` divisor, so the shell estimates
+ * the group rather than describing the handful of samples that happen to have
+ * been measured — both exactly as the two-dimensional outline does.
+ * @param points - The group's samples, in the cloud's own units.
+ * @param options - See {@link ConfidenceEllipsoidOptions}.
+ * @returns The shell, or `null` when the group is too small to have a shape, when its coordinates are so large that squaring them overflows, or when the requested size is unbounded.
+ */
+export declare function confidenceEllipsoid(points: readonly Vector3[], options?: ConfidenceEllipsoidOptions): ConfidenceEllipsoid | null;
+/**
+ * How far out a size asks the shell to be drawn.
+ *
+ * Public because a legend has to be able to say what a share came to without
+ * building a shell to ask it.
+ * @param size - How large the shell was asked to be.
+ * @returns The number of standard deviations.
+ */
+export declare function ellipsoidStandardDeviations(size: EllipseSize): number;
+//# sourceMappingURL=confidenceEllipsoid.d.ts.map

package/lib/scatter3d/core/confidenceEllipsoid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"confidenceEllipsoid.d.ts","sourceRoot":"","sources":["../../../src/scatter3d/core/confidenceEllipsoid.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,yCAAyC,CAAC;AAG3E,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAIhD,iDAAiD;AACjD,MAAM,WAAW,mBAAmB;IAClC,yDAAyD;IACzD,MAAM,EAAE,OAAO,CAAC;IAChB;;;;OAIG;IACH,IAAI,EAAE,SAAS,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;IAC3C,2CAA2C;IAC3C,KAAK,EAAE,MAAM,CAAC;CACf;AAED,wDAAwD;AACxD,MAAM,WAAW,0BAA0B;IACzC;;;;OAIG;IACH,IAAI,CAAC,EAAE,WAAW,CAAC;IACnB;;;;;OAKG;IACH,aAAa,CAAC,EAAE,MAAM,CAAC;CACxB;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,mBAAmB,CACjC,MAAM,EAAE,SAAS,OAAO,EAAE,EAC1B,OAAO,GAAE,0BAA+B,GACvC,mBAAmB,GAAG,IAAI,CA0C5B;AAED;;;;;;;GAOG;AACH,wBAAgB,2BAA2B,CAAC,IAAI,EAAE,WAAW,GAAG,MAAM,CAGrE"}

package/lib/scatter3d/core/confidenceEllipsoid.js

@@ -0,0 +1,128 @@
+/**
+ * The shell that summarises where a group of samples sits in the cloud.
+ *
+ * It is the outline of `scatter/core` with one dimension more, and it is
+ * described the same way — a centre and the axes of the group's own spread,
+ * in the data's units — so that turning it into something an SVG can draw
+ * stays a separate step. Here that step matters more than it does on a map:
+ * a shell has no `<ellipse>` to fall back on and has to be tessellated,
+ * projected and painted back to front.
+ */
+import { standardDeviationsForCoverage3 } from "./ellipsoidCoverage.js";
+import { symmetricEigen3 } from "./symmetricEigen.js";
+/**
+ * The shell that summarises a group of samples.
+ *
+ * Points that are not finite are skipped rather than poisoning the average,
+ * and the spread is measured with the `n - 1` divisor, so the shell estimates
+ * the group rather than describing the handful of samples that happen to have
+ * been measured — both exactly as the two-dimensional outline does.
+ * @param points - The group's samples, in the cloud's own units.
+ * @param options - See {@link ConfidenceEllipsoidOptions}.
+ * @returns The shell, or `null` when the group is too small to have a shape, when its coordinates are so large that squaring them overflows, or when the requested size is unbounded.
+ */
+export function confidenceEllipsoid(points, options = {}) {
+    const { size = DEFAULT_SIZE, minimumPoints = 4 } = options;
+    const floor = Math.max(4, minimumPoints);
+    let count = 0;
+    let sumX = 0;
+    let sumY = 0;
+    let sumZ = 0;
+    for (const point of points) {
+        if (!isFinitePoint(point))
+            continue;
+        count++;
+        sumX += point[0];
+        sumY += point[1];
+        sumZ += point[2];
+    }
+    if (count < floor)
+        return null;
+    const center = [sumX / count, sumY / count, sumZ / count];
+    const covariance = spreadAbout(points, center, count);
+    if (covariance === null)
+        return null;
+    const radius = ellipsoidStandardDeviations(size);
+    if (!Number.isFinite(radius))
+        return null;
+    const { values, vectors } = symmetricEigen3(covariance);
+    const axes = [];
+    for (let i = 0; i < 3; i++) {
+        const length = radius * Math.sqrt(Math.max(0, values[i]));
+        if (!Number.isFinite(length))
+            return null;
+        const direction = vectors[i];
+        axes.push([
+            direction[0] * length,
+            direction[1] * length,
+            direction[2] * length,
+        ]);
+    }
+    return {
+        center,
+        axes: [axes[0], axes[1], axes[2]],
+        count,
+    };
+}
+/**
+ * How far out a size asks the shell to be drawn.
+ *
+ * Public because a legend has to be able to say what a share came to without
+ * building a shell to ask it.
+ * @param size - How large the shell was asked to be.
+ * @returns The number of standard deviations.
+ */
+export function ellipsoidStandardDeviations(size) {
+    if (size.kind === 'standardDeviations')
+        return size.standardDeviations;
+    return standardDeviationsForCoverage3(size.probability);
+}
+function spreadAbout(points, center, count) {
+    let xx = 0;
+    let xy = 0;
+    let xz = 0;
+    let yy = 0;
+    let yz = 0;
+    let zz = 0;
+    for (const point of points) {
+        if (!isFinitePoint(point))
+            continue;
+        const dx = point[0] - center[0];
+        const dy = point[1] - center[1];
+        const dz = point[2] - center[2];
+        xx += dx * dx;
+        xy += dx * dy;
+        xz += dx * dz;
+        yy += dy * dy;
+        yz += dy * dz;
+        zz += dz * dz;
+    }
+    const divisor = count - 1;
+    const covariance = {
+        xx: xx / divisor,
+        xy: xy / divisor,
+        xz: xz / divisor,
+        yy: yy / divisor,
+        yz: yz / divisor,
+        zz: zz / divisor,
+    };
+    if (!Number.isFinite(center[0]) ||
+        !Number.isFinite(center[1]) ||
+        !Number.isFinite(center[2]) ||
+        !Number.isFinite(covariance.xx) ||
+        !Number.isFinite(covariance.yy) ||
+        !Number.isFinite(covariance.zz) ||
+        !Number.isFinite(covariance.xy) ||
+        !Number.isFinite(covariance.xz) ||
+        !Number.isFinite(covariance.yz)) {
+        return null;
+    }
+    return covariance;
+}
+function isFinitePoint(point) {
+    return (Number.isFinite(point[0]) &&
+        Number.isFinite(point[1]) &&
+        Number.isFinite(point[2]));
+}
+const DEFAULT_SIZE = { kind: 'coverage', probability: 0.95 };
+//# sourceMappingURL=confidenceEllipsoid.js.map

package/lib/scatter3d/core/confidenceEllipsoid.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"confidenceEllipsoid.js","sourceRoot":"","sources":["../../../src/scatter3d/core/confidenceEllipsoid.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAIH,OAAO,EAAE,8BAA8B,EAAE,MAAM,wBAAwB,CAAC;AAGxE,OAAO,EAAE,eAAe,EAAE,MAAM,qBAAqB,CAAC;AAiCtD;;;;;;;;;;GAUG;AACH,MAAM,UAAU,mBAAmB,CACjC,MAA0B,EAC1B,UAAsC,EAAE;IAExC,MAAM,EAAE,IAAI,GAAG,YAAY,EAAE,aAAa,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAC3D,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,aAAa,CAAC,CAAC;IAEzC,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC;YAAE,SAAS;QACpC,KAAK,EAAE,CAAC;QACR,IAAI,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;QACjB,IAAI,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;QACjB,IAAI,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IACnB,CAAC;IACD,IAAI,KAAK,GAAG,KAAK;QAAE,OAAO,IAAI,CAAC;IAE/B,MAAM,MAAM,GAAY,CAAC,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,CAAC,CAAC;IACnE,MAAM,UAAU,GAAG,WAAW,CAAC,MAAM,EAAE,MAAM,EAAE,KAAK,CAAC,CAAC;IACtD,IAAI,UAAU,KAAK,IAAI;QAAE,OAAO,IAAI,CAAC;IAErC,MAAM,MAAM,GAAG,2BAA2B,CAAC,IAAI,CAAC,CAAC;IACjD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;QAAE,OAAO,IAAI,CAAC;IAE1C,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,eAAe,CAAC,UAAU,CAAC,CAAC;IACxD,MAAM,IAAI,GAAc,EAAE,CAAC;IAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAW,CAAC,CAAC,CAAC;QACpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YAAE,OAAO,IAAI,CAAC;QAC1C,MAAM,SAAS,GAAG,OAAO,CAAC,CAAC,CAAY,CAAC;QACxC,IAAI,CAAC,IAAI,CAAC;YACR,SAAS,CAAC,CAAC,CAAC,GAAG,MAAM;YACrB,SAAS,CAAC,CAAC,CAAC,GAAG,MAAM;YACrB,SAAS,CAAC,CAAC,CAAC,GAAG,MAAM;SACtB,CAAC,CAAC;IACL,CAAC;IAED,OAAO;QACL,MAAM;QACN,IAAI,EAAE,CAAC,IAAI,CAAC,CAAC,CAAY,EAAE,IAAI,CAAC,CAAC,CAAY,EAAE,IAAI,CAAC,CAAC,CAAY,CAAC;QAClE,KAAK;KACN,CAAC;AACJ,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,2BAA2B,CAAC,IAAiB;IAC3D,IAAI,IAAI,CAAC,IAAI,KAAK,oBAAoB;QAAE,OAAO,IAAI,CAAC,kBAAkB,CAAC;IACvE,OAAO,8BAA8B,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;AAC1D,CAAC;AAED,SAAS,WAAW,CAClB,MAA0B,EAC1B,MAAe,EACf,KAAa;IAEb,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC;YAAE,SAAS;QACpC,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAChC,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAChC,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAChC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QACd,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QACd,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QACd,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QACd,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QACd,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;IAChB,CAAC;IACD,MAAM,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IAC1B,MAAM,UAAU,GAAqB;QACnC,EAAE,EAAE,EAAE,GAAG,OAAO;QAChB,EAAE,EAAE,EAAE,GAAG,OAAO;QAChB,EAAE,EAAE,EAAE,GAAG,OAAO;QAChB,EAAE,EAAE,EAAE,GAAG,OAAO;QAChB,EAAE,EAAE,EAAE,GAAG,OAAO;QAChB,EAAE,EAAE,EAAE,GAAG,OAAO;KACjB,CAAC;IACF,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,EAAE,CAAC;QAC/B,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,EAAE,CAAC;QAC/B,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,EAAE,CAAC;QAC/B,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,EAAE,CAAC;QAC/B,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,EAAE,CAAC;QAC/B,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC,EAAE,CAAC,EAC/B,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,UAAU,CAAC;AACpB,CAAC;AAED,SAAS,aAAa,CAAC,KAAc;IACnC,OAAO,CACL,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAC1B,CAAC;AACJ,CAAC;AAED,MAAM,YAAY,GAAgB,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,EAAE,CAAC"}

package/lib/scatter3d/core/cubeFrame.d.ts

@@ -0,0 +1,40 @@
+/**
+ * The box the cloud sits in.
+ *
+ * A cloud of dots on an empty background has no orientation at all: turn it
+ * and nothing tells the reader it turned, and a dot near the top could be near
+ * or far. Twelve faint edges and three labelled arms fix that for the cost of
+ * fifteen lines, and they are the reason the projection can stay orthographic
+ * — the frame carries the depth that perspective would otherwise have to, and
+ * unlike perspective it distorts no distance.
+ */
+import type { Vector3 } from './orbitCamera.ts';
+/** The corner the three labelled arms are drawn from. */
+export declare const CUBE_ORIGIN: Vector3;
+/**
+ * The twelve edges of the cube, each as its two corners.
+ *
+ * Every edge, not just the three at the back: which edges are hidden changes
+ * as the box turns, and drawing all twelve faintly reads as a wireframe, while
+ * drawing three reads as an axis system that keeps rearranging itself.
+ * @returns The edges. The order is stable but carries no meaning.
+ */
+export declare function cubeEdges(): Array<[Vector3, Vector3]>;
+/** One labelled arm of the frame. */
+export interface CubeArm {
+    /** Which of the three axes it stands for. */
+    axis: 'x' | 'y' | 'z';
+    /** The far end of the arm. */
+    end: Vector3;
+    /** Where its name is written, just beyond the end. */
+    label: Vector3;
+}
+/**
+ * The three arms, from the origin corner along each axis.
+ *
+ * The label sits a little past the end of its arm so that it is outside the
+ * cloud whatever the box is turned to; at the end itself it would fall on the
+ * frame from about a third of the angles.
+ */
+export declare const CUBE_ARMS: readonly CubeArm[];
+//# sourceMappingURL=cubeFrame.d.ts.map

package/lib/scatter3d/core/cubeFrame.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cubeFrame.d.ts","sourceRoot":"","sources":["../../../src/scatter3d/core/cubeFrame.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAEhD,yDAAyD;AACzD,eAAO,MAAM,WAAW,EAAE,OAAsB,CAAC;AAEjD;;;;;;;GAOG;AACH,wBAAgB,SAAS,IAAI,KAAK,CAAC,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAwBrD;AAED,qCAAqC;AACrC,MAAM,WAAW,OAAO;IACtB,6CAA6C;IAC7C,IAAI,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IACtB,8BAA8B;IAC9B,GAAG,EAAE,OAAO,CAAC;IACb,sDAAsD;IACtD,KAAK,EAAE,OAAO,CAAC;CAChB;AAED;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,EAAE,SAAS,OAAO,EAIvC,CAAC"}

package/lib/scatter3d/core/cubeFrame.js

@@ -0,0 +1,58 @@
+/**
+ * The box the cloud sits in.
+ *
+ * A cloud of dots on an empty background has no orientation at all: turn it
+ * and nothing tells the reader it turned, and a dot near the top could be near
+ * or far. Twelve faint edges and three labelled arms fix that for the cost of
+ * fifteen lines, and they are the reason the projection can stay orthographic
+ * — the frame carries the depth that perspective would otherwise have to, and
+ * unlike perspective it distorts no distance.
+ */
+/** The corner the three labelled arms are drawn from. */
+export const CUBE_ORIGIN = [-1, -1, -1];
+/**
+ * The twelve edges of the cube, each as its two corners.
+ *
+ * Every edge, not just the three at the back: which edges are hidden changes
+ * as the box turns, and drawing all twelve faintly reads as a wireframe, while
+ * drawing three reads as an axis system that keeps rearranging itself.
+ * @returns The edges. The order is stable but carries no meaning.
+ */
+export function cubeEdges() {
+    const corners = [];
+    for (const x of [-1, 1]) {
+        for (const y of [-1, 1]) {
+            for (const z of [-1, 1])
+                corners.push([x, y, z]);
+        }
+    }
+    const edges = [];
+    for (let i = 0; i < corners.length; i++) {
+        for (let j = i + 1; j < corners.length; j++) {
+            const a = corners[i];
+            const b = corners[j];
+            // Two corners are joined by an edge exactly when they differ in one
+            // coordinate; differing in two is a face diagonal and in three is the
+            // body diagonal, and drawing either would turn the frame into a mesh.
+            const differing = (a[0] === b[0] ? 0 : 1) +
+                (a[1] === b[1] ? 0 : 1) +
+                (a[2] === b[2] ? 0 : 1);
+            if (differing === 1)
+                edges.push([a, b]);
+        }
+    }
+    return edges;
+}
+/**
+ * The three arms, from the origin corner along each axis.
+ *
+ * The label sits a little past the end of its arm so that it is outside the
+ * cloud whatever the box is turned to; at the end itself it would fall on the
+ * frame from about a third of the angles.
+ */
+export const CUBE_ARMS = [
+    { axis: 'x', end: [1, -1, -1], label: [1.18, -1, -1] },
+    { axis: 'y', end: [-1, 1, -1], label: [-1, 1.18, -1] },
+    { axis: 'z', end: [-1, -1, 1], label: [-1, -1, 1.18] },
+];
+//# sourceMappingURL=cubeFrame.js.map

package/lib/scatter3d/core/cubeFrame.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cubeFrame.js","sourceRoot":"","sources":["../../../src/scatter3d/core/cubeFrame.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAIH,yDAAyD;AACzD,MAAM,CAAC,MAAM,WAAW,GAAY,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAEjD;;;;;;;GAOG;AACH,MAAM,UAAU,SAAS;IACvB,MAAM,OAAO,GAAc,EAAE,CAAC;IAC9B,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;QACxB,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;YACxB,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACnD,CAAC;IACH,CAAC;IAED,MAAM,KAAK,GAA8B,EAAE,CAAC;IAC5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAY,CAAC;YAChC,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAY,CAAC;YAChC,oEAAoE;YACpE,sEAAsE;YACtE,sEAAsE;YACtE,MAAM,SAAS,GACb,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACvB,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACvB,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1B,IAAI,SAAS,KAAK,CAAC;gBAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC1C,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAYD;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,SAAS,GAAuB;IAC3C,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE;IACtD,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,EAAE;IACtD,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE;CACvD,CAAC"}

package/lib/scatter3d/core/cubePoints.d.ts

@@ -0,0 +1,35 @@
+/**
+ * Three runs of numbers, stretched into the cube the cloud is drawn in.
+ *
+ * A map draws its two axes at whatever scale each one needs, and the tick
+ * labels are the reader's record of that. A cloud has no such record: it is
+ * one solid seen from an angle, so three axes drawn at three scales would make
+ * a group lean in a direction that is an artefact of the scaling and not
+ * anything in the data. Every axis is therefore stretched to the same cube,
+ * the frame is labelled with which axis is which, and the numbers stay in the
+ * card the pointer raises.
+ *
+ * The shells are fitted in this cube for the same reason the map fits its
+ * outlines in pixels: a shell has to sit on the cloud the reader is looking
+ * at, not on a cloud in units nothing draws.
+ */
+import type { Vector3 } from './orbitCamera.ts';
+/** How {@link cubePoints} stretches the runs it is given. */
+export interface CubePointsOptions {
+    /**
+     * Share of each axis' span left empty at both ends, so the cloud does not
+     * touch the frame.
+     * @default 0.05
+     */
+    padding?: number;
+}
+/**
+ * Three runs of numbers, each stretched to `[-1, 1]`.
+ * @param x - Where every sample sits along the axis that runs left to right.
+ * @param y - Along the axis that runs bottom to top, in the same order.
+ * @param z - Along the axis that runs away from the reader.
+ * @param options - See {@link CubePointsOptions}.
+ * @returns One point per sample, as many as the shortest run holds. An axis with no spread at all puts every sample at the middle of it rather than dividing by zero, and a value that is not finite stays that way so nothing draws it.
+ */
+export declare function cubePoints(x: ArrayLike<number>, y: ArrayLike<number>, z: ArrayLike<number>, options?: CubePointsOptions): Vector3[];
+//# sourceMappingURL=cubePoints.d.ts.map

package/lib/scatter3d/core/cubePoints.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cubePoints.d.ts","sourceRoot":"","sources":["../../../src/scatter3d/core/cubePoints.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAOH,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,kBAAkB,CAAC;AAEhD,6DAA6D;AAC7D,MAAM,WAAW,iBAAiB;IAChC;;;;OAIG;IACH,OAAO,CAAC,EAAE,MAAM,CAAC;CAClB;AAED;;;;;;;GAOG;AACH,wBAAgB,UAAU,CACxB,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EACpB,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EACpB,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EACpB,OAAO,GAAE,iBAAsB,GAC9B,OAAO,EAAE,CAiBX"}

package/lib/scatter3d/core/cubePoints.js

@@ -0,0 +1,50 @@
+/**
+ * Three runs of numbers, stretched into the cube the cloud is drawn in.
+ *
+ * A map draws its two axes at whatever scale each one needs, and the tick
+ * labels are the reader's record of that. A cloud has no such record: it is
+ * one solid seen from an angle, so three axes drawn at three scales would make
+ * a group lean in a direction that is an artefact of the scaling and not
+ * anything in the data. Every axis is therefore stretched to the same cube,
+ * the frame is labelled with which axis is which, and the numbers stay in the
+ * card the pointer raises.
+ *
+ * The shells are fitted in this cube for the same reason the map fits its
+ * outlines in pixels: a shell has to sit on the cloud the reader is looking
+ * at, not on a cloud in units nothing draws.
+ */
+import { chartPadExtent, chartValuesExtent, } from "../../chart/core/chartExtent.js";
+/**
+ * Three runs of numbers, each stretched to `[-1, 1]`.
+ * @param x - Where every sample sits along the axis that runs left to right.
+ * @param y - Along the axis that runs bottom to top, in the same order.
+ * @param z - Along the axis that runs away from the reader.
+ * @param options - See {@link CubePointsOptions}.
+ * @returns One point per sample, as many as the shortest run holds. An axis with no spread at all puts every sample at the middle of it rather than dividing by zero, and a value that is not finite stays that way so nothing draws it.
+ */
+export function cubePoints(x, y, z, options = {}) {
+    const { padding = 0.05 } = options;
+    const toCube = [x, y, z].map((values) => axisMapper(values, padding));
+    const across = toCube[0];
+    const up = toCube[1];
+    const away = toCube[2];
+    const count = Math.min(x.length, y.length, z.length);
+    const points = new Array(count);
+    for (let index = 0; index < count; index++) {
+        points[index] = [
+            across(x[index]),
+            up(y[index]),
+            away(z[index]),
+        ];
+    }
+    return points;
+}
+function axisMapper(values, padding) {
+    const { min, max } = chartPadExtent(chartValuesExtent(values), padding);
+    const span = max - min;
+    if (!Number.isFinite(span) || span <= 0) {
+        return (value) => (Number.isFinite(value) ? 0 : Number.NaN);
+    }
+    return (value) => (2 * (value - min)) / span - 1;
+}
+//# sourceMappingURL=cubePoints.js.map

package/lib/scatter3d/core/cubePoints.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cubePoints.js","sourceRoot":"","sources":["../../../src/scatter3d/core/cubePoints.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EACL,cAAc,EACd,iBAAiB,GAClB,MAAM,iCAAiC,CAAC;AAczC;;;;;;;GAOG;AACH,MAAM,UAAU,UAAU,CACxB,CAAoB,EACpB,CAAoB,EACpB,CAAoB,EACpB,UAA6B,EAAE;IAE/B,MAAM,EAAE,OAAO,GAAG,IAAI,EAAE,GAAG,OAAO,CAAC;IACnC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,UAAU,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC,CAAC;IACtE,MAAM,MAAM,GAAG,MAAM,CAAC,CAAC,CAAW,CAAC;IACnC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAW,CAAC;IAC/B,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAW,CAAC;IAEjC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;IACrD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAU,KAAK,CAAC,CAAC;IACzC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QAC3C,MAAM,CAAC,KAAK,CAAC,GAAG;YACd,MAAM,CAAC,CAAC,CAAC,KAAK,CAAW,CAAC;YAC1B,EAAE,CAAC,CAAC,CAAC,KAAK,CAAW,CAAC;YACtB,IAAI,CAAC,CAAC,CAAC,KAAK,CAAW,CAAC;SACzB,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAID,SAAS,UAAU,CAAC,MAAyB,EAAE,OAAe;IAC5D,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,cAAc,CAAC,iBAAiB,CAAC,MAAM,CAAC,EAAE,OAAO,CAAC,CAAC;IACxE,MAAM,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC;IACvB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE,CAAC;QACxC,OAAO,CAAC,KAAa,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACtE,CAAC;IACD,OAAO,CAAC,KAAa,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;AAC3D,CAAC"}

package/lib/scatter3d/core/ellipsoidCoverage.d.ts

@@ -0,0 +1,32 @@
+/**
+ * How far out a group's shell is drawn, and how much of the group it then
+ * holds.
+ *
+ * The two-dimensional pair in `scatter/core` cannot answer this, and the
+ * difference is not a rounding: a shell has to hold a sample in three
+ * directions at once, so the same two standard deviations that keep 86 per
+ * cent of a flat cloud keep only 74 per cent of a solid one. Quoting the flat
+ * figure beside a cloud is the usual way a legend comes to promise more than
+ * its shells contain, which is why the map and the cloud each ask their own
+ * dimension.
+ */
+/**
+ * The share of a three-dimensional normal cloud inside the shell drawn at a
+ * number of standard deviations.
+ * @param standardDeviations - How far out the shell is drawn.
+ * @returns The share, from 0 to 1. Zero at or below zero, since a shell of no size holds nothing.
+ */
+export declare function coverageForStandardDeviations3(standardDeviations: number): number;
+/**
+ * How far out to draw the shell that holds a share of a three-dimensional
+ * normal cloud.
+ *
+ * The inverse of {@link coverageForStandardDeviations3}. There is no closed
+ * form for it, so it is bisected — the share grows strictly with the radius,
+ * so fifty halvings settle it far below the width of the line the shell is
+ * drawn with, and the search runs once per group rather than once per face.
+ * @param coverage - The share to enclose, from 0 to 1.
+ * @returns The number of standard deviations. Zero at or below zero, and infinite at or above one, because no finite shell holds every point of a cloud with unbounded tails.
+ */
+export declare function standardDeviationsForCoverage3(coverage: number): number;
+//# sourceMappingURL=ellipsoidCoverage.d.ts.map

package/lib/scatter3d/core/ellipsoidCoverage.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoidCoverage.d.ts","sourceRoot":"","sources":["../../../src/scatter3d/core/ellipsoidCoverage.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH;;;;;GAKG;AACH,wBAAgB,8BAA8B,CAC5C,kBAAkB,EAAE,MAAM,GACzB,MAAM,CAUR;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,8BAA8B,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM,CAiBvE"}

package/lib/scatter3d/core/ellipsoidCoverage.js

@@ -0,0 +1,83 @@
+/**
+ * How far out a group's shell is drawn, and how much of the group it then
+ * holds.
+ *
+ * The two-dimensional pair in `scatter/core` cannot answer this, and the
+ * difference is not a rounding: a shell has to hold a sample in three
+ * directions at once, so the same two standard deviations that keep 86 per
+ * cent of a flat cloud keep only 74 per cent of a solid one. Quoting the flat
+ * figure beside a cloud is the usual way a legend comes to promise more than
+ * its shells contain, which is why the map and the cloud each ask their own
+ * dimension.
+ */
+/**
+ * The share of a three-dimensional normal cloud inside the shell drawn at a
+ * number of standard deviations.
+ * @param standardDeviations - How far out the shell is drawn.
+ * @returns The share, from 0 to 1. Zero at or below zero, since a shell of no size holds nothing.
+ */
+export function coverageForStandardDeviations3(standardDeviations) {
+    if (!Number.isFinite(standardDeviations)) {
+        return standardDeviations > 0 ? 1 : 0;
+    }
+    if (standardDeviations <= 0)
+        return 0;
+    const radius = standardDeviations;
+    const inside = erf(radius / Math.SQRT2);
+    const shell = radius * Math.sqrt(2 / Math.PI) * Math.exp(-(radius * radius) / 2);
+    return Math.min(1, Math.max(0, inside - shell));
+}
+/**
+ * How far out to draw the shell that holds a share of a three-dimensional
+ * normal cloud.
+ *
+ * The inverse of {@link coverageForStandardDeviations3}. There is no closed
+ * form for it, so it is bisected — the share grows strictly with the radius,
+ * so fifty halvings settle it far below the width of the line the shell is
+ * drawn with, and the search runs once per group rather than once per face.
+ * @param coverage - The share to enclose, from 0 to 1.
+ * @returns The number of standard deviations. Zero at or below zero, and infinite at or above one, because no finite shell holds every point of a cloud with unbounded tails.
+ */
+export function standardDeviationsForCoverage3(coverage) {
+    if (!Number.isFinite(coverage) || coverage <= 0)
+        return 0;
+    if (coverage >= 1)
+        return Number.POSITIVE_INFINITY;
+    let low = 0;
+    // Nine standard deviations hold every share a control can ask for: the
+    // remainder there is under 1e-17, which no `number` can tell from zero.
+    let high = 9;
+    for (let step = 0; step < 50; step++) {
+        const middle = (low + high) / 2;
+        if (coverageForStandardDeviations3(middle) < coverage) {
+            low = middle;
+        }
+        else {
+            high = middle;
+        }
+    }
+    return (low + high) / 2;
+}
+/**
+ * The error function, to about a part in ten million.
+ *
+ * Abramowitz and Stegun 7.1.26. That is coarser than the library's usual
+ * standard and is deliberately so: it decides the radius of a shell drawn a
+ * few hundred pixels across, where a part in ten million is a part in ten
+ * thousand of one pixel, and the alternative is a series that costs more to
+ * read than the whole module around it.
+ * @param x - The argument.
+ * @returns The error function at `x`.
+ */
+function erf(x) {
+    const sign = x < 0 ? -1 : 1;
+    const absolute = Math.abs(x);
+    const t = 1 / (1 + 0.327_591_1 * absolute);
+    const series = t *
+        (0.254_829_592 +
+            t *
+                (-0.284_496_736 +
+                    t * (1.421_413_741 + t * (-1.453_152_027 + t * 1.061_405_429))));
+    return sign * (1 - series * Math.exp(-absolute * absolute));
+}
+//# sourceMappingURL=ellipsoidCoverage.js.map

package/lib/scatter3d/core/ellipsoidCoverage.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoidCoverage.js","sourceRoot":"","sources":["../../../src/scatter3d/core/ellipsoidCoverage.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH;;;;;GAKG;AACH,MAAM,UAAU,8BAA8B,CAC5C,kBAA0B;IAE1B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,kBAAkB,CAAC,EAAE,CAAC;QACzC,OAAO,kBAAkB,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxC,CAAC;IACD,IAAI,kBAAkB,IAAI,CAAC;QAAE,OAAO,CAAC,CAAC;IACtC,MAAM,MAAM,GAAG,kBAAkB,CAAC;IAClC,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;IACxC,MAAM,KAAK,GACT,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;IACrE,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC;AAClD,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,8BAA8B,CAAC,QAAgB;IAC7D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC;QAAE,OAAO,CAAC,CAAC;IAC1D,IAAI,QAAQ,IAAI,CAAC;QAAE,OAAO,MAAM,CAAC,iBAAiB,CAAC;IAEnD,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,uEAAuE;IACvE,wEAAwE;IACxE,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,EAAE,EAAE,IAAI,EAAE,EAAE,CAAC;QACrC,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAChC,IAAI,8BAA8B,CAAC,MAAM,CAAC,GAAG,QAAQ,EAAE,CAAC;YACtD,GAAG,GAAG,MAAM,CAAC;QACf,CAAC;aAAM,CAAC;YACN,IAAI,GAAG,MAAM,CAAC;QAChB,CAAC;IACH,CAAC;IACD,OAAO,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;AAC1B,CAAC;AAED;;;;;;;;;;GAUG;AACH,SAAS,GAAG,CAAC,CAAS;IACpB,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5B,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,WAAW,GAAG,QAAQ,CAAC,CAAC;IAC3C,MAAM,MAAM,GACV,CAAC;QACD,CAAC,aAAa;YACZ,CAAC;gBACC,CAAC,CAAC,aAAa;oBACb,CAAC,GAAG,CAAC,aAAa,GAAG,CAAC,GAAG,CAAC,CAAC,aAAa,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;IACzE,OAAO,IAAI,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC;AAC9D,CAAC"}