@math.gl/geospatial 5.0.0-alpha.17 → 5.0.0-alpha.19

package/dist/projection/analysis.d.ts

@@ -0,0 +1,72 @@
+import type { ProjectionContext, ProjectionPlugin, ProjectionPoint } from "./experimental/types.js";
+/** Closed geographic rectangle in radians, without antimeridian wrapping. */
+export type ProjectionDomain = {
+    readonly west: number;
+    readonly east: number;
+    readonly south: number;
+    readonly north: number;
+};
+/** Physical projected-axis units per input radian. No datum or CRS axis/unit stages. */
+export type ProjectionJacobian = {
+    dxDLongitude: number;
+    dxDLatitude: number;
+    dyDLongitude: number;
+    dyDLatitude: number;
+};
+export type ProjectionFactors = ProjectionJacobian & {
+    meridionalScale: number;
+    parallelScale: number;
+    arealScale: number;
+    meridianConvergence: number;
+    meridianParallelAngle: number;
+    angularDistortion: number;
+    maximumScale: number;
+    minimumScale: number;
+};
+export type ProjectionAnalysisOptions = {
+    projection: ProjectionPlugin;
+    context: ProjectionContext;
+    /** Ellipsoid used to measure ground distances, independently of the projection kernel.
+     * Axes are in the kernel's physical output units. Defaults to context geometry.
+     * Changes factors only; projected coordinates and raw Jacobians are unchanged.
+     */
+    groundEllipsoid?: Pick<ProjectionContext, 'semiMajorAxis' | 'eccentricitySquared'>;
+    /** Application-qualified region, not a promise of a globally valid algorithm. */
+    domain: ProjectionDomain;
+    /** Central difference step in radians. Default 1e-4; stencil uses +/-2 steps. */
+    step?: number;
+    /** Agreement required between fourth-order stencils at step and half-step. */
+    derivativeTolerance?: number;
+};
+export declare function createProjectionJacobian(): ProjectionJacobian;
+export declare function createProjectionFactors(): ProjectionFactors;
+/** Optional single-projection analysis. Allocate scratch once; reuse all method outputs.
+ * false leaves outputs untouched. Plugin errors, nonfinite values, domain edges,
+ * singularities and unresolved derivatives are observable failures. Recursive plugin
+ * calls fail while scratch is in use; output setters can safely reenter afterwards.
+ */
+export declare class ProjectionAnalysis {
+    readonly domain: ProjectionDomain;
+    private readonly implementation;
+    private readonly a;
+    private readonly groundSemiMajorAxis;
+    private readonly groundEs;
+    private readonly step;
+    private readonly tolerance;
+    private readonly point;
+    private readonly derivative;
+    private busy;
+    private longitude;
+    private latitude;
+    private axis;
+    private readonly numerical;
+    constructor(options: ProjectionAnalysisOptions);
+    contains(longitude: number, latitude: number): boolean;
+    projectTo(longitude: number, latitude: number, result: ProjectionPoint): boolean;
+    unprojectTo(x: number, y: number, result: ProjectionPoint): boolean;
+    jacobian(longitude: number, latitude: number, result: ProjectionJacobian): boolean;
+    factors(longitude: number, latitude: number, result: ProjectionFactors): boolean;
+    private sample;
+    private calculate;
+}
+//# sourceMappingURL=analysis.d.ts.map

package/dist/projection/analysis.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"analysis.d.ts","sourceRoot":"","sources":["../../src/projection/analysis.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EACV,iBAAiB,EAEjB,gBAAgB,EAChB,eAAe,EAChB,gCAA6B;AAE9B,6EAA6E;AAC7E,MAAM,MAAM,gBAAgB,GAAG;IAC7B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACxB,CAAC;AACF,wFAAwF;AACxF,MAAM,MAAM,kBAAkB,GAAG;IAC/B,YAAY,EAAE,MAAM,CAAC;IACrB,WAAW,EAAE,MAAM,CAAC;IACpB,YAAY,EAAE,MAAM,CAAC;IACrB,WAAW,EAAE,MAAM,CAAC;CACrB,CAAC;AACF,MAAM,MAAM,iBAAiB,GAAG,kBAAkB,GAAG;IACnD,eAAe,EAAE,MAAM,CAAC;IACxB,aAAa,EAAE,MAAM,CAAC;IACtB,UAAU,EAAE,MAAM,CAAC;IACnB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,qBAAqB,EAAE,MAAM,CAAC;IAC9B,iBAAiB,EAAE,MAAM,CAAC;IAC1B,YAAY,EAAE,MAAM,CAAC;IACrB,YAAY,EAAE,MAAM,CAAC;CACtB,CAAC;AACF,MAAM,MAAM,yBAAyB,GAAG;IACtC,UAAU,EAAE,gBAAgB,CAAC;IAC7B,OAAO,EAAE,iBAAiB,CAAC;IAC3B;;;OAGG;IACH,eAAe,CAAC,EAAE,IAAI,CAAC,iBAAiB,EAAE,eAAe,GAAG,qBAAqB,CAAC,CAAC;IACnF,iFAAiF;IACjF,MAAM,EAAE,gBAAgB,CAAC;IACzB,iFAAiF;IACjF,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,8EAA8E;IAC9E,mBAAmB,CAAC,EAAE,MAAM,CAAC;CAC9B,CAAC;AACF,wBAAgB,wBAAwB,IAAI,kBAAkB,CAE7D;AACD,wBAAgB,uBAAuB,IAAI,iBAAiB,CAY3D;AAED;;;;GAIG;AACH,qBAAa,kBAAkB;IAC7B,QAAQ,CAAC,MAAM,EAAE,gBAAgB,CAAC;IAClC,OAAO,CAAC,QAAQ,CAAC,cAAc,CAA2B;IAC1D,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAS;IAC3B,OAAO,CAAC,QAAQ,CAAC,mBAAmB,CAAS;IAC7C,OAAO,CAAC,QAAQ,CAAC,QAAQ,CAAS;IAClC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAS;IAC9B,OAAO,CAAC,QAAQ,CAAC,SAAS,CAAS;IACnC,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAuC;IAE7D,OAAO,CAAC,QAAQ,CAAC,UAAU,CAAgB;IAC3C,OAAO,CAAC,IAAI,CAAS;IACrB,OAAO,CAAC,SAAS,CAAK;IACtB,OAAO,CAAC,QAAQ,CAAK;IACrB,OAAO,CAAC,IAAI,CAAK;IACjB,OAAO,CAAC,QAAQ,CAAC,SAAS,CAYvB;gBACS,OAAO,EAAE,yBAAyB;IA+C9C,QAAQ,CAAC,SAAS,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,OAAO;IAWtD,SAAS,CAAC,SAAS,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,eAAe,GAAG,OAAO;IAehF,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,eAAe,GAAG,OAAO;IAuBnE,QAAQ,CAAC,SAAS,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,kBAAkB,GAAG,OAAO;IAUlF,OAAO,CAAC,SAAS,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,iBAAiB,GAAG,OAAO;IA+ChF,OAAO,CAAC,MAAM;IAYd,OAAO,CAAC,SAAS;CAiBlB"}

package/dist/projection/analysis.js

