Merge pull request #492 from neutrinoceros/numpy2/copy_semantics

unyt/array.py

@@ -10,6 +10,7 @@
 import sys
 import warnings
 from functools import lru_cache
+from importlib.metadata import version
 from numbers import Number as numeric_type
 
 import numpy as np
@@ -102,6 +103,7 @@
     true_divide,
     trunc,
 )
+from packaging.version import Version
 from sympy import Rational
 
 from unyt._on_demand_imports import _astropy, _dask, _pint
@@ -146,6 +148,8 @@
     ("unyt_array.from_astropy", "unyt_array.to_astropy"): ["astropy"],
 }
 
+_COPY_IF_NEEDED = None if Version(version("numpy")) >= Version("2.0.0dev0") else False
+
 # This is partially adapted from the following SO thread
 # https://stackoverflow.com/questions/41668588/regex-to-match-scientific-notation
 _NUMB_PATTERN = (
@@ -1985,7 +1989,7 @@ def __array_ufunc__(self, ufunc, method, *inputs, **kwargs):
                     "added to unyt_array." % (str(ufunc), len(inputs))
                 )
         if unit is None:
-            out_arr = np.array(out_arr, copy=False)
+            out_arr = np.array(out_arr, copy=_COPY_IF_NEEDED)
         elif ufunc in (modf, divmod_):
             out_arr = tuple(ret_class(o, unit) for o in out_arr)
         elif out_arr.shape == ():