Exemple #1
0
def test_behavior2():
    x = Symbol('x')
    p = Wild('p')

    e = Rational(6)
    assert e.match(2*p) == {p: 3}

    e = 3*x + 3 + 6/x
    a = Wild('a', exclude = [x])
    assert e.expand().match(a*x**2 + a*x + 2*a) == None
    assert e.expand().match(p*x**2 + p*x + 2*p) == {p: 3/x}
Exemple #2
0
def test_behavior2():
    x = Symbol('x')
    p = Wild('p')

    e = Rational(6)
    assert e.match(2 * p) == {p: 3}

    e = 3 * x + 3 + 6 / x
    a = Wild('a', exclude=[x])
    assert e.expand().match(a * x**2 + a * x + 2 * a) == None
    assert e.expand().match(p * x**2 + p * x + 2 * p) == {p: 3 / x}
Exemple #3
0
T = sum(get_ccsd_t_operators())
T_2 = sum(get_ccsd_t_operators())
T_t = sum(get_ccsd_t_operators(ast_symb="t(t)"))
T_t_2 = sum(get_ccsd_t_operators(ast_symb="t(t)"))
L = sum(get_ccsd_lambda_operators())
L_t = sum(get_ccsd_lambda_operators(ast_symb="l(t)"))

tilde_t_eq = Rational(1, 1)
tilde_t_eq += eval_equation(-L_t * T_t)
tilde_t_eq += eval_equation(L_t * T)
tilde_t_eq += eval_equation(-L_t * T_t * T)
tilde_t_eq += eval_equation(Rational(1, 2) * L_t * T_t * T_t_2)
tilde_t_eq += eval_equation(Rational(1, 2) * L_t * T * T_2)

tilde_t_eq = tilde_t_eq.expand()
tilde_t_eq = evaluate_deltas(tilde_t_eq)
tilde_t_eq = substitute_dummies(tilde_t_eq, **sub_kwargs)

tilde_eq = Rational(1, 1)
tilde_eq += eval_equation(-L * T)
tilde_eq += eval_equation(L * T_t)
tilde_eq += eval_equation(-L * T * T_t)
tilde_eq += eval_equation(Rational(1, 2) * L * T * T_2)
tilde_eq += eval_equation(Rational(1, 2) * L * T_t * T_t_2)

tilde_eq = tilde_eq.expand()
tilde_eq = evaluate_deltas(tilde_eq)
tilde_eq = substitute_dummies(tilde_eq, **sub_kwargs)

print(latex(tilde_t_eq))