def test_quadratic_phase(self, tol): """Test the quadratic phase symplectic transform.""" s = 0.543 S = quadratic_phase(s) # apply to a coherent state. P[x, p] -> [x, p+sx] alpha = 0.23+0.12j out = S @ np.array([alpha.real, alpha.imag])*np.sqrt(2*hbar) expected = np.array([alpha.real, alpha.imag+s*alpha.real])*np.sqrt(2*hbar) assert out == pytest.approx(expected, abs=tol)
def test_quadratic_phase(self): """Test the quadratic phase symplectic transform.""" self.logTestName() s = 0.543 S = quadratic_phase(s) # apply to a coherent state. P[x, p] -> [x, p+sx] alpha = 0.23+0.12j out = S @ np.array([alpha.real, alpha.imag])*np.sqrt(2*hbar) expected = np.array([alpha.real, alpha.imag+s*alpha.real])*np.sqrt(2*hbar) self.assertAllAlmostEqual(out, expected, delta=self.tol)