diff --git a/pint/numpy_func.py b/pint/numpy_func.py
index 75a4de6..24bb8d5 100644
--- a/pint/numpy_func.py
+++ b/pint/numpy_func.py
@@ -51,7 +51,8 @@ def _get_first_input_units(args, kwargs={}):
 def convert_arg(arg, pre_calc_units):
     """Convert quantities and sequences of quantities to pre_calc_units and strip units.
 
-    Helper function for convert_to_consistent_units.
+    Helper function for convert_to_consistent_units. pre_calc_units must be given as a pint
+    Unit or None.
     """
     if pre_calc_units is not None:
         if _is_quantity(arg):
@@ -183,7 +184,30 @@ def implements(numpy_func_string, func_type):
 
 
 def implement_func(func_type, func_str, input_units=None, output_unit=None):
-    """Add default-behavior NumPy function/ufunc to the handled list."""
+    """Add default-behavior NumPy function/ufunc to the handled list.
+
+    Parameters
+    ----------
+    func_type : str
+        "function" for NumPy functions, "ufunc" for NumPy ufuncs
+    func_str : str
+        String representing the name of the NumPy function/ufunc to add
+    input_units : pint.Unit or str or None
+        Parameter to control how the function downcasts to magnitudes of arguments. If
+        `pint.Unit`, converts all args and kwargs to this unit before downcasting to
+        magnitude. If "all_consistent", converts all args and kwargs to the unit of the
+        first Quantity in args and kwargs before downcasting to magnitude. If some
+        other string, the string is parsed as a unit, and all args and kwargs are
+        converted to that unit. If None, units are stripped without conversion.
+    output_unit : pint.Unit or str or None
+        Parameter to control the unit of the output. If `pint.Unit`, output is wrapped
+        with that unit. If "match_input", output is wrapped with the unit of the first
+        Quantity in args and kwargs. If a string representing a unit operation defined
+        in `get_op_output_unit`, output is wrapped by the unit determined by
+        `get_op_output_unit`. If some other string, the string is parsed as a unit,
+        which becomes the unit of the output. If None, the bare magnitude is returned.
+
+    """
     # If NumPy is not available, do not attempt implement that which does not exist
     if np is None:
         return
@@ -199,10 +223,16 @@ def implement_func(func_type, func_str, input_units=None, output_unit=None):
                 *args, pre_calc_units=first_input_units, **kwargs
             )
         else:
+            if isinstance(input_units, str):
+                # Conversion requires Unit, not str
+                pre_calc_units = first_input_units._REGISTRY.parse_units(input_units)
+            else:
+                pre_calc_units = input_units
+
             # Match all input args/kwargs to input_units, or if input_units is None,
             # simply strip units
             stripped_args, stripped_kwargs = convert_to_consistent_units(
-                *args, pre_calc_units=input_units, **kwargs
+                *args, pre_calc_units=pre_calc_units, **kwargs
             )
 
         # Determine result through base numpy function on stripped arguments
@@ -675,17 +705,31 @@ def implement_atleast_nd(func_str):
 for func_str in ["atleast_1d", "atleast_2d", "atleast_3d"]:
     implement_atleast_nd(func_str)
 
-# Handle single-argument consistent unit functions
-for func_str in ["block", "hstack", "vstack", "dstack", "column_stack"]:
-    implement_func(
-        "function", func_str, input_units="all_consistent", output_unit="match_input"
-    )
 
 # Handle cumulative products (which must be dimensionless for consistent units across
 # output array)
+def implement_single_dimensionless_argument_func(func_str):
+    # If NumPy is not available, do not attempt implement that which does not exist
+    if np is None:
+        return
+
+    func = getattr(np, func_str)
+
+    @implements(func_str, "function")
+    def implementation(a, *args, **kwargs):
+        (a_stripped,), _ = convert_to_consistent_units(
+            a, pre_calc_units=a._REGISTRY.parse_units("dimensionless")
+        )
+        return a._REGISTRY.Quantity(func(a_stripped, *args, **kwargs))
+
+
 for func_str in ["cumprod", "cumproduct", "nancumprod"]:
