Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add support for scalar arguments to xp.where #78

Merged
merged 6 commits into from
Feb 3, 2025
Merged
Show file tree
Hide file tree
Changes from 3 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 10 additions & 3 deletions array_api_strict/_searching_functions.py
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,7 @@

from ._array_object import Array
from ._dtypes import _result_type, _real_numeric_dtypes
from ._flags import requires_data_dependent_shapes, requires_api_version
from ._flags import requires_data_dependent_shapes, requires_api_version, get_array_api_strict_flags

from typing import TYPE_CHECKING
if TYPE_CHECKING:
Expand Down Expand Up @@ -72,17 +72,24 @@ def searchsorted(
# x1 must be 1-D, but NumPy already requires this.
return Array._new(np.searchsorted(x1._array, x2._array, side=side, sorter=sorter), device=x1.device)

def where(condition: Array, x1: Array, x2: Array, /) -> Array:
def where(condition: Array, x1: bool | int | float | Array, x2: bool | int | float | Array, /) -> Array:
"""
Array API compatible wrapper for :py:func:`np.where <numpy.where>`.

See its docstring for more information.
"""
if get_array_api_strict_flags()['api_version'] > '2023.12':
if isinstance(x1, (bool, float, int)):
ev-br marked this conversation as resolved.
Show resolved Hide resolved
x1 = Array._new(np.asarray(x1), device=condition.device)

if isinstance(x2, (bool, float, int)):
x2 = Array._new(np.asarray(x2), device=condition.device)

# Call result type here just to raise on disallowed type combinations
_result_type(x1.dtype, x2.dtype)

if len({a.device for a in (condition, x1, x2)}) > 1:
raise ValueError("where inputs must all be on the same device")
raise ValueError("Inputs to `where` must all use the same device")

x1, x2 = Array._normalize_two_args(x1, x2)
return Array._new(np.where(condition._array, x1._array, x2._array), device=x1.device)
26 changes: 26 additions & 0 deletions array_api_strict/tests/test_searching_functions.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,26 @@
import pytest

import array_api_strict as xp

from array_api_strict import ArrayAPIStrictFlags
from array_api_strict._flags import draft_version


def test_where_with_scalars():
x = xp.asarray([1, 2, 3, 1])

# Versions up to and including 2023.12 don't support scalar arguments
with pytest.raises(AttributeError, match="object has no attribute 'dtype'"):
xp.where(x == 1, 42, 44)

# Versions after 2023.12 support scalar arguments
with (pytest.warns(
UserWarning,
match="The 2024.12 version of the array API specification is in draft status"
),
ArrayAPIStrictFlags(api_version=draft_version),
):
x_where = xp.where(x == 1, 42, 44)

expected = xp.asarray([42, 44, 44, 42])
assert xp.all(x_where == expected)
Loading