Fixed bug in regridding.regrid_from_weights() where zeros were bein…

…g returned if the input arrays were not contiguous.

regridding/_regrid/_regrid_from_weights.py

@@ -70,7 +70,7 @@ def regrid_from_weights(
  if values_output.shape != shape_output:
  raise ValueError(f"")
  values_output.fill(0)
- values_output_copy = values_output
+
  ndim_output = len(shape_output)
  axis_output = _util._normalize_axis(axis_output, ndim=ndim_output)
 
@@ -83,17 +83,26 @@ def regrid_from_weights(
  values_input = np.moveaxis(values_input, axis_input, axis_input_numba)
  values_output = np.moveaxis(values_output, axis_output, axis_output_numba)
 
+ shape_output_tmp = values_output.shape
+
  weights = numba.typed.List(weights.reshape(-1))
  values_input = values_input.reshape(-1, *shape_input_numba)
  values_output = values_output.reshape(-1, *shape_output_numba)
 
+ values_input = np.ascontiguousarray(values_input)
+ values_output = np.ascontiguousarray(values_output)
+
  _regrid_from_weights(
  weights=weights,
  values_input=values_input,
  values_output=values_output,
  )
 
- return values_output_copy
+ values_output = values_output.reshape(*shape_output_tmp)
+
+ values_output = np.moveaxis(values_output, axis_output_numba, axis_output)
+
+ return values_output
 
 
 @numba.njit()