+    implement_single_dimensionless_argument_func(func_str)
+
+# Handle single-argument consistent unit functions
+for func_str in ["block", "hstack", "vstack", "dstack", "column_stack"]:
     implement_func(
-        "function", func_str, input_units="dimensionless", output_unit="match_input"
+        "function", func_str, input_units="all_consistent", output_unit="match_input"
     )
 
 # Handle functions that ignore units on input and output

array function implementation for np.cumprod / np.nancumprod fails with axis
Passing a `axis` argument to `cumprod` / `nancumprod` will make them fail. I tracked it down to https://github.com/hgrecco/pint/blob/4c3114c72bdae1e3c5c9f839659169eeb06071cd/pint/numpy_func.py#L62 where `pre_calc_units` can be a string (`"dimensionless"`), but a `Unit` object was expected. The string was created when `implement_func` was called for `cumprod` and `nancumprod` (by the way, `cumprod` also gets implemented with `set_units_ufuncs`). We can fix this by passing `None` as `input_units` instead, but if feasible I'd prefer `"dimensionless"` since that is way easier to read and understand. For that, we might need to preprocess the `input_units` parameter in `implement_func`.
diff --git a/pint/testsuite/test_numpy.py b/pint/testsuite/test_numpy.py
index 8fd39d8..399886d 100644
--- a/pint/testsuite/test_numpy.py
+++ b/pint/testsuite/test_numpy.py
@@ -290,6 +290,9 @@ class TestNumpyMathematicalFunctions(TestNumpyMethods):
         self.assertRaises(DimensionalityError, np.cumproduct, self.q)
         self.assertQuantityEqual(np.cumprod(self.q / self.ureg.m), [1, 2, 6, 24])
         self.assertQuantityEqual(np.cumproduct(self.q / self.ureg.m), [1, 2, 6, 24])
+        self.assertQuantityEqual(
+            np.cumprod(self.q / self.ureg.m, axis=1), [[1, 2], [3, 12]]
+        )
 
     @helpers.requires_array_function_protocol()
     def test_nancumprod_numpy_func(self):

pip install -e .[test]
parse_log_pytest
requirements.txt
pytest pytest-cov pytest-mpl
3.7
pytest --no-header -rA --tb=line --color=no -p no:cacheprovider -W ignore::DeprecationWarning
certifi @ file:///croot/certifi_1671487769961/work/certifi
coverage==7.2.7
cycler==0.11.0
exceptiongroup==1.2.2
fonttools==4.38.0
importlib-metadata==6.7.0
importlib-resources==5.12.0
iniconfig==2.0.0
Jinja2==3.1.6
kiwisolver==1.4.5
MarkupSafe==2.1.5
matplotlib==3.5.3
numpy==1.21.6
packaging==24.0
Pillow==9.5.0
-e git+https://github.com/hgrecco/pint.git@4c3114c72bdae1e3c5c9f839659169eeb06071cd#egg=Pint
pluggy==1.2.0
pyparsing==3.1.4
pytest==7.4.4
pytest-cov==4.1.0
pytest-mpl==0.17.0
python-dateutil==2.9.0.post0
six==1.17.0
tomli==2.0.1
typing_extensions==4.7.1
zipp==3.15.0

name: pint
channels:
  - defaults
  - https://repo.anaconda.com/pkgs/main
  - https://repo.anaconda.com/pkgs/r
  - conda-forge