@@ -0,0 +1,221 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original domain enforcement, numerical differentiation and cartographic factors.
+import { NumericalDerivative } from "./numerical-derivative.js";
+export function createProjectionJacobian() {
+    return { dxDLongitude: 0, dxDLatitude: 0, dyDLongitude: 0, dyDLatitude: 0 };
+}
+export function createProjectionFactors() {
+    return {
+        ...createProjectionJacobian(),
+        meridionalScale: 0,
+        parallelScale: 0,
+        arealScale: 0,
+        meridianConvergence: 0,
+        meridianParallelAngle: 0,
+        angularDistortion: 0,
+        maximumScale: 0,
+        minimumScale: 0
+    };
+}
+/** Optional single-projection analysis. Allocate scratch once; reuse all method outputs.
+ * false leaves outputs untouched. Plugin errors, nonfinite values, domain edges,
+ * singularities and unresolved derivatives are observable failures. Recursive plugin
+ * calls fail while scratch is in use; output setters can safely reenter afterwards.
+ */
+export class ProjectionAnalysis {
+    constructor(options) {
+        this.point = { x: 0, y: 0, z: 0 };
+        // Internal row-major Jacobian scratch; public outputs retain named fields.
+        this.derivative = [0, 0, 0, 0];
+        this.busy = false;
+        this.longitude = 0;
+        this.latitude = 0;
+        this.axis = 0;
+        this.numerical = new NumericalDerivative((offset, output) => {
+            if (!this.sample(this.longitude + (this.axis === 0 ? offset : 0), this.latitude + (this.axis === 1 ? offset : 0)))
+                return false;
+            output[0] = this.point.x;
+            output[1] = this.point.y;
+            output[2] = 0;
+            return true;
+        });
+        const { semiMajorAxis: a } = options.context;
+        const ground = options.groundEllipsoid ?? options.context;
+        const { west, east, south, north } = options.domain;
+        const step = options.step ?? 1e-4, tolerance = options.derivativeTolerance ?? 1e-6;
+        if (!([options.context, ground].every(({ semiMajorAxis, eccentricitySquared }) => Number.isFinite(semiMajorAxis) &&
+            semiMajorAxis > 0 &&
+            Number.isFinite(eccentricitySquared) &&
+            eccentricitySquared >= 0 &&
+            eccentricitySquared < 1) &&
+            Number.isFinite(west) &&
+            Number.isFinite(east) &&
+            west < east &&
+            west >= -Math.PI &&
+            east <= Math.PI &&
+            Number.isFinite(south) &&
+            Number.isFinite(north) &&
+            south < north &&
+            south >= -Math.PI / 2 &&
+            north <= Math.PI / 2 &&
+            Number.isFinite(step) &&
+            step > 0 &&
+            step < 0.1 &&
+            Number.isFinite(tolerance) &&
+            tolerance > 0 &&
+            tolerance < 1)) {
+            throw new Error('Invalid projection analysis geometry, domain or derivative options');
+        }
+        this.domain = Object.freeze({ west, east, south, north });
+        this.implementation = options.projection.create(options.context);
+        if (!this.implementation.forwardInPlace || !this.implementation.inverseInPlace) {
+            throw new Error('Projection analysis requires mutable horizontal projection hooks');
+        }
+        this.a = a;
+        this.groundSemiMajorAxis = ground.semiMajorAxis;
+        this.groundEs = ground.eccentricitySquared;
+        this.step = step;
+        this.tolerance = tolerance;
+    }
+    contains(longitude, latitude) {
+        const { west, east, south, north } = this.domain;
+        return (Number.isFinite(longitude) &&
+            Number.isFinite(latitude) &&
+            longitude >= west &&
+            longitude <= east &&
+            latitude >= south &&
+            latitude <= north);
+    }
+    projectTo(longitude, latitude, result) {
+        if (this.busy || !this.contains(longitude, latitude))
+            return false;
+        this.busy = true;
+        let x, y;
+        try {
+            if (!this.sample(longitude, latitude))
+                return false;
+            x = this.point.x;
+            y = this.point.y;
+        }
+        finally {
+            this.busy = false;
+        }
+        result.x = x;
+        result.y = y;
+        return true;
+    }
+    unprojectTo(x, y, result) {
+        if (this.busy || !Number.isFinite(x) || !Number.isFinite(y))
+            return false;
+        this.busy = true;
+        let longitude, latitude;
+        try {
+            this.point.x = x;
+            this.point.y = y;
+            this.point.z = 0;
+            try {
+                this.implementation.inverseInPlace(this.point);
+            }
+            catch {
+                return false;
+            }
+            longitude = this.point.x;
+            latitude = this.point.y;
+            if (!this.contains(longitude, latitude))
+                return false;
+        }
+        finally {
+            this.busy = false;
+        }
+        result.x = longitude;
+        result.y = latitude;
+        return true;
+    }
+    jacobian(longitude, latitude, result) {
+        if (!this.calculate(longitude, latitude))
+            return false;
+        // Snapshot before writing caller-owned outputs, whose setters can reenter.
+        const [xx, xy, yx, yy] = this.derivative;
+        result.dxDLongitude = xx;
+        result.dxDLatitude = xy;
+        result.dyDLongitude = yx;
+        result.dyDLatitude = yy;
+        return true;
+    }
+    factors(longitude, latitude, result) {
+        if (!this.calculate(longitude, latitude))
+            return false;
+        // Snapshot before writing caller-owned outputs, whose setters can reenter.
+        const [xx, xy, yx, yy] = this.derivative;
+        const a = this.groundSemiMajorAxis;
+        const es = this.groundEs;
+        const d = 1 - es * Math.sin(latitude) ** 2;
+        const parallelRadius = (a * Math.cos(latitude)) / Math.sqrt(d);
+        const meridionalRadius = (a * (1 - es)) / (d * Math.sqrt(d));
+        const ex = xx / parallelRadius, ey = yx / parallelRadius;
+        const nx = xy / meridionalRadius, ny = yy / meridionalRadius;
+        const k = Math.hypot(ex, ey), h = Math.hypot(nx, ny), area = ex * ny - ey * nx;
+        // Singular values of the local east/north map; retain a signed area determinant.
+        const sum = Math.hypot(ex + ny, ey - nx), difference = Math.hypot(ex - ny, ey + nx);
+        const maximum = (sum + difference) / 2;
+        const minimum = maximum > 0 ? Math.min(maximum, Math.abs(area) / maximum) : 0;
+        if (!(Number.isFinite(maximum) &&
+            minimum > maximum * 1e-12 &&
+            Number.isFinite(h) &&
+            Number.isFinite(k)))
+            return false;
+        const angle = Math.atan2(Math.abs(area), Math.abs(ex * nx + ey * ny));
+        const convergence = -Math.atan2(xy, yy);
+        const distortion = 2 * Math.asin(Math.min(1, (maximum - minimum) / (maximum + minimum)));
+        result.dxDLongitude = xx;
+        result.dxDLatitude = xy;
+        result.dyDLongitude = yx;
+        result.dyDLatitude = yy;
+        result.meridionalScale = h;
+        result.parallelScale = k;
+        result.arealScale = area;
+        result.meridianConvergence = convergence;
+        result.meridianParallelAngle = angle;
+        result.angularDistortion = distortion;
+        result.maximumScale = maximum;
+        result.minimumScale = minimum;
+        return true;
+    }
+    sample(longitude, latitude) {
+        if (!this.contains(longitude, latitude))
+            return false;
+        this.point.x = longitude;
+        this.point.y = latitude;
+        this.point.z = 0;
+        try {
+            this.implementation.forwardInPlace(this.point);
+        }
+        catch {
+            return false;
+        }
+        return Number.isFinite(this.point.x) && Number.isFinite(this.point.y);
+    }
+    calculate(longitude, latitude) {
+        if (this.busy)
+            return false;
+        this.busy = true;
+        try {
+            this.longitude = longitude;
+            this.latitude = latitude;
+            for (let axis = 0; axis < 2; axis++) {
+                this.axis = axis;
+                if (!this.numerical.evaluate(this.step, this.a, this.tolerance))
+                    return false;
+                this.derivative[axis] = this.numerical.result[0];
+                this.derivative[axis + 2] = this.numerical.result[1];
+            }
+            return true;
+        }
+        finally {
+            this.busy = false;
+        }
+    }
+}
+//# sourceMappingURL=analysis.js.map

