xsuite 0.40.2__cp312-cp312-macosx_11_0_arm64.whl → 0.41.1__cp312-cp312-macosx_11_0_arm64.whl

xsuite/lib/non_tracking_kernels.c

@@ -5755,156 +5755,156 @@ if (_sin_rot_s > -2.) {
 //            #endif                                                     //only_for_context cpu_openmp
         }
 
-#ifndef XOBJ_TYPEDEF_RandomExponentialData
-#define XOBJ_TYPEDEF_RandomExponentialData
-typedef   struct RandomExponentialData_s * RandomExponentialData;
- static inline RandomExponentialData RandomExponentialData_getp(RandomExponentialData restrict  obj){
+#ifndef XOBJ_TYPEDEF_RandomNormalData
+#define XOBJ_TYPEDEF_RandomNormalData
+typedef   struct RandomNormalData_s * RandomNormalData;
+ static inline RandomNormalData RandomNormalData_getp(RandomNormalData restrict  obj){
   int64_t offset=0;
-  return (RandomExponentialData)(( char*) obj+offset);
+  return (RandomNormalData)(( char*) obj+offset);
 }
- static inline uint8_t RandomExponentialData_get__dummy(const RandomExponentialData restrict  obj){
+ static inline uint8_t RandomNormalData_get__dummy(const RandomNormalData restrict  obj){
   int64_t offset=0;
   return *(( uint8_t*) obj+offset);
 }
- static inline void RandomExponentialData_set__dummy(RandomExponentialData restrict  obj, uint8_t value){
+ static inline void RandomNormalData_set__dummy(RandomNormalData restrict  obj, uint8_t value){
   int64_t offset=0;
   *(( uint8_t*) obj+offset)=value;
 }
- static inline  uint8_t* RandomExponentialData_getp__dummy(RandomExponentialData restrict  obj){
+ static inline  uint8_t* RandomNormalData_getp__dummy(RandomNormalData restrict  obj){
   int64_t offset=0;
   return ( uint8_t*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__sin_rot_s(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__sin_rot_s(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=8;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__sin_rot_s(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__sin_rot_s(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=8;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__sin_rot_s(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__sin_rot_s(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=8;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__cos_rot_s(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__cos_rot_s(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=16;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__cos_rot_s(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__cos_rot_s(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=16;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__cos_rot_s(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__cos_rot_s(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=16;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__shift_x(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__shift_x(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=24;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__shift_x(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__shift_x(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=24;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__shift_x(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__shift_x(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=24;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__shift_y(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__shift_y(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=32;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__shift_y(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__shift_y(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=32;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__shift_y(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__shift_y(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=32;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__shift_s(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__shift_s(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=40;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__shift_s(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__shift_s(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=40;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__shift_s(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__shift_s(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=40;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__rot_x_rad(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__rot_x_rad(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=48;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__rot_x_rad(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__rot_x_rad(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=48;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__rot_x_rad(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__rot_x_rad(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=48;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__rot_y_rad(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__rot_y_rad(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=56;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__rot_y_rad(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__rot_y_rad(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=56;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__rot_y_rad(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__rot_y_rad(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=56;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get__rot_s_rad_no_frame(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get__rot_s_rad_no_frame(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=64;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set__rot_s_rad_no_frame(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set__rot_s_rad_no_frame(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=64;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp__rot_s_rad_no_frame(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp__rot_s_rad_no_frame(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=64;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomExponentialData_get_rot_shift_anchor(const RandomExponentialData restrict  obj){
+ static inline double RandomNormalData_get_rot_shift_anchor(const RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=72;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomExponentialData_set_rot_shift_anchor(RandomExponentialData restrict  obj, double value){
+ static inline void RandomNormalData_set_rot_shift_anchor(RandomNormalData restrict  obj, double value){
   int64_t offset=0;
   offset+=72;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomExponentialData_getp_rot_shift_anchor(RandomExponentialData restrict  obj){
+ static inline  double* RandomNormalData_getp_rot_shift_anchor(RandomNormalData restrict  obj){
   int64_t offset=0;
   offset+=72;
   return ( double*)(( char*) obj+offset);
@@ -6119,11 +6119,11 @@ void Drift_single_particle(LocalParticle* part, double length){
 
 #endif // XSTUITE_TRACK_FLAGS_H
 
-#include <random/random_src/exponential.h>
+#include <random/random_src/normal.h>
 
              
-            void sample_exp(
-               RandomExponentialData el,
+            void sample_gauss(
+               RandomNormalData el,
 
                              ParticlesData particles,
     double* samples,  int64_t n_samples_per_seed, 
@@ -6171,7 +6171,7 @@ void Drift_single_particle(LocalParticle* part, double length){
                 if (part_id<part_capacity){
                     Particles_to_LocalParticle(particles, &lpart, part_id, end_id);
                     if (check_is_active(&lpart)>0){
-              RandomExponential_sample(el, &lpart, samples, n_samples_per_seed);
+              RandomNormal_sample(el, &lpart, samples, n_samples_per_seed);
 
                     }
                     if (check_is_active(&lpart)>0 && flag_increment_at_element){
@@ -6190,219 +6190,159 @@ void Drift_single_particle(LocalParticle* part, double length){
 //            #endif                                                     //only_for_context cpu_openmp
         }
 
-#ifndef XOBJ_TYPEDEF_RandomRutherfordData
-#define XOBJ_TYPEDEF_RandomRutherfordData
-typedef   struct RandomRutherfordData_s * RandomRutherfordData;
- static inline RandomRutherfordData RandomRutherfordData_getp(RandomRutherfordData restrict  obj){
+#ifndef XOBJ_TYPEDEF_RandomExponentialData
+#define XOBJ_TYPEDEF_RandomExponentialData
+typedef   struct RandomExponentialData_s * RandomExponentialData;
+ static inline RandomExponentialData RandomExponentialData_getp(RandomExponentialData restrict  obj){
   int64_t offset=0;
-  return (RandomRutherfordData)(( char*) obj+offset);
+  return (RandomExponentialData)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get_lower_val(const RandomRutherfordData restrict  obj){
+ static inline uint8_t RandomExponentialData_get__dummy(const RandomExponentialData restrict  obj){
   int64_t offset=0;
-  return *( double*)(( char*) obj+offset);
+  return *(( uint8_t*) obj+offset);
 }
- static inline void RandomRutherfordData_set_lower_val(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__dummy(RandomExponentialData restrict  obj, uint8_t value){
   int64_t offset=0;
-  *( double*)(( char*) obj+offset)=value;
+  *(( uint8_t*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp_lower_val(RandomRutherfordData restrict  obj){
+ static inline  uint8_t* RandomExponentialData_getp__dummy(RandomExponentialData restrict  obj){
   int64_t offset=0;
-  return ( double*)(( char*) obj+offset);
+  return ( uint8_t*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get_upper_val(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__sin_rot_s(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=8;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set_upper_val(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__sin_rot_s(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=8;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp_upper_val(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__sin_rot_s(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=8;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get_A(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__cos_rot_s(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=16;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set_A(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__cos_rot_s(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=16;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp_A(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__cos_rot_s(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=16;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get_B(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__shift_x(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=24;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set_B(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__shift_x(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=24;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp_B(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__shift_x(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=24;
   return ( double*)(( char*) obj+offset);
 }
- static inline int8_t RandomRutherfordData_get_Newton_iterations(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__shift_y(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=32;
-  return *(( int8_t*) obj+offset);
+  return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set_Newton_iterations(RandomRutherfordData restrict  obj, int8_t value){
+ static inline void RandomExponentialData_set__shift_y(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=32;
-  *(( int8_t*) obj+offset)=value;
+  *( double*)(( char*) obj+offset)=value;
 }
- static inline  int8_t* RandomRutherfordData_getp_Newton_iterations(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__shift_y(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=32;
-  return ( int8_t*)(( char*) obj+offset);
+  return ( double*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get__sin_rot_s(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__shift_s(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=40;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set__sin_rot_s(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__shift_s(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=40;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp__sin_rot_s(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__shift_s(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=40;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get__cos_rot_s(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__rot_x_rad(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=48;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set__cos_rot_s(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__rot_x_rad(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=48;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp__cos_rot_s(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__rot_x_rad(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=48;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get__shift_x(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__rot_y_rad(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=56;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set__shift_x(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__rot_y_rad(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=56;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp__shift_x(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__rot_y_rad(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=56;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get__shift_y(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get__rot_s_rad_no_frame(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=64;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set__shift_y(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set__rot_s_rad_no_frame(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=64;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp__shift_y(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp__rot_s_rad_no_frame(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=64;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomRutherfordData_get__shift_s(const RandomRutherfordData restrict  obj){
+ static inline double RandomExponentialData_get_rot_shift_anchor(const RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=72;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomRutherfordData_set__shift_s(RandomRutherfordData restrict  obj, double value){
+ static inline void RandomExponentialData_set_rot_shift_anchor(RandomExponentialData restrict  obj, double value){
   int64_t offset=0;
   offset+=72;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomRutherfordData_getp__shift_s(RandomRutherfordData restrict  obj){
+ static inline  double* RandomExponentialData_getp_rot_shift_anchor(RandomExponentialData restrict  obj){
   int64_t offset=0;
   offset+=72;
   return ( double*)(( char*) obj+offset);
-}
- static inline double RandomRutherfordData_get__rot_x_rad(const RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=80;
-  return *( double*)(( char*) obj+offset);
-}
- static inline void RandomRutherfordData_set__rot_x_rad(RandomRutherfordData restrict  obj, double value){
-  int64_t offset=0;
-  offset+=80;
-  *( double*)(( char*) obj+offset)=value;
-}
- static inline  double* RandomRutherfordData_getp__rot_x_rad(RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=80;
-  return ( double*)(( char*) obj+offset);
-}
- static inline double RandomRutherfordData_get__rot_y_rad(const RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=88;
-  return *( double*)(( char*) obj+offset);
-}
- static inline void RandomRutherfordData_set__rot_y_rad(RandomRutherfordData restrict  obj, double value){
-  int64_t offset=0;
-  offset+=88;
-  *( double*)(( char*) obj+offset)=value;
-}
- static inline  double* RandomRutherfordData_getp__rot_y_rad(RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=88;
-  return ( double*)(( char*) obj+offset);
-}
- static inline double RandomRutherfordData_get__rot_s_rad_no_frame(const RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=96;
-  return *( double*)(( char*) obj+offset);
-}
- static inline void RandomRutherfordData_set__rot_s_rad_no_frame(RandomRutherfordData restrict  obj, double value){
-  int64_t offset=0;
-  offset+=96;
-  *( double*)(( char*) obj+offset)=value;
-}
- static inline  double* RandomRutherfordData_getp__rot_s_rad_no_frame(RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=96;
-  return ( double*)(( char*) obj+offset);
-}
- static inline double RandomRutherfordData_get_rot_shift_anchor(const RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=104;
-  return *( double*)(( char*) obj+offset);
-}
- static inline void RandomRutherfordData_set_rot_shift_anchor(RandomRutherfordData restrict  obj, double value){
-  int64_t offset=0;
-  offset+=104;
-  *( double*)(( char*) obj+offset)=value;
-}
- static inline  double* RandomRutherfordData_getp_rot_shift_anchor(RandomRutherfordData restrict  obj){
-  int64_t offset=0;
-  offset+=104;
-  return ( double*)(( char*) obj+offset);
 }
 #endif
 // copyright ############################### //
@@ -6614,11 +6554,11 @@ void Drift_single_particle(LocalParticle* part, double length){
 
 #endif // XSTUITE_TRACK_FLAGS_H
 
-#include <random/random_src/rutherford.h>
+#include <random/random_src/exponential.h>
 
              
-            void sample_ruth(
-               RandomRutherfordData el,
+            void sample_exp(
+               RandomExponentialData el,
 
                              ParticlesData particles,
     double* samples,  int64_t n_samples_per_seed, 
@@ -6666,7 +6606,7 @@ void Drift_single_particle(LocalParticle* part, double length){
                 if (part_id<part_capacity){
                     Particles_to_LocalParticle(particles, &lpart, part_id, end_id);
                     if (check_is_active(&lpart)>0){
-              RandomRutherford_sample(el, &lpart, samples, n_samples_per_seed);
+              RandomExponential_sample(el, &lpart, samples, n_samples_per_seed);
 
                     }
                     if (check_is_active(&lpart)>0 && flag_increment_at_element){
@@ -6685,159 +6625,219 @@ void Drift_single_particle(LocalParticle* part, double length){
 //            #endif                                                     //only_for_context cpu_openmp
         }
 
-#ifndef XOBJ_TYPEDEF_RandomNormalData
-#define XOBJ_TYPEDEF_RandomNormalData
-typedef   struct RandomNormalData_s * RandomNormalData;
- static inline RandomNormalData RandomNormalData_getp(RandomNormalData restrict  obj){
+#ifndef XOBJ_TYPEDEF_RandomRutherfordData
+#define XOBJ_TYPEDEF_RandomRutherfordData
+typedef   struct RandomRutherfordData_s * RandomRutherfordData;
+ static inline RandomRutherfordData RandomRutherfordData_getp(RandomRutherfordData restrict  obj){
   int64_t offset=0;
-  return (RandomNormalData)(( char*) obj+offset);
+  return (RandomRutherfordData)(( char*) obj+offset);
 }
- static inline uint8_t RandomNormalData_get__dummy(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get_lower_val(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
-  return *(( uint8_t*) obj+offset);
+  return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__dummy(RandomNormalData restrict  obj, uint8_t value){
+ static inline void RandomRutherfordData_set_lower_val(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
-  *(( uint8_t*) obj+offset)=value;
+  *( double*)(( char*) obj+offset)=value;
 }
- static inline  uint8_t* RandomNormalData_getp__dummy(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp_lower_val(RandomRutherfordData restrict  obj){
   int64_t offset=0;
-  return ( uint8_t*)(( char*) obj+offset);
+  return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__sin_rot_s(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get_upper_val(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=8;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__sin_rot_s(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set_upper_val(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=8;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__sin_rot_s(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp_upper_val(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=8;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__cos_rot_s(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get_A(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=16;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__cos_rot_s(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set_A(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=16;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__cos_rot_s(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp_A(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=16;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__shift_x(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get_B(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=24;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__shift_x(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set_B(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=24;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__shift_x(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp_B(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=24;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__shift_y(const RandomNormalData restrict  obj){
+ static inline int8_t RandomRutherfordData_get_Newton_iterations(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=32;
-  return *( double*)(( char*) obj+offset);
+  return *(( int8_t*) obj+offset);
 }
- static inline void RandomNormalData_set__shift_y(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set_Newton_iterations(RandomRutherfordData restrict  obj, int8_t value){
   int64_t offset=0;
   offset+=32;
-  *( double*)(( char*) obj+offset)=value;
+  *(( int8_t*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__shift_y(RandomNormalData restrict  obj){
+ static inline  int8_t* RandomRutherfordData_getp_Newton_iterations(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=32;
-  return ( double*)(( char*) obj+offset);
+  return ( int8_t*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__shift_s(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get__sin_rot_s(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=40;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__shift_s(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set__sin_rot_s(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=40;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__shift_s(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp__sin_rot_s(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=40;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__rot_x_rad(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get__cos_rot_s(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=48;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__rot_x_rad(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set__cos_rot_s(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=48;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__rot_x_rad(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp__cos_rot_s(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=48;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__rot_y_rad(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get__shift_x(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=56;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__rot_y_rad(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set__shift_x(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=56;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__rot_y_rad(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp__shift_x(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=56;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get__rot_s_rad_no_frame(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get__shift_y(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=64;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set__rot_s_rad_no_frame(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set__shift_y(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=64;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp__rot_s_rad_no_frame(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp__shift_y(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=64;
   return ( double*)(( char*) obj+offset);
 }
- static inline double RandomNormalData_get_rot_shift_anchor(const RandomNormalData restrict  obj){
+ static inline double RandomRutherfordData_get__shift_s(const RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=72;
   return *( double*)(( char*) obj+offset);
 }
- static inline void RandomNormalData_set_rot_shift_anchor(RandomNormalData restrict  obj, double value){
+ static inline void RandomRutherfordData_set__shift_s(RandomRutherfordData restrict  obj, double value){
   int64_t offset=0;
   offset+=72;
   *( double*)(( char*) obj+offset)=value;
 }
- static inline  double* RandomNormalData_getp_rot_shift_anchor(RandomNormalData restrict  obj){
+ static inline  double* RandomRutherfordData_getp__shift_s(RandomRutherfordData restrict  obj){
   int64_t offset=0;
   offset+=72;
   return ( double*)(( char*) obj+offset);
+}
+ static inline double RandomRutherfordData_get__rot_x_rad(const RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=80;
+  return *( double*)(( char*) obj+offset);
+}
+ static inline void RandomRutherfordData_set__rot_x_rad(RandomRutherfordData restrict  obj, double value){
+  int64_t offset=0;
+  offset+=80;
+  *( double*)(( char*) obj+offset)=value;
+}
+ static inline  double* RandomRutherfordData_getp__rot_x_rad(RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=80;
+  return ( double*)(( char*) obj+offset);
+}
+ static inline double RandomRutherfordData_get__rot_y_rad(const RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=88;
+  return *( double*)(( char*) obj+offset);
+}
+ static inline void RandomRutherfordData_set__rot_y_rad(RandomRutherfordData restrict  obj, double value){
+  int64_t offset=0;
+  offset+=88;
+  *( double*)(( char*) obj+offset)=value;
+}
+ static inline  double* RandomRutherfordData_getp__rot_y_rad(RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=88;
+  return ( double*)(( char*) obj+offset);
+}
+ static inline double RandomRutherfordData_get__rot_s_rad_no_frame(const RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=96;
+  return *( double*)(( char*) obj+offset);
+}
+ static inline void RandomRutherfordData_set__rot_s_rad_no_frame(RandomRutherfordData restrict  obj, double value){
+  int64_t offset=0;
+  offset+=96;
+  *( double*)(( char*) obj+offset)=value;
+}
+ static inline  double* RandomRutherfordData_getp__rot_s_rad_no_frame(RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=96;
+  return ( double*)(( char*) obj+offset);
+}
+ static inline double RandomRutherfordData_get_rot_shift_anchor(const RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=104;
+  return *( double*)(( char*) obj+offset);
+}
+ static inline void RandomRutherfordData_set_rot_shift_anchor(RandomRutherfordData restrict  obj, double value){
+  int64_t offset=0;
+  offset+=104;
+  *( double*)(( char*) obj+offset)=value;
+}
+ static inline  double* RandomRutherfordData_getp_rot_shift_anchor(RandomRutherfordData restrict  obj){
+  int64_t offset=0;
+  offset+=104;
+  return ( double*)(( char*) obj+offset);
 }
 #endif
 // copyright ############################### //
@@ -7049,11 +7049,11 @@ void Drift_single_particle(LocalParticle* part, double length){
 
 #endif // XSTUITE_TRACK_FLAGS_H
 
-#include <random/random_src/normal.h>
+#include <random/random_src/rutherford.h>
 
              
-            void sample_gauss(
-               RandomNormalData el,
+            void sample_ruth(
+               RandomRutherfordData el,
 
                              ParticlesData particles,
     double* samples,  int64_t n_samples_per_seed, 
@@ -7101,7 +7101,7 @@ void Drift_single_particle(LocalParticle* part, double length){
                 if (part_id<part_capacity){
                     Particles_to_LocalParticle(particles, &lpart, part_id, end_id);
                     if (check_is_active(&lpart)>0){
-              RandomNormal_sample(el, &lpart, samples, n_samples_per_seed);
+              RandomRutherford_sample(el, &lpart, samples, n_samples_per_seed);
 
                     }
                     if (check_is_active(&lpart)>0 && flag_increment_at_element){