@thi.ng/tensors 0.9.0 → 0.10.4

package/README.md

@@ -250,7 +250,7 @@ For Node.js REPL:
 const ten = await import("@thi.ng/tensors");
 ```
 
-Package sizes (brotli'd, pre-treeshake): ESM: 10.06 KB
+Package sizes (brotli'd, pre-treeshake): ESM: 11.49 KB
 
 ## Dependencies
 

package/index.d.ts

@@ -50,6 +50,7 @@ export * from "./rand-distrib.js";
 export * from "./range.js";
 export * from "./relu.js";
 export * from "./relun.js";
+export * from "./sample.js";
 export * from "./select.js";
 export * from "./set.js";
 export * from "./setn.js";

package/index.js

@@ -50,6 +50,7 @@ export * from "./rand-distrib.js";
 export * from "./range.js";
 export * from "./relu.js";
 export * from "./relun.js";
+export * from "./sample.js";
 export * from "./select.js";
 export * from "./set.js";
 export * from "./setn.js";

package/package.json

@@ -1,343 +1,346 @@
 {
-  "name": "@thi.ng/tensors",
-  "version": "0.9.0",
-  "description": "0D/1D/2D/3D/4D tensors with extensible polymorphic operations and customizable storage",
-  "type": "module",
-  "module": "./index.js",
-  "typings": "./index.d.ts",
-  "sideEffects": false,
-  "repository": {
-    "type": "git",
-    "url": "https://github.com/thi-ng/umbrella.git"
-  },
-  "homepage": "https://thi.ng/tensors",
-  "funding": [
-    {
-      "type": "github",
-      "url": "https://github.com/sponsors/postspectacular"
-    },
-    {
-      "type": "patreon",
-      "url": "https://patreon.com/thing_umbrella"
-    },
-    {
-      "type": "liberapay",
-      "url": "https://liberapay.com/thi.ng"
-    }
-  ],
-  "author": "Karsten Schmidt (https://thi.ng)",
-  "license": "Apache-2.0",
-  "scripts": {
-    "build": "yarn build:esbuild && yarn build:decl",
-    "build:decl": "tsc --declaration --emitDeclarationOnly",
-    "build:esbuild": "esbuild --format=esm --platform=neutral --target=es2022 --tsconfig=tsconfig.json --outdir=. src/**/*.ts",
-    "clean": "bun ../../tools/src/clean-package.ts",
-    "doc": "typedoc --options ../../typedoc.json --out doc src/index.ts",
-    "doc:readme": "bun ../../tools/src/module-stats.ts && bun ../../tools/src/readme.ts",
-    "pub": "yarn npm publish --access public",
-    "test": "bun test",
-    "tool:tangle": "../../node_modules/.bin/tangle src/**/*.ts"
-  },
-  "dependencies": {
-    "@thi.ng/api": "^8.12.1",
-    "@thi.ng/arrays": "^2.13.10",
-    "@thi.ng/checks": "^3.7.17",
-    "@thi.ng/equiv": "^2.1.91",
-    "@thi.ng/errors": "^2.5.41",
-    "@thi.ng/math": "^5.11.35",
-    "@thi.ng/random": "^4.1.26",
-    "@thi.ng/strings": "^3.9.21",
-    "@thi.ng/vectors": "^8.6.5"
-  },
-  "devDependencies": {
-    "esbuild": "^0.25.8",
-    "typedoc": "^0.28.9",
-    "typescript": "^5.9.2"
-  },
-  "keywords": [
-    "1d",
-    "2d",
-    "3d",
-    "4d",
-    "activation",
-    "algebra",
-    "array",
-    "blur",
-    "broadcast",
-    "conversion",
-    "convolution",
-    "data-oriented",
-    "datastructure",
-    "decomposition",
-    "detect",
-    "diagonal",
-    "dotproduct",
-    "edge",
-    "equality",
-    "filter",
-    "gaussian",
-    "gradient",
-    "interval",
-    "math",
-    "matrix",
-    "memory-mapped",
-    "nd",
-    "polymorphic",
-    "pool",
-    "presets",
-    "random",
-    "smoothstep",
-    "step",
-    "svd",
-    "tensor",
-    "typedarray",
-    "typescript",
-    "vector"
-  ],
-  "publishConfig": {
-    "access": "public"
-  },
-  "browser": {
-    "process": false,
-    "setTimeout": false
-  },
-  "engines": {
-    "node": ">=18"
-  },
-  "files": [
-    "./*.js",
-    "./*.d.ts"
-  ],
-  "exports": {
-    ".": {
-      "default": "./index.js"
-    },
-    "./abs": {
-      "default": "./abs.js"
-    },
-    "./add": {
-      "default": "./add.js"
-    },
-    "./addn": {
-      "default": "./addn.js"
-    },
-    "./api": {
-      "default": "./api.js"
-    },
-    "./apply-kernel": {
-      "default": "./apply-kernel.js"
-    },
-    "./broadcast": {
-      "default": "./broadcast.js"
-    },
-    "./clamp": {
-      "default": "./clamp.js"
-    },
-    "./clampn": {
-      "default": "./clampn.js"
-    },
-    "./convert": {
-      "default": "./convert.js"
-    },
-    "./convolve": {
-      "default": "./convolve.js"
-    },
-    "./cos": {
-      "default": "./cos.js"
-    },
-    "./defopn": {
-      "default": "./defopn.js"
-    },
-    "./defoprt": {
-      "default": "./defoprt.js"
-    },
-    "./defoprtt": {
-      "default": "./defoprtt.js"
-    },
-    "./defopt": {
-      "default": "./defopt.js"
-    },
-    "./defoptn": {
-      "default": "./defoptn.js"
-    },
-    "./defoptnn": {
-      "default": "./defoptnn.js"
-    },
-    "./defoptt": {
-      "default": "./defoptt.js"
-    },
-    "./defopttt": {
-      "default": "./defopttt.js"
-    },
-    "./diagonal": {
-      "default": "./diagonal.js"
-    },
-    "./div": {
-      "default": "./div.js"
-    },
-    "./divn": {
-      "default": "./divn.js"
-    },
-    "./dot": {
-      "default": "./dot.js"
-    },
-    "./errors": {
-      "default": "./errors.js"
-    },
-    "./exp": {
-      "default": "./exp.js"
-    },
-    "./exp2": {
-      "default": "./exp2.js"
-    },
-    "./filtered-indices": {
-      "default": "./filtered-indices.js"
-    },
-    "./format": {
-      "default": "./format.js"
-    },
-    "./identity": {
-      "default": "./identity.js"
-    },
-    "./integrate": {
-      "default": "./integrate.js"
-    },
-    "./kernels": {
-      "default": "./kernels.js"
-    },
-    "./log": {
-      "default": "./log.js"
-    },
-    "./log2": {
-      "default": "./log2.js"
-    },
-    "./mag": {
-      "default": "./mag.js"
-    },
-    "./magsq": {
-      "default": "./magsq.js"
-    },
-    "./max": {
-      "default": "./max.js"
-    },
-    "./maxn": {
-      "default": "./maxn.js"
-    },
-    "./mean": {
-      "default": "./mean.js"
-    },
-    "./min": {
-      "default": "./min.js"
-    },
-    "./minn": {
-      "default": "./minn.js"
-    },
-    "./mul": {
-      "default": "./mul.js"
-    },
-    "./mulm": {
-      "default": "./mulm.js"
-    },
-    "./muln": {
-      "default": "./muln.js"
-    },
-    "./mulv": {
-      "default": "./mulv.js"
-    },
-    "./normalize": {
-      "default": "./normalize.js"
-    },
-    "./pow": {
-      "default": "./pow.js"
-    },
-    "./pown": {
-      "default": "./pown.js"
-    },
-    "./product": {
-      "default": "./product.js"
-    },
-    "./rand-distrib": {
-      "default": "./rand-distrib.js"
-    },
-    "./range": {
-      "default": "./range.js"
-    },
-    "./relu": {
-      "default": "./relu.js"
-    },
-    "./relun": {
-      "default": "./relun.js"
-    },
-    "./select": {
-      "default": "./select.js"
-    },
-    "./set": {
-      "default": "./set.js"
-    },
-    "./setn": {
-      "default": "./setn.js"
-    },
-    "./sigmoid": {
-      "default": "./sigmoid.js"
-    },
-    "./sin": {
-      "default": "./sin.js"
-    },
-    "./smoothstep": {
-      "default": "./smoothstep.js"
-    },
-    "./smoothstepn": {
-      "default": "./smoothstepn.js"
-    },
-    "./softmax": {
-      "default": "./softmax.js"
-    },
-    "./softplus": {
-      "default": "./softplus.js"
-    },
-    "./sqrt": {
-      "default": "./sqrt.js"
-    },
-    "./step": {
-      "default": "./step.js"
-    },
-    "./stepn": {
-      "default": "./stepn.js"
-    },
-    "./storage": {
-      "default": "./storage.js"
-    },
-    "./sub": {
-      "default": "./sub.js"
-    },
-    "./subn": {
-      "default": "./subn.js"
-    },
-    "./sum": {
-      "default": "./sum.js"
-    },
-    "./svd": {
-      "default": "./svd.js"
-    },
-    "./swap": {
-      "default": "./swap.js"
-    },
-    "./tan": {
-      "default": "./tan.js"
-    },
-    "./tanh": {
-      "default": "./tanh.js"
-    },
-    "./tensor": {
-      "default": "./tensor.js"
-    },
-    "./top": {
-      "default": "./top.js"
-    }
-  },
-  "thi.ng": {
-    "status": "alpha",
-    "year": 2018
-  },
-  "gitHead": "b133e905b546bff0d244826684269cdbb2fcfa25\n"
-}
+	"name": "@thi.ng/tensors",
+	"version": "0.10.4",
+	"description": "0D/1D/2D/3D/4D tensors with extensible polymorphic operations and customizable storage",
+	"type": "module",
+	"module": "./index.js",
+	"typings": "./index.d.ts",
+	"sideEffects": false,
+	"repository": {
+		"type": "git",
+		"url": "https://github.com/thi-ng/umbrella.git"
+	},
+	"homepage": "https://thi.ng/tensors",
+	"funding": [
+		{
+			"type": "github",
+			"url": "https://github.com/sponsors/postspectacular"
+		},
+		{
+			"type": "patreon",
+			"url": "https://patreon.com/thing_umbrella"
+		},
+		{
+			"type": "liberapay",
+			"url": "https://liberapay.com/thi.ng"
+		}
+	],
+	"author": "Karsten Schmidt (https://thi.ng)",
+	"license": "Apache-2.0",
+	"scripts": {
+		"build": "yarn build:esbuild && yarn build:decl",
+		"build:decl": "tsc --declaration --emitDeclarationOnly",
+		"build:esbuild": "esbuild --format=esm --platform=neutral --target=es2022 --tsconfig=tsconfig.json --outdir=. src/**/*.ts",
+		"clean": "bun ../../tools/src/clean-package.ts",
+		"doc": "typedoc --options ../../typedoc.json --out doc src/index.ts",
+		"doc:readme": "bun ../../tools/src/module-stats.ts && bun ../../tools/src/readme.ts",
+		"pub": "npm publish --access public",
+		"test": "bun test",
+		"tool:tangle": "../../node_modules/.bin/tangle src/**/*.ts"
+	},
+	"dependencies": {
+		"@thi.ng/api": "^8.12.6",
+		"@thi.ng/arrays": "^2.13.15",
+		"@thi.ng/checks": "^3.7.22",
+		"@thi.ng/equiv": "^2.1.96",
+		"@thi.ng/errors": "^2.5.46",
+		"@thi.ng/math": "^5.13.3",
+		"@thi.ng/random": "^4.1.31",
+		"@thi.ng/strings": "^3.9.26",
+		"@thi.ng/vectors": "^8.6.10"
+	},
+	"devDependencies": {
+		"esbuild": "^0.25.11",
+		"typedoc": "^0.28.14",
+		"typescript": "^5.9.3"
+	},
+	"keywords": [
+		"1d",
+		"2d",
+		"3d",
+		"4d",
+		"activation",
+		"algebra",
+		"array",
+		"blur",
+		"broadcast",
+		"conversion",
+		"convolution",
+		"data-oriented",
+		"datastructure",
+		"decomposition",
+		"detect",
+		"diagonal",
+		"dotproduct",
+		"edge",
+		"equality",
+		"filter",
+		"gaussian",
+		"gradient",
+		"interval",
+		"math",
+		"matrix",
+		"memory-mapped",
+		"nd",
+		"polymorphic",
+		"pool",
+		"presets",
+		"random",
+		"smoothstep",
+		"step",
+		"svd",
+		"tensor",
+		"typedarray",
+		"typescript",
+		"vector"
+	],
+	"publishConfig": {
+		"access": "public"
+	},
+	"browser": {
+		"process": false,
+		"setTimeout": false
+	},
+	"engines": {
+		"node": ">=18"
+	},
+	"files": [
+		"./*.js",
+		"./*.d.ts"
+	],
+	"exports": {
+		".": {
+			"default": "./index.js"
+		},
+		"./abs": {
+			"default": "./abs.js"
+		},
+		"./add": {
+			"default": "./add.js"
+		},
+		"./addn": {
+			"default": "./addn.js"
+		},
+		"./api": {
+			"default": "./api.js"
+		},
+		"./apply-kernel": {
+			"default": "./apply-kernel.js"
+		},
+		"./broadcast": {
+			"default": "./broadcast.js"
+		},
+		"./clamp": {
+			"default": "./clamp.js"
+		},
+		"./clampn": {
+			"default": "./clampn.js"
+		},
+		"./convert": {
+			"default": "./convert.js"
+		},
+		"./convolve": {
+			"default": "./convolve.js"
+		},
+		"./cos": {
+			"default": "./cos.js"
+		},
+		"./defopn": {
+			"default": "./defopn.js"
+		},
+		"./defoprt": {
+			"default": "./defoprt.js"
+		},
+		"./defoprtt": {
+			"default": "./defoprtt.js"
+		},
+		"./defopt": {
+			"default": "./defopt.js"
+		},
+		"./defoptn": {
+			"default": "./defoptn.js"
+		},
+		"./defoptnn": {
+			"default": "./defoptnn.js"
+		},
+		"./defoptt": {
+			"default": "./defoptt.js"
+		},
+		"./defopttt": {
+			"default": "./defopttt.js"
+		},
+		"./diagonal": {
+			"default": "./diagonal.js"
+		},
+		"./div": {
+			"default": "./div.js"
+		},
+		"./divn": {
+			"default": "./divn.js"
+		},
+		"./dot": {
+			"default": "./dot.js"
+		},
+		"./errors": {
+			"default": "./errors.js"
+		},
+		"./exp": {
+			"default": "./exp.js"
+		},
+		"./exp2": {
+			"default": "./exp2.js"
+		},
+		"./filtered-indices": {
+			"default": "./filtered-indices.js"
+		},
+		"./format": {
+			"default": "./format.js"
+		},
+		"./identity": {
+			"default": "./identity.js"
+		},
+		"./integrate": {
+			"default": "./integrate.js"
+		},
+		"./kernels": {
+			"default": "./kernels.js"
+		},
+		"./log": {
+			"default": "./log.js"
+		},
+		"./log2": {
+			"default": "./log2.js"
+		},
+		"./mag": {
+			"default": "./mag.js"
+		},
+		"./magsq": {
+			"default": "./magsq.js"
+		},
+		"./max": {
+			"default": "./max.js"
+		},
+		"./maxn": {
+			"default": "./maxn.js"
+		},
+		"./mean": {
+			"default": "./mean.js"
+		},
+		"./min": {
+			"default": "./min.js"
+		},
+		"./minn": {
+			"default": "./minn.js"
+		},
+		"./mul": {
+			"default": "./mul.js"
+		},
+		"./mulm": {
+			"default": "./mulm.js"
+		},
+		"./muln": {
+			"default": "./muln.js"
+		},
+		"./mulv": {
+			"default": "./mulv.js"
+		},
+		"./normalize": {
+			"default": "./normalize.js"
+		},
+		"./pow": {
+			"default": "./pow.js"
+		},
+		"./pown": {
+			"default": "./pown.js"
+		},
+		"./product": {
+			"default": "./product.js"
+		},
+		"./rand-distrib": {
+			"default": "./rand-distrib.js"
+		},
+		"./range": {
+			"default": "./range.js"
+		},
+		"./relu": {
+			"default": "./relu.js"
+		},
+		"./relun": {
+			"default": "./relun.js"
+		},
+		"./sample": {
+			"default": "./sample.js"
+		},
+		"./select": {
+			"default": "./select.js"
+		},
+		"./set": {
+			"default": "./set.js"
+		},
+		"./setn": {
+			"default": "./setn.js"
+		},
+		"./sigmoid": {
+			"default": "./sigmoid.js"
+		},
+		"./sin": {
+			"default": "./sin.js"
+		},
+		"./smoothstep": {
+			"default": "./smoothstep.js"
+		},
+		"./smoothstepn": {
+			"default": "./smoothstepn.js"
+		},
+		"./softmax": {
+			"default": "./softmax.js"
+		},
+		"./softplus": {
+			"default": "./softplus.js"
+		},
+		"./sqrt": {
+			"default": "./sqrt.js"
+		},
+		"./step": {
+			"default": "./step.js"
+		},
+		"./stepn": {
+			"default": "./stepn.js"
+		},
+		"./storage": {
+			"default": "./storage.js"
+		},
+		"./sub": {
+			"default": "./sub.js"
+		},
+		"./subn": {
+			"default": "./subn.js"
+		},
+		"./sum": {
+			"default": "./sum.js"
+		},
+		"./svd": {
+			"default": "./svd.js"
+		},
+		"./swap": {
+			"default": "./swap.js"
+		},
+		"./tan": {
+			"default": "./tan.js"
+		},
+		"./tanh": {
+			"default": "./tanh.js"
+		},
+		"./tensor": {
+			"default": "./tensor.js"
+		},
+		"./top": {
+			"default": "./top.js"
+		}
+	},
+	"thi.ng": {
+		"status": "alpha",
+		"year": 2018
+	},
+	"gitHead": "136a5e5ef0b69e82329db00d806c3c4e8f1aa063\n"
+}

package/sample.d.ts

@@ -0,0 +1,27 @@
+import type { Fn2, Fn3, Fn4, Fn5, FnN, FnN2 } from "@thi.ng/api";
+import type { ITensor } from "./api.js";
+import { type Tensor1 } from "./tensor.js";
+type BoundaryFn = FnN2;
+type Sampler1 = Fn2<ITensor, number, number>;
+type Sampler2 = Fn3<ITensor, number, number, number>;
+type Sampler3 = Fn4<ITensor, number, number, number, number>;
+type Sampler4 = Fn5<ITensor, number, number, number, number, number>;
+type BoundaryType = "clamp" | "mirror" | "wrap" | "zero";
+type SamplerType = "nearest" | "linear" | "cubic" | "lanczos";
+interface SamplerKernel {
+    radius: number;
+    weight: FnN;
+}
+export declare const SAMPLE_NEAREST: SamplerKernel;
+export declare const SAMPLE_LINEAR: SamplerKernel;
+export declare const SAMPLE_CUBIC: SamplerKernel;
+export declare const SAMPLE_LANCZOS: (r?: number) => SamplerKernel;
+export declare const defSampler1: (kernel: SamplerKernel | SamplerType, boundary?: BoundaryFn | BoundaryType) => Sampler1;
+export declare const defSampler2: (kernel: SamplerType | SamplerKernel | [SamplerType | SamplerKernel, SamplerType | SamplerKernel], boundary?: BoundaryType | BoundaryFn | [BoundaryType | BoundaryFn, BoundaryType | BoundaryFn]) => Sampler2;
+export declare const defSampler3: (kernel: SamplerType | SamplerKernel | [SamplerType | SamplerKernel, SamplerType | SamplerKernel, SamplerType | SamplerKernel], boundary?: BoundaryType | BoundaryFn | [BoundaryType | BoundaryFn, BoundaryType | BoundaryFn, BoundaryType | BoundaryFn]) => Sampler3;
+export declare const defSampler4: (kernel: SamplerType | SamplerKernel | [SamplerType | SamplerKernel, SamplerType | SamplerKernel, SamplerType | SamplerKernel, SamplerType | SamplerKernel], boundary?: BoundaryType | BoundaryFn | [BoundaryType | BoundaryFn, BoundaryType | BoundaryFn, BoundaryType | BoundaryFn, BoundaryType | BoundaryFn]) => Sampler4;
+export declare const resample1: (out: Tensor1, a: Tensor1, sampler: Sampler1) => Tensor1<number>;
+export declare const resample2: (out: ITensor, a: ITensor, [samplerX, samplerY]: Sampler1[]) => ITensor<number>;
+export declare const resample3: (out: ITensor, a: ITensor, [samplerX, samplerY, samplerZ]: Sampler1[]) => ITensor<number>;
+export {};
+//# sourceMappingURL=sample.d.ts.map

package/sample.js

@@ -0,0 +1,290 @@
+import { isArray } from "@thi.ng/checks/is-array";
+import { isString } from "@thi.ng/checks/is-string";
+import { PI } from "@thi.ng/math/api";
+import { clamp } from "@thi.ng/math/interval";
+import { tensor } from "./tensor.js";
+const { abs, floor, sin } = Math;
+const __boundaryClamp = (i, max) => clamp(i, 0, max - 1);
+const __boundaryWrap = (i, max) => (i % max + max) % max;
+const __boundaryMirror = (i, max) => {
+  if (max === 1) return 0;
+  const period = 2 * (max - 1);
+  i = (i % period + period) % period;
+  return i < max ? i : period - i;
+};
+const __boundaryZero = (i, max) => i >= 0 || i < max ? i : -1;
+const BOUNDARY_MODES = {
+  clamp: __boundaryClamp,
+  wrap: __boundaryWrap,
+  mirror: __boundaryMirror,
+  zero: __boundaryZero
+};
+const SAMPLE_NEAREST = {
+  radius: 1,
+  weight: (t) => t >= -0.5 && t < 0.5 ? 1 : 0
+};
+const SAMPLE_LINEAR = {
+  radius: 1,
+  weight: (t) => {
+    t = abs(t);
+    return t < 1 ? 1 - t : 0;
+  }
+};
+const SAMPLE_CUBIC = {
+  radius: 2,
+  weight: (t) => {
+    t = abs(t);
+    const a = -0.5;
+    return t < 2 ? t <= 1 ? ((a + 2) * t - (a + 3)) * t * t + 1 : ((a * t - 5 * a) * t + 8 * a) * t - 4 * a : 0;
+  }
+};
+const SAMPLE_LANCZOS = (r = 2) => ({
+  radius: r,
+  weight: (t) => {
+    if (t === 0) return 1;
+    if (t > -r && t < r) {
+      t *= PI;
+      return r * sin(t) * sin(t / r) / (t * t);
+    }
+    return 0;
+  }
+});
+const SAMPLER_TYPES = {
+  nearest: SAMPLE_NEAREST,
+  linear: SAMPLE_LINEAR,
+  cubic: SAMPLE_CUBIC,
+  lanczos: SAMPLE_LANCZOS(2)
+};
+const defSampler1 = (kernel, boundary = "clamp") => {
+  const index = __resolveBoundary(boundary);
+  const { radius, weight } = __resolveKernel(kernel);
+  return ({ data, shape: [sx], stride: [tx], offset }, x) => {
+    const c = floor(x);
+    let sum = 0;
+    let wsum = 0;
+    for (let k = -radius; k <= radius; k++) {
+      const $x = c + k;
+      const w = weight(x - $x);
+      if (w === 0) continue;
+      const i = index($x, sx);
+      if (i < 0) {
+        wsum += w;
+        continue;
+      }
+      sum += data[offset + i * tx] * w;
+      wsum += w;
+    }
+    return wsum !== 0 ? sum / wsum : 0;
+  };
+};
+const defSampler2 = (kernel, boundary = "clamp") => {
+  const [indexX, indexY] = (isArray(boundary) ? boundary : [boundary, boundary]).map(__resolveBoundary);
+  const [kernelX, kernelY] = isArray(kernel) ? kernel : [kernel, kernel];
+  const { radius: rx, weight: weightX } = __resolveKernel(kernelX);
+  const { radius: ry, weight: weightY } = __resolveKernel(kernelY);
+  const sizeY = 2 * ry + 1;
+  const cacheY = new Array(sizeY);
+  return ({ data, shape: [sx, sy], stride: [tx, ty], offset }, x, y) => {
+    const cx = floor(x);
+    const cy = floor(y);
+    let sum = 0;
+    let wsum = 0;
+    for (let i = 0; i < sizeY; i++) {
+      cacheY[i] = weightY(y - (cy + i - ry));
+    }
+    for (let kx = -rx; kx <= rx; kx++) {
+      const $x = cx + kx;
+      const ix = indexX($x, sx);
+      const wx = weightX(x - $x);
+      const idx = offset + ix * tx;
+      for (let ky = -ry; ky <= ry; ky++) {
+        const wy = cacheY[ky + ry];
+        const weight = wx * wy;
+        if (weight === 0) continue;
+        const iy = indexY(cy + ky, sy);
+        if (ix < 0 || iy < 0) {
+          wsum += weight;
+          continue;
+        }
+        sum += data[idx + iy * ty] * weight;
+        wsum += weight;
+      }
+    }
+    return wsum !== 0 ? sum / wsum : 0;
+  };
+};
+const defSampler3 = (kernel, boundary = "clamp") => {
+  const [indexX, indexY, indexZ] = (isArray(boundary) ? boundary : [boundary, boundary, boundary]).map(__resolveBoundary);
+  const [kernelX, kernelY, kernelZ] = (isArray(kernel) ? kernel : [kernel, kernel, kernel]).map(__resolveKernel);
+  const { radius: rx, weight: weightX } = kernelX;
+  const { radius: ry, weight: weightY } = kernelY;
+  const { radius: rz, weight: weightZ } = kernelZ;
+  const sizeY = 2 * ry + 1;
+  const sizeZ = 2 * rz + 1;
+  const cacheY = new Array(sizeY);
+  const cacheZ = new Array(sizeZ);
+  return ({ data, shape: [sx, sy, sz], stride: [tx, ty, tz], offset }, x, y, z) => {
+    const cx = floor(x);
+    const cy = floor(y);
+    const cz = floor(z);
+    let sum = 0;
+    let wsum = 0;
+    for (let i = 0; i < sizeY; i++) cacheY[i] = weightY(y - (cy + i - ry));
+    for (let i = 0; i < sizeZ; i++) cacheZ[i] = weightZ(z - (cz + i - rz));
+    for (let kx = -rx; kx <= rx; kx++) {
+      const $x = cx + kx;
+      const ix = indexX($x, sx);
+      const wx = weightX(x - $x);
+      const idxX = offset + ix * tx;
+      for (let ky = -ry; ky <= ry; ky++) {
+        const iy = indexY(cy + ky, sy);
+        const wy = cacheY[ky + ry];
+        const idxXY = idxX + iy * ty;
+        for (let kz = -rz; kz <= rz; kz++) {
+          const wz = cacheZ[kz + rz];
+          const weight = wx * wy * wz;
+          if (weight === 0) continue;
+          const iz = indexZ(cz + kz, sz);
+          if (ix < 0 || iy < 0 || iz < 0) {
+            wsum += weight;
+            continue;
+          }
+          sum += data[idxXY + iz * tz] * weight;
+          wsum += weight;
+        }
+      }
+    }
+    return wsum !== 0 ? sum / wsum : 0;
+  };
+};
+const defSampler4 = (kernel, boundary = "clamp") => {
+  const [indexX, indexY, indexZ, indexW] = (isArray(boundary) ? boundary : [boundary, boundary, boundary, boundary]).map(__resolveBoundary);
+  const [kernelX, kernelY, kernelZ, kernelW] = (isArray(kernel) ? kernel : [kernel, kernel, kernel, kernel]).map(__resolveKernel);
+  const { radius: rx, weight: weightX } = kernelX;
+  const { radius: ry, weight: weightY } = kernelY;
+  const { radius: rz, weight: weightZ } = kernelZ;
+  const { radius: rw, weight: weightW } = kernelW;
+  const sizeY = 2 * ry + 1;
+  const sizeZ = 2 * rz + 1;
+  const sizeW = 2 * rw + 1;
+  const cacheY = new Array(sizeY);
+  const cacheZ = new Array(sizeZ);
+  const cacheW = new Array(sizeW);
+  return ({ data, shape: [sx, sy, sz, sw], stride: [tx, ty, tz, tw], offset }, x, y, z, w) => {
+    const cx = floor(x);
+    const cy = floor(y);
+    const cz = floor(z);
+    const cw = floor(w);
+    let sum = 0;
+    let wsum = 0;
+    for (let i = 0; i < sizeY; i++) cacheY[i] = weightY(y - (cy + i - ry));
+    for (let i = 0; i < sizeZ; i++) cacheZ[i] = weightZ(z - (cz + i - rz));
+    for (let i = 0; i < sizeW; i++) cacheW[i] = weightW(w - (cw + i - rw));
+    for (let kx = -rx; kx <= rx; kx++) {
+      const $x = cx + kx;
+      const ix = indexX($x, sx);
+      const wx = weightX(x - $x);
+      const idxX = offset + ix * tx;
+      for (let ky = -ry; ky <= ry; ky++) {
+        const iy = indexY(cy + ky, sy);
+        const wy = cacheY[ky + ry];
+        const idxXY = idxX + iy * ty;
+        for (let kz = -rz; kz <= rz; kz++) {
+          const wz = cacheZ[kz + rz];
+          const iz = indexZ(cz + kz, sz);
+          const idxXYZ = idxXY + iz * tz;
+          for (let kw = -rw; kw <= rw; kw++) {
+            const ww = cacheW[kw + rw];
+            const weight = wx * wy * wz * ww;
+            if (weight === 0) continue;
+            const iw = indexW(cw + kw, sw);
+            if (ix < 0 || iy < 0 || iz < 0 || iw < 0) {
+              wsum += weight;
+              continue;
+            }
+            sum += data[idxXYZ + iw * tw] * weight;
+            wsum += w;
+          }
+        }
+      }
+    }
+    return wsum !== 0 ? sum / wsum : 0;
+  };
+};
+const __resolveBoundary = (mode) => isString(mode) ? BOUNDARY_MODES[mode] : mode ?? __boundaryClamp;
+const __resolveKernel = (kernel) => isString(kernel) ? SAMPLER_TYPES[kernel] : kernel;
+const resample1 = (out, a, sampler) => {
+  const {
+    data: odata,
+    shape: [sxo],
+    stride: [txo],
+    offset: oo
+  } = out;
+  const {
+    shape: [sxa]
+  } = a;
+  const scale = sxa / sxo;
+  for (let i = 0; i < sxo; i++) {
+    odata[oo + i * txo] = sampler(a, (i + 0.5) * scale - 0.5);
+  }
+  return out;
+};
+const resample2 = (out, a, [samplerX, samplerY]) => {
+  const {
+    shape: [_, syo]
+  } = out;
+  const {
+    shape: [sxa]
+  } = a;
+  const tmp = tensor("num", [sxa, syo]);
+  for (let x = 0; x < sxa; x++) {
+    resample1(tmp.pick([x, -1]), a.pick([x, -1]), samplerX);
+  }
+  for (let y = 0; y < syo; y++) {
+    resample1(out.pick([-1, y]), tmp.pick([-1, y]), samplerY);
+  }
+  return out;
+};
+const resample3 = (out, a, [samplerX, samplerY, samplerZ]) => {
+  const {
+    shape: [_, syo, szo]
+  } = out;
+  const {
+    shape: [sxa, sya]
+  } = a;
+  const tmpX = tensor("num", [sxa, sya, szo]);
+  for (let i = 0; i < sxa; i++) {
+    const src = a.pick([i]);
+    const dest = tmpX.pick([i]);
+    for (let j = 0; j < sya; j++) {
+      resample1(dest.pick([j]), src.pick([j]), samplerX);
+    }
+  }
+  const tmpY = tensor("num", [sxa, syo, szo]);
+  for (let i = 0; i < sxa; i++) {
+    const src = tmpX.pick([i]);
+    const dest = tmpY.pick([i]);
+    for (let j = 0; j < szo; j++) {
+      resample1(dest.pick([-1, j]), src.pick([-1, j]), samplerY);
+    }
+  }
+  for (let i = 0; i < syo; i++) {
+    for (let j = 0; j < szo; j++) {
+      resample1(out.pick([-1, i, j]), tmpY.pick([-1, i, j]), samplerZ);
+    }
+  }
+  return out;
+};
+export {
+  SAMPLE_CUBIC,
+  SAMPLE_LANCZOS,
+  SAMPLE_LINEAR,
+  SAMPLE_NEAREST,
+  defSampler1,
+  defSampler2,
+  defSampler3,
+  defSampler4,
+  resample1,
+  resample2,
+  resample3
+};

package/CHANGELOG.md

@@ -1,146 +0,0 @@
-# Change Log
-
-- **Last updated**: 2025-08-16T18:12:02Z
-- **Generator**: [thi.ng/monopub](https://thi.ng/monopub)
-
-All notable changes to this project will be documented in this file.
-Only versions published since **2022-01-01** are listed here.
-Please consult the Git history for older version information.
-See [Conventional Commits](https://conventionalcommits.org/) for commit guidelines.
-
-**Note:** Unlisted _patch_ versions only involve non-code or otherwise excluded changes
-and/or version bumps of transitive dependencies.
-
-## [0.9.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.9.0) (2025-08-16)
-
-#### 🚀 Features
-
-- add 0-dim (wrapped scalar) tensor support ([cd22ee6](https://github.com/thi-ng/umbrella/commit/cd22ee6))
-  - add `Tensor0`, `Shape0`
-  - update HOF tensor op generators
-  - add/update tests
-
-## [0.8.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.8.0) (2025-07-10)
-
-#### 🚀 Features
-
-- add `step()`/`smoothStep()` variations ([b433462](https://github.com/thi-ng/umbrella/commit/b433462))
-- add ITensor.crop() ([77e294d](https://github.com/thi-ng/umbrella/commit/77e294d))
-
-## [0.7.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.7.0) (2025-07-09)
-
-#### 🚀 Features
-
-- add `integrate()` tensor op ([fe39a4a](https://github.com/thi-ng/umbrella/commit/fe39a4a))
-- add `TensorLike` and `asTensor()` ([f76d9a6](https://github.com/thi-ng/umbrella/commit/f76d9a6))
-
-#### ♻️ Refactoring
-
-- migrate `asTensor()` to convert.ts ([fcb83c5](https://github.com/thi-ng/umbrella/commit/fcb83c5))
-
-## [0.6.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.6.0) (2025-06-05)
-
-#### 🚀 Features
-
-- update defOpRT()/defOpRTT() ([8099abc](https://github.com/thi-ng/umbrella/commit/8099abc))
-  - add optional `complete` fn to produce/transform final result (reducer like)
-- add mean() tensor op ([2f15d21](https://github.com/thi-ng/umbrella/commit/2f15d21))
-
-## [0.5.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.5.0) (2025-05-28)
-
-#### 🚀 Features
-
-- add convolve() (1D/2D/3D versions) ([bc9a6eb](https://github.com/thi-ng/umbrella/commit/bc9a6eb))
-- add `constant()`, refactor `ones()` ([8cdd68c](https://github.com/thi-ng/umbrella/commit/8cdd68c))
-- add convolution kernel presets ([e553f06](https://github.com/thi-ng/umbrella/commit/e553f06))
-  - add 1D/2D/3D sobel kernels
-  - add 2D box blur
-  - add 2D gaussian blur
-- add/update/rename convolution kernels ([025c7ba](https://github.com/thi-ng/umbrella/commit/025c7ba))
-- add `applyKernel()`, add/update kernel presets ([b8c2a7a](https://github.com/thi-ng/umbrella/commit/b8c2a7a))
-- add `fromFloatBuffer()` conversion ([d78a0d3](https://github.com/thi-ng/umbrella/commit/d78a0d3))
-
-## [0.4.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.4.0) (2025-05-12)
-
-#### 🚀 Features
-
-- add release() handling for tensors ([da06312](https://github.com/thi-ng/umbrella/commit/da06312))
-- add data type support for range() options ([90d02bc](https://github.com/thi-ng/umbrella/commit/90d02bc))
-
-## [0.3.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.3.0) (2025-05-08)
-
-#### 🚀 Features
-
-- add swap(), update defOpRT/TT() fns ([5d2dd2a](https://github.com/thi-ng/umbrella/commit/5d2dd2a))
-- add range() tensor factory ([a8a6365](https://github.com/thi-ng/umbrella/commit/a8a6365))
-- add ITensor.position() & impls ([595764d](https://github.com/thi-ng/umbrella/commit/595764d))
-- add filteredIndices() & presets ([edd8983](https://github.com/thi-ng/umbrella/commit/edd8983))
-  - add nonZeroIndices()
-  - add negativeIndices(), positiveIndices()
-
-#### 🩹 Bug fixes
-
-- update ITensor.broadcast() return type ([f7b2e7a](https://github.com/thi-ng/umbrella/commit/f7b2e7a))
-
-#### ♻️ Refactoring
-
-- update range() impl (w/o dependencies) ([ce89d71](https://github.com/thi-ng/umbrella/commit/ce89d71))
-
-## [0.2.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.2.0) (2025-05-02)
-
-#### 🚀 Features
-
-- add `zeroes()`/`ones()` ([eb6f82d](https://github.com/thi-ng/umbrella/commit/eb6f82d))
-- add singular value decomp `svd()` ([8902157](https://github.com/thi-ng/umbrella/commit/8902157))
-- add `diagonal()` and `trace()` fns ([47a0e73](https://github.com/thi-ng/umbrella/commit/47a0e73))
-- add broadcast(), add ITensor.broadcast() ([b9b2dfc](https://github.com/thi-ng/umbrella/commit/b9b2dfc))
-- update all tensor ops & op generators ([47b3e4e](https://github.com/thi-ng/umbrella/commit/47b3e4e))
-  - add broadcasting support where possible
-  - update all `defOpXXX` HOF generators
-  - replace & remove various `TensorOpXXX` function types
-  - update all tensor ops, only keep one version per op
-  - update doc strings
-
-#### 🩹 Bug fixes
-
-- add missing `randDistrib4()` ([eb2ded2](https://github.com/thi-ng/umbrella/commit/eb2ded2))
-
-#### ♻️ Refactoring
-
-- update internal local var handling in various tensor ops ([c02c2a6](https://github.com/thi-ng/umbrella/commit/c02c2a6))
-- internal updates dot(), identity(), select() ([c0ed334](https://github.com/thi-ng/umbrella/commit/c0ed334))
-- update top() return type ([dcab223](https://github.com/thi-ng/umbrella/commit/dcab223))
-
-## [0.1.0](https://github.com/thi-ng/umbrella/tree/@thi.ng/tensors@0.1.0) (2025-04-30)
-
-#### 🚀 Features
-
-- initial import (old code, needs update) ([dc4f2b4](https://github.com/thi-ng/umbrella/commit/dc4f2b4))
-- major updates & refactoring/cleanup ([49249a9](https://github.com/thi-ng/umbrella/commit/49249a9))
-- update & restructure defOpXX() fns ([a96ecf8](https://github.com/thi-ng/umbrella/commit/a96ecf8))
-- update TensorOpts and tensor factory fn, add docs ([21b8d48](https://github.com/thi-ng/umbrella/commit/21b8d48))
-- update/replace tensor op generators, add op types ([2a00949](https://github.com/thi-ng/umbrella/commit/2a00949))
-- add first set of common math ops ([00c02c1](https://github.com/thi-ng/umbrella/commit/00c02c1))
-- add matrix-vector multiply, add custom error type ([a005276](https://github.com/thi-ng/umbrella/commit/a005276))
-- add `identity()` factory fn ([1f2e914](https://github.com/thi-ng/umbrella/commit/1f2e914))
-- add defOpRT(), add product/sum reducers ([19b6a9a](https://github.com/thi-ng/umbrella/commit/19b6a9a))
-- add exp, log, softMax ops, update pkg exports ([465b0d2](https://github.com/thi-ng/umbrella/commit/465b0d2))
-- add various math ops ([9e385ad](https://github.com/thi-ng/umbrella/commit/9e385ad))
-- add/update various math ops, types, op factories ([a64468e](https://github.com/thi-ng/umbrella/commit/a64468e))
-- add select(), argMin, argMax() fns ([fa508e1](https://github.com/thi-ng/umbrella/commit/fa508e1))
-- add `print()` helper ([6c2e347](https://github.com/thi-ng/umbrella/commit/6c2e347))
-  - rename existing fn => `debug()`
-- add softPlus() op ([1d85630](https://github.com/thi-ng/umbrella/commit/1d85630))
-- add 4D tensor support & ops ([e3c9a0f](https://github.com/thi-ng/umbrella/commit/e3c9a0f))
-- add exp2(), log2() ops ([18544a4](https://github.com/thi-ng/umbrella/commit/18544a4))
-- update tensor(), add tensorFromArray(), add types/generics ([204734d](https://github.com/thi-ng/umbrella/commit/204734d))
-
-#### 🩹 Bug fixes
-
-- fix tensor scalar ops ([a26309d](https://github.com/thi-ng/umbrella/commit/a26309d))
-
-#### ♻️ Refactoring
-
-- minor update mulM() ([8bec8d1](https://github.com/thi-ng/umbrella/commit/8bec8d1))
-- deduplicate tensor impls ([5fac485](https://github.com/thi-ng/umbrella/commit/5fac485))
-- minor internal update arg types ([be56f5a](https://github.com/thi-ng/umbrella/commit/be56f5a))