package/dist/projection/analysis.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"analysis.js","sourceRoot":"","sources":["../../src/projection/analysis.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,qGAAqG;AACrG,OAAO,EAAC,mBAAmB,EAAC,kCAA+B;AA+C3D,MAAM,UAAU,wBAAwB;IACtC,OAAO,EAAC,YAAY,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAC,CAAC;AAC5E,CAAC;AACD,MAAM,UAAU,uBAAuB;IACrC,OAAO;QACL,GAAG,wBAAwB,EAAE;QAC7B,eAAe,EAAE,CAAC;QAClB,aAAa,EAAE,CAAC;QAChB,UAAU,EAAE,CAAC;QACb,mBAAmB,EAAE,CAAC;QACtB,qBAAqB,EAAE,CAAC;QACxB,iBAAiB,EAAE,CAAC;QACpB,YAAY,EAAE,CAAC;QACf,YAAY,EAAE,CAAC;KAChB,CAAC;AACJ,CAAC;AAED;;;;GAIG;AACH,MAAM,OAAO,kBAAkB;IA4B7B,YAAY,OAAkC;QApB7B,UAAK,GAAoB,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAC,CAAC;QAC7D,2EAA2E;QAC1D,eAAU,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACnC,SAAI,GAAG,KAAK,CAAC;QACb,cAAS,GAAG,CAAC,CAAC;QACd,aAAQ,GAAG,CAAC,CAAC;QACb,SAAI,GAAG,CAAC,CAAC;QACA,cAAS,GAAG,IAAI,mBAAmB,CAAC,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE;YACtE,IACE,CAAC,IAAI,CAAC,MAAM,CACV,IAAI,CAAC,SAAS,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAC/C,IAAI,CAAC,QAAQ,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAC/C;gBAED,OAAO,KAAK,CAAC;YACf,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;YACzB,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;YACzB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACd,OAAO,IAAI,CAAC;QACd,CAAC,CAAC,CAAC;QAED,MAAM,EAAC,aAAa,EAAE,CAAC,EAAC,GAAG,OAAO,CAAC,OAAO,CAAC;QAC3C,MAAM,MAAM,GAAG,OAAO,CAAC,eAAe,IAAI,OAAO,CAAC,OAAO,CAAC;QAC1D,MAAM,EAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAC,GAAG,OAAO,CAAC,MAAM,CAAC;QAClD,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,IAAI,IAAI,EAC/B,SAAS,GAAG,OAAO,CAAC,mBAAmB,IAAI,IAAI,CAAC;QAClD,IACE,CAAC,CACC,CAAC,OAAO,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC,KAAK,CAC7B,CAAC,EAAC,aAAa,EAAE,mBAAmB,EAAC,EAAE,EAAE,CACvC,MAAM,CAAC,QAAQ,CAAC,aAAa,CAAC;YAC9B,aAAa,GAAG,CAAC;YACjB,MAAM,CAAC,QAAQ,CAAC,mBAAmB,CAAC;YACpC,mBAAmB,IAAI,CAAC;YACxB,mBAAmB,GAAG,CAAC,CAC1B;YACD,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;YACrB,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;YACrB,IAAI,GAAG,IAAI;YACX,IAAI,IAAI,CAAC,IAAI,CAAC,EAAE;YAChB,IAAI,IAAI,IAAI,CAAC,EAAE;YACf,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;YACtB,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;YACtB,KAAK,GAAG,KAAK;YACb,KAAK,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC;YACrB,KAAK,IAAI,IAAI,CAAC,EAAE,GAAG,CAAC;YACpB,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;YACrB,IAAI,GAAG,CAAC;YACR,IAAI,GAAG,GAAG;YACV,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC;YAC1B,SAAS,GAAG,CAAC;YACb,SAAS,GAAG,CAAC,CACd,EACD,CAAC;YACD,MAAM,IAAI,KAAK,CAAC,oEAAoE,CAAC,CAAC;QACxF,CAAC;QACD,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,EAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAC,CAAC,CAAC;QACxD,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC,UAAU,CAAC,MAAM,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QACjE,IAAI,CAAC,IAAI,CAAC,cAAc,CAAC,cAAc,IAAI,CAAC,IAAI,CAAC,cAAc,CAAC,cAAc,EAAE,CAAC;YAC/E,MAAM,IAAI,KAAK,CAAC,kEAAkE,CAAC,CAAC;QACtF,CAAC;QACD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,mBAAmB,GAAG,MAAM,CAAC,aAAa,CAAC;QAChD,IAAI,CAAC,QAAQ,GAAG,MAAM,CAAC,mBAAmB,CAAC;QAC3C,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;IAC7B,CAAC;IACD,QAAQ,CAAC,SAAiB,EAAE,QAAgB;QAC1C,MAAM,EAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QAC/C,OAAO,CACL,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC;YAC1B,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;YACzB,SAAS,IAAI,IAAI;YACjB,SAAS,IAAI,IAAI;YACjB,QAAQ,IAAI,KAAK;YACjB,QAAQ,IAAI,KAAK,CAClB,CAAC;IACJ,CAAC;IACD,SAAS,CAAC,SAAiB,EAAE,QAAgB,EAAE,MAAuB;QACpE,IAAI,IAAI,CAAC,IAAI,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,SAAS,EAAE,QAAQ,CAAC;YAAE,OAAO,KAAK,CAAC;QACnE,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAS,EAAE,CAAS,CAAC;QACzB,IAAI,CAAC;YACH,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,QAAQ,CAAC;gBAAE,OAAO,KAAK,CAAC;YACpD,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;YACjB,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QACnB,CAAC;gBAAS,CAAC;YACT,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC;QACpB,CAAC;QACD,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,OAAO,IAAI,CAAC;IACd,CAAC;IACD,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,MAAuB;QACvD,IAAI,IAAI,CAAC,IAAI,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,OAAO,KAAK,CAAC;QAC1E,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,SAAiB,EAAE,QAAgB,CAAC;QACxC,IAAI,CAAC;YACH,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;YACjB,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;YACjB,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;YACjB,IAAI,CAAC;gBACH,IAAI,CAAC,cAAc,CAAC,cAAe,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YAClD,CAAC;YAAC,MAAM,CAAC;gBACP,OAAO,KAAK,CAAC;YACf,CAAC;YACD,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;YACzB,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;YACxB,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,SAAS,EAAE,QAAQ,CAAC;gBAAE,OAAO,KAAK,CAAC;QACxD,CAAC;gBAAS,CAAC;YACT,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC;QACpB,CAAC;QACD,MAAM,CAAC,CAAC,GAAG,SAAS,CAAC;QACrB,MAAM,CAAC,CAAC,GAAG,QAAQ,CAAC;QACpB,OAAO,IAAI,CAAC;IACd,CAAC;IACD,QAAQ,CAAC,SAAiB,EAAE,QAAgB,EAAE,MAA0B;QACtE,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,SAAS,EAAE,QAAQ,CAAC;YAAE,OAAO,KAAK,CAAC;QACvD,2EAA2E;QAC3E,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;QACzC,MAAM,CAAC,YAAY,GAAG,EAAE,CAAC;QACzB,MAAM,CAAC,WAAW,GAAG,EAAE,CAAC;QACxB,MAAM,CAAC,YAAY,GAAG,EAAE,CAAC;QACzB,MAAM,CAAC,WAAW,GAAG,EAAE,CAAC;QACxB,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,CAAC,SAAiB,EAAE,QAAgB,EAAE,MAAyB;QACpE,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,SAAS,EAAE,QAAQ,CAAC;YAAE,OAAO,KAAK,CAAC;QACvD,2EAA2E;QAC3E,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;QACzC,MAAM,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC;QACzB,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QAC3C,MAAM,cAAc,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/D,MAAM,gBAAgB,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7D,MAAM,EAAE,GAAG,EAAE,GAAG,cAAc,EAC5B,EAAE,GAAG,EAAE,GAAG,cAAc,CAAC;QAC3B,MAAM,EAAE,GAAG,EAAE,GAAG,gBAAgB,EAC9B,EAAE,GAAG,EAAE,GAAG,gBAAgB,CAAC;QAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,EAC1B,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,EACtB,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC3B,iFAAiF;QACjF,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,EACtC,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;QAC5C,MAAM,OAAO,GAAG,CAAC,GAAG,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC;QACvC,MAAM,OAAO,GAAG,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9E,IACE,CAAC,CACC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;YACxB,OAAO,GAAG,OAAO,GAAG,KAAK;YACzB,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YAClB,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CACnB;YAED,OAAO,KAAK,CAAC;QACf,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QACtE,MAAM,WAAW,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACxC,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACzF,MAAM,CAAC,YAAY,GAAG,EAAE,CAAC;QACzB,MAAM,CAAC,WAAW,GAAG,EAAE,CAAC;QACxB,MAAM,CAAC,YAAY,GAAG,EAAE,CAAC;QACzB,MAAM,CAAC,WAAW,GAAG,EAAE,CAAC;QACxB,MAAM,CAAC,eAAe,GAAG,CAAC,CAAC;QAC3B,MAAM,CAAC,aAAa,GAAG,CAAC,CAAC;QACzB,MAAM,CAAC,UAAU,GAAG,IAAI,CAAC;QACzB,MAAM,CAAC,mBAAmB,GAAG,WAAW,CAAC;QACzC,MAAM,CAAC,qBAAqB,GAAG,KAAK,CAAC;QACrC,MAAM,CAAC,iBAAiB,GAAG,UAAU,CAAC;QACtC,MAAM,CAAC,YAAY,GAAG,OAAO,CAAC;QAC9B,MAAM,CAAC,YAAY,GAAG,OAAO,CAAC;QAC9B,OAAO,IAAI,CAAC;IACd,CAAC;IACO,MAAM,CAAC,SAAiB,EAAE,QAAgB;QAChD,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,SAAS,EAAE,QAAQ,CAAC;YAAE,OAAO,KAAK,CAAC;QACtD,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,SAAS,CAAC;QACzB,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,QAAQ,CAAC;QACxB,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC;YACH,IAAI,CAAC,cAAc,CAAC,cAAe,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAClD,CAAC;QAAC,MAAM,CAAC;YACP,OAAO,KAAK,CAAC;QACf,CAAC;QACD,OAAO,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACO,SAAS,CAAC,SAAiB,EAAE,QAAgB;QACnD,IAAI,IAAI,CAAC,IAAI;YAAE,OAAO,KAAK,CAAC;QAC5B,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC;YACH,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;YAC3B,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;YACzB,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;gBACpC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;gBACjB,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC;oBAAE,OAAO,KAAK,CAAC;gBAC9E,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;gBACjD,IAAI,CAAC,UAAU,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;YACvD,CAAC;YACD,OAAO,IAAI,CAAC;QACd,CAAC;gBAAS,CAAC;YACT,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC;QACpB,CAAC;IACH,CAAC;CACF"}

