docs: def pauli_expansion(
      val: Any,
      *,
      default: Union[value.LinearDict[str], TDefault]
          = RaiseTypeErrorIfNotProvided,
      atol: float = 1e-9
  ) -> Union[value.LinearDict[str], TDefault]:
      """Returns coefficients of the expansion of val in the Pauli basis.
  
      Args:
          val: The value whose Pauli expansion is to returned.
          default: Determines what happens when `val` does not have methods that
              allow Pauli expansion to be obtained (see below). If set, the value
              is returned in that case. Otherwise, TypeError is raised.
          atol: Ignore coefficients whose absolute value is smaller than this.
  
      Returns:
          If `val` has a _pauli_expansion_ method, then its result is returned.
          Otherwise, if `val` has a small unitary then that unitary is expanded
          in the Pauli basis and coefficients are returned. Otherwise, if default
          is set to None or other value then default is returned. Otherwise,
          TypeError is raised.
  
      Raises:
          TypeError if `val` has none of the methods necessary to obtain its Pauli
          expansion and no default value has been provided.
      """
      method = getattr(val, '_pauli_expansion_', None)
      expansion = NotImplemented if method is None else method()
  
      if expansion is not NotImplemented:
          return expansion.clean(atol=atol)
  
      matrix = unitary(val, default=None)
      if matrix is None:
          if default is RaiseTypeErrorIfNotProvided:
              raise TypeError('No Pauli expansion for object {} of type {}'
                      .format(val, type(val)))
          return default
  
      num_qubits = matrix.shape[0].bit_length() - 1
      basis = operator_spaces.kron_bases(operator_spaces.PAULI_BASIS,
                                         repeat=num_qubits)
  
      expansion = operator_spaces.expand_matrix_in_orthogonal_basis(matrix, basis)
      return expansion.clean(atol=atol)