sync 2e7068faf55e

Browse files

Files changed (7) hide show

README.md +57 -0
build/webgpu/bench.json +53 -0
build/webgpu/manifest.json +80 -0
build/webgpu/metadata.json +19 -0
build/webgpu/test.json +209 -0
build/webgpu/unary-scalar.wgsl.jinja +41 -0
build/webgpu/unary-vec4.wgsl.jinja +34 -0

README.md CHANGED Viewed

@@ -1,3 +1,60 @@
 ---
 license: apache-2.0
 ---

 ---
+library_name: kernels
 license: apache-2.0
+tags:
+- kernel
+- webgpu
+- wgsl
 ---
+# ai.onnx.Asinh
+`ai.onnx`  ·  standard ONNX operator  ·  ONNX opset ≥ 9
+## Description
+Computes the hyperbolic arcsine (`asinh`) of each element in the input tensor. The output has the same shape as the input.
+See the [ONNX `Asinh` spec](https://onnx.ai/onnx/operators/onnx__Asinh.html) for the reference semantics.
+## Inputs
+| Name | Bind key | Logical dtype | Rank | Shape | Description | Presence |
+| --- | --- | --- | --- | --- | --- | --- |
+| `input` | `x` | `T` | — | — | Values to map elementwise through inverse hyperbolic sine. | required |
+## Outputs
+| Name | Bind key | Logical dtype | Rank | Shape | Description | Presence |
+| --- | --- | --- | --- | --- | --- | --- |
+| `output` | `y` | `T` | same as `input` | same as `input` | The hyperbolic arcsine of each input element, same shape as the input. | required |
+## Type constraints
+| Variable | Allowed dtypes |
+| --- | --- |
+| `T` | `float32`, `float16` |
+## Files
+- [`metadata.json`](build/webgpu/metadata.json) — kernel metadata (id, digests, provenance)
+- [`manifest.json`](build/webgpu/manifest.json) — the op contract (source of truth)
+- [`test.json`](build/webgpu/test.json) — correctness cases
+- [`bench.json`](build/webgpu/bench.json) — benchmark + tuning cases
+- [`unary-scalar.wgsl.jinja`](build/webgpu/unary-scalar.wgsl.jinja)
+- [`unary-vec4.wgsl.jinja`](build/webgpu/unary-vec4.wgsl.jinja)
+## Use with `@huggingface/kernels`
+The loader derives every required output's shape and logical dtype from the manifest contract and this call.
+It then allocates the result tensors automatically.
+The `version: 1` option selects the published kernel contract; it is independent of any operator opset, contrib `since_version`, or model version.
+Replace each `*Data` placeholder with a typed array containing the corresponding input data.
+```js
+import { getKernel } from "@huggingface/kernels";
+const kernel = await getKernel("webgpu-kernels/ai.onnx.Asinh", { version: 1 });
+const { y } = await kernel({ x: { data: xData, shape: [3] } });
+```

build/webgpu/bench.json ADDED Viewed

	@@ -0,0 +1,53 @@

+{
+  "op": "ai.onnx.Asinh",
+  "cases": [
+    {
+      "name": "1m_f32",
+      "inputs": { "x": { "dtype": "float32", "shape": [1048576] } },
+      "outputs": { "y": { "dtype": "float32", "shape": [1048576] } }
+    },
+    {
+      "name": "asinh-f32-1m-vec4",
+      "preset": "smoke",
+      "vars": { "dtype": "float32", "count": 1048576 },
+      "inputs": {
+        "x": { "shape": [1048576], "dtype": "float32", "dist": "uniform", "seed": 911, "scale": 20, "offset": -10 }
+      },
+      "outputs": { "y": { "shape": [1048576], "dtype": "float32" } },
+      "bench": {
+        "primary": true,
+        "metrics": [{ "type": "bandwidth", "value": "args.count * dtypeBytes(args.dtype) * 2" }]
+      }
+    },
+    {
+      "name": "asinh-f32-1m-minus-1-scalar-fallback",
+      "preset": "smoke",
+      "vars": { "dtype": "float32", "count": 1048575 },
+      "inputs": {
+        "x": { "shape": [1048575], "dtype": "float32", "dist": "uniform", "seed": 912, "scale": 20, "offset": -10 }
+      },
+      "outputs": { "y": { "shape": [1048575], "dtype": "float32" } },
+      "bench": { "metrics": [{ "type": "bandwidth", "value": "args.count * dtypeBytes(args.dtype) * 2" }] }
+    },
+    {
+      "name": "asinh-f16-1m-vec4",
+      "preset": "smoke",
+      "vars": { "dtype": "float16", "count": 1048576 },
+      "inputs": {
+        "x": { "shape": [1048576], "dtype": "float16", "dist": "uniform", "seed": 913, "scale": 20, "offset": -10 }
+      },
+      "outputs": { "y": { "shape": [1048576], "dtype": "float16" } },
+      "bench": { "metrics": [{ "type": "bandwidth", "value": "args.count * dtypeBytes(args.dtype) * 2" }] }
+    },
+    {
+      "name": "asinh-f16-1m-minus-1-scalar-fallback",
+      "preset": "smoke",
+      "vars": { "dtype": "float16", "count": 1048575 },
+      "inputs": {
+        "x": { "shape": [1048575], "dtype": "float16", "dist": "uniform", "seed": 914, "scale": 20, "offset": -10 }
+      },
+      "outputs": { "y": { "shape": [1048575], "dtype": "float16" } },
+      "bench": { "metrics": [{ "type": "bandwidth", "value": "args.count * dtypeBytes(args.dtype) * 2" }] }
+    }
+  ]
+}

build/webgpu/manifest.json ADDED Viewed

	@@ -0,0 +1,80 @@

+{
+  "domain": "ai.onnx",
+  "name": "Asinh",
+  "sinceVersion": 9,
+  "description": "Computes the hyperbolic arcsine (`asinh`) of each element in the input tensor. The output has the same shape as the input.",
+  "inputs": [
+    { "role": "input", "dtype": "T", "description": "Values to map elementwise through inverse hyperbolic sine." }
+  ],
+  "outputs": [
+    {
+      "role": "output",
+      "dtype": "T",
+      "rank": "ranks.input",
+      "description": "The hyperbolic arcsine of each input element, same shape as the input.",
+      "shape": "shapes.input"
+    }
+  ],
+  "typeConstraints": { "T": ["float32", "float16"] },
+  "args": {
+    "x": { "kind": "tensor", "semantic": "input", "role": "input" },
+    "y": { "kind": "tensor", "semantic": "output", "role": "output" }
+  },
+  "tunables": { "WORKGROUP_SIZE": 256 },
+  "variants": [
+    {
+      "id": "same_layout_vec4",
+      "when": ["numel(shapes.x) > 0", "numel(shapes.x) % 4 == 0", "numel(shapes.x) == numel(shapes.y)", "f16Ok(dtypes.T)"],
+      "constants": {
+        "scalar": "dtypes.T",
+        "usesF16": "dtypes.T == \"f16\"",
+        "vectorScalar": "\"vec4<\" ~ dtypes.T ~ \">\""
+      },
+      "passes": [
+        {
+          "id": "main",
+          "name": "Asinh.vec4",
+          "source": { "shader": "unary-vec4.wgsl.jinja", "inputs": { "op": "\"asinh\"" } },
+          "bindings": [
+            { "name": "x", "arg": "x", "buffer": { "type": "read-only-storage" }, "elementType": "$vectorScalar" },
+            { "name": "y", "arg": "y", "buffer": { "type": "storage" }, "elementType": "$vectorScalar" },
+            {
+              "name": "params",
+              "semantic": "kernel.params",
+              "buffer": { "type": "uniform" },
+              "struct": {
+                "name": "Params",
+                "fields": [{ "name": "count", "type": "u32", "value": "numel(shapes.y) / 4" }]
+              }
+            }
+          ],
+          "dispatch": { "threads": "numel(shapes.y) / 4", "workgroupSize": "tunables.WORKGROUP_SIZE" }
+        }
+      ],
+      "priority": 20
+    },
+    {
+      "id": "elementwise",
+      "when": ["numel(shapes.x) == numel(shapes.y)", "f16Ok(dtypes.T)"],
+      "constants": { "scalar": "dtypes.T", "usesF16": "dtypes.T == \"f16\"" },
+      "passes": [
+        {
+          "id": "main",
+          "name": "Asinh",
+          "source": { "shader": "unary-scalar.wgsl.jinja", "inputs": { "op": "\"asinh\"", "itemsPerInvocation": 4 } },
+          "bindings": [
+            { "name": "x", "arg": "x", "buffer": { "type": "read-only-storage" }, "elementType": "$scalar" },
+            { "name": "y", "arg": "y", "buffer": { "type": "storage" }, "elementType": "$scalar" },
+            {
+              "name": "params",
+              "semantic": "kernel.params",
+              "buffer": { "type": "uniform" },
+              "struct": { "name": "Params", "fields": [{ "name": "count", "type": "u32", "value": "numel(shapes.y)" }] }
+            }
+          ],
+          "dispatch": { "threads": "ceilDiv(numel(shapes.y), 4)", "workgroupSize": "tunables.WORKGROUP_SIZE" }
+        }
+      ]
+    }
+  ]
+}

build/webgpu/metadata.json ADDED Viewed

	@@ -0,0 +1,19 @@

+{
+  "name": "ai.onnx.Asinh",
+  "id": "_ai_onnx_asinh_webgpu_a503438",
+  "version": 1,
+  "license": "Apache-2.0",
+  "backend": { "type": "webgpu" },
+  "digest": {
+    "algorithm": "sha256",
+    "files": {
+      "bench.json": "UW3ydjAwtJ6/U0sBoKtbFqtbwYgZfSSu7LuQs6VApTE=",
+      "manifest.json": "/3bRWF5ZfvXDNsa8jL5H5Q47aMtB735wo54xo5avPxM=",
+      "test.json": "3PxMi6OpSAIMhT+AZG7i1kw5jmD1uKj7ru2+snR8mBE=",
+      "unary-scalar.wgsl.jinja": "Zr4SHvEJ1w76n2vKN8vj5IrnQyIISQpUSzlsNriOAnQ=",
+      "unary-vec4.wgsl.jinja": "WSYxLj9wsnlDWUSeLwPHDn1fc0XFtaSLZxQgiaIxg98="
+    }
+  },
+  "provenance": { "kernel": { "sha": "2e7068faf55e7f43df740015f6d1ee49391a41c5", "dirty": false } },
+  "webgpu": { "manifestSpec": "1.0", "specialized": true, "opPath": "ops/ai.onnx.Asinh" }
+}

build/webgpu/test.json ADDED Viewed

	@@ -0,0 +1,209 @@

+{
+  "op": "ai.onnx.Asinh",
+  "cases": [
+    {
+      "name": "f32_values",
+      "inputs": {
+        "x": {
+          "dtype": "float32",
+          "shape": [6],
+          "data": { "kind": "values", "values": [-4.0, -1.0, 0.0, 0.5, 1.0, 4.0] }
+        }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [6], "tolerance": 0.000001 } }
+    },
+    {
+      "name": "f32_subnormal_identity_tail_gpu_gap",
+      "skipGpu": {
+        "category": "permanent",
+        "reason": "Portable WGSL floating-point semantics do not guarantee preservation of the subnormal values required by this fixture. Backend evidence: WebGPU/Metal flushes subnormals to zero (f32 and f16); the kernel cannot preserve denormal inputs/outputs bit-exactly."
+      },
+      "provenance": {
+        "source": "onnxruntime/test/providers/cpu/math/element_wise_ops_test.cc",
+        "test": "MathOpTest.Asinh",
+        "notes": "For tiny finite inputs asinh(x) rounds back to x in float32; zero-flushing erases the signed tail."
+      },
+      "inputs": {
+        "x": {
+          "dtype": "float32",
+          "shape": [4],
+          "data": { "kind": "values", "values": [-1e-39, -1e-40, 1e-40, 1e-39] }
+        }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [4], "tolerance": 0 } }
+    },
+    {
+      "name": "f32_subnormal_identity_tail_scalar_gpu_gap",
+      "skipGpu": {
+        "category": "permanent",
+        "reason": "Portable WGSL floating-point semantics do not guarantee preservation of the subnormal values required by this fixture. Backend evidence: WebGPU/Metal flushes subnormals to zero (f32 and f16); the kernel cannot preserve denormal inputs/outputs bit-exactly."
+      },
+      "provenance": {
+        "source": "onnxruntime/test/providers/cpu/math/element_wise_ops_test.cc",
+        "test": "MathOpTest.Asinh",
+        "notes": "Scalar-path companion: subnormal inputs are valid finite Asinh outputs."
+      },
+      "inputs": {
+        "x": { "dtype": "float32", "shape": [3], "data": { "kind": "values", "values": [-1e-40, 0.0, 1e-40] } }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [3], "tolerance": 0 } }
+    },
+    {
+      "name": "f16_values",
+      "inputs": {
+        "x": {
+          "dtype": "float16",
+          "shape": [2, 3],
+          "data": { "kind": "values", "values": [-4.0, -1.0, 0.0, 0.5, 1.0, 4.0] }
+        }
+      },
+      "outputs": { "y": { "dtype": "float16", "shape": [2, 3] } },
+      "tolerance": 0.002
+    },
+    {
+      "name": "large_in_domain_args_no_overflow",
+      "inputs": {
+        "x": {
+          "dtype": "float32",
+          "shape": [6],
+          "data": {
+            "kind": "values",
+            "values": [10000000000.0, 100000000000000000000.0, -100000000000000000000.0, 1e+30, -1e+30, 3e+38]
+          }
+        }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [6], "tolerance": 0.001 } }
+    },
+    {
+      "name": "f32_large_magnitude_accuracy_gpu_gap",
+      "skipGpu": {
+        "category": "todo",
+        "reason": "The current f32 implementation differs from the CPU reference by about one ULP at large finite magnitudes. A higher-accuracy f32 algorithm or software-extended precision could close this implementable gap."
+      },
+      "provenance": {
+        "source": "onnxruntime/test/providers/cpu/math/element_wise_ops_test.cc",
+        "test": "MathOpTest.Asinh",
+        "notes": "Large finite magnitudes exercise Asinh's overflow-avoidance region; CPU and ORT agree within a tighter bound than the current WGSL builtin path."
+      },
+      "inputs": {
+        "x": {
+          "dtype": "float32",
+          "shape": [6],
+          "data": {
+            "kind": "values",
+            "values": [10000000000.0, 100000000000000000000.0, -100000000000000000000.0, 1e+30, -1e+30, 3e+38]
+          }
+        }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [6], "tolerance": 0.000001, "relTolerance": 0 } }
+    },
+    {
+      "name": "ort_float_small_range",
+      "provenance": {
+        "source": "onnxruntime/test/providers/cpu/math/element_wise_ops_test.cc",
+        "test": "MathOpTest.Asinh"
+      },
+      "inputs": {
+        "x": { "dtype": "float32", "shape": [5], "data": { "kind": "values", "values": [-1.0, -0.5, 0.0, 0.5, 1.0] } }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [5], "tolerance": 0.00001 } }
+    },
+    {
+      "name": "ort_nonfinite_values",
+      "provenance": {
+        "source": "onnxruntime/test/providers/cpu/math/element_wise_ops_test.cc",
+        "test": "MathOpTest.Asinh",
+        "notes": "Extends ORT's small-range case with signed infinities, signed zero, and NaN."
+      },
+      "inputs": {
+        "x": {
+          "dtype": "float32",
+          "shape": [7],
+          "data": { "kind": "values", "values": ["-Infinity", -1.0, 0.0, 0.0, 1.0, "Infinity", "NaN"] }
+        }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [7], "tolerance": 0.00001, "allowNaN": true } }
+    },
+    {
+      "name": "onnx_backend_example",
+      "provenance": { "source": "cmake/external/onnx/onnx/backend/test/data/node/test_asinh_example" },
+      "inputs": { "x": { "dtype": "float32", "shape": [3], "data": { "kind": "values", "values": [-1.0, 0.0, 1.0] } } },
+      "outputs": { "y": { "dtype": "float32", "shape": [3], "tolerance": 0.000001 } }
+    },
+    {
+      "name": "onnx_backend_asinh",
+      "provenance": { "source": "cmake/external/onnx/onnx/backend/test/data/node/test_asinh" },
+      "inputs": {
+        "x": {
+          "dtype": "float32",
+          "shape": [3, 4, 5],
+          "data": {
+            "kind": "values",
+            "values": [1.764052391052246, 0.40015721321105957, 0.978738009929657, 2.2408931255340576, 1.8675580024719238, -0.9772778749465942, 0.9500884413719177, -0.15135720372200012, -0.10321885347366333, 0.4105985164642334, 0.14404356479644775, 1.4542734622955322, 0.7610377073287964, 0.12167501449584961, 0.44386324286460876, 0.3336743414402008, 1.4940791130065918, -0.2051582634449005, 0.3130677044391632, -0.8540957570075989, -2.5529897212982178, 0.653618574142456, 0.8644362092018127, -0.7421650290489197, 2.269754648208618, -1.4543657302856445, 0.04575851559638977, -0.18718385696411133, 1.5327792167663574, 1.4693588018417358, 0.154947429895401, 0.37816253304481506, -0.8877857327461243, -1.980796456336975, -0.34791216254234314, 0.15634897351264954, 1.2302906513214111, 1.202379822731018, -0.38732680678367615, -0.302302747964859, -1.0485529899597168, -1.420017957687378, -1.7062702178955078, 1.950775384902954, -0.5096521973609924, -0.4380742907524109, -1.2527953386306763, 0.7774903774261475, -1.6138978004455566, -0.21274028718471527, -0.8954665660858154, 0.38690251111984253, -0.5108051300048828, -1.18063223361969, -0.02818222902715206, 0.4283318817615509, 0.06651721894741058, 0.30247190594673157, -0.6343221068382263, -0.3627411723136902]
+          }
+        }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [3, 4, 5], "tolerance": 0.00001 } }
+    },
+    {
+      "name": "onnx_backend_asinh_example",
+      "provenance": { "source": "cmake/external/onnx/onnx/backend/test/data/node/test_asinh_example" },
+      "inputs": { "x": { "dtype": "float32", "shape": [3], "data": { "kind": "values", "values": [-1.0, 0.0, 1.0] } } },
+      "outputs": { "y": { "dtype": "float32", "shape": [3], "tolerance": 0.00001 } }
+    },
+    {
+      "name": "vec4_f16_lanes",
+      "inputs": {
+        "x": {
+          "dtype": "float16",
+          "shape": [16],
+          "data": {
+            "kind": "values",
+            "values": [-8.0, -4.0, -2.0, -1.0, -0.5, -0.25, 0.0, 0.25, 0.5, 1.0, 2.0, 4.0, 8.0, 16.0, 50.0, 100.0]
+          }
+        }
+      },
+      "outputs": { "y": { "dtype": "float16", "shape": [16], "tolerance": 0.001, "relTolerance": 0.002 } }
+    },
+    {
+      "name": "vec4_f32_large_magnitude_guard",
+      "inputs": {
+        "x": {
+          "dtype": "float32",
+          "shape": [4],
+          "data": {
+            "kind": "values",
+            "values": [1000000000000000000.0, -1000000000000000000.0, 100000000000000000000.0, -100000000000000000000.0]
+          }
+        }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [4] } }
+    },
+    {
+      "name": "empty_input_zero_dim",
+      "inputs": { "x": { "dtype": "float32", "shape": [0], "data": { "kind": "values", "values": [] } } },
+      "outputs": { "y": { "dtype": "float32", "shape": [0], "tolerance": 0 } }
+    },
+    {
+      "name": "f16_unaligned_scalar_fallback_tail",
+      "inputs": {
+        "x": {
+          "dtype": "float16",
+          "shape": [7],
+          "data": { "kind": "values", "values": [-100.0, -8.0, -0.75, 0.0, 0.75, 8.0, 100.0] }
+        }
+      },
+      "outputs": { "y": { "dtype": "float16", "shape": [7], "tolerance": 0.005, "relTolerance": 0.005 } }
+    },
+    {
+      "name": "f32_scalar_tail_4095",
+      "provenance": {
+        "notes": "Compact sibling for the f32 Asinh scalar-fallback benchmark; count is intentionally not divisible by four."
+      },
+      "inputs": {
+        "x": { "dtype": "float32", "shape": [4095], "data": { "kind": "linspace", "start": -10.0, "end": 10.0 } }
+      },
+      "outputs": { "y": { "dtype": "float32", "shape": [4095], "tolerance": 0.00001, "relTolerance": 0.00001 } }
+    }
+  ]
+}

