From 7fcf23d1795e530014cea61288c435114e8cc886 Mon Sep 17 00:00:00 2001 From: Angel Ferran Pousa Date: Wed, 15 Nov 2023 10:44:48 +0100 Subject: [PATCH] Remove phase jump correction --- wake_t/physics_models/laser/envelope_solver.py | 10 ---------- .../laser/envelope_solver_non_centered.py | 10 ---------- 2 files changed, 20 deletions(-) diff --git a/wake_t/physics_models/laser/envelope_solver.py b/wake_t/physics_models/laser/envelope_solver.py index d460e48..1bc4ab8 100644 --- a/wake_t/physics_models/laser/envelope_solver.py +++ b/wake_t/physics_models/laser/envelope_solver.py @@ -102,16 +102,6 @@ def evolve_envelope( d_theta1 = phases[j + 1] - phases[j] d_theta2 = phases[j + 2] - phases[j + 1] - # Prevent phase jumps bigger than 1.5*pi. - if d_theta1 < -1.5 * np.pi: - d_theta1 += 2 * np.pi - if d_theta2 < -1.5 * np.pi: - d_theta2 += 2 * np.pi - if d_theta1 > 1.5 * np.pi: - d_theta1 -= 2 * np.pi - if d_theta2 > 1.5 * np.pi: - d_theta2 -= 2 * np.pi - # Calculate D factor [Eq. (6)]. D_jkn = (1.5 * d_theta1 - 0.5 * d_theta2) * inv_dz diff --git a/wake_t/physics_models/laser/envelope_solver_non_centered.py b/wake_t/physics_models/laser/envelope_solver_non_centered.py index 5940190..d1892a3 100644 --- a/wake_t/physics_models/laser/envelope_solver_non_centered.py +++ b/wake_t/physics_models/laser/envelope_solver_non_centered.py @@ -102,16 +102,6 @@ def evolve_envelope_non_centered( d_theta1 = phases[j + 1] - phases[j] d_theta2 = phases[j + 2] - phases[j + 1] - # Prevent phase jumps bigger than 1.5*pi. - if d_theta1 < -1.5 * np.pi: - d_theta1 += 2 * np.pi - if d_theta2 < -1.5 * np.pi: - d_theta2 += 2 * np.pi - if d_theta1 > 1.5 * np.pi: - d_theta1 -= 2 * np.pi - if d_theta2 > 1.5 * np.pi: - d_theta2 -= 2 * np.pi - # Calculate D factor [Eq. (6)]. D_jkn = (1.5 * d_theta1 - 0.5 * d_theta2) * inv_dz