docs: def check_semidefinite_positiveness(A):
      """Check if ``A`` is a semi-definite positive matrix.
  
      Args:
          A (array_like): Matrix.
  
      Returns:
          bool: ``True`` if ``A`` is definite positive; ``False`` otherwise.
      """
      B = empty_like(A)
      B[:] = A
      B[diag_indices_from(B)] += sqrt(finfo(float).eps)
      try:
          cholesky(B)
      except LinAlgError:
          return False
      return True