From 5f704f79b63dd9b418f5dad7683d77d9f80bad44 Mon Sep 17 00:00:00 2001
From: oscarxblanco <oblanco@cells.es>
Date: Fri, 20 Dec 2024 12:46:56 +0100
Subject: [PATCH] branchless if

---
 atintegrators/VariableThinMPolePass.c | 13 +++++++------
 1 file changed, 7 insertions(+), 6 deletions(-)

diff --git a/atintegrators/VariableThinMPolePass.c b/atintegrators/VariableThinMPolePass.c
index 40e768cf8..7347a464a 100644
--- a/atintegrators/VariableThinMPolePass.c
+++ b/atintegrators/VariableThinMPolePass.c
@@ -146,9 +146,9 @@ double get_pol(
         freq = elem->Frequency;
         ph = elem->Phase;
         sinval = sin(TWOPI * freq * t + ph);
-        // branchless if, to get
-        // the whole sin if whole_sin_limit is -1 or less
-        // the positive side of sin uf whole_sin_limit is 0.
+        // branchless if
+        // sinval >= wholelimit  -> sinval
+        // else                  -> 0
         ampt = ampt*sinval*(sinval >= whole_sin_limit);
         return ampt;
     case 1:
@@ -239,9 +239,10 @@ void VariableThinMPolePass(
     };
 
     // offset the time when applying the sin function
-    if (mode == 0){
-      time_in_this_mode = t0 * turn;
-    }
+    // branchless if
+    // mode == 0 -> time_in_this_mode = t0 *turn
+    // else      -> 0
+    time_in_this_mode = t0 * turn * (mode == 0);
 
     /* cycle over all particles */
     for (c = 0; c < num_particles; c++) {