react-cheminfo 0.10.0 → 0.12.0

package/lib/router/ui/useTabRoute.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"useTabRoute.js","sourceRoot":"","sources":["../../../src/router/ui/useTabRoute.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,WAAW,EAAE,SAAS,EAAE,OAAO,EAAE,oBAAoB,EAAE,MAAM,OAAO,CAAC;AAE9E,OAAO,EACL,cAAc,EACd,SAAS,EACT,gBAAgB,GACjB,MAAM,oBAAoB,CAAC;AAG5B;;;;;;;;GAQG;AACH,MAAM,UAAU,WAAW,CACzB,MAAsB;IAEtB,MAAM,SAAS,GAAG,WAAW,CAC3B,CAAC,QAAoB,EAAE,EAAE,CAAC,gBAAgB,CAAC,MAAM,EAAE,QAAQ,CAAC,EAC5D,CAAC,MAAM,CAAC,CACT,CAAC;IACF,MAAM,OAAO,GAAG,oBAAoB,CAClC,SAAS,EACT,cAAc,EACd,cAAc,CACf,CAAC;IACF,SAAS,CAAC,GAAG,EAAE;QACb,SAAS,CAAC,MAAM,CAAC,CAAC;IACpB,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC;IACb,OAAO,OAAO,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,EAAE,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC,CAAC;AACjE,CAAC"}

package/lib/scatter/core/confidenceEllipse.d.ts

@@ -5,6 +5,8 @@
  * can draw is `projectEllipse`, which has to be a separate step because the two
  * axes rarely carry the same number of pixels per unit.
  */
