From b1218aa1e061c632ccc64eafe6968c7a302707d7 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?R=C3=A9mi=20Delaporte-Mathurin?= <40028739+RemDelaporteMathurin@users.noreply.github.com> Date: Mon, 15 Jan 2024 23:54:06 +0000 Subject: [PATCH] fixed formulation test --- test/unit/test_mobile.py | 21 +++++++++++++++++---- 1 file changed, 17 insertions(+), 4 deletions(-) diff --git a/test/unit/test_mobile.py b/test/unit/test_mobile.py index 90c2fdab3..9f9128b35 100644 --- a/test/unit/test_mobile.py +++ b/test/unit/test_mobile.py @@ -159,13 +159,26 @@ def test_with_traps_transient(self): c_0_n = my_mobile.previous_solution expected_form = ((c_0 - c_0_n) / self.dt.value) * v * self.my_mesh.dx(1) expected_form += f.dot(D * f.grad(c_0), f.grad(v)) * self.my_mesh.dx(1) + form_trapping_expected = 0 for trap in my_traps.traps: - expected_form += ( - ((trap.solution - trap.previous_solution) / self.dt.value) + form_trapping_expected += ( + ( + -trap.k_0 + * f.exp(-trap.E_k / festim.k_B / self.my_temp.T) + * c_0 + * (trap.density[0] - trap.solution) + ) * v - * self.my_mesh.dx + * self.my_mesh.dx(1) ) - + form_trapping_expected += ( + trap.p_0 + * f.exp(-trap.E_p / festim.k_B / self.my_temp.T) + * trap.solution + * v + * self.my_mesh.dx(1) + ) + expected_form += -form_trapping_expected print("expected F:") print(expected_form) print("produced F:")