package/dist/projection/bulk.cjs

@@ -0,0 +1,203 @@
+var __defProp = Object.defineProperty;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
+var __getOwnPropNames = Object.getOwnPropertyNames;
+var __hasOwnProp = Object.prototype.hasOwnProperty;
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
+
+// dist/projection/bulk.js
+var bulk_exports = {};
+__export(bulk_exports, {
+  ProjectionBuffer: () => ProjectionBuffer
+});
+module.exports = __toCommonJS(bulk_exports);
+function createScratch(dimension) {
+  return { input: new Float64Array(dimension), output: new Float64Array(dimension) };
+}
+function isBuffer(value) {
+  return value instanceof Float32Array || value instanceof Float64Array;
+}
+function shared(buffer) {
+  return typeof SharedArrayBuffer !== "undefined" && buffer instanceof SharedArrayBuffer;
+}
+function overlap(a, aStart, aEnd, b, bStart, bEnd) {
+  return (a.buffer === b.buffer || shared(a.buffer) && shared(b.buffer)) && aStart < bEnd && bStart < aEnd;
+}
+function capacity(length, offset, stride, width) {
+  return length < offset + width ? 0 : Math.floor((length - offset - width) / stride) + 1;
+}
+var ProjectionBuffer = class {
+  constructor(options) {
+    var _a, _b;
+    this.depth = 0;
+    const dimension = options.dimension ?? 2;
+    const inputOffset = options.inputOffset ?? 0, outputOffset = options.outputOffset ?? 0;
+    if (!Number.isSafeInteger(dimension) || dimension < 2 || !Number.isSafeInteger(inputOffset) || inputOffset < 0 || !Number.isSafeInteger(outputOffset) || outputOffset < 0 || options.inputStride !== void 0 && (!Number.isSafeInteger(options.inputStride) || options.inputStride < 1) || options.outputStride !== void 0 && (!Number.isSafeInteger(options.outputStride) || options.outputStride < 1) || typeof ((_a = options.projection) == null ? void 0 : _a.projectToSync) !== "function" || typeof ((_b = options.projection) == null ? void 0 : _b.unprojectToSync) !== "function") {
+      throw new Error("Invalid projection buffer transform or layout");
+    }
+    this.dimension = dimension;
+    this.projection = options.projection;
+    this.inputOffset = inputOffset;
+    this.outputOffset = outputOffset;
+    this.inputStride = options.inputStride;
+    this.outputStride = options.outputStride;
+    this.scratch = [createScratch(dimension)];
+  }
+  projectFlatTo(input, output, count, start = 0, epochs) {
+    return this.flat(input, output, count, start, false, epochs);
+  }
+  unprojectFlatTo(input, output, count, start = 0, epochs) {
+    return this.flat(input, output, count, start, true, epochs);
+  }
+  projectColumnsTo(input, output, count, start = 0, epochs) {
+    return this.columns(input, output, count, start, false, epochs);
+  }
+  unprojectColumnsTo(input, output, count, start = 0, epochs) {
+    return this.columns(input, output, count, start, true, epochs);
+  }
+  range(available, count, start) {
+    const records = count ?? Math.max(0, available - start);
+    if (!Number.isSafeInteger(start) || start < 0 || !Number.isSafeInteger(records) || records < 0 || !Number.isSafeInteger(start + records) || start + records > available)
+      throw new Error("Invalid projection buffer record range");
+    return records;
+  }
+  epochStorage(epochs, start, count) {
+    if (epochs === void 0)
+      return;
+    if (typeof epochs === "number") {
+      if (!Number.isFinite(epochs))
+        throw new Error("Coordinate epoch must be finite");
+    } else if (!isBuffer(epochs) || epochs.length < start + count)
+      throw new Error("Epoch buffer requires one value per addressed record");
+  }
+  aliases(input, output, start, count, inputStride, outputStride, width) {
+    if (!count)
+      return;
+    const ia = input.byteOffset + (this.inputOffset + start * inputStride) * input.BYTES_PER_ELEMENT;
+    const ib = ia + ((count - 1) * inputStride + width) * input.BYTES_PER_ELEMENT;
+    const oa = output.byteOffset + (this.outputOffset + start * outputStride) * output.BYTES_PER_ELEMENT;
+    const ob = oa + ((count - 1) * outputStride + width) * output.BYTES_PER_ELEMENT;
+    const exact = input.buffer === output.buffer && input.constructor === output.constructor && ia === oa && inputStride === outputStride;
+    if (overlap(input, ia, ib, output, oa, ob) && !exact)
+      throw new Error("Projection input/output ranges overlap without an exact in-place mapping");
+  }
+  epochAlias(epochs, output, start, count, stride, width) {
+    if (typeof epochs !== "object" || !epochs || !count)
+      return;
+    const a = epochs.byteOffset + start * epochs.BYTES_PER_ELEMENT;
+    const b = a + count * epochs.BYTES_PER_ELEMENT;
+    const c = output.byteOffset + (this.outputOffset + start * stride) * output.BYTES_PER_ELEMENT;
+    const d = c + ((count - 1) * stride + width) * output.BYTES_PER_ELEMENT;
+    if (overlap(epochs, a, b, output, c, d))
+      throw new Error("Epoch and output buffers must not overlap");
+  }
+  acquire() {
+    var _a;
+    const depth = this.depth++;
+    (_a = this.scratch)[depth] || (_a[depth] = createScratch(this.dimension));
+    return this.scratch[depth];
+  }
+  transform(scratch, inverse, epoch) {
+    if (epoch !== void 0 && !Number.isFinite(epoch))
+      throw new Error("Coordinate epoch must be finite");
+    if (inverse)
+      this.projection.unprojectToSync(scratch.input, scratch.output, epoch);
+    else
+      this.projection.projectToSync(scratch.input, scratch.output, epoch);
+  }
+  floatRange(result, float32) {
+    for (let i = 0; i < this.dimension; i++) {
+      const value = result[i];
+      if (i < 3 && !Number.isFinite(value) || float32 && Number.isFinite(value) && Math.abs(value) > 34028234663852886e22)
+        throw new Error("Projection output is nonfinite or exceeds Float32 range");
+    }
+  }
+  flat(input, output, count, start, inverse, epochs) {
+    if (!isBuffer(input) || !isBuffer(output))
+      throw new Error("Projection buffers require Float32Array or Float64Array");
+    const inputStride = this.inputStride ?? this.dimension, outputStride = this.outputStride ?? this.dimension;
+    if (inputStride < this.dimension || outputStride < this.dimension)
+      throw new Error("Interleaved stride must contain every ordinate");
+    const records = this.range(capacity(input.length, this.inputOffset, inputStride, this.dimension), count, start);
+    if (start + records > capacity(output.length, this.outputOffset, outputStride, this.dimension))
+      throw new Error("Projection output buffer is too short");
+    this.epochStorage(epochs, start, records);
+    this.aliases(input, output, start, records, inputStride, outputStride, this.dimension);
+    this.epochAlias(epochs, output, start, records, outputStride, this.dimension);
+    const scratch = this.acquire();
+    try {
+      const epochBuffer = typeof epochs === "number" ? void 0 : epochs;
+      const float32 = output instanceof Float32Array;
+      for (let record = start; record < start + records; record++) {
+        const source = this.inputOffset + record * inputStride, target = this.outputOffset + record * outputStride;
+        for (let i = 0; i < this.dimension; i++)
+          scratch.input[i] = input[source + i];
+        this.transform(scratch, inverse, epochBuffer ? epochBuffer[record] : epochs);
+        this.floatRange(scratch.output, float32);
+        for (let i = 0; i < this.dimension; i++)
+          output[target + i] = scratch.output[i];
+      }
+      return output;
+    } finally {
+      this.depth--;
+    }
+  }
+  columns(input, output, count, start, inverse, epochs) {
+    if (input.length !== this.dimension || output.length !== this.dimension)
+      throw new Error("Projection columns must match dimension");
+    const inputStride = this.inputStride ?? 1, outputStride = this.outputStride ?? 1;
+    let available = Number.MAX_SAFE_INTEGER;
+    for (let i = 0; i < this.dimension; i++) {
+      if (!isBuffer(input[i]) || !isBuffer(output[i]))
+        throw new Error("Projection columns require floating typed arrays");
+      available = Math.min(available, capacity(input[i].length, this.inputOffset, inputStride, 1));
+    }
+    const records = this.range(available, count, start);
+    this.epochStorage(epochs, start, records);
+    for (let i = 0; i < this.dimension; i++) {
+      if (start + records > capacity(output[i].length, this.outputOffset, outputStride, 1))
+        throw new Error("Projection output column is too short");
+      this.epochAlias(epochs, output[i], start, records, outputStride, 1);
+      for (let j = 0; j < this.dimension; j++)
+        this.aliases(input[i], output[j], start, records, inputStride, outputStride, 1);
+      for (let j = 0; j < i && records; j++) {
+        const a = output[i].byteOffset + (this.outputOffset + start * outputStride) * output[i].BYTES_PER_ELEMENT;
+        const b = output[j].byteOffset + (this.outputOffset + start * outputStride) * output[j].BYTES_PER_ELEMENT;
+        if (overlap(output[i], a, a + ((records - 1) * outputStride + 1) * output[i].BYTES_PER_ELEMENT, output[j], b, b + ((records - 1) * outputStride + 1) * output[j].BYTES_PER_ELEMENT))
+          throw new Error("Projection output columns overlap");
+      }
+    }
+    const scratch = this.acquire();
+    try {
+      const epochBuffer = typeof epochs === "number" ? void 0 : epochs;
+      for (let record = start; record < start + records; record++) {
+        const source = this.inputOffset + record * inputStride, target = this.outputOffset + record * outputStride;
+        for (let i = 0; i < this.dimension; i++)
+          scratch.input[i] = input[i][source];
+        this.transform(scratch, inverse, epochBuffer ? epochBuffer[record] : epochs);
+        for (let i = 0; i < this.dimension; i++) {
+          const value = scratch.output[i];
+          if (i < 3 && !Number.isFinite(value) || output[i] instanceof Float32Array && Number.isFinite(value) && Math.abs(value) > 34028234663852886e22)
+            throw new Error("Projection output is nonfinite or exceeds Float32 range");
+        }
+        for (let i = 0; i < this.dimension; i++)
+          output[i][target] = scratch.output[i];
+      }
+      return output;
+    } finally {
+      this.depth--;
+    }
+  }
+};
+//# sourceMappingURL=bulk.cjs.map

