def test_consistency(self): a = array([[0., 1], [-1, 0]]) assert_array_almost_equal(expm(a), expm2(a)) assert_array_almost_equal(expm(a), expm3(a)) a = array([[1j, 1], [-1, -2j]]) assert_array_almost_equal(expm(a), expm2(a)) assert_array_almost_equal(expm(a), expm3(a))
def test_consistency(self): a = array([[0.,1],[-1,0]]) assert_array_almost_equal(expm(a), expm2(a)) assert_array_almost_equal(expm(a), expm3(a)) a = array([[1j,1],[-1,-2j]]) assert_array_almost_equal(expm(a), expm2(a)) assert_array_almost_equal(expm(a), expm3(a))
def test_consistency(self): with warnings.catch_warnings(): warnings.simplefilter("ignore", DeprecationWarning) a = array([[0.,1],[-1,0]]) assert_array_almost_equal(expm(a), expm2(a)) assert_array_almost_equal(expm(a), expm3(a)) a = array([[1j,1],[-1,-2j]]) assert_array_almost_equal(expm(a), expm2(a)) assert_array_almost_equal(expm(a), expm3(a))
def test_consistency(self): a = array([[0., 1], [-1, 0]]) b = array([[1j, 1], [-1, -2j]]) with suppress_warnings() as sup: sup.filter(DeprecationWarning, "`expm.` is deprecated") assert_array_almost_equal(expm(a), expm2(a)) assert_array_almost_equal(expm(a), expm3(a)) assert_array_almost_equal(expm(b), expm2(b)) assert_array_almost_equal(expm(b), expm3(b))
def test_zero(self): with warnings.catch_warnings(): warnings.simplefilter("ignore", DeprecationWarning) a = array([[0.,0],[0,0]]) assert_array_almost_equal(expm(a),[[1,0],[0,1]]) assert_array_almost_equal(expm2(a),[[1,0],[0,1]]) assert_array_almost_equal(expm3(a),[[1,0],[0,1]])
def test_zero(self): a = array([[0., 0], [0, 0]]) assert_array_almost_equal(expm(a), [[1, 0], [0, 1]]) with suppress_warnings() as sup: sup.filter(DeprecationWarning, "`expm.` is deprecated") assert_array_almost_equal(expm2(a), [[1, 0], [0, 1]]) assert_array_almost_equal(expm3(a), [[1, 0], [0, 1]])
def test_zero(self): a = array([[0.,0],[0,0]]) assert_array_almost_equal(expm(a),[[1,0],[0,1]]) assert_array_almost_equal(expm2(a),[[1,0],[0,1]]) assert_array_almost_equal(expm3(a),[[1,0],[0,1]])
def evolve_x1(x0,t): return x0 + expect(p[0] * HI * t + .5j * HI * (alpha1 - HI) * (expm2(-2j*H*t) - 1), u_p0)
def test_npmatrix(self): with warnings.catch_warnings(): warnings.simplefilter("ignore", DeprecationWarning) a = matrix([[3., 0], [0, -3.]]) assert_array_almost_equal(expm(a), expm2(a))
def test_zero(self): a = array([[0., 0], [0, 0]]) assert_array_almost_equal(expm(a), [[1, 0], [0, 1]]) assert_array_almost_equal(expm2(a), [[1, 0], [0, 1]]) assert_array_almost_equal(expm3(a), [[1, 0], [0, 1]])
def test_npmatrix(self): a = matrix([[3., 0], [0, -3.]]) with suppress_warnings() as sup: sup.filter(DeprecationWarning, "`expm2` is deprecated") assert_array_almost_equal(expm(a), expm2(a))
def test_npmatrix(self): with warnings.catch_warnings(): warnings.simplefilter("ignore", DeprecationWarning) a = matrix([[3.,0],[0,-3.]]) assert_array_almost_equal(expm(a), expm2(a))
def propagator(H,t): return expm2(-1j*H*t) #expm2: compute the matrix exponential using eigenvalue decomposition