arraybridge 0.3.6__tar.gz → 0.3.7__tar.gz

{arraybridge-0.3.6 → arraybridge-0.3.7}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: arraybridge
-Version: 0.3.6
+Version: 0.3.7
 Summary: Unified API for NumPy, CuPy, PyTorch, TensorFlow, JAX, and pyclesperanto with automatic memory type conversion
 Project-URL: Homepage, https://github.com/OpenHCSDev/arraybridge
 Project-URL: Documentation, https://arraybridge.readthedocs.io
@@ -143,3 +143,17 @@ pip install "arraybridge[cupy]"
 ```
 
 Documentation: <https://arraybridge.readthedocs.io>
+
+### Native array geometry
+
+`MemoryType.reshape(array, shape)` and `MemoryType.broadcast_to(array, shape)`
+operate in the array's framework and device. `MemoryType.ones_like(reference,
+shape=..., dtype=...)` allocates on the reference's device; omitting `shape`
+retains its shape. Reshape and broadcast preserve views where the backend can
+represent them. Broadcasting is not a promise of a writable, independent buffer.
+
+These operations never silently project pixels to NumPy. Pyclesperanto cannot
+provide native changed-shape reshape or broadcasting, so those requests raise
+`NotImplementedError`. Its native ones allocation supports one to three
+dimensions and its supported integer/float32 dtypes; boolean allocation is
+explicitly unsupported. Existing conversion and disk boundaries remain separate.

{arraybridge-0.3.6 → arraybridge-0.3.7}/README.md

@@ -71,3 +71,17 @@ pip install "arraybridge[cupy]"
 ```
 
 Documentation: <https://arraybridge.readthedocs.io>