package/dist/projection/bulk.cjs.map

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../src/projection/bulk.ts"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAkCA;;;;;SAAS,cAAc,WAAiB;AACtC,SAAO,EAAC,OAAO,IAAI,aAAa,SAAS,GAAG,QAAQ,IAAI,aAAa,SAAS,EAAC;AACjF;AACA,SAAS,SAAS,OAAc;AAC9B,SAAO,iBAAiB,gBAAgB,iBAAiB;AAC3D;AACA,SAAS,OAAO,QAAuB;AACrC,SAAO,OAAO,sBAAsB,eAAe,kBAAkB;AACvE;AACA,SAAS,QACP,GACA,QACA,MACA,GACA,QACA,MAAY;AAEZ,UACG,EAAE,WAAW,EAAE,UAAW,OAAO,EAAE,MAAM,KAAK,OAAO,EAAE,MAAM,MAC9D,SAAS,QACT,SAAS;AAEb;AACA,SAAS,SAAS,QAAgB,QAAgB,QAAgB,OAAa;AAC7E,SAAO,SAAS,SAAS,QAAQ,IAAI,KAAK,OAAO,SAAS,SAAS,SAAS,MAAM,IAAI;AACxF;AAQM,IAAO,mBAAP,MAAuB;EAS3B,YAAY,SAAgC;AA1C9C;AAyCU,SAAA,QAAQ;AAEd,UAAM,YAAY,QAAQ,aAAa;AACvC,UAAM,cAAc,QAAQ,eAAe,GACzC,eAAe,QAAQ,gBAAgB;AACzC,QACE,CAAC,OAAO,cAAc,SAAS,KAC/B,YAAY,KACZ,CAAC,OAAO,cAAc,WAAW,KACjC,cAAc,KACd,CAAC,OAAO,cAAc,YAAY,KAClC,eAAe,KACd,QAAQ,gBAAgB,WACtB,CAAC,OAAO,cAAc,QAAQ,WAAW,KAAK,QAAQ,cAAc,MACtE,QAAQ,iBAAiB,WACvB,CAAC,OAAO,cAAc,QAAQ,YAAY,KAAK,QAAQ,eAAe,MACzE,SAAO,aAAQ,eAAR,mBAAoB,mBAAkB,cAC7C,SAAO,aAAQ,eAAR,mBAAoB,qBAAoB,YAC/C;AACA,YAAM,IAAI,MAAM,+CAA+C;IACjE;AACA,SAAK,YAAY;AACjB,SAAK,aAAa,QAAQ;AAC1B,SAAK,cAAc;AACnB,SAAK,eAAe;AACpB,SAAK,cAAc,QAAQ;AAC3B,SAAK,eAAe,QAAQ;AAC5B,SAAK,UAAU,CAAC,cAAc,SAAS,CAAC;EAC1C;EACA,cACE,OACA,QACA,OACA,QAAQ,GACR,QAAuB;AAEvB,WAAO,KAAK,KAAK,OAAO,QAAQ,OAAO,OAAO,OAAO,MAAM;EAC7D;EACA,gBACE,OACA,QACA,OACA,QAAQ,GACR,QAAuB;AAEvB,WAAO,KAAK,KAAK,OAAO,QAAQ,OAAO,OAAO,MAAM,MAAM;EAC5D;EACA,iBACE,OACA,QACA,OACA,QAAQ,GACR,QAAuB;AAEvB,WAAO,KAAK,QAAQ,OAAO,QAAQ,OAAO,OAAO,OAAO,MAAM;EAChE;EACA,mBACE,OACA,QACA,OACA,QAAQ,GACR,QAAuB;AAEvB,WAAO,KAAK,QAAQ,OAAO,QAAQ,OAAO,OAAO,MAAM,MAAM;EAC/D;EACQ,MAAM,WAAmB,OAA2B,OAAa;AACvE,UAAM,UAAU,SAAS,KAAK,IAAI,GAAG,YAAY,KAAK;AACtD,QACE,CAAC,OAAO,cAAc,KAAK,KAC3B,QAAQ,KACR,CAAC,OAAO,cAAc,OAAO,KAC7B,UAAU,KACV,CAAC,OAAO,cAAc,QAAQ,OAAO,KACrC,QAAQ,UAAU;AAElB,YAAM,IAAI,MAAM,wCAAwC;AAC1D,WAAO;EACT;EACQ,aAAa,QAAoC,OAAe,OAAa;AACnF,QAAI,WAAW;AAAW;AAC1B,QAAI,OAAO,WAAW,UAAU;AAC9B,UAAI,CAAC,OAAO,SAAS,MAAM;AAAG,cAAM,IAAI,MAAM,iCAAiC;IACjF,WAAW,CAAC,SAAS,MAAM,KAAK,OAAO,SAAS,QAAQ;AACtD,YAAM,IAAI,MAAM,sDAAsD;EAC1E;EACQ,QACN,OACA,QACA,OACA,OACA,aACA,cACA,OAAa;AAEb,QAAI,CAAC;AAAO;AACZ,UAAM,KACJ,MAAM,cAAc,KAAK,cAAc,QAAQ,eAAe,MAAM;AACtE,UAAM,KAAK,OAAO,QAAQ,KAAK,cAAc,SAAS,MAAM;AAC5D,UAAM,KACJ,OAAO,cAAc,KAAK,eAAe,QAAQ,gBAAgB,OAAO;AAC1E,UAAM,KAAK,OAAO,QAAQ,KAAK,eAAe,SAAS,OAAO;AAC9D,UAAM,QACJ,MAAM,WAAW,OAAO,UACxB,MAAM,gBAAgB,OAAO,eAC7B,OAAO,MACP,gBAAgB;AAClB,QAAI,QAAQ,OAAO,IAAI,IAAI,QAAQ,IAAI,EAAE,KAAK,CAAC;AAC7C,YAAM,IAAI,MAAM,0EAA0E;EAC9F;EACQ,WACN,QACA,QACA,OACA,OACA,QACA,OAAa;AAEb,QAAI,OAAO,WAAW,YAAY,CAAC,UAAU,CAAC;AAAO;AACrD,UAAM,IAAI,OAAO,aAAa,QAAQ,OAAO;AAC7C,UAAM,IAAI,IAAI,QAAQ,OAAO;AAC7B,UAAM,IAAI,OAAO,cAAc,KAAK,eAAe,QAAQ,UAAU,OAAO;AAC5E,UAAM,IAAI,MAAM,QAAQ,KAAK,SAAS,SAAS,OAAO;AACtD,QAAI,QAAQ,QAAQ,GAAG,GAAG,QAAQ,GAAG,CAAC;AACpC,YAAM,IAAI,MAAM,2CAA2C;EAC/D;EACQ,UAAO;;AACb,UAAM,QAAQ,KAAK;AACnB,KAAA,KAAA,KAAK,SAAQ,KAAK,MAAA,GAAL,KAAK,IAAM,cAAc,KAAK,SAAS;AACpD,WAAO,KAAK,QAAQ,KAAK;EAC3B;EACQ,UAAU,SAAkB,SAAkB,OAAyB;AAC7E,QAAI,UAAU,UAAa,CAAC,OAAO,SAAS,KAAK;AAC/C,YAAM,IAAI,MAAM,iCAAiC;AACnD,QAAI;AAAS,WAAK,WAAW,gBAAgB,QAAQ,OAAO,QAAQ,QAAQ,KAAK;;AAC5E,WAAK,WAAW,cAAc,QAAQ,OAAO,QAAQ,QAAQ,KAAK;EACzE;EACQ,WAAW,QAAsB,SAAgB;AAEvD,aAAS,IAAI,GAAG,IAAI,KAAK,WAAW,KAAK;AACvC,YAAM,QAAQ,OAAO,CAAC;AACtB,UACG,IAAI,KAAK,CAAC,OAAO,SAAS,KAAK,KAC/B,WAAW,OAAO,SAAS,KAAK,KAAK,KAAK,IAAI,KAAK,IAAI;AAExD,cAAM,IAAI,MAAM,yDAAyD;IAC7E;EACF;EACQ,KACN,OACA,QACA,OACA,OACA,SACA,QAAuB;AAEvB,QAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS,MAAM;AACtC,YAAM,IAAI,MAAM,yDAAyD;AAC3E,UAAM,cAAc,KAAK,eAAe,KAAK,WAC3C,eAAe,KAAK,gBAAgB,KAAK;AAC3C,QAAI,cAAc,KAAK,aAAa,eAAe,KAAK;AACtD,YAAM,IAAI,MAAM,gDAAgD;AAClE,UAAM,UAAU,KAAK,MACnB,SAAS,MAAM,QAAQ,KAAK,aAAa,aAAa,KAAK,SAAS,GACpE,OACA,KAAK;AAEP,QAAI,QAAQ,UAAU,SAAS,OAAO,QAAQ,KAAK,cAAc,cAAc,KAAK,SAAS;AAC3F,YAAM,IAAI,MAAM,uCAAuC;AACzD,SAAK,aAAa,QAAQ,OAAO,OAAO;AACxC,SAAK,QAAQ,OAAO,QAAQ,OAAO,SAAS,aAAa,cAAc,KAAK,SAAS;AACrF,SAAK,WAAW,QAAQ,QAAQ,OAAO,SAAS,cAAc,KAAK,SAAS;AAC5E,UAAM,UAAU,KAAK,QAAO;AAC5B,QAAI;AACF,YAAM,cAAc,OAAO,WAAW,WAAW,SAAY;AAC7D,YAAM,UAAU,kBAAkB;AAClC,eAAS,SAAS,OAAO,SAAS,QAAQ,SAAS,UAAU;AAC3D,cAAM,SAAS,KAAK,cAAc,SAAS,aACzC,SAAS,KAAK,eAAe,SAAS;AACxC,iBAAS,IAAI,GAAG,IAAI,KAAK,WAAW;AAAK,kBAAQ,MAAM,CAAC,IAAI,MAAM,SAAS,CAAC;AAC5E,aAAK,UACH,SACA,SACA,cAAc,YAAY,MAAM,IAAK,MAA6B;AAEpE,aAAK,WAAW,QAAQ,QAAQ,OAAO;AACvC,iBAAS,IAAI,GAAG,IAAI,KAAK,WAAW;AAAK,iBAAO,SAAS,CAAC,IAAI,QAAQ,OAAO,CAAC;MAChF;AACA,aAAO;IACT;AACE,WAAK;IACP;EACF;EACQ,QACN,OACA,QACA,OACA,OACA,SACA,QAAuB;AAEvB,QAAI,MAAM,WAAW,KAAK,aAAa,OAAO,WAAW,KAAK;AAC5D,YAAM,IAAI,MAAM,yCAAyC;AAC3D,UAAM,cAAc,KAAK,eAAe,GACtC,eAAe,KAAK,gBAAgB;AACtC,QAAI,YAAY,OAAO;AACvB,aAAS,IAAI,GAAG,IAAI,KAAK,WAAW,KAAK;AACvC,UAAI,CAAC,SAAS,MAAM,CAAC,CAAC,KAAK,CAAC,SAAS,OAAO,CAAC,CAAC;AAC5C,cAAM,IAAI,MAAM,kDAAkD;AACpE,kBAAY,KAAK,IAAI,WAAW,SAAS,MAAM,CAAC,EAAE,QAAQ,KAAK,aAAa,aAAa,CAAC,CAAC;IAC7F;AACA,UAAM,UAAU,KAAK,MAAM,WAAW,OAAO,KAAK;AAClD,SAAK,aAAa,QAAQ,OAAO,OAAO;AACxC,aAAS,IAAI,GAAG,IAAI,KAAK,WAAW,KAAK;AACvC,UAAI,QAAQ,UAAU,SAAS,OAAO,CAAC,EAAE,QAAQ,KAAK,cAAc,cAAc,CAAC;AACjF,cAAM,IAAI,MAAM,uCAAuC;AACzD,WAAK,WAAW,QAAQ,OAAO,CAAC,GAAG,OAAO,SAAS,cAAc,CAAC;AAClE,eAAS,IAAI,GAAG,IAAI,KAAK,WAAW;AAClC,aAAK,QAAQ,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,OAAO,SAAS,aAAa,cAAc,CAAC;AAChF,eAAS,IAAI,GAAG,IAAI,KAAK,SAAS,KAAK;AACrC,cAAM,IACJ,OAAO,CAAC,EAAE,cACT,KAAK,eAAe,QAAQ,gBAAgB,OAAO,CAAC,EAAE;AACzD,cAAM,IACJ,OAAO,CAAC,EAAE,cACT,KAAK,eAAe,QAAQ,gBAAgB,OAAO,CAAC,EAAE;AACzD,YACE,QACE,OAAO,CAAC,GACR,GACA,MAAM,UAAU,KAAK,eAAe,KAAK,OAAO,CAAC,EAAE,mBACnD,OAAO,CAAC,GACR,GACA,MAAM,UAAU,KAAK,eAAe,KAAK,OAAO,CAAC,EAAE,iBAAiB;AAGtE,gBAAM,IAAI,MAAM,mCAAmC;MACvD;IACF;AACA,UAAM,UAAU,KAAK,QAAO;AAC5B,QAAI;AACF,YAAM,cAAc,OAAO,WAAW,WAAW,SAAY;AAC7D,eAAS,SAAS,OAAO,SAAS,QAAQ,SAAS,UAAU;AAC3D,cAAM,SAAS,KAAK,cAAc,SAAS,aACzC,SAAS,KAAK,eAAe,SAAS;AACxC,iBAAS,IAAI,GAAG,IAAI,KAAK,WAAW;AAAK,kBAAQ,MAAM,CAAC,IAAI,MAAM,CAAC,EAAE,MAAM;AAC3E,aAAK,UACH,SACA,SACA,cAAc,YAAY,MAAM,IAAK,MAA6B;AAEpE,iBAAS,IAAI,GAAG,IAAI,KAAK,WAAW,KAAK;AACvC,gBAAM,QAAQ,QAAQ,OAAO,CAAC;AAC9B,cACG,IAAI,KAAK,CAAC,OAAO,SAAS,KAAK,KAC/B,OAAO,CAAC,aAAa,gBACpB,OAAO,SAAS,KAAK,KACrB,KAAK,IAAI,KAAK,IAAI;AAEpB,kBAAM,IAAI,MAAM,yDAAyD;QAC7E;AACA,iBAAS,IAAI,GAAG,IAAI,KAAK,WAAW;AAAK,iBAAO,CAAC,EAAE,MAAM,IAAI,QAAQ,OAAO,CAAC;MAC/E;AACA,aAAO;IACT;AACE,WAAK;IACP;EACF;;",
+  "names": []
+}

