""" use the symetry of lineshape function in the exciton indices """
#D = CumulantExpr(expr)._leading_index(a)
#expr = D._getExpr()
A = Anorm.rewrite(gg)
norm = CumulantExpr(A)
""" use option large=T to calculate evaluate in T --> oo """
norm = norm.evaluate(large=T)
""" use the symetry of lineshape function in the exciton indices """
#D = CumulantExpr(expr)._leading_index(a)
#expr = D._getExpr()
expr = (expr-norm).simplify()
if verbatim:
print("Cumulant: ")
print(" ")
print(expr)
print(" ")
print("Foerster rate")
expr = expr.subs(e,d).simplify()
a1 = Symbol('a1')
b1 = Symbol('b1')
expr = expr.subs({d:b1,c:a1,b:a1,a:b1})
expr = expr.subs({a1:a,b1:b}).simplify()
print(expr)
A = A.rewrite(gg)
expr = CumulantExpr(A)
""" use option large=T to calculate evaluate in T --> oo """
expr = expr.evaluate(large=T)
""" use the symetry of lineshape function in the exciton indices """
#D = CumulantExpr(expr)._leading_index(a)
#expr = D._getExpr()
A = Anorm.rewrite(gg)
norm = CumulantExpr(A)
""" use option large=T to calculate evaluate in T --> oo """
norm = norm.evaluate(large=T)
""" use the symetry of lineshape function in the exciton indices """
#D = CumulantExpr(expr)._leading_index(a)
#expr = D._getExpr()
expr = (expr-norm).simplify()
if verbatim:
print("Cumulant: ")
print(" ")
print(expr)
print(" ")
print("In simplified form: ")
print(" ")
expr = expr.subs(c,b)
expr = expr.subs(d,b).simplify()
print(expr)
expr = expr.evaluate(large=T)
""" use the symetry of lineshape function in the exciton indices """
#D = CumulantExpr(expr)._leading_index(a)
#expr = D._getExpr()
A = Anorm.rewrite(gg)
norm = CumulantExpr(A)
""" use option large=T to calculate evaluate in T --> oo """
norm = norm.evaluate(large=T)
""" use the symetry of lineshape function in the exciton indices """
#D = CumulantExpr(expr)._leading_index(a)
#expr = D._getExpr()
expr = (expr - norm).simplify()
if verbatim:
print("Cumulant: ")
print(" ")
print(expr)
print(" ")
print("Foerster rate")
expr = expr.subs(e, d).simplify()
a1 = Symbol('a1')
b1 = Symbol('b1')
expr = expr.subs({d: b1, c: a1, b: a1, a: b1})
expr = expr.subs({a1: a, b1: b}).simplify()
print(expr)