build/webgpu/unary-scalar.wgsl.jinja ADDED Viewed

	@@ -0,0 +1,41 @@

+{% macro flat_tail_open() %}
+@compute @workgroup_size({{ tunables.WORKGROUP_SIZE }})
+fn main(@builtin(global_invocation_id) gid: vec3<u32>, @builtin(num_workgroups) nwg: vec3<u32>) {
+  // 2D-folded flat index: gid.y carries the high bits when the element count exceeds the
+  // maxComputeWorkgroupsPerDimension limit (the dispatch caps x and spills the rest into y).
+  let invocation = gid.x + gid.y * nwg.x * {{ tunables.WORKGROUP_SIZE }}u;
+  // Tail-safe scalar x4 keeps vector-like dispatch density without requiring
+  // the logical tensor length (or its storage binding) to be vec4 aligned.
+  let begin = invocation * {{ source.itemsPerInvocation }}u;
+  let end = min(begin + {{ source.itemsPerInvocation }}u, params.count);
+  for (var i = begin; i < end; i = i + 1u) {
+{%- endmacro %}
+{% macro flat_tail_close() %}
+  }
+{% endmacro %}
+// Scalar unary fallback. Each branch retains the operation's numeric hardening,
+// including Payne-Hanek trigonometric range reduction and NaN/overflow guards.
+{% if usesF16 %}
+enable f16;
+{% endif %}
+{{ env.wgsl.resourceDeclarations }}
+{% macro emit_asinh_safe() %}
+fn asinh_safe(x: f32) -> f32 {
+  // The log(x + sqrt(x*x + 1)) identity overflows in the x*x term for very
+  // large |x|, returning Inf instead of the true ~log(2|x|). Route large
+  // magnitudes through the algebraically-equivalent, overflow-free form
+  // (sqrt(x*x+1) == |x| to full f32 precision long before this point). The
+  // builtin path handles the normal range.
+  let a = abs(x);
+  if (a > 1.0e16) {
+    return sign(x) * (log(a) + 0.6931471805599453);
+  }
+  return asinh(x);
+}{% endmacro %}
+{{ emit_asinh_safe() }}
+{{ flat_tail_open() }}
+  let v = f32(x[i]);
+  y[i] = {{ scalar }}(asinh_safe(v));
+{{ flat_tail_close() -}}
+}