package/dist/projection/bulk.d.ts

@@ -0,0 +1,51 @@
+import type { ProjectionCoordinate, ProjectionOutput } from "./experimental/scalar-output.js";
+import type { ProjectionArray } from "./experimental/typescript-projection.js";
+import type { PipelineEpochs } from "./experimental/projection-pipeline.js";
+/** ProjectionTransform, Projection or ProjectionPipeline with synchronous reusable outputs.
+ * Lazy projections must have been preloaded. Epochs are interpreted by the supplied transform.
+ */
+export type BulkProjection = {
+    projectToSync<T extends ProjectionOutput>(coordinate: ProjectionCoordinate, result: T, epoch?: number): T;
+    unprojectToSync<T extends ProjectionOutput>(coordinate: ProjectionCoordinate, result: T, epoch?: number): T;
+};
+export type ProjectionBufferOptions = {
+    projection: BulkProjection;
+    dimension?: number;
+    /** Element offsets relative to the supplied typed-array views. Default zero. */
+    inputOffset?: number;
+    outputOffset?: number;
+    /** Elements per record; defaults to dimension for interleaved, one for columns. */
+    inputStride?: number;
+    outputStride?: number;
+};
+/** Capture a layout once, reuse scratch and output buffers for each batch/chunk.
+ * Validation errors write nothing. Coordinate failures preserve the failed record and
+ * commit preceding records. Input/output storage must be disjoint or exactly mapped
+ * in place; shifted overlap is rejected before transformation. All offsets/counts
+ * refer to records/elements of the supplied views, never to the backing buffer.
+ */
+export declare class ProjectionBuffer {
+    readonly dimension: number;
+    private readonly projection;
+    private readonly inputOffset;
+    private readonly outputOffset;
+    private readonly inputStride?;
+    private readonly outputStride?;
+    private readonly scratch;
+    private depth;
+    constructor(options: ProjectionBufferOptions);
+    projectFlatTo<T extends ProjectionArray>(input: ProjectionArray, output: T, count?: number, start?: number, epochs?: PipelineEpochs): T;
+    unprojectFlatTo<T extends ProjectionArray>(input: ProjectionArray, output: T, count?: number, start?: number, epochs?: PipelineEpochs): T;
+    projectColumnsTo<T extends readonly ProjectionArray[]>(input: readonly ProjectionArray[], output: T, count?: number, start?: number, epochs?: PipelineEpochs): T;
+    unprojectColumnsTo<T extends readonly ProjectionArray[]>(input: readonly ProjectionArray[], output: T, count?: number, start?: number, epochs?: PipelineEpochs): T;
+    private range;
+    private epochStorage;
+    private aliases;
+    private epochAlias;
+    private acquire;
+    private transform;
+    private floatRange;
+    private flat;
+    private columns;
+}
+//# sourceMappingURL=bulk.d.ts.map