+
+### Native array geometry
+
+`MemoryType.reshape(array, shape)` and `MemoryType.broadcast_to(array, shape)`
+operate in the array's framework and device. `MemoryType.ones_like(reference,
+shape=..., dtype=...)` allocates on the reference's device; omitting `shape`
+retains its shape. Reshape and broadcast preserve views where the backend can
+represent them. Broadcasting is not a promise of a writable, independent buffer.
+
+These operations never silently project pixels to NumPy. Pyclesperanto cannot
+provide native changed-shape reshape or broadcasting, so those requests raise
+`NotImplementedError`. Its native ones allocation supports one to three
+dimensions and its supported integer/float32 dtypes; boolean allocation is
+explicitly unsupported. Existing conversion and disk boundaries remain separate.

{arraybridge-0.3.6 → arraybridge-0.3.7}/docs/source/index.rst

@@ -20,6 +20,8 @@ applications.
    examples/index
    contributing
    ci-cd
+   release_035_validation
+   release_036_validation
 
 Boundary rule
 -------------

{arraybridge-0.3.6 → arraybridge-0.3.7}/docs/source/release_035_validation.rst

@@ -1,5 +1,5 @@
 ArrayBridge 0.3.5 release preparation
-===================================
+=====================================
 
 Parent integration owner; base ea3f2a4cc91c4810d12343f58f85c1195e1a41e6.
 The reviewed callable-native dtype default and durable decorator-context

{arraybridge-0.3.6 → arraybridge-0.3.7}/docs/source/release_036_validation.rst

@@ -1,5 +1,5 @@
 ArrayBridge typed thread-local release correction
-================================================
+=================================================
 
 Parent integration owner. Actual failed original 0.3.5 publisher:
 https://github.com/OpenHCSDev/ArrayBridge/actions/runs/36946585391.

{arraybridge-0.3.6 → arraybridge-0.3.7}/src/arraybridge/__init__.py

@@ -10,7 +10,7 @@ attribute access instead of at import time. Declaration-only consumers
 NumPy/numcodecs import cost.
 """
 
-__version__ = "0.3.6"
+__version__ = "0.3.7"
 
 _LAZY_EXPORTS: dict[str, str] = {
     "MemoryType": ".types",

arraybridge-0.3.7/src/arraybridge/array_operations.py

@@ -0,0 +1,450 @@
+"""Native array operations carried by MemoryType declarations."""
+
+from abc import ABC, abstractmethod
+from collections.abc import Sequence
+from typing import Any
+
+import numpy as np
+
+_SCALING_RANGES: dict[str, float | tuple[float, float]] = {
+    "uint8": 255.0,
+    "uint16": 65535.0,
+    "uint32": 4294967295.0,
+    "int16": (65535.0, 32768.0),
+    "int32": (4294967295.0, 2147483648.0),
+}
+
+
+def _dtype_name(dtype: Any) -> str:
+    declared_name = getattr(dtype, "name", None)
+    if declared_name is not None:
+        return str(declared_name)
+    return getattr(dtype, "__name__", str(dtype).rsplit(".", maxsplit=1)[-1])
+
+
+def _numpy_dtype_name(dtype: Any) -> str:
+    return str(np.dtype(dtype).name)
+
+
+def _torch_dtype_name(dtype: Any) -> str:
+    return str(dtype).rsplit(".", maxsplit=1)[-1]
+
+
+def _tensorflow_dtype_name(dtype: Any) -> str:
+    numpy_dtype = getattr(dtype, "as_numpy_dtype", dtype)
+    return str(np.dtype(numpy_dtype).name)
+
+
+def _scaled_values(result: Any, result_min: Any, result_max: Any, target_dtype: Any) -> Any:
+    normalized = (result - result_min) / (result_max - result_min)
+    range_info = _SCALING_RANGES.get(_dtype_name(target_dtype))
+    if range_info is None:
+        return normalized
+    if isinstance(range_info, tuple):
+        scale, offset = range_info
+        return normalized * scale - offset
+    return normalized * range_info
+
+
+def _clamp_bounds(target_dtype: Any) -> tuple[float, float] | None:
+    range_info = _SCALING_RANGES.get(_dtype_name(target_dtype))
+    if range_info is None:
+        return None
+    if isinstance(range_info, tuple):
+        scale, offset = range_info
+        return -offset, scale - offset - 128
+    return 0, range_info
+
+
+def _mapped_dtype(target_dtype: Any, module: Any) -> Any:
+    try:
+        dtype_name = np.dtype(target_dtype).name
+    except TypeError as error:
+        raise TypeError(f"Unsupported target dtype {target_dtype!r}") from error
+    mapped = getattr(module, dtype_name, None)
+    if mapped is None:
+        module_name = getattr(module, "__name__", type(module).__name__)
+        raise TypeError(f"{module_name} does not expose dtype {dtype_name}")
+    return mapped
+
+
+class ArrayOperations(ABC):
+    """Native array semantics selected by the existing MemoryType declaration."""
+
+    @staticmethod
+    @abstractmethod
+    def to_numpy(data: Any, module: Any) -> Any:
+        """Project pixels to the explicit host boundary."""
+
+    @staticmethod
+    @abstractmethod
+    def from_numpy(data: Any, module: Any, device_id: int) -> Any:
+        """Admit host pixels on the selected framework device."""
+
+    @staticmethod
+    @abstractmethod
+    def scale_dtype(result: Any, target_dtype: Any, module: Any) -> Any:
+        """Apply the framework's existing intensity conversion law."""
+
+    @staticmethod
+    def dtype_name(dtype: Any) -> str:
+        return _numpy_dtype_name(dtype)
+
+    @staticmethod
+    def stack(values: Sequence[Any], module: Any) -> Any:
+        return module.stack(tuple(values), axis=0)
+
+    @staticmethod
+    def cast(data: Any, dtype: Any, module: Any) -> Any:
+        return data.astype(dtype, copy=False)
+
+    @staticmethod
+    def logical_and(left: Any, right: Any, module: Any) -> Any:
+        return module.logical_and(left, right)
+
+    @staticmethod
+    def reshape(data: Any, shape: tuple[int, ...], module: Any) -> Any:
+        return data.reshape(shape)
+
+    @staticmethod
+    def broadcast_to(data: Any, shape: tuple[int, ...], module: Any) -> Any:
+        return module.broadcast_to(data, shape)
+
+    @staticmethod
+    def ones_like(reference: Any, shape: tuple[int, ...], dtype: Any, module: Any) -> Any:
+        return module.ones(shape, dtype=dtype)
+
+
+class NumpyArrayOperations(ArrayOperations):
+    """Numpy native operation leaves."""
+
+    @staticmethod
+    def to_numpy(data: Any, module: Any) -> Any:
+        del module
+        return data
+
+    @staticmethod
+    def from_numpy(data: Any, module: Any, device_id: int) -> Any:
+        del module, device_id
+        return data
+
+    @staticmethod
+    def scale_dtype(result: Any, target_dtype: Any, module: Any) -> Any:
+        if not hasattr(result, "dtype"):
+            return result
+        if not (
+            module.issubdtype(result.dtype, module.floating)
+            and module.issubdtype(target_dtype, module.integer)
+        ):
+            return result.astype(target_dtype)
+        result_min = result.min()
+        result_max = result.max()
+        if result_max <= result_min:
+            return result.astype(target_dtype)
+        scaled = _scaled_values(result, result_min, result_max, target_dtype)
+        bounds = _clamp_bounds(target_dtype)
+        if bounds is not None:
+            scaled = module.clip(scaled, *bounds)
+        return scaled.astype(target_dtype)
+
+
+class CupyArrayOperations(ArrayOperations):
+    """Cupy native operation leaves."""
+
+    @staticmethod
+    def to_numpy(data: Any, module: Any) -> Any:
+        del module
+        return data.get()
+
+    @staticmethod
+    def from_numpy(data: Any, module: Any, device_id: int) -> Any:
+        del device_id
+        return module.array(data)
+
+    @staticmethod
+    def scale_dtype(result: Any, target_dtype: Any, module: Any) -> Any:
+        if not hasattr(result, "dtype"):
+            return result
+        if not (
+            module.issubdtype(result.dtype, module.floating)
+            and not module.issubdtype(target_dtype, module.floating)
+        ):
+            return result.astype(target_dtype)
+        result_min = module.min(result)
+        result_max = module.max(result)
+        if result_max <= result_min:
+            return result.astype(target_dtype)
+        scaled = _scaled_values(result, result_min, result_max, target_dtype)
+        bounds = _clamp_bounds(target_dtype)
+        if bounds is not None:
+            scaled = module.clip(scaled, *bounds)
+        return scaled.astype(target_dtype)
+
+
+class TorchArrayOperations(ArrayOperations):
+    """Torch native operation leaves."""
+
+    @staticmethod
+    def to_numpy(data: Any, module: Any) -> Any:
+        del module
+        return data.cpu().numpy()
+
+    @staticmethod
+    def from_numpy(data: Any, module: Any, device_id: int) -> Any:
+        host_data = (
+            np.ascontiguousarray(data)
+            if any(stride < 0 for stride in getattr(data, "strides", ()))
+            else data
+        )
+        return module.from_numpy(host_data).to(f"cuda:{device_id}")
+
+    @staticmethod
+    def scale_dtype(result: Any, target_dtype: Any, module: Any) -> Any:
+        if not hasattr(result, "dtype"):
+            return result
+        mapped = _mapped_dtype(target_dtype, module)
+        floats = (module.float16, module.float32, module.float64)
+        if not (result.dtype in floats and np.issubdtype(np.dtype(target_dtype), np.integer)):
+            return result.to(mapped)
+        result_min = result.min()
+        result_max = result.max()
+        if result_max <= result_min:
+            return result.to(mapped)
+        scaled = _scaled_values(result, result_min, result_max, target_dtype)
+        bounds = _clamp_bounds(target_dtype)
+        if bounds is not None:
+            scaled = module.clamp(scaled, min=bounds[0], max=bounds[1])
+        return scaled.to(mapped)
+
+    @staticmethod
+    def stack(values: Sequence[Any], module: Any) -> Any:
+        return module.stack(tuple(values), dim=0)
+
+    @staticmethod
+    def cast(data: Any, dtype: Any, module: Any) -> Any:
+        return data.to(dtype=_mapped_dtype(dtype, module))
+
+    @staticmethod
+    def dtype_name(dtype: Any) -> str:
+        return str(dtype).rsplit(".", maxsplit=1)[-1]
+
+    @staticmethod
+    def broadcast_to(data: Any, shape: tuple[int, ...], module: Any) -> Any:
+        return data.expand(shape)
+
+    @staticmethod
+    def ones_like(reference: Any, shape: tuple[int, ...], dtype: Any, module: Any) -> Any:
+        return module.ones(shape, dtype=_mapped_dtype(dtype, module), device=reference.device)
+
+
+class TensorflowArrayOperations(ArrayOperations):
+    """Tensorflow native operation leaves."""
+
+    @staticmethod
+    def to_numpy(data: Any, module: Any) -> Any:
+        del module
+        return data.numpy()
+
+    @staticmethod
+    def from_numpy(data: Any, module: Any, device_id: int) -> Any:
+        del device_id
+        return module.convert_to_tensor(data)
+
+    @staticmethod
+    def scale_dtype(result: Any, target_dtype: Any, module: Any) -> Any:
+        if not hasattr(result, "dtype"):
+            return result
+        mapped = _mapped_dtype(target_dtype, module)
+        floats = (module.float16, module.float32, module.float64)
+        if not (result.dtype in floats and np.issubdtype(np.dtype(target_dtype), np.integer)):
+            return module.cast(result, mapped)
+        result_min = module.reduce_min(result)
+        result_max = module.reduce_max(result)
+        if result_max <= result_min:
+            return module.cast(result, mapped)
+        scaled = _scaled_values(result, result_min, result_max, target_dtype)
+        bounds = _clamp_bounds(target_dtype)
+        if bounds is not None:
+            scaled = module.clip_by_value(scaled, *bounds)
+        return module.cast(scaled, mapped)
+
+    @staticmethod
+    def cast(data: Any, dtype: Any, module: Any) -> Any:
+        return module.cast(data, _mapped_dtype(dtype, module))
+
+    @staticmethod
+    def reshape(data: Any, shape: tuple[int, ...], module: Any) -> Any:
+        return module.reshape(data, shape)
+
+    @staticmethod
+    def ones_like(reference: Any, shape: tuple[int, ...], dtype: Any, module: Any) -> Any:
+        with module.device(reference.device):
+            return module.ones(shape, dtype=_mapped_dtype(dtype, module))
+
+
+class JaxArrayOperations(ArrayOperations):
+    """Jax native operation leaves."""
+
+    @staticmethod
+    def to_numpy(data: Any, module: Any) -> Any:
+        del module
+        return np.asarray(data)
+
+    @staticmethod
+    def from_numpy(data: Any, module: Any, device_id: int) -> Any:
+        devices = tuple(device for device in module.devices() if device.platform == "gpu")
+        return module.device_put(data, devices[device_id])
+
+    @staticmethod
+    def scale_dtype(result: Any, target_dtype: Any, module: Any) -> Any:
+        if not hasattr(result, "dtype"):
+            return result
+        if np.dtype(target_dtype) == np.dtype(np.float64):
+            x64_enabled = getattr(module.config, "x64_enabled", None)
+            if x64_enabled is None:
+                x64_enabled = module.config.read("jax_enable_x64")
+            if not x64_enabled:
+                raise ValueError(
+                    "JAX float64 output requires x64 mode; set JAX_ENABLE_X64=true before import"
+                )
+        jnp = module.numpy
+        mapped = _mapped_dtype(target_dtype, jnp)
+        floats = (jnp.float16, jnp.float32, jnp.float64)
+        if not (result.dtype in floats and np.issubdtype(np.dtype(target_dtype), np.integer)):
+            return result.astype(mapped)
+        result_min = jnp.min(result)
+        result_max = jnp.max(result)
+        if result_max <= result_min:
+            return result.astype(mapped)
+        scaled = _scaled_values(result, result_min, result_max, target_dtype)
+        bounds = _clamp_bounds(target_dtype)
+        if bounds is not None:
+            scaled = jnp.clip(scaled, *bounds)
+        return scaled.astype(mapped)
+
+    @staticmethod
+    def stack(values: Sequence[Any], module: Any) -> Any:
+        return module.numpy.stack(tuple(values), axis=0)
+
+    @staticmethod
+    def cast(data: Any, dtype: Any, module: Any) -> Any:
+        return data.astype(_mapped_dtype(dtype, module.numpy))
+
+    @staticmethod
+    def logical_and(left: Any, right: Any, module: Any) -> Any:
+        return module.numpy.logical_and(left, right)
+
+    @staticmethod
+    def broadcast_to(data: Any, shape: tuple[int, ...], module: Any) -> Any:
+        return module.numpy.broadcast_to(data, shape)
+
+    @staticmethod
+    def ones_like(reference: Any, shape: tuple[int, ...], dtype: Any, module: Any) -> Any:
+        return module.numpy.ones(shape, dtype=dtype, device=reference.device)
+
+
+class PyclesperantoArrayOperations(ArrayOperations):
+    """Pyclesperanto native operation leaves."""
+
+    @staticmethod
+    def to_numpy(data: Any, module: Any) -> Any:
+        return module.pull(data)
+
+    @staticmethod
+    def from_numpy(data: Any, module: Any, device_id: int) -> Any:
+        del device_id
+        return module.push(data)
+
+    @staticmethod
+    def scale_dtype(result: Any, target_dtype: Any, module: Any) -> Any:
+        if not hasattr(result, "dtype"):
+            return result
+        target_is_int = np.issubdtype(np.dtype(target_dtype), np.integer)
+        if not (np.issubdtype(result.dtype, np.floating) and target_is_int):
+            return module.push(module.pull(result).astype(target_dtype))
+        result_min = float(module.minimum_of_all_pixels(result))
+        result_max = float(module.maximum_of_all_pixels(result))
+        if result_max <= result_min:
+            return module.push(module.pull(result).astype(target_dtype))
+        normalized = module.subtract_image_from_scalar(result, scalar=result_min)
+        normalized = module.multiply_image_and_scalar(
+            normalized,
+            scalar=1.0 / (result_max - result_min),
+        )
+        range_info = _SCALING_RANGES.get(_dtype_name(target_dtype))
+        if isinstance(range_info, tuple):
+            scale, offset = range_info
+            scaled = module.multiply_image_and_scalar(normalized, scalar=scale)
+            scaled = module.subtract_image_from_scalar(scaled, scalar=offset)
+        elif range_info is not None:
+            scaled = module.multiply_image_and_scalar(normalized, scalar=range_info)
+        else:
+            scaled = normalized
+        host_values = module.pull(scaled)
+        bounds = _clamp_bounds(target_dtype)
+        if bounds is not None:
+            host_values = np.clip(host_values, *bounds)
+        return module.push(host_values.astype(target_dtype))
+
+    @staticmethod
+    def stack(values: Sequence[Any], module: Any) -> Any:
+        if not values:
+            raise ValueError("Cannot stack an empty pyclesperanto sequence")
+        if len(values) == 1:
+            source = values[0]
+            result = module.create((1, *source.shape), dtype=source.dtype)
+            return module.copy_slice(source, result, 0)
+        result = values[0]
+        for value in values[1:]:
+            result = module.concatenate_along_z(result, value)
+        return result
+
+    @staticmethod
+    def cast(data: Any, dtype: Any, module: Any) -> Any:
+        return module.push(module.pull(data).astype(dtype, copy=False))
+
+    @staticmethod
+    def logical_and(left: Any, right: Any, module: Any) -> Any:
+        return module.push(np.logical_and(module.pull(left), module.pull(right)))
+
+    @staticmethod
+    def reshape(data: Any, shape: tuple[int, ...], module: Any) -> Any:
+        if tuple(data.shape) == shape:
+            return data
+        raise NotImplementedError(
+            "pyclesperanto does not provide native array reshape; its reshape downloads pixels"
+        )
+
+    @staticmethod
+    def broadcast_to(data: Any, shape: tuple[int, ...], module: Any) -> Any:
+        if tuple(data.shape) == shape:
+            return data
+        raise NotImplementedError("pyclesperanto does not provide native array broadcasting")
+
+    @staticmethod
+    def ones_like(reference: Any, shape: tuple[int, ...], dtype: Any, module: Any) -> Any:
+        if np.dtype(dtype).name not in {
+            "float32",
+            "int8",
+            "int16",
+            "int32",
+            "uint8",
+            "uint16",
+            "uint32",
+        }:
+            raise NotImplementedError(f"pyclesperanto cannot allocate native dtype {dtype!r}")
+        if not 1 <= len(shape) <= 3:
+            raise NotImplementedError(
+                "pyclesperanto native allocation requires one to three dimensions"
+            )
+        result = module.create(shape, dtype=dtype, device=reference.device)
+        module.set(result, scalar=1, device=reference.device)
+        return result
+
+
+NUMPY_OPERATIONS = NumpyArrayOperations()
+CUPY_OPERATIONS = CupyArrayOperations()
+TORCH_OPERATIONS = TorchArrayOperations()
+TENSORFLOW_OPERATIONS = TensorflowArrayOperations()
+JAX_OPERATIONS = JaxArrayOperations()
+PYCLESPERANTO_OPERATIONS = PyclesperantoArrayOperations()

{arraybridge-0.3.6 → arraybridge-0.3.7}/src/arraybridge/types.py

@@ -12,7 +12,7 @@ import importlib.util
 import logging
 import os
 import sys
-from collections.abc import Callable, Iterator, Mapping, MutableMapping
+from collections.abc import Callable, Iterator, Mapping, MutableMapping, Sequence
 from contextlib import AbstractContextManager, contextmanager, nullcontext
 from dataclasses import dataclass
 from enum import Enum
@@ -834,11 +834,44 @@ class MemoryType(_MemoryTypeFields, Enum):
         """Cast one array through this framework member's operation leaf."""
 
         framework = module if module is not None else self.import_module()
-        device_id = self.device_id_of(data, framework)
-        scope = nullcontext() if device_id is None else self.device_scope(device_id, framework)
-        with scope:
+        with self._array_device_scope(data, framework):
             return self._operations.cast(data, dtype, framework)
 
+    def _array_device_scope(self, reference: Any, module: Any) -> AbstractContextManager[None]:
+        """Derive every native operation scope from its actual input array."""
+        device_id = self.device_id_of(reference, module)
+        return nullcontext() if device_id is None else self.device_scope(device_id, module)
+
+    def reshape(self, data: Any, shape: Sequence[int], module: Any | None = None) -> Any:
+        """Reshape on the input device without an implicit host projection."""
+        framework = module if module is not None else self.import_module()
+        with self._array_device_scope(data, framework):
+            return self._operations.reshape(data, tuple(shape), framework)
+
+    def broadcast_to(self, data: Any, shape: Sequence[int], module: Any | None = None) -> Any:
+        """Broadcast on the input device without an implicit host projection."""
+        framework = module if module is not None else self.import_module()
+        with self._array_device_scope(data, framework):
+            return self._operations.broadcast_to(data, tuple(shape), framework)
+
+    def ones_like(
+        self,
+        reference: Any,
+        *,
+        shape: Sequence[int] | None = None,
+        dtype: Any = bool,
+        module: Any | None = None,
+    ) -> Any:
+        """Allocate ones on the supplied reference's framework-local device."""
+        framework = module if module is not None else self.import_module()
+        with self._array_device_scope(reference, framework):
+            return self._operations.ones_like(
+                reference,
+                tuple(reference.shape if shape is None else shape),
+                dtype,
+                framework,
+            )
+
     def logical_and(
         self,
         left: Any,
@@ -848,9 +881,7 @@ class MemoryType(_MemoryTypeFields, Enum):
         """Intersect two arrays through this framework member's operation leaf."""
 
         framework = module if module is not None else self.import_module()
-        device_id = self.device_id_of(left, framework)
-        scope = nullcontext() if device_id is None else self.device_scope(device_id, framework)
-        with scope:
+        with self._array_device_scope(left, framework):
             return self._operations.logical_and(left, right, framework)
 
     def available_device_ids(self, module: Any | None = None) -> tuple[int, ...]:

{arraybridge-0.3.6 → arraybridge-0.3.7}/tests/test_callable_dtype_default.py

@@ -1,7 +1,7 @@
 """Callable defaults use the existing typed dtype policy and real conversion."""
 
-from dataclasses import dataclass
 import inspect
+from dataclasses import dataclass
 
 import numpy as np
 import pytest
@@ -10,6 +10,8 @@ from arraybridge.decorators import (
     DtypeConversion,
     DtypeConversionConfig,
     PreserveInputDtypeConfig,
+)
+from arraybridge.decorators import (
     numpy as numpy_func,
 )
 

arraybridge-0.3.7/tests/test_native_geometry.py

@@ -0,0 +1,65 @@
+"""Native geometry must preserve pixels, sharing and the reference device."""
+
+import numpy as np
+import pytest
+
+from arraybridge.array_operations import PYCLESPERANTO_OPERATIONS
+from arraybridge.types import MemoryType
+
+
+def test_numpy_views_and_independent_reference_allocation():
+    source = np.arange(12, dtype=np.float32).reshape(3, 4)[:, ::-1]
+    shaped = MemoryType.NUMPY.reshape(source, (3, 1, 4))
+    broadcast = MemoryType.NUMPY.broadcast_to(shaped, (3, 2, 4))
+    assert np.shares_memory(shaped, source)
+    assert np.shares_memory(broadcast, source)
+    np.testing.assert_array_equal(broadcast, np.broadcast_to(source[:, None], (3, 2, 4)))
+    ones = MemoryType.NUMPY.ones_like(source, shape=(2, 4), dtype=bool)
+    assert ones.dtype == np.bool_
+    assert not np.shares_memory(ones, source)
+    np.testing.assert_array_equal(ones, np.ones((2, 4), dtype=bool))
+
+
+def test_numpy_geometry_rejects_invalid_shapes():
+    source = np.arange(6).reshape(2, 3)
+    with pytest.raises(ValueError):
+        MemoryType.NUMPY.reshape(source, (7,))
+    with pytest.raises(ValueError):
+        MemoryType.NUMPY.broadcast_to(source, (2, 4))
+
+
+def test_cupy_geometry_and_allocation_stay_native():
+    cp = pytest.importorskip("cupy")
+    if not MemoryType.CUPY.available_device_ids(cp):
+        pytest.skip("No CUDA device")
+    source = cp.asarray(np.arange(12, dtype=np.float32).reshape(3, 4))[:, ::-1]
+    shaped = MemoryType.CUPY.reshape(source, (3, 1, 4))
+    broadcast = MemoryType.CUPY.broadcast_to(shaped, (3, 2, 4))
+    assert shaped.data.ptr == source.data.ptr
+    assert broadcast.data.ptr == source.data.ptr
+    assert broadcast.device.id == source.device.id
+    np.testing.assert_array_equal(
+        broadcast.get(), np.broadcast_to(source.get()[:, None], (3, 2, 4))
+    )
+    ones = MemoryType.CUPY.ones_like(source, shape=(2, 4), dtype=bool)
+    assert ones.device.id == source.device.id
+    assert ones.dtype == cp.bool_
+    np.testing.assert_array_equal(ones.get(), np.ones((2, 4), dtype=bool))
+
+
+def test_pycles_geometry_rejects_host_fallback():
+    class DeviceOnlyArray:
+        shape = (2, 3)
+
+        def __array__(self):
+            pytest.fail("Native geometry must never download pixels")
+
+    source = DeviceOnlyArray()
+    assert PYCLESPERANTO_OPERATIONS.reshape(source, (2, 3), None) is source
+    assert PYCLESPERANTO_OPERATIONS.broadcast_to(source, (2, 3), None) is source
+    with pytest.raises(NotImplementedError, match="reshape"):
+        PYCLESPERANTO_OPERATIONS.reshape(source, (1, 2, 3), None)
+    with pytest.raises(NotImplementedError, match="broadcasting"):
+        PYCLESPERANTO_OPERATIONS.broadcast_to(source, (4, 2, 3), None)
+    with pytest.raises(NotImplementedError, match="dtype"):
+        PYCLESPERANTO_OPERATIONS.ones_like(source, (2, 3), bool, None)