build/webgpu/unary-vec4.wgsl.jinja ADDED Viewed

	@@ -0,0 +1,34 @@

+// Loads and stores vec4<T> (128 bits) while retaining scalar per-component
+// arithmetic, including per-component helper calls for guard-heavy operations.
+{% if usesF16 %}
+enable f16;
+{% endif %}
+{{ env.wgsl.resourceDeclarations }}
+{% macro emit_asinh_safe() %}
+fn asinh_safe(x: f32) -> f32 {
+  // The log(x + sqrt(x*x + 1)) identity overflows in the x*x term for very
+  // large |x|, returning Inf instead of the true ~log(2|x|). Route large
+  // magnitudes through the algebraically-equivalent, overflow-free form
+  // (sqrt(x*x+1) == |x| to full f32 precision long before this point). The
+  // builtin path handles the normal range.
+  let a = abs(x);
+  if (a > 1.0e16) {
+    return sign(x) * (log(a) + 0.6931471805599453);
+  }
+  return asinh(x);
+}{% endmacro %}
+{{ emit_asinh_safe() }}
+@compute @workgroup_size({{ tunables.WORKGROUP_SIZE }})
+fn main(@builtin(global_invocation_id) gid: vec3<u32>, @builtin(num_workgroups) nwg: vec3<u32>) {
+  // 2D-folded flat index: gid.y carries the high bits when the element count exceeds the
+  // maxComputeWorkgroupsPerDimension limit (the dispatch caps x and spills the rest into y).
+  let i = gid.x + gid.y * nwg.x * {{ tunables.WORKGROUP_SIZE }}u;
+  if (i >= params.count) {
+    return;
+  }
+  let xv = x[i];
+  let fv = vec4<f32>(xv);
+  y[i] = {{ vectorScalar }}(vec4<f32>(asinh_safe(fv.x), asinh_safe(fv.y), asinh_safe(fv.z), asinh_safe(fv.w)));
+}