code: def phase2t(self, psi):
          
          def _deg(rads, domain_lower_limit):
              # Convert rads to degrees in [0, 360) domain
              degs = degrees(rads % (2*pi))
  
              # Convert to [domain_lower_limit, domain_lower_limit + 360) domain
              k = domain_lower_limit // 360
              degs += k * 360
              if degs < domain_lower_limit:
                  degs += 360
              return degs
  
          if self.delta > 0:
              degs = _deg(psi, domain_lower_limit=self.theta)
          else:
              degs = _deg(psi, domain_lower_limit=self.theta)
          return (degs - self.theta)/self.delta
docs: Given phase -pi < psi <= pi,
          returns the t value such that
          exp(1j*psi) = self.u1transform(self.point(t)).