def test_E_elektroni_ydin(self):
""" testataan elektronit-ytimet vuotovaikutusenergia """
energiat.get_V_hartree(
self.V_hartree2, self.elektroni_tiheys, tolerance=0.001)
e_elektroni_ydin = energiat.E_elektroni_ydin(
self.ydin_tiheys, self.V_hartree2)
self.assertAlmostEqual(e_elektroni_ydin, -157.07796, places=4)
def test_E_elektroni_elektroni(self):
""" testataan elektronien itseisvuorovaikutusenergia """
energiat.get_V_hartree(
self.V_hartree1, self.elektroni_tiheys, tolerance=0.001)
e_elektroni_elektroni = energiat.E_elektroni_elektroni(
self.elektroni_tiheys, self.V_hartree1)
self.assertAlmostEqual(e_elektroni_elektroni, 57.2074, places=4)
def test_E_tot(self):
""" testataan kokonaisenergia """
energiat.get_V_hartree(
self.V_hartree3, self.elektroni_tiheys, tolerance=0.001)
e_tot = energiat.E_T(self.elektroni_tiheys)+\
energiat.E_vaihtokorrelaatio(self.elektroni_tiheys) + \
energiat.E_elektroni_elektroni(
self.elektroni_tiheys, self.V_hartree3) + \
energiat.E_elektroni_ydin(
self.ydin_tiheys, self.V_hartree3)
self.assertAlmostEqual(e_tot, 1177.64253, places=4)
