code: def densities_close(rho0: Density, rho1: Density,
                      tolerance: float = TOLERANCE) -> bool:
      
      return vectors_close(rho0.vec, rho1.vec, tolerance)
docs: Returns True if densities are almost identical.
  
      Closeness is measured with the metric Fubini-Study angle.