Fixed error for zero beam sizes (mostly with zero current slices) for…

… the calculation of the space charge field.
ZeugAusHH · Jun 6, 2024 · d35b727 · d35b727
1 parent 9e8fb80
commit d35b727
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 0 deletions.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -2,6 +2,7 @@
 
 ### [4.6.6] - 2024/06/05
 - fixed: Running non 1-4-1 simulations (aka quiet loading) with zero current (e.g. step profile) causes NaN events in the particle phase for the calculation of the shot noise. It ripples through the simulation producing more NaN in other parameters.
+- fixed: Zero current slices caused a division by zero error in the calculation of the space charge field.
 
 ### [4.6.6] - 2024/04/18
 - fixed: changing elements via the lattice-namelist methods works again after ignoring any sequence.

diff --git a/src/Core/EFieldSolver.cpp b/src/Core/EFieldSolver.cpp
@@ -87,6 +87,10 @@ void EFieldSolver::longRange(Beam *beam, double gamma0, double aw) {
     for (int i=0; i<nsize;i++){
         work1[i]=beam->current[i];
         work2[i]=beam->getSize(i);
+        if (work2[i] <= 0) {  // check for zero current slices which have also no beam size and therefore will cause a NaN
+            work1[i] = 0;
+            work2[i] = 1;
+        }
     }
 
     // gathering information on full current and beam size profile