diff --git a/python/damask/_orientation.py b/python/damask/_orientation.py
index b28ba2aae..359d6c0a4 100644
--- a/python/damask/_orientation.py
+++ b/python/damask/_orientation.py
@@ -534,21 +534,23 @@ class Orientation(Rotation,Crystal):
 
         r_ = s[:,np.newaxis,...].misorientation(o[np.newaxis,:,...]) # type: ignore[index]
         _r = ~r_
+        shp = r_.shape[2:]
 
         forward = r_.in_disorientation_FZ
         reverse = _r.in_disorientation_FZ
         ok  = forward | reverse
-        ok &= (np.cumsum(ok.reshape((-1,)+blend),axis=0) == 1).reshape(ok.shape)
+        ok &= (np.cumsum(ok.reshape((-1,*shp)),axis=0) == 1).reshape(ok.shape)
         r = np.where(np.any((ok&forward)[...,np.newaxis],axis=(0,1),keepdims=True),
                      r_.quaternion,
                      _r.quaternion)
         loc  = np.where(ok)
         sort = 0 if len(loc) == 2 else np.lexsort(loc[:1:-1])
-        quat = r[ok][sort].reshape(blend+(4,))
+
+        quat = r[ok][sort].reshape((*shp,4))
 
         return (
                 (self.copy(rotation=quat),
-                 (np.vstack(loc[:2]).T)[sort].reshape(blend+(2,)))
+                 (np.vstack(loc[:2]).T)[sort].reshape((*shp,2)))
                 if return_operators else
                 self.copy(rotation=quat)
                )
diff --git a/python/damask/_rotation.py b/python/damask/_rotation.py
index f4f6a65f3..c955b2da0 100644
--- a/python/damask/_rotation.py
+++ b/python/damask/_rotation.py
@@ -616,6 +616,7 @@ class Rotation:
         ----------
         other : damask.Rotation
             Rotation to which the misorientation is computed.
+            Compatible innermost dimensions will blend.
 
         Returns
         -------
@@ -623,7 +624,7 @@ class Rotation:
             Misorientation.
 
         """
-        return other/self
+        return ~(self*~other)
 
 
     ################################################################################################