package/dist/projection/bulk.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"bulk.d.ts","sourceRoot":"","sources":["../../src/projection/bulk.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,oBAAoB,EAAE,gBAAgB,EAAC,wCAAqC;AACzF,OAAO,KAAK,EAAC,eAAe,EAAC,gDAA6C;AAC1E,OAAO,KAAK,EAAC,cAAc,EAAC,8CAA2C;AAEvE;;GAEG;AACH,MAAM,MAAM,cAAc,GAAG;IAC3B,aAAa,CAAC,CAAC,SAAS,gBAAgB,EACtC,UAAU,EAAE,oBAAoB,EAChC,MAAM,EAAE,CAAC,EACT,KAAK,CAAC,EAAE,MAAM,GACb,CAAC,CAAC;IACL,eAAe,CAAC,CAAC,SAAS,gBAAgB,EACxC,UAAU,EAAE,oBAAoB,EAChC,MAAM,EAAE,CAAC,EACT,KAAK,CAAC,EAAE,MAAM,GACb,CAAC,CAAC;CACN,CAAC;AACF,MAAM,MAAM,uBAAuB,GAAG;IACpC,UAAU,EAAE,cAAc,CAAC;IAC3B,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,gFAAgF;IAChF,WAAW,CAAC,EAAE,MAAM,CAAC;IACrB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,mFAAmF;IACnF,WAAW,CAAC,EAAE,MAAM,CAAC;IACrB,YAAY,CAAC,EAAE,MAAM,CAAC;CACvB,CAAC;AA6BF;;;;;GAKG;AACH,qBAAa,gBAAgB;IAC3B,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,OAAO,CAAC,QAAQ,CAAC,UAAU,CAAiB;IAC5C,OAAO,CAAC,QAAQ,CAAC,WAAW,CAAS;IACrC,OAAO,CAAC,QAAQ,CAAC,YAAY,CAAS;IACtC,OAAO,CAAC,QAAQ,CAAC,WAAW,CAAC,CAAS;IACtC,OAAO,CAAC,QAAQ,CAAC,YAAY,CAAC,CAAS;IACvC,OAAO,CAAC,QAAQ,CAAC,OAAO,CAAY;IACpC,OAAO,CAAC,KAAK,CAAK;gBACN,OAAO,EAAE,uBAAuB;IA4B5C,aAAa,CAAC,CAAC,SAAS,eAAe,EACrC,KAAK,EAAE,eAAe,EACtB,MAAM,EAAE,CAAC,EACT,KAAK,CAAC,EAAE,MAAM,EACd,KAAK,SAAI,EACT,MAAM,CAAC,EAAE,cAAc,GACtB,CAAC;IAGJ,eAAe,CAAC,CAAC,SAAS,eAAe,EACvC,KAAK,EAAE,eAAe,EACtB,MAAM,EAAE,CAAC,EACT,KAAK,CAAC,EAAE,MAAM,EACd,KAAK,SAAI,EACT,MAAM,CAAC,EAAE,cAAc,GACtB,CAAC;IAGJ,gBAAgB,CAAC,CAAC,SAAS,SAAS,eAAe,EAAE,EACnD,KAAK,EAAE,SAAS,eAAe,EAAE,EACjC,MAAM,EAAE,CAAC,EACT,KAAK,CAAC,EAAE,MAAM,EACd,KAAK,SAAI,EACT,MAAM,CAAC,EAAE,cAAc,GACtB,CAAC;IAGJ,kBAAkB,CAAC,CAAC,SAAS,SAAS,eAAe,EAAE,EACrD,KAAK,EAAE,SAAS,eAAe,EAAE,EACjC,MAAM,EAAE,CAAC,EACT,KAAK,CAAC,EAAE,MAAM,EACd,KAAK,SAAI,EACT,MAAM,CAAC,EAAE,cAAc,GACtB,CAAC;IAGJ,OAAO,CAAC,KAAK;IAab,OAAO,CAAC,YAAY;IAOpB,OAAO,CAAC,OAAO;IAwBf,OAAO,CAAC,UAAU;IAgBlB,OAAO,CAAC,OAAO;IAKf,OAAO,CAAC,SAAS;IAMjB,OAAO,CAAC,UAAU;IAWlB,OAAO,CAAC,IAAI;IA6CZ,OAAO,CAAC,OAAO;CA2EhB"}