dependencies:
  - _libgcc_mutex=0.1=main
  - _openmp_mutex=5.1=1_gnu
  - ca-certificates=2025.2.25=h06a4308_0
  - certifi=2022.12.7=py37h06a4308_0
  - ld_impl_linux-64=2.40=h12ee557_0
  - libffi=3.4.4=h6a678d5_1
  - libgcc-ng=11.2.0=h1234567_1
  - libgomp=11.2.0=h1234567_1
  - libstdcxx-ng=11.2.0=h1234567_1
  - ncurses=6.4=h6a678d5_0
  - openssl=1.1.1w=h7f8727e_0
  - pip=22.3.1=py37h06a4308_0
  - python=3.7.16=h7a1cb2a_0
  - readline=8.2=h5eee18b_0
  - setuptools=65.6.3=py37h06a4308_0
  - sqlite=3.45.3=h5eee18b_0
  - tk=8.6.14=h39e8969_0
  - wheel=0.38.4=py37h06a4308_0
  - xz=5.6.4=h5eee18b_1
  - zlib=1.2.13=h5eee18b_1
  - pip:
      - coverage==7.2.7
      - cycler==0.11.0
      - exceptiongroup==1.2.2
      - fonttools==4.38.0
      - importlib-metadata==6.7.0
      - importlib-resources==5.12.0
      - iniconfig==2.0.0
      - jinja2==3.1.6
      - kiwisolver==1.4.5
      - markupsafe==2.1.5
      - matplotlib==3.5.3
      - numpy==1.21.6
      - packaging==24.0
      - pillow==9.5.0
      - pluggy==1.2.0
      - pyparsing==3.1.4
      - pytest==7.4.4
      - pytest-cov==4.1.0
      - pytest-mpl==0.17.0
      - python-dateutil==2.9.0.post0
      - six==1.17.0
      - tomli==2.0.1
      - typing-extensions==4.7.1
      - zipp==3.15.0
prefix: /opt/conda/envs/pint

pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_cumprod_numpy_func
pint/testsuite/test_numpy.py::TestNumpyArrayCreation::test_empty_like pint/testsuite/test_numpy.py::TestNumpyArrayCreation::test_full_like pint/testsuite/test_numpy.py::TestNumpyArrayCreation::test_ones_like pint/testsuite/test_numpy.py::TestNumpyArrayCreation::test_zeros_like pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_append pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_astype pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_atleast_1d pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_atleast_2d pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_atleast_3d pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_block pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_broadcast_to pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_column_stack pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_concatentate pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_dstack pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_expand_dims pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_flat pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_flatten pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_flip_numpy_func pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_hstack pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_item pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_moveaxis pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_ravel pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_ravel_numpy_func pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_reshape pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_rollaxis pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_squeeze pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_stack pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_swapaxes pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_transpose pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_transpose_numpy_func pint/testsuite/test_numpy.py::TestNumpyArrayManipulation::test_vstack pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_addition_with_incompatible_scalar pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_addition_with_scalar pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_cross pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_cumprod pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_diff pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_dot pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_dot_numpy_func pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_ediff1d pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_einsum pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_fix pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_gradient pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_nancumprod_numpy_func pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_nansum_numpy_func pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_power pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_prod pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_sum pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_sum_numpy_func pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_trapz pint/testsuite/test_numpy.py::TestNumpyMathematicalFunctions::test_unwrap pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_alen_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_argmax pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_argmax_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_argmin pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_argmin_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_argsort pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_argsort_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_average_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_clip pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_clip_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_comparisons pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_compress pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_compress_nep18 pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_conj pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_copy_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_copyto pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_count_nonzero_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_cumprod pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_cumsum pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_cumsum_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_diagonal pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_diagonal_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_equal pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_fabs pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_fill pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_getitem pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_insert pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_interp_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_isclose_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_isin pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_iterable pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_iterator pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_max pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_max_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_max_with_axis_arg pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_max_with_initial_arg pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_maximum pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_mean pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_mean_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_median_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_meshgrid_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_min pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_min_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_min_with_axis_arg pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_min_with_initial_arg pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_minimum pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nan_to_num_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanargmax_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanargmin_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nancumsum_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanmax pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanmean_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanmedian_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanmin pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanpercentile pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanstd_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nanvar_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_ndim_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nonzero pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_nonzero_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_percentile pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_pickle pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_prod pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_ptp pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_ptp_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_put pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_repeat pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_result_type_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_reversible_op pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_rot90 pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_round pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_round_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_searchsorted pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_searchsorted_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_setitem pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_shape pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_shape_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_sort pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_sort_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_std pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_std_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_take pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_tile pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_tolist pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_trace pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_trim_zeros_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_var pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_var_numpy_func pint/testsuite/test_numpy.py::TestNumpyUnclassified::test_where
pint/numpy_func.py
pint/numpy_func.py