diff --git a/docs/reference.md b/docs/reference.md
index d8d6abc..92d92d0 100644
--- a/docs/reference.md
+++ b/docs/reference.md
@@ -26,6 +26,8 @@ Univariate
     vrgksuv_ensemble
     crps_quantile
 
+    dssuv_ensemble
+
 Multivariate
 --------------------
 
@@ -47,6 +49,8 @@ Multivariate
     owgksmv_ensemble
     vrgksmv_ensemble
 
+    dssmv_ensemble
+
 Parametric distributions forecasts
 ====================================
 .. autosummary::
diff --git a/scoringrules/core/utils.py b/scoringrules/core/utils.py
index 22b9f23..6aa1c13 100644
--- a/scoringrules/core/utils.py
+++ b/scoringrules/core/utils.py
@@ -45,7 +45,10 @@ def lazy_gufunc_wrapper_uv(func):
     def wrapper(*args):
         out = np.empty_like(args[0])
         func(*args, out)
-        return out
+        if out.shape == ():
+            return out.item()
+        else:
+            return out
 
     return wrapper
 
@@ -62,6 +65,9 @@ def lazy_gufunc_wrapper_mv(func):
     def wrapper(*args):
         out = np.empty_like(args[0][..., 0])
         func(*args, out)
-        return out
+        if out.shape == ():
+            return out.item()
+        else:
+            return out
 
     return wrapper