save

Author: StaticObserver

setups/srpic/1d_polar_cap/pgen.hpp

@@ -82,7 +82,8 @@ namespace user {
     const real_t  drift_u_1, drift_u_2;
     const real_t  j0;
     InitFields<D> init_flds;
-    cdfTable cdf;
+    // const real_t e_min, gamma_emit, gamma_rad, gamma_pc, gamma_min, curvR;
+    // cdfTable cdf;
 
     inline PGen(const SimulationParams& p, const Metadomain<S, M>& m)
       : arch::ProblemGenerator<S, M>(p)
@@ -95,8 +96,14 @@ namespace user {
       , drift_u_1 { p.template get<real_t>("setup.drift_u_1") }
       , drift_u_2 { p.template get<real_t>("setup.drift_u_2") }
       , j0 { p.template get<real_t>("setup.j0") }
-      , init_flds { b0, TWO * FOUR * constant::PI * b0 * Omega * SQR(skindepth0) / larmor0 }
-      , cdf { "cdf_table.txt", "inverse_cdf_table.txt" } {}
+      , init_flds { b0, TWO * FOUR * constant::PI * b0 * Omega * SQR(skindepth0) / larmor0 } {}
+      // , e_min { p.template get<real_t>("setup.e_min") }
+      // , gamma_emit { p.template get<real_t>("setup.gamma_emit") }
+      // , gamma_rad { p.template get<real_t>("setup.gamma_rad") }
+      // , gamma_pc { p.template get<real_t>("setup.gamma_pc") }
+      // , gamma_min { p.template get<real_t>("setup.gamma_min") }
+      // , curvR { p.template get<real_t>("setup.curvR") }
+      // , cdf { "cdf_table.txt", "inverse_cdf_table.txt" } {}
 
     inline PGen() {}
 
@@ -137,9 +144,6 @@ namespace user {
   //       ONE + j0);
   //     }
 
-    void CustomPostStep(std::size_t, long double, Domain<S, M>& domain) {
-
-    }
   };
 
 } // namespace user

src/kernels/QED_process.hpp

@@ -362,7 +362,8 @@ namespace kernel::QED{
         const size_t                              n_steps;
 
         Inline auto Rho(const real_t x) const -> real_t{
-            return rho0 * (ONE + 0.8 * x / L);
+            // return rho0 * (ONE + 0.8 * x / L);
+            return rho0;
         }
 
 

src/kernels/tests/QED_process.cpp

@@ -115,8 +115,8 @@ auto main(int argc, char* argv[]) -> int {
                                                                        photon, 
                                                                        e_min, 
                                                                        gamma_emit, 
-                                                                       rho, 
                                                                        100,
+                                                                       rho, 
                                                                        N_phs,
                                                                        offsets,
                                                                        random_pool,