+import type { EllipseCovariance } from './ellipseAxes.ts';
+import type { ScatterGroupMoments } from './scatterGroupMoments.ts';
 /** A point of a cloud, or of an outline drawn around one, in data units. */
 export interface EllipsePoint {
     /** Position along the horizontal axis. */
@@ -35,15 +37,8 @@ export interface EllipseStandardDeviationSize {
  * one into the other by hand.
  */
 export type EllipseSize = EllipseCoverageSize | EllipseStandardDeviationSize;
-/** How a cloud is spread, as the three distinct entries of its covariance. */
-export interface EllipseCovariance {
-    /** The spread along the horizontal axis. */
-    xx: number;
-    /** How the two axes move together; zero when they are unrelated. */
-    xy: number;
-    /** The spread along the vertical axis. */
-    yy: number;
-}
+/** What a group is outlined at when the caller says nothing: the share holding 95 % of it. */
+export declare const DEFAULT_ELLIPSE_SIZE: EllipseSize;
 /** A group outline, in data units. */
 export interface ConfidenceEllipse {
     /** Centre of the outline, which is the group's average. */
@@ -89,19 +84,30 @@ export interface ConfidenceEllipseOptions {
  * @returns The outline, or `null` when the group is too small to have a shape, or when its coordinates are so large that squaring them overflows, or when the requested size is unbounded.
  */
 export declare function confidenceEllipse(points: readonly EllipsePoint[], options?: ConfidenceEllipseOptions): ConfidenceEllipse | null;
+/** Which outline {@link momentsEllipse} cuts out of a block of moments. */
+interface MomentsEllipseOptions {
+    /** Which group. */
+    group: number;
+    /** The axis the outline's `x` is read along. */
+    xAxis: number;
+    /** The axis its `y` is read along. */
+    yAxis: number;
+    /** How far out to draw, in standard deviations. */
+    standardDeviations: number;
+    /** How many rows the group needs before it has a shape; below two is read as two. */
+    minimumPoints: number;
+}
 /**
- * The semi-axes and the rotation that a covariance describes.
+ * One group's outline for one pair of axes, cut out of measured moments.
  *
- * The shape of an outline is the eigen-decomposition of a symmetric two by two
- * matrix, which has a closed form and so never needs an iterative solver. It
- * is public because the same three numbers are what turns any covariance into
- * something drawable — one measured elsewhere, or one already mapped into
- * pixels.
- * @param covariance - How the cloud is spread.
- * @param standardDeviations - How far out to measure. Defaults to `1`.
- * @returns The two semi-axis lengths and the rotation of the major one, in radians, in `(-π/2, π/2]`. A collinear cloud gives `ry` of exactly `0`, and one with no spread at all gives an angle of `0` because it has no direction to report.
+ * The two-by-two block of a covariance is a covariance in its own right, so a
+ * pair grid measuring every axis at once and a map measuring one pair give
+ * the same outline for the same samples.
+ * @param moments - The block, from `scatterGroupMoments`.
+ * @param options - See {@link MomentsEllipseOptions}.
+ * @returns The outline, or `null` when the group is too small to have a shape, when either axis or the group is outside the block, or when the size or the coordinates are not finite.
  */
-export declare function ellipseAxes(covariance: EllipseCovariance, standardDeviations?: number): Pick<ConfidenceEllipse, 'rx' | 'ry' | 'angle'>;
+export declare function momentsEllipse(moments: ScatterGroupMoments, options: MomentsEllipseOptions): ConfidenceEllipse | null;
 /**
  * How far out an outline of a given size is drawn.
  *
@@ -114,4 +120,5 @@ export declare function ellipseAxes(covariance: EllipseCovariance, standardDevia
  * @returns The number of standard deviations, never below zero.
  */
 export declare function ellipseStandardDeviations(size: EllipseSize): number;
+export {};
 //# sourceMappingURL=confidenceEllipse.d.ts.map

package/lib/scatter/core/confidenceEllipse.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"confidenceEllipse.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/confidenceEllipse.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAgBH,4EAA4E;AAC5E,MAAM,WAAW,YAAY;IAC3B,0CAA0C;IAC1C,CAAC,EAAE,MAAM,CAAC;IACV,wCAAwC;IACxC,CAAC,EAAE,MAAM,CAAC;CACX;AAED,qDAAqD;AACrD,MAAM,WAAW,mBAAmB;IAClC,gEAAgE;IAChE,IAAI,EAAE,UAAU,CAAC;IACjB,sDAAsD;IACtD,WAAW,EAAE,MAAM,CAAC;CACrB;AAED,4DAA4D;AAC5D,MAAM,WAAW,4BAA4B;IAC3C,gEAAgE;IAChE,IAAI,EAAE,oBAAoB,CAAC;IAC3B,mDAAmD;IACnD,kBAAkB,EAAE,MAAM,CAAC;CAC5B;AAED;;;;;;;GAOG;AACH,MAAM,MAAM,WAAW,GAAG,mBAAmB,GAAG,4BAA4B,CAAC;AAE7E,8EAA8E;AAC9E,MAAM,WAAW,iBAAiB;IAChC,4CAA4C;IAC5C,EAAE,EAAE,MAAM,CAAC;IACX,oEAAoE;IACpE,EAAE,EAAE,MAAM,CAAC;IACX,0CAA0C;IAC1C,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,sCAAsC;AACtC,MAAM,WAAW,iBAAiB;IAChC,2DAA2D;IAC3D,EAAE,EAAE,MAAM,CAAC;IACX,2DAA2D;IAC3D,EAAE,EAAE,MAAM,CAAC;IACX,8BAA8B;IAC9B,EAAE,EAAE,MAAM,CAAC;IACX,uEAAuE;IACvE,EAAE,EAAE,MAAM,CAAC;IACX,kFAAkF;IAClF,KAAK,EAAE,MAAM,CAAC;IACd,4EAA4E;IAC5E,KAAK,EAAE,MAAM,CAAC;IACd,0EAA0E;IAC1E,UAAU,EAAE,iBAAiB,CAAC;CAC/B;AAED,sDAAsD;AACtD,MAAM,WAAW,wBAAwB;IACvC;;;OAGG;IACH,IAAI,CAAC,EAAE,WAAW,CAAC;IACnB;;;;;OAKG;IACH,aAAa,CAAC,EAAE,MAAM,CAAC;CACxB;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,SAAS,YAAY,EAAE,EAC/B,OAAO,GAAE,wBAA6B,GACrC,iBAAiB,GAAG,IAAI,CAmD1B;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,WAAW,CACzB,UAAU,EAAE,iBAAiB,EAC7B,kBAAkB,SAAI,GACrB,IAAI,CAAC,iBAAiB,EAAE,IAAI,GAAG,IAAI,GAAG,OAAO,CAAC,CAahD;AAiBD;;;;;;;;;;GAUG;AACH,wBAAgB,yBAAyB,CAAC,IAAI,EAAE,WAAW,GAAG,MAAM,CAKnE"}
+{"version":3,"file":"confidenceEllipse.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/confidenceEllipse.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAIH,OAAO,KAAK,EAAE,iBAAiB,EAAE,MAAM,kBAAkB,CAAC;AAG1D,OAAO,KAAK,EAAE,mBAAmB,EAAE,MAAM,0BAA0B,CAAC;AAGpE,4EAA4E;AAC5E,MAAM,WAAW,YAAY;IAC3B,0CAA0C;IAC1C,CAAC,EAAE,MAAM,CAAC;IACV,wCAAwC;IACxC,CAAC,EAAE,MAAM,CAAC;CACX;AAED,qDAAqD;AACrD,MAAM,WAAW,mBAAmB;IAClC,gEAAgE;IAChE,IAAI,EAAE,UAAU,CAAC;IACjB,sDAAsD;IACtD,WAAW,EAAE,MAAM,CAAC;CACrB;AAED,4DAA4D;AAC5D,MAAM,WAAW,4BAA4B;IAC3C,gEAAgE;IAChE,IAAI,EAAE,oBAAoB,CAAC;IAC3B,mDAAmD;IACnD,kBAAkB,EAAE,MAAM,CAAC;CAC5B;AAED;;;;;;;GAOG;AACH,MAAM,MAAM,WAAW,GAAG,mBAAmB,GAAG,4BAA4B,CAAC;AAE7E,8FAA8F;AAC9F,eAAO,MAAM,oBAAoB,EAAE,WAGlC,CAAC;AAEF,sCAAsC;AACtC,MAAM,WAAW,iBAAiB;IAChC,2DAA2D;IAC3D,EAAE,EAAE,MAAM,CAAC;IACX,2DAA2D;IAC3D,EAAE,EAAE,MAAM,CAAC;IACX,8BAA8B;IAC9B,EAAE,EAAE,MAAM,CAAC;IACX,uEAAuE;IACvE,EAAE,EAAE,MAAM,CAAC;IACX,kFAAkF;IAClF,KAAK,EAAE,MAAM,CAAC;IACd,4EAA4E;IAC5E,KAAK,EAAE,MAAM,CAAC;IACd,0EAA0E;IAC1E,UAAU,EAAE,iBAAiB,CAAC;CAC/B;AAED,sDAAsD;AACtD,MAAM,WAAW,wBAAwB;IACvC;;;OAGG;IACH,IAAI,CAAC,EAAE,WAAW,CAAC;IACnB;;;;;OAKG;IACH,aAAa,CAAC,EAAE,MAAM,CAAC;CACxB;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,SAAS,YAAY,EAAE,EAC/B,OAAO,GAAE,wBAA6B,GACrC,iBAAiB,GAAG,IAAI,CAU1B;AAED,2EAA2E;AAC3E,UAAU,qBAAqB;IAC7B,mBAAmB;IACnB,KAAK,EAAE,MAAM,CAAC;IACd,gDAAgD;IAChD,KAAK,EAAE,MAAM,CAAC;IACd,sCAAsC;IACtC,KAAK,EAAE,MAAM,CAAC;IACd,mDAAmD;IACnD,kBAAkB,EAAE,MAAM,CAAC;IAC3B,qFAAqF;IACrF,aAAa,EAAE,MAAM,CAAC;CACvB;AAED;;;;;;;;;GASG;AACH,wBAAgB,cAAc,CAC5B,OAAO,EAAE,mBAAmB,EAC5B,OAAO,EAAE,qBAAqB,GAC7B,iBAAiB,GAAG,IAAI,CAqB1B;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,yBAAyB,CAAC,IAAI,EAAE,WAAW,GAAG,MAAM,CAKnE"}

package/lib/scatter/core/confidenceEllipse.js

@@ -5,18 +5,14 @@
  * can draw is `projectEllipse`, which has to be a separate step because the two
  * axes rarely carry the same number of pixels per unit.
  */
+import { ellipseAxes } from "./ellipseAxes.js";
 import { standardDeviationsForCoverage } from "./ellipseCoverage.js";
-/**
- * A smaller eigenvalue below this share of the larger one is read as zero.
- *
- * A group whose spread across the major axis is under a millionth of its
- * spread along it is collinear for any purpose a plot has, and the rounding of
- * the covariance alone reaches that far: the nine points on a straight line
- * used in the tests leave a residue of about 5e-17 where the exact answer is
- * 0. Without the floor those points would draw a hairline ellipse on one
- * machine and a segment on the next.
- */
-const COLLINEAR_FLOOR = 1e-12;
+import { scatterGroupMoments } from "./scatterGroupMoments.js";
+/** What a group is outlined at when the caller says nothing: the share holding 95 % of it. */
+export const DEFAULT_ELLIPSE_SIZE = {
+    kind: 'coverage',
+    probability: 0.95,
+};
 /**
  * The outline that summarises a group of samples.
  *
@@ -30,88 +26,52 @@ const COLLINEAR_FLOOR = 1e-12;
  * @returns The outline, or `null` when the group is too small to have a shape, or when its coordinates are so large that squaring them overflows, or when the requested size is unbounded.
  */
 export function confidenceEllipse(points, options = {}) {
-    const { size = { kind: 'coverage', probability: 0.95 }, minimumPoints = 3 } = options;
-    let count = 0;
-    let sumX = 0;
-    let sumY = 0;
-    for (const point of points) {
-        if (!Number.isFinite(point.x) || !Number.isFinite(point.y))
-            continue;
-        count++;
-        sumX += point.x;
-        sumY += point.y;
-    }
+    const { size = DEFAULT_ELLIPSE_SIZE, minimumPoints = 3 } = options;
+    const moments = scatterGroupMoments(pointMatrix(points), null, 1, 2);
+    return momentsEllipse(moments, {
+        group: 0,
+        xAxis: 0,
+        yAxis: 1,
+        standardDeviations: ellipseStandardDeviations(size),
+        minimumPoints,
+    });
+}
+/**
+ * One group's outline for one pair of axes, cut out of measured moments.
+ *
+ * The two-by-two block of a covariance is a covariance in its own right, so a
+ * pair grid measuring every axis at once and a map measuring one pair give
+ * the same outline for the same samples.
+ * @param moments - The block, from `scatterGroupMoments`.
+ * @param options - See {@link MomentsEllipseOptions}.
+ * @returns The outline, or `null` when the group is too small to have a shape, when either axis or the group is outside the block, or when the size or the coordinates are not finite.
+ */
+export function momentsEllipse(moments, options) {
+    const { axes, groups, counts, means, covariances } = moments;
+    const { group, xAxis, yAxis, standardDeviations, minimumPoints } = options;
+    if (group < 0 || group >= groups)
+        return null;
+    if (xAxis < 0 || yAxis < 0 || xAxis >= axes || yAxis >= axes)
+        return null;
+    const count = counts[group] ?? 0;
     if (count < Math.max(2, minimumPoints))
         return null;
-    const cx = sumX / count;
-    const cy = sumY / count;
-    let spreadX = 0;
-    let together = 0;
-    let spreadY = 0;
-    for (const point of points) {
-        if (!Number.isFinite(point.x) || !Number.isFinite(point.y))
-            continue;
-        const dx = point.x - cx;
-        const dy = point.y - cy;
-        spreadX += dx * dx;
-        together += dx * dy;
-        spreadY += dy * dy;
-    }
-    const divisor = count - 1;
+    const base = group * axes;
+    const block = group * axes * axes;
     const covariance = {
-        xx: spreadX / divisor,
-        xy: together / divisor,
-        yy: spreadY / divisor,
+        xx: covariances[block + xAxis * axes + xAxis] ?? 0,
+        xy: covariances[block + xAxis * axes + yAxis] ?? 0,
+        yy: covariances[block + yAxis * axes + yAxis] ?? 0,
     };
-    if (!Number.isFinite(cx) ||
-        !Number.isFinite(cy) ||
-        !Number.isFinite(covariance.xx) ||
-        !Number.isFinite(covariance.xy) ||
-        !Number.isFinite(covariance.yy)) {
-        return null;
-    }
-    const { rx, ry, angle } = ellipseAxes(covariance, ellipseStandardDeviations(size));
+    const { rx, ry, angle } = ellipseAxes(covariance, standardDeviations);
+    const cx = means[base + xAxis] ?? 0;
+    const cy = means[base + yAxis] ?? 0;
     if (!Number.isFinite(rx) || !Number.isFinite(ry))
         return null;
+    if (!Number.isFinite(cx) || !Number.isFinite(cy))
+        return null;
     return { cx, cy, rx, ry, angle, count, covariance };
 }
-/**
- * The semi-axes and the rotation that a covariance describes.
- *
- * The shape of an outline is the eigen-decomposition of a symmetric two by two
- * matrix, which has a closed form and so never needs an iterative solver. It
- * is public because the same three numbers are what turns any covariance into
- * something drawable — one measured elsewhere, or one already mapped into
- * pixels.
- * @param covariance - How the cloud is spread.
- * @param standardDeviations - How far out to measure. Defaults to `1`.
- * @returns The two semi-axis lengths and the rotation of the major one, in radians, in `(-π/2, π/2]`. A collinear cloud gives `ry` of exactly `0`, and one with no spread at all gives an angle of `0` because it has no direction to report.
- */
-export function ellipseAxes(covariance, standardDeviations = 1) {
-    const { xx, xy, yy } = covariance;
-    const middle = (xx + yy) / 2;
-    const half = (xx - yy) / 2;
-    const radius = Math.hypot(half, xy);
-    const major = middle + radius;
-    let minor = middle - radius;
-    if (!(minor > major * COLLINEAR_FLOOR))
-        minor = 0;
-    return {
-        rx: standardDeviations * Math.sqrt(major),
-        ry: standardDeviations * Math.sqrt(minor),
-        angle: majorAxisAngle(xx, xy, yy, major),
-    };
-}
-function majorAxisAngle(xx, xy, yy, major) {
-    // With no cross term the axes already are the eigenvectors, and asking
-    // `atan2` would only turn the exact answer into a rounded one.
-    if (xy === 0)
-        return xx >= yy ? 0 : Math.PI / 2;
-    // `major - xx` is never negative, so this lands in [0, π] and a negative
-    // cross term arrives half a turn away from the direction actually wanted.
-    const angle = Math.atan2(major - xx, xy);
-    return angle > Math.PI / 2 ? angle - Math.PI : angle;
-}
 /**
  * How far out an outline of a given size is drawn.
  *
@@ -129,4 +89,22 @@ export function ellipseStandardDeviations(size) {
     }
     return standardDeviationsForCoverage(size.probability);
 }
+/*
+ * A run of points read as a two-column matrix, so a group handed as objects
+ * is measured by the same scan as one handed as a score matrix.
+ */
+function pointMatrix(points) {
+    return {
+        rows: points.length,
+        columns: 2,
+        get(rowIndex, columnIndex) {
+            const point = points[rowIndex];
+            if (point === undefined)
+                return Number.NaN;
+            if (columnIndex === 0)
+                return point.x;
+            return columnIndex === 1 ? point.y : Number.NaN;
+        },
+    };
+}
 //# sourceMappingURL=confidenceEllipse.js.map

package/lib/scatter/core/confidenceEllipse.js.map

@@ -1 +1 @@
-{"version":3,"file":"confidenceEllipse.js","sourceRoot":"","sources":["../../../src/scatter/core/confidenceEllipse.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,EAAE,6BAA6B,EAAE,MAAM,sBAAsB,CAAC;AAErE;;;;;;;;;GASG;AACH,MAAM,eAAe,GAAG,KAAK,CAAC;AAgF9B;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAA+B,EAC/B,UAAoC,EAAE;IAEtC,MAAM,EAAE,IAAI,GAAG,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,EAAE,EAAE,aAAa,GAAG,CAAC,EAAE,GACzE,OAAO,CAAC;IAEV,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC;YAAE,SAAS;QACrE,KAAK,EAAE,CAAC;QACR,IAAI,IAAI,KAAK,CAAC,CAAC,CAAC;QAChB,IAAI,IAAI,KAAK,CAAC,CAAC,CAAC;IAClB,CAAC;IACD,IAAI,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,aAAa,CAAC;QAAE,OAAO,IAAI,CAAC;IAEpD,MAAM,EAAE,GAAG,IAAI,GAAG,KAAK,CAAC;IACxB,MAAM,EAAE,GAAG,IAAI,GAAG,KAAK,CAAC;IACxB,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC;YAAE,SAAS;QACrE,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;QACxB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;QACxB,OAAO,IAAI,EAAE,GAAG,EAAE,CAAC;QACnB,QAAQ,IAAI,EAAE,GAAG,EAAE,CAAC;QACpB,OAAO,IAAI,EAAE,GAAG,EAAE,CAAC;IACrB,CAAC;IAED,MAAM,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IAC1B,MAAM,UAAU,GAAsB;QACpC,EAAE,EAAE,OAAO,GAAG,OAAO;QACrB,EAAE,EAAE,QAAQ,GAAG,OAAO;QACtB,EAAE,EAAE,OAAO,GAAG,OAAO;KACtB,CAAC;IACF,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACpB,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACpB,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;IAED,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,GAAG,WAAW,CACnC,UAAU,EACV,yBAAyB,CAAC,IAAI,CAAC,CAChC,CAAC;IACF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QAAE,OAAO,IAAI,CAAC;IAC9D,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,EAAE,CAAC;AACtD,CAAC;AAED;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,WAAW,CACzB,UAA6B,EAC7B,kBAAkB,GAAG,CAAC;IAEtB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,UAAU,CAAC;IAClC,MAAM,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;IAC3B,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,EAAE,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,MAAM,GAAG,MAAM,CAAC;IAC9B,IAAI,KAAK,GAAG,MAAM,GAAG,MAAM,CAAC;IAC5B,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,GAAG,eAAe,CAAC;QAAE,KAAK,GAAG,CAAC,CAAC;IAClD,OAAO;QACL,EAAE,EAAE,kBAAkB,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC;QACzC,EAAE,EAAE,kBAAkB,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC;QACzC,KAAK,EAAE,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC;KACzC,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CACrB,EAAU,EACV,EAAU,EACV,EAAU,EACV,KAAa;IAEb,uEAAuE;IACvE,+DAA+D;IAC/D,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IAChD,yEAAyE;IACzE,0EAA0E;IAC1E,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC;IACzC,OAAO,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;AACvD,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,yBAAyB,CAAC,IAAiB;IACzD,IAAI,IAAI,CAAC,IAAI,KAAK,oBAAoB,EAAE,CAAC;QACvC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,kBAAkB,EAAE,CAAC,CAAC,CAAC;IAC9C,CAAC;IACD,OAAO,6BAA6B,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;AACzD,CAAC"}
+{"version":3,"file":"confidenceEllipse.js","sourceRoot":"","sources":["../../../src/scatter/core/confidenceEllipse.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAKH,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAE,6BAA6B,EAAE,MAAM,sBAAsB,CAAC;AAErE,OAAO,EAAE,mBAAmB,EAAE,MAAM,0BAA0B,CAAC;AAoC/D,8FAA8F;AAC9F,MAAM,CAAC,MAAM,oBAAoB,GAAgB;IAC/C,IAAI,EAAE,UAAU;IAChB,WAAW,EAAE,IAAI;CAClB,CAAC;AAoCF;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAA+B,EAC/B,UAAoC,EAAE;IAEtC,MAAM,EAAE,IAAI,GAAG,oBAAoB,EAAE,aAAa,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IACnE,MAAM,OAAO,GAAG,mBAAmB,CAAC,WAAW,CAAC,MAAM,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACrE,OAAO,cAAc,CAAC,OAAO,EAAE;QAC7B,KAAK,EAAE,CAAC;QACR,KAAK,EAAE,CAAC;QACR,KAAK,EAAE,CAAC;QACR,kBAAkB,EAAE,yBAAyB,CAAC,IAAI,CAAC;QACnD,aAAa;KACd,CAAC,CAAC;AACL,CAAC;AAgBD;;;;;;;;;GASG;AACH,MAAM,UAAU,cAAc,CAC5B,OAA4B,EAC5B,OAA8B;IAE9B,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,WAAW,EAAE,GAAG,OAAO,CAAC;IAC7D,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,kBAAkB,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC3E,IAAI,KAAK,GAAG,CAAC,IAAI,KAAK,IAAI,MAAM;QAAE,OAAO,IAAI,CAAC;IAC9C,IAAI,KAAK,GAAG,CAAC,IAAI,KAAK,GAAG,CAAC,IAAI,KAAK,IAAI,IAAI,IAAI,KAAK,IAAI,IAAI;QAAE,OAAO,IAAI,CAAC;IAC1E,MAAM,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IACjC,IAAI,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,aAAa,CAAC;QAAE,OAAO,IAAI,CAAC;IAEpD,MAAM,IAAI,GAAG,KAAK,GAAG,IAAI,CAAC;IAC1B,MAAM,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;IAClC,MAAM,UAAU,GAAsB;QACpC,EAAE,EAAE,WAAW,CAAC,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;QAClD,EAAE,EAAE,WAAW,CAAC,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;QAClD,EAAE,EAAE,WAAW,CAAC,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;KACnD,CAAC;IACF,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,GAAG,WAAW,CAAC,UAAU,EAAE,kBAAkB,CAAC,CAAC;IACtE,MAAM,EAAE,GAAG,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;IACpC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QAAE,OAAO,IAAI,CAAC;IAC9D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QAAE,OAAO,IAAI,CAAC;IAC9D,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,EAAE,CAAC;AACtD,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,yBAAyB,CAAC,IAAiB;IACzD,IAAI,IAAI,CAAC,IAAI,KAAK,oBAAoB,EAAE,CAAC;QACvC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,kBAAkB,EAAE,CAAC,CAAC,CAAC;IAC9C,CAAC;IACD,OAAO,6BAA6B,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;AACzD,CAAC;AAED;;;GAGG;AACH,SAAS,WAAW,CAAC,MAA+B;IAClD,OAAO;QACL,IAAI,EAAE,MAAM,CAAC,MAAM;QACnB,OAAO,EAAE,CAAC;QACV,GAAG,CAAC,QAAgB,EAAE,WAAmB;YACvC,MAAM,KAAK,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC;YAC/B,IAAI,KAAK,KAAK,SAAS;gBAAE,OAAO,MAAM,CAAC,GAAG,CAAC;YAC3C,IAAI,WAAW,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC,CAAC,CAAC;YACtC,OAAO,WAAW,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;QAClD,CAAC;KACF,CAAC;AACJ,CAAC"}

package/lib/scatter/core/ellipseAxes.d.ts

@@ -0,0 +1,40 @@
+/**
+ * The shape a covariance describes: two semi-axes and the turn of the longer
+ * one, from the closed form of a symmetric two by two eigen-decomposition.
+ *
+ * It is its own module because a map outline, a pair-grid cell and the
+ * silhouette of a shell in a cloud all end in this same decomposition, and a
+ * second copy of it is a second collinear floor to forget.
+ */
+/** How a cloud is spread, as the three distinct entries of its covariance. */
+export interface EllipseCovariance {
+    /** The spread along the horizontal axis. */
+    xx: number;
+    /** How the two axes move together; zero when they are unrelated. */
+    xy: number;
+    /** The spread along the vertical axis. */
+    yy: number;
+}
+/** The drawable shape of a covariance. */
+export interface EllipseAxes {
+    /** Semi-major axis length. */
+    rx: number;
+    /** Semi-minor axis length; exactly `0` when the spread is collinear. */
+    ry: number;
+    /** Rotation of the major axis, in radians counter-clockwise, in `(-π/2, π/2]`. */
+    angle: number;
+}
+/**
+ * The semi-axes and the rotation that a covariance describes.
+ *
+ * The shape of an outline is the eigen-decomposition of a symmetric two by two
+ * matrix, which has a closed form and so never needs an iterative solver. It
+ * is public because the same three numbers are what turns any covariance into
+ * something drawable — one measured elsewhere, or one already mapped into
+ * pixels.
+ * @param covariance - How the cloud is spread.
+ * @param standardDeviations - How far out to measure. Defaults to `1`.
+ * @returns The two semi-axis lengths and the rotation of the major one, in radians, in `(-π/2, π/2]`. A collinear cloud gives `ry` of exactly `0`, and one with no spread at all gives an angle of `0` because it has no direction to report.
+ */
+export declare function ellipseAxes(covariance: EllipseCovariance, standardDeviations?: number): EllipseAxes;
+//# sourceMappingURL=ellipseAxes.d.ts.map

package/lib/scatter/core/ellipseAxes.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipseAxes.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/ellipseAxes.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAcH,8EAA8E;AAC9E,MAAM,WAAW,iBAAiB;IAChC,4CAA4C;IAC5C,EAAE,EAAE,MAAM,CAAC;IACX,oEAAoE;IACpE,EAAE,EAAE,MAAM,CAAC;IACX,0CAA0C;IAC1C,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,0CAA0C;AAC1C,MAAM,WAAW,WAAW;IAC1B,8BAA8B;IAC9B,EAAE,EAAE,MAAM,CAAC;IACX,wEAAwE;IACxE,EAAE,EAAE,MAAM,CAAC;IACX,kFAAkF;IAClF,KAAK,EAAE,MAAM,CAAC;CACf;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,WAAW,CACzB,UAAU,EAAE,iBAAiB,EAC7B,kBAAkB,SAAI,GACrB,WAAW,CAab"}

package/lib/scatter/core/ellipseAxes.js

@@ -0,0 +1,57 @@
+/**
+ * The shape a covariance describes: two semi-axes and the turn of the longer
+ * one, from the closed form of a symmetric two by two eigen-decomposition.
+ *
+ * It is its own module because a map outline, a pair-grid cell and the
+ * silhouette of a shell in a cloud all end in this same decomposition, and a
+ * second copy of it is a second collinear floor to forget.
+ */
+/**
+ * A smaller eigenvalue below this share of the larger one is read as zero.
+ *
+ * A group whose spread across the major axis is under a millionth of its
+ * spread along it is collinear for any purpose a plot has, and the rounding of
+ * the covariance alone reaches that far: the nine points on a straight line
+ * used in the tests leave a residue of about 5e-17 where the exact answer is
+ * 0. Without the floor those points would draw a hairline ellipse on one
+ * machine and a segment on the next.
+ */
+const COLLINEAR_FLOOR = 1e-12;
+/**
+ * The semi-axes and the rotation that a covariance describes.
+ *
+ * The shape of an outline is the eigen-decomposition of a symmetric two by two
+ * matrix, which has a closed form and so never needs an iterative solver. It
+ * is public because the same three numbers are what turns any covariance into
+ * something drawable — one measured elsewhere, or one already mapped into
+ * pixels.
+ * @param covariance - How the cloud is spread.
+ * @param standardDeviations - How far out to measure. Defaults to `1`.
+ * @returns The two semi-axis lengths and the rotation of the major one, in radians, in `(-π/2, π/2]`. A collinear cloud gives `ry` of exactly `0`, and one with no spread at all gives an angle of `0` because it has no direction to report.
+ */
+export function ellipseAxes(covariance, standardDeviations = 1) {
+    const { xx, xy, yy } = covariance;
+    const middle = (xx + yy) / 2;
+    const half = (xx - yy) / 2;
+    const radius = Math.hypot(half, xy);
+    const major = middle + radius;
+    let minor = middle - radius;
+    if (!(minor > major * COLLINEAR_FLOOR))
+        minor = 0;
+    return {
+        rx: standardDeviations * Math.sqrt(major),
+        ry: standardDeviations * Math.sqrt(minor),
+        angle: majorAxisAngle(xx, xy, yy, major),
+    };
+}
+function majorAxisAngle(xx, xy, yy, major) {
+    // With no cross term the axes already are the eigenvectors, and asking
+    // `atan2` would only turn the exact answer into a rounded one.
+    if (xy === 0)
+        return xx >= yy ? 0 : Math.PI / 2;
+    // `major - xx` is never negative, so this lands in [0, π] and a negative
+    // cross term arrives half a turn away from the direction actually wanted.
+    const angle = Math.atan2(major - xx, xy);
+    return angle > Math.PI / 2 ? angle - Math.PI : angle;
+}
+//# sourceMappingURL=ellipseAxes.js.map

package/lib/scatter/core/ellipseAxes.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipseAxes.js","sourceRoot":"","sources":["../../../src/scatter/core/ellipseAxes.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH;;;;;;;;;GASG;AACH,MAAM,eAAe,GAAG,KAAK,CAAC;AAsB9B;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,WAAW,CACzB,UAA6B,EAC7B,kBAAkB,GAAG,CAAC;IAEtB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,UAAU,CAAC;IAClC,MAAM,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;IAC3B,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,EAAE,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,MAAM,GAAG,MAAM,CAAC;IAC9B,IAAI,KAAK,GAAG,MAAM,GAAG,MAAM,CAAC;IAC5B,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,GAAG,eAAe,CAAC;QAAE,KAAK,GAAG,CAAC,CAAC;IAClD,OAAO;QACL,EAAE,EAAE,kBAAkB,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC;QACzC,EAAE,EAAE,kBAAkB,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC;QACzC,KAAK,EAAE,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC;KACzC,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CACrB,EAAU,EACV,EAAU,EACV,EAAU,EACV,KAAa;IAEb,uEAAuE;IACvE,+DAA+D;IAC/D,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IAChD,yEAAyE;IACzE,0EAA0E;IAC1E,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC;IACzC,OAAO,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;AACvD,CAAC"}

package/lib/scatter/core/ellipseProjection.js

@@ -6,7 +6,7 @@
  * is drawn again at every zoom and every resize, and only this half has to be
  * redone.
  */
-import { ellipseAxes } from "./confidenceEllipse.js";
+import { ellipseAxes } from "./ellipseAxes.js";
 const HALF_TURN_DEGREES = 180;
 const DEFAULT_STEPS = 64;
 const MINIMUM_STEPS = 3;

package/lib/scatter/core/ellipseProjection.js.map

@@ -1 +1 @@
-{"version":3,"file":"ellipseProjection.js","sourceRoot":"","sources":["../../../src/scatter/core/ellipseProjection.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAKH,OAAO,EAAE,WAAW,EAAE,MAAM,wBAAwB,CAAC;AAErD,MAAM,iBAAiB,GAAG,GAAG,CAAC;AAC9B,MAAM,aAAa,GAAG,EAAE,CAAC;AACzB,MAAM,aAAa,GAAG,CAAC,CAAC;AAgBxB;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,UAAU,cAAc,CAC5B,OAA0B,EAC1B,CAAa,EACb,CAAa;IAEb,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACvC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,OAAO,CAAC,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;IACtC,MAAM,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;IAEvC,yEAAyE;IACzE,6CAA6C;IAC7C,MAAM,EAAE,GAAG,KAAK,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI,GAAG,IAAI,CAAC;IAC1D,MAAM,EAAE,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC;IAC5C,MAAM,EAAE,GAAG,KAAK,GAAG,IAAI,GAAG,IAAI,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC;IAE1D,0EAA0E;IAC1E,0EAA0E;IAC1E,oEAAoE;IACpE,gCAAgC;IAChC,MAAM,OAAO,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;IACpC,MAAM,SAAS,GAAG,WAAW,CAAC;QAC5B,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,EAAE;QAC5B,EAAE,EAAE,OAAO,GAAG,EAAE;QAChB,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,EAAE;KAC7B,CAAC,CAAC;IAEH,MAAM,MAAM,GAAiB;QAC3B,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM;QACpC,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM;QACpC,EAAE,EAAE,SAAS,CAAC,EAAE;QAChB,EAAE,EAAE,SAAS,CAAC,EAAE;QAChB,YAAY,EAAE,CAAC,SAAS,CAAC,KAAK,GAAG,iBAAiB,CAAC,GAAG,IAAI,CAAC,EAAE;KAC9D,CAAC;IACF,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,YAAY,CAAC,EACrC,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,cAAc,CAC5B,OAA0B,EAC1B,KAAK,GAAG,aAAa;IAErB,MAAM,KAAK,GAAG,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAClC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,aAAa,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QAC5C,CAAC,CAAC,aAAa,CAAC;IAClB,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACvC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACrC,MAAM,OAAO,GAAmB,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC;IACjD,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QAC3C,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;QAC3C,MAAM,KAAK,GAAG,OAAO,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC1C,MAAM,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC3C,OAAO,CAAC,KAAK,CAAC,GAAG;YACf,CAAC,EAAE,OAAO,CAAC,EAAE,GAAG,KAAK,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI;YAC9C,CAAC,EAAE,OAAO,CAAC,EAAE,GAAG,KAAK,GAAG,IAAI,GAAG,MAAM,GAAG,MAAM;SAC/C,CAAC;IACJ,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC"}
+{"version":3,"file":"ellipseProjection.js","sourceRoot":"","sources":["../../../src/scatter/core/ellipseProjection.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAKH,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAE/C,MAAM,iBAAiB,GAAG,GAAG,CAAC;AAC9B,MAAM,aAAa,GAAG,EAAE,CAAC;AACzB,MAAM,aAAa,GAAG,CAAC,CAAC;AAgBxB;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,UAAU,cAAc,CAC5B,OAA0B,EAC1B,CAAa,EACb,CAAa;IAEb,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACvC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,OAAO,CAAC,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;IACtC,MAAM,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;IAEvC,yEAAyE;IACzE,6CAA6C;IAC7C,MAAM,EAAE,GAAG,KAAK,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI,GAAG,IAAI,CAAC;IAC1D,MAAM,EAAE,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC;IAC5C,MAAM,EAAE,GAAG,KAAK,GAAG,IAAI,GAAG,IAAI,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC;IAE1D,0EAA0E;IAC1E,0EAA0E;IAC1E,oEAAoE;IACpE,gCAAgC;IAChC,MAAM,OAAO,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;IACpC,MAAM,SAAS,GAAG,WAAW,CAAC;QAC5B,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,EAAE;QAC5B,EAAE,EAAE,OAAO,GAAG,EAAE;QAChB,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,EAAE;KAC7B,CAAC,CAAC;IAEH,MAAM,MAAM,GAAiB;QAC3B,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM;QACpC,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM;QACpC,EAAE,EAAE,SAAS,CAAC,EAAE;QAChB,EAAE,EAAE,SAAS,CAAC,EAAE;QAChB,YAAY,EAAE,CAAC,SAAS,CAAC,KAAK,GAAG,iBAAiB,CAAC,GAAG,IAAI,CAAC,EAAE;KAC9D,CAAC;IACF,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,YAAY,CAAC,EACrC,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,cAAc,CAC5B,OAA0B,EAC1B,KAAK,GAAG,aAAa;IAErB,MAAM,KAAK,GAAG,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAClC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,aAAa,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QAC5C,CAAC,CAAC,aAAa,CAAC;IAClB,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACvC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACrC,MAAM,OAAO,GAAmB,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC;IACjD,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QAC3C,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;QAC3C,MAAM,KAAK,GAAG,OAAO,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC1C,MAAM,MAAM,GAAG,OAAO,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC3C,OAAO,CAAC,KAAK,CAAC,GAAG;YACf,CAAC,EAAE,OAAO,CAAC,EAAE,GAAG,KAAK,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI;YAC9C,CAAC,EAAE,OAAO,CAAC,EAAE,GAAG,KAAK,GAAG,IAAI,GAAG,MAAM,GAAG,MAAM;SAC/C,CAAC;IACJ,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC"}

package/lib/scatter/core/index.d.ts

@@ -9,14 +9,19 @@
  * both — which is why the outline is projected rather than handed to an SVG
  * `<ellipse>` with a data-space angle.
  */
-export type { ConfidenceEllipse, ConfidenceEllipseOptions, EllipseCovariance, EllipseCoverageSize, EllipsePoint, EllipseSize, EllipseStandardDeviationSize, } from './confidenceEllipse.ts';
-export { confidenceEllipse, ellipseAxes, ellipseStandardDeviations, } from './confidenceEllipse.ts';
+export type { ConfidenceEllipse, ConfidenceEllipseOptions, EllipseCoverageSize, EllipsePoint, EllipseSize, EllipseStandardDeviationSize, } from './confidenceEllipse.ts';
+export { confidenceEllipse, ellipseStandardDeviations, } from './confidenceEllipse.ts';
+export type { EllipseAxes, EllipseCovariance } from './ellipseAxes.ts';
+export { ellipseAxes } from './ellipseAxes.ts';
 export { coverageForStandardDeviations, standardDeviationsForCoverage, } from './ellipseCoverage.ts';
 export type { PixelEllipse } from './ellipseProjection.ts';
 export { ellipsePolygon, projectEllipse } from './ellipseProjection.ts';
 export type { LassoPath } from './lassoPath.ts';
 export { appendLassoPoint, createLassoPath, lassoPathData, resetLassoPath, } from './lassoPath.ts';
 export { pointInPolygon, pointsInPolygon, polygonBounds } from './polygon.ts';
+export type { ScatterGroupMoments } from './scatterGroupMoments.ts';
+export type { ScatterGroupSpread, ScatterGroupSpreadOptions, } from './scatterGroupSpread.ts';
+export { scatterGroupSpread, scatterPairEllipse, scatterSkippedGroups, } from './scatterGroupSpread.ts';
 export type { ScatterSelectionMode } from './scatterSelection.ts';
 export { mergeScatterSelection, scatterSelectedIndices, scatterSelectionMask, } from './scatterSelection.ts';
 export type { PixelBounds, ScreenPoints } from './screenPoints.ts';

package/lib/scatter/core/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH,YAAY,EACV,iBAAiB,EACjB,wBAAwB,EACxB,iBAAiB,EACjB,mBAAmB,EACnB,YAAY,EACZ,WAAW,EACX,4BAA4B,GAC7B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,iBAAiB,EACjB,WAAW,EACX,yBAAyB,GAC1B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,6BAA6B,EAC7B,6BAA6B,GAC9B,MAAM,sBAAsB,CAAC;AAC9B,YAAY,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AACxE,YAAY,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAChD,OAAO,EACL,gBAAgB,EAChB,eAAe,EACf,aAAa,EACb,cAAc,GACf,MAAM,gBAAgB,CAAC;AACxB,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,aAAa,EAAE,MAAM,cAAc,CAAC;AAC9E,YAAY,EAAE,oBAAoB,EAAE,MAAM,uBAAuB,CAAC;AAClE,OAAO,EACL,qBAAqB,EACrB,sBAAsB,EACtB,oBAAoB,GACrB,MAAM,uBAAuB,CAAC;AAC/B,YAAY,EAAE,WAAW,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACnE,OAAO,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,MAAM,mBAAmB,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH,YAAY,EACV,iBAAiB,EACjB,wBAAwB,EACxB,mBAAmB,EACnB,YAAY,EACZ,WAAW,EACX,4BAA4B,GAC7B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,iBAAiB,EACjB,yBAAyB,GAC1B,MAAM,wBAAwB,CAAC;AAChC,YAAY,EAAE,WAAW,EAAE,iBAAiB,EAAE,MAAM,kBAAkB,CAAC;AACvE,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EACL,6BAA6B,EAC7B,6BAA6B,GAC9B,MAAM,sBAAsB,CAAC;AAC9B,YAAY,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AACxE,YAAY,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAChD,OAAO,EACL,gBAAgB,EAChB,eAAe,EACf,aAAa,EACb,cAAc,GACf,MAAM,gBAAgB,CAAC;AACxB,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,aAAa,EAAE,MAAM,cAAc,CAAC;AAC9E,YAAY,EAAE,mBAAmB,EAAE,MAAM,0BAA0B,CAAC;AACpE,YAAY,EACV,kBAAkB,EAClB,yBAAyB,GAC1B,MAAM,yBAAyB,CAAC;AACjC,OAAO,EACL,kBAAkB,EAClB,kBAAkB,EAClB,oBAAoB,GACrB,MAAM,yBAAyB,CAAC;AACjC,YAAY,EAAE,oBAAoB,EAAE,MAAM,uBAAuB,CAAC;AAClE,OAAO,EACL,qBAAqB,EACrB,sBAAsB,EACtB,oBAAoB,GACrB,MAAM,uBAAuB,CAAC;AAC/B,YAAY,EAAE,WAAW,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AACnE,OAAO,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,MAAM,mBAAmB,CAAC"}

package/lib/scatter/core/index.js

@@ -9,11 +9,13 @@
  * both — which is why the outline is projected rather than handed to an SVG
  * `<ellipse>` with a data-space angle.
  */
-export { confidenceEllipse, ellipseAxes, ellipseStandardDeviations, } from "./confidenceEllipse.js";
+export { confidenceEllipse, ellipseStandardDeviations, } from "./confidenceEllipse.js";
+export { ellipseAxes } from "./ellipseAxes.js";
 export { coverageForStandardDeviations, standardDeviationsForCoverage, } from "./ellipseCoverage.js";
 export { ellipsePolygon, projectEllipse } from "./ellipseProjection.js";
 export { appendLassoPoint, createLassoPath, lassoPathData, resetLassoPath, } from "./lassoPath.js";
 export { pointInPolygon, pointsInPolygon, polygonBounds } from "./polygon.js";
+export { scatterGroupSpread, scatterPairEllipse, scatterSkippedGroups, } from "./scatterGroupSpread.js";
 export { mergeScatterSelection, scatterSelectedIndices, scatterSelectionMask, } from "./scatterSelection.js";
 export { nearestPointIndex, pointsWithinBounds } from "./screenPoints.js";
 //# sourceMappingURL=index.js.map

package/lib/scatter/core/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/scatter/core/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAWH,OAAO,EACL,iBAAiB,EACjB,WAAW,EACX,yBAAyB,GAC1B,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,6BAA6B,EAC7B,6BAA6B,GAC9B,MAAM,sBAAsB,CAAC;AAE9B,OAAO,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AAExE,OAAO,EACL,gBAAgB,EAChB,eAAe,EACf,aAAa,EACb,cAAc,GACf,MAAM,gBAAgB,CAAC;AACxB,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,aAAa,EAAE,MAAM,cAAc,CAAC;AAE9E,OAAO,EACL,qBAAqB,EACrB,sBAAsB,EACtB,oBAAoB,GACrB,MAAM,uBAAuB,CAAC;AAE/B,OAAO,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,MAAM,mBAAmB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/scatter/core/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAUH,OAAO,EACL,iBAAiB,EACjB,yBAAyB,GAC1B,MAAM,wBAAwB,CAAC;AAEhC,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EACL,6BAA6B,EAC7B,6BAA6B,GAC9B,MAAM,sBAAsB,CAAC;AAE9B,OAAO,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AAExE,OAAO,EACL,gBAAgB,EAChB,eAAe,EACf,aAAa,EACb,cAAc,GACf,MAAM,gBAAgB,CAAC;AACxB,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,aAAa,EAAE,MAAM,cAAc,CAAC;AAM9E,OAAO,EACL,kBAAkB,EAClB,kBAAkB,EAClB,oBAAoB,GACrB,MAAM,yBAAyB,CAAC;AAEjC,OAAO,EACL,qBAAqB,EACrB,sBAAsB,EACtB,oBAAoB,GACrB,MAAM,uBAAuB,CAAC;AAE/B,OAAO,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,MAAM,mBAAmB,CAAC"}

package/lib/scatter/core/lassoPath.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"lassoPath.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/lassoPath.ts"],"names":[],"mappings":"AAEA;;;;;;;;GAQG;AACH,MAAM,WAAW,SAAS;IACxB,8EAA8E;IAC9E,EAAE,EAAE,YAAY,CAAC;IACjB,4EAA4E;IAC5E,EAAE,EAAE,YAAY,CAAC;IACjB,wCAAwC;IACxC,MAAM,EAAE,MAAM,CAAC;CAChB;AAED;;;;GAIG;AACH,wBAAgB,eAAe,CAAC,QAAQ,SAAmB,GAAG,SAAS,CAKtE;AAED;;;;;;GAMG;AACH,wBAAgB,cAAc,CAAC,IAAI,EAAE,SAAS,GAAG,IAAI,CAEpD;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,gBAAgB,CAC9B,IAAI,EAAE,SAAS,EACf,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,WAAW,EAAE,MAAM,GAClB,OAAO,CAwBT;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,aAAa,CAAC,IAAI,EAAE,SAAS,EAAE,KAAK,EAAE,OAAO,GAAG,MAAM,CAWrE"}
+{"version":3,"file":"lassoPath.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/lassoPath.ts"],"names":[],"mappings":"AAIA;;;;;;;;GAQG;AACH,MAAM,WAAW,SAAS;IACxB,8EAA8E;IAC9E,EAAE,EAAE,YAAY,CAAC;IACjB,4EAA4E;IAC5E,EAAE,EAAE,YAAY,CAAC;IACjB,wCAAwC;IACxC,MAAM,EAAE,MAAM,CAAC;CAChB;AAED;;;;GAIG;AACH,wBAAgB,eAAe,CAAC,QAAQ,SAAmB,GAAG,SAAS,CAKtE;AAED;;;;;;GAMG;AACH,wBAAgB,cAAc,CAAC,IAAI,EAAE,SAAS,GAAG,IAAI,CAEpD;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,gBAAgB,CAC9B,IAAI,EAAE,SAAS,EACf,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,WAAW,EAAE,MAAM,GAClB,OAAO,CAwBT;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,aAAa,CAAC,IAAI,EAAE,SAAS,EAAE,KAAK,EAAE,OAAO,GAAG,MAAM,CAWrE"}

package/lib/scatter/core/lassoPath.js

@@ -1,3 +1,4 @@
+import { chartRoundPixel } from "../../chart/core/chartScale.js";
 const DEFAULT_CAPACITY = 256;
 /**
  * An empty path, with room to grow.
@@ -94,8 +95,6 @@ function growLassoPath(path) {
     path.ys = ys;
 }
 function roundToTenth(value) {
-    if (!Number.isFinite(value))
-        return 0;
-    return Math.round(value * 10) / 10;
+    return Number.isFinite(value) ? chartRoundPixel(value, 1) : 0;
 }
 //# sourceMappingURL=lassoPath.js.map

package/lib/scatter/core/lassoPath.js.map

@@ -1 +1 @@
-{"version":3,"file":"lassoPath.js","sourceRoot":"","sources":["../../../src/scatter/core/lassoPath.ts"],"names":[],"mappings":"AAAA,MAAM,gBAAgB,GAAG,GAAG,CAAC;AAoB7B;;;;GAIG;AACH,MAAM,UAAU,eAAe,CAAC,QAAQ,GAAG,gBAAgB;IACzD,MAAM,IAAI,GAAG,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;QACpC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QACnC,CAAC,CAAC,gBAAgB,CAAC;IACrB,OAAO,EAAE,EAAE,EAAE,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;AAC/E,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,cAAc,CAAC,IAAe;IAC5C,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;AAClB,CAAC;AAED;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,gBAAgB,CAC9B,IAAe,EACf,CAAS,EACT,CAAS,EACT,WAAmB;IAEnB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,OAAO,KAAK,CAAC;IAE7D,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACvC,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACvC,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;YAC/C,MAAM,IAAI,GACR,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,IAAI,WAAW,GAAG,CAAC;gBAC7C,CAAC,CAAC,WAAW,GAAG,WAAW;gBAC3B,CAAC,CAAC,CAAC,CAAC;YACR,MAAM,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;YACzB,MAAM,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;YACzB,IAAI,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI;gBAAE,OAAO,KAAK,CAAC;QAC7D,CAAC;IACH,CAAC;IAED,IAAI,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,EAAE,CAAC,MAAM,IAAI,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC;QACnE,aAAa,CAAC,IAAI,CAAC,CAAC;IACtB,CAAC;IACD,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACzB,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACzB,IAAI,CAAC,MAAM,IAAI,CAAC,CAAC;IACjB,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,aAAa,CAAC,IAAe,EAAE,KAAc;IAC3D,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IACpE,IAAI,IAAI,GAAG,EAAE,CAAC;IACd,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QAC3C,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;QACzB,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YAAE,MAAM;QAC9C,IAAI,IAAI,GAAG,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,YAAY,CAAC,CAAC,CAAC,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;IAC7E,CAAC;IACD,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,EAAE,CAAC;IACjC,OAAO,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC;AACnC,CAAC;AAED,SAAS,aAAa,CAAC,IAAe;IACpC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IAC3E,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACtC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChB,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChB,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACb,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;AACf,CAAC;AAED,SAAS,YAAY,CAAC,KAAa;IACjC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtC,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;AACrC,CAAC"}
+{"version":3,"file":"lassoPath.js","sourceRoot":"","sources":["../../../src/scatter/core/lassoPath.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,eAAe,EAAE,MAAM,gCAAgC,CAAC;AAEjE,MAAM,gBAAgB,GAAG,GAAG,CAAC;AAoB7B;;;;GAIG;AACH,MAAM,UAAU,eAAe,CAAC,QAAQ,GAAG,gBAAgB;IACzD,MAAM,IAAI,GAAG,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;QACpC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QACnC,CAAC,CAAC,gBAAgB,CAAC;IACrB,OAAO,EAAE,EAAE,EAAE,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;AAC/E,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,cAAc,CAAC,IAAe;IAC5C,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;AAClB,CAAC;AAED;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,gBAAgB,CAC9B,IAAe,EACf,CAAS,EACT,CAAS,EACT,WAAmB;IAEnB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,OAAO,KAAK,CAAC;IAE7D,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpB,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACvC,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACvC,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;YAC/C,MAAM,IAAI,GACR,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,IAAI,WAAW,GAAG,CAAC;gBAC7C,CAAC,CAAC,WAAW,GAAG,WAAW;gBAC3B,CAAC,CAAC,CAAC,CAAC;YACR,MAAM,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;YACzB,MAAM,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;YACzB,IAAI,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI;gBAAE,OAAO,KAAK,CAAC;QAC7D,CAAC;IACH,CAAC;IAED,IAAI,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,EAAE,CAAC,MAAM,IAAI,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC;QACnE,aAAa,CAAC,IAAI,CAAC,CAAC;IACtB,CAAC;IACD,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACzB,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACzB,IAAI,CAAC,MAAM,IAAI,CAAC,CAAC;IACjB,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,aAAa,CAAC,IAAe,EAAE,KAAc;IAC3D,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IACpE,IAAI,IAAI,GAAG,EAAE,CAAC;IACd,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QAC3C,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;QACzB,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YAAE,MAAM;QAC9C,IAAI,IAAI,GAAG,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,YAAY,CAAC,CAAC,CAAC,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;IAC7E,CAAC;IACD,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,EAAE,CAAC;IACjC,OAAO,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC;AACnC,CAAC;AAED,SAAS,aAAa,CAAC,IAAe;IACpC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IAC3E,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACtC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChB,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChB,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACb,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;AACf,CAAC;AAED,SAAS,YAAY,CAAC,KAAa;IACjC,OAAO,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,eAAe,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAChE,CAAC"}

package/lib/scatter/core/scatterGroupMoments.d.ts

@@ -0,0 +1,42 @@
+/**
+ * The one scan every group outline starts from: how many rows a group holds,
+ * where its middle is, and how it varies between every pair of axes.
+ *
+ * A map outline, a cell of a pair grid and a shell in a cloud all ask this of
+ * the same samples, and three copies of the loop are three places for the
+ * divisor or the treatment of a missing value to drift apart.
+ */
+import type { MatrixLike } from '../../chart/core/matrix.ts';
+/** How every group is spread over a run of axes. */
+export interface ScatterGroupMoments {
+    /** How many leading axes each block covers. */
+    axes: number;
+    /** How many groups it holds. */
+    groups: number;
+    /** How many rows went into each group, in group order. */
+    counts: Int32Array;
+    /** Each group's average along each axis, at `group * axes + axis`. */
+    means: Float64Array;
+    /**
+     * How each group varies between two axes, at `(group * axes + a) * axes + b`
+     * and again at the transposed position, so a reader reads either way round.
+     */
+    covariances: Float64Array;
+}
+/**
+ * Measure every group over the leading axes, in two passes over the scores
+ * read where they stand.
+ *
+ * A row missing a value in any one of the axes is left out of its group
+ * altogether, so every pair of axes describes the same samples. The spread
+ * uses the `n - 1` divisor, which is what makes it an estimate of the group
+ * rather than a statement about the handful of samples that happen to have
+ * been measured.
+ * @param scores - The coordinates: one row per sample, one column per axis.
+ * @param groupOf - Which group each row belongs to, or `null` when every row is one group.
+ * @param groups - How many groups there are.
+ * @param axes - How many leading axes to measure over.
+ * @returns The moments. A group of fewer than two rows keeps a zero spread.
+ */
+export declare function scatterGroupMoments(scores: MatrixLike, groupOf: ArrayLike<number> | null, groups: number, axes: number): ScatterGroupMoments;
+//# sourceMappingURL=scatterGroupMoments.d.ts.map

package/lib/scatter/core/scatterGroupMoments.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"scatterGroupMoments.d.ts","sourceRoot":"","sources":["../../../src/scatter/core/scatterGroupMoments.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAGH,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,4BAA4B,CAAC;AAE7D,oDAAoD;AACpD,MAAM,WAAW,mBAAmB;IAClC,+CAA+C;IAC/C,IAAI,EAAE,MAAM,CAAC;IACb,gCAAgC;IAChC,MAAM,EAAE,MAAM,CAAC;IACf,0DAA0D;IAC1D,MAAM,EAAE,UAAU,CAAC;IACnB,sEAAsE;IACtE,KAAK,EAAE,YAAY,CAAC;IACpB;;;OAGG;IACH,WAAW,EAAE,YAAY,CAAC;CAC3B;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,mBAAmB,CACjC,MAAM,EAAE,UAAU,EAClB,OAAO,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,IAAI,EACjC,MAAM,EAAE,MAAM,EACd,IAAI,EAAE,MAAM,GACX,mBAAmB,CA8DrB"}