docs: def orthogonal(*args) -> bool:
      """Determine if a set of arrays are orthogonal.
  
      Parameters
      ----------
      args : array-likes or array shapes
  
      Returns
      -------
      bool
          Array orthogonality condition.
      """
      for i, arg in enumerate(args):
          if hasattr(arg, "shape"):
              args[i] = arg.shape
      for s in zip(*args):
          if np.product(s) != max(s):
              return False
      return True