[WIP]Implementation of Bremsstrahlung collisions

Source/Particles/Collision/BinaryCollision/BinaryCollision.H

@@ -11,6 +11,7 @@
 #include "Particles/Collision/BinaryCollision/Coulomb/ComputeTemperature.H"
 #include "Particles/Collision/BinaryCollision/DSMC/DSMCFunc.H"
 #include "Particles/Collision/BinaryCollision/NuclearFusion/NuclearFusionFunc.H"
+#include "Particles/Collision/BinaryCollision/Bremsstrahlung/BremsstrahlungFunc.H"
 #include "Particles/Collision/BinaryCollision/ParticleCreationFunc.H"
 #include "Particles/Collision/BinaryCollision/ShuffleFisherYates.H"
 #include "Particles/Collision/CollisionBase.H"
@@ -363,6 +364,10 @@ public:
             amrex::Gpu::DeviceVector<amrex::ParticleReal> pair_reaction_weight(n_total_pairs);
             amrex::ParticleReal* AMREX_RESTRICT p_pair_reaction_weight =
                                                                pair_reaction_weight.dataPtr();
+            // Extra data needed when products are created
+            int const n_product_data = (binary_collision_functor.m_need_product_data ? n_total_pairs : 0);
+            amrex::Gpu::DeviceVector<amrex::ParticleReal> product_data(n_product_data);
+            amrex::ParticleReal* AMREX_RESTRICT p_product_data = product_data.dataPtr();
             /*
               End of calculations only required when creating product particles
             */
@@ -458,7 +463,7 @@ public:
 
                     // Call the function in order to perform collisions
                     // If there are product species, mask, p_pair_indices_1/2, and
-                    // p_pair_reaction_weight are filled here
+                    // p_pair_reaction_weight and p_product_data are filled here
                     binary_collision_functor(
                         cell_start_1, cell_half_1,
                         cell_half_1, cell_stop_1,
@@ -467,7 +472,7 @@ public:
                         n1, n1, T1, T1,
                         q1, q1, m1, m1, dt, dV*volume_factor(i_cell), coll_idx,
                         cell_start_pair, p_mask, p_pair_indices_1, p_pair_indices_2,
-                        p_pair_reaction_weight, engine);
+                        p_pair_reaction_weight, p_product_data, engine);
                 }
             );
 
@@ -481,7 +486,7 @@ public:
                                                     products_mass, p_mask, products_np,
                                                     copy_species1, copy_species2,
                                                     p_pair_indices_1, p_pair_indices_2,
-                                                    p_pair_reaction_weight);
+                                                    p_pair_reaction_weight, p_product_data);
 
             for (int i = 0; i < n_product_species; i++)
             {
@@ -589,6 +594,10 @@ public:
             amrex::Gpu::DeviceVector<amrex::ParticleReal> pair_reaction_weight(n_total_pairs);
             amrex::ParticleReal* AMREX_RESTRICT p_pair_reaction_weight =
                                                                pair_reaction_weight.dataPtr();
+            // Extra data needed when products are created
+            int const n_product_data = (binary_collision_functor.m_need_product_data ? n_total_pairs : 0);
+            amrex::Gpu::DeviceVector<amrex::ParticleReal> product_data(n_product_data);
+            amrex::ParticleReal* AMREX_RESTRICT p_product_data = product_data.dataPtr();
             /*
               End of calculations only required when creating product particles
             */
@@ -717,15 +726,15 @@ public:
 
                     // Call the function in order to perform collisions
                     // If there are product species, p_mask, p_pair_indices_1/2, and
-                    // p_pair_reaction_weight are filled here
+                    // p_pair_reaction_weight and p_product_data are filled here
                     binary_collision_functor(
                         cell_start_1, cell_stop_1, cell_start_2, cell_stop_2,
                         indices_1, indices_2,
                         soa_1, soa_2, get_position_1, get_position_2,
                         n1, n2, T1, T2,
                         q1, q2, m1, m2, dt, dV*volume_factor(i_cell), coll_idx,
                         cell_start_pair, p_mask, p_pair_indices_1, p_pair_indices_2,
-                        p_pair_reaction_weight, engine);
+                        p_pair_reaction_weight, p_product_data, engine);
                 }
             );
 
@@ -739,7 +748,7 @@ public:
                                                     products_mass, p_mask, products_np,
                                                     copy_species1, copy_species2,
                                                     p_pair_indices_1, p_pair_indices_2,
-                                                    p_pair_reaction_weight);
+                                                    p_pair_reaction_weight, p_product_data);
 
             for (int i = 0; i < n_product_species; i++)
             {

Source/Particles/Collision/BinaryCollision/BinaryCollisionUtils.H

@@ -19,6 +19,7 @@ enum struct CollisionType { DeuteriumTritiumToNeutronHeliumFusion,
                             DeuteriumDeuteriumToNeutronHeliumFusion,
                             DeuteriumHeliumToProtonHeliumFusion,
                             ProtonBoronToAlphasFusion,
+                            Bremsstrahlung,
                             DSMC,
                             PairwiseCoulomb,
                             Undefined };

Source/Particles/Collision/BinaryCollision/BinaryCollisionUtils.cpp

@@ -31,6 +31,9 @@ namespace BinaryCollisionUtils{
             const NuclearFusionType fusion_type = get_nuclear_fusion_type(collision_name, mypc);
             return nuclear_fusion_type_to_collision_type(fusion_type);
         }
+        else if (type == "bremsstrahlung") {
+            return CollisionType::Bremsstrahlung;
+        }
         else if (type == "dsmc") {
             return CollisionType::DSMC;
         }