Check for subclass of unyt_quantity instead of exact matching type.

unyt/array.py

@@ -1997,7 +1997,7 @@ def __array_ufunc__(self, ufunc, method, *inputs, **kwargs):
         elif out_arr.size == 1:
             out_arr = unyt_array(np.asarray(out_arr), unit)
         else:
-            if ret_class is unyt_quantity:
+            if issubclass(ret_class, unyt_quantity):
                 # This happens if you do ndarray * unyt_quantity.
                 # Explicitly casting to unyt_array avoids creating a
                 # unyt_quantity with size > 1