def test_normal_order():
a = BosonOp('a')
b = BosonOp('b')
c = FermionOp('c')
d = FermionOp('d')
assert normal_order(a * Dagger(a)) == Dagger(a) * a
assert normal_order(Dagger(a) * a) == Dagger(a) * a
assert normal_order(a * Dagger(a)**2) == Dagger(a)**2 * a
assert normal_order(c * Dagger(c)) == -Dagger(c) * c
assert normal_order(Dagger(c) * c) == Dagger(c) * c
assert normal_order(c * Dagger(c)**2) == Dagger(c)**2 * c
def test_normal_ordered_form():
a = BosonOp('a')
b = BosonOp('b')
c = FermionOp('c')
d = FermionOp('d')
assert normal_ordered_form(Dagger(a) * a) == Dagger(a) * a
assert normal_ordered_form(a * Dagger(a)) == 1 + Dagger(a) * a
assert normal_ordered_form(a ** 2 * Dagger(a)) == \
2 * a + Dagger(a) * a ** 2
assert normal_ordered_form(a ** 3 * Dagger(a)) == \
3 * a ** 2 + Dagger(a) * a ** 3
assert normal_ordered_form(Dagger(c) * c) == Dagger(c) * c
assert normal_ordered_form(c * Dagger(c)) == 1 - Dagger(c) * c
assert normal_ordered_form(c**2 * Dagger(c)) == Dagger(c) * c**2
assert normal_ordered_form(c ** 3 * Dagger(c)) == \
c ** 2 - Dagger(c) * c ** 3
def test_normal_order():
a = BosonOp("a")
c = FermionOp("c")
assert normal_order(a * Dagger(a)) == Dagger(a) * a
assert normal_order(Dagger(a) * a) == Dagger(a) * a
assert normal_order(a * Dagger(a)**2) == Dagger(a)**2 * a
assert normal_order(c * Dagger(c)) == -Dagger(c) * c
assert normal_order(Dagger(c) * c) == Dagger(c) * c
assert normal_order(c * Dagger(c)**2) == Dagger(c)**2 * c
def test_fermion_states():
c = FermionOp("c")
# Fock states
assert (FermionFockBra(0) * FermionFockKet(1)).doit() == 0
assert (FermionFockBra(1) * FermionFockKet(1)).doit() == 1
assert qapply(c * FermionFockKet(1)) == FermionFockKet(0)
assert qapply(c * FermionFockKet(0)) == 0
assert qapply(Dagger(c) * FermionFockKet(0)) == FermionFockKet(1)
assert qapply(Dagger(c) * FermionFockKet(1)) == 0
def test_normal_ordered_form():
a = BosonOp("a")
c = FermionOp("c")
assert normal_ordered_form(Dagger(a) * a) == Dagger(a) * a
assert normal_ordered_form(a * Dagger(a)) == 1 + Dagger(a) * a
assert normal_ordered_form(a**2 * Dagger(a)) == 2 * a + Dagger(a) * a**2
assert normal_ordered_form(a**3 * Dagger(a)) == 3 * a**2 + Dagger(a) * a**3
assert normal_ordered_form(Dagger(c) * c) == Dagger(c) * c
assert normal_ordered_form(c * Dagger(c)) == 1 - Dagger(c) * c
assert normal_ordered_form(c**2 * Dagger(c)) == Dagger(c) * c**2
assert normal_ordered_form(c**3 * Dagger(c)) == c**2 - Dagger(c) * c**3
def test_fermionoperator():
c = FermionOp('c')
d = FermionOp('d')
assert isinstance(c, FermionOp)
assert isinstance(Dagger(c), FermionOp)
assert c.is_annihilation
assert not Dagger(c).is_annihilation
assert FermionOp("c") == FermionOp("c")
assert FermionOp("c") != FermionOp("d")
assert FermionOp("c", True) != FermionOp("c", False)
assert AntiCommutator(c, Dagger(c)).doit() == 1
assert AntiCommutator(c, Dagger(d)).doit() == c * Dagger(d) + Dagger(d) * c