Debugging

Author: EZoni

Source/Particles/Collision/BinaryCollision/BinaryCollision.H

@@ -181,7 +181,9 @@ public:
         auto& species2 = mypc->GetParticleContainerFromName(m_species_names[1]);
 
         // In case of particle creation, create the necessary vectors
+        amrex::Print() << "Starting binary collisions between " << m_species_names[0] << " and " << m_species_names[1] << std::endl;
         const int n_product_species = m_product_species.size();
+        amrex::Print() << "n_product_species = " << n_product_species << std::endl;
         amrex::Vector<WarpXParticleContainer*> product_species_vector;
         amrex::Vector<SmartCopyFactory> copy_factory_species1;
         amrex::Vector<SmartCopyFactory> copy_factory_species2;
@@ -192,6 +194,14 @@ public:
             auto& product = mypc->GetParticleContainerFromName(m_product_species[i]);
             product.defineAllParticleTiles();
             product_species_vector.push_back(&product);
+            amrex::Print() << "BinaryCollision::doCollisions: preparing to create " << m_product_species[i] << std::endl;
+            amrex::Print() << "product.NumRealComps() = " << product.NumRealComps() << std::endl;
+            amrex::Print() << "product.NumRuntimeRealComps() = " << product.NumRuntimeRealComps() << std::endl;
+            amrex::Print() << "product.NStructReal = " << product.NStructReal << std::endl;
+            amrex::Print() << "product.NArrayReal = " << product.NArrayReal << std::endl;
+            for (auto const & comp : product.GetRealSoANames()) {
+                amrex::Print() << "product.GetRealSoANames(): " << comp << std::endl;
+            }
             // Although the copy factories are not explicitly reused past this point, we need to
             // store them in vectors so that the data that they own, which is used by the smart
             // copy functors, does not go out of scope at the end of this for loop.

Source/Particles/Collision/BinaryCollision/DSMC/SplitAndScatterFunc.H

@@ -286,6 +286,7 @@ public:
                 // create a copy of the first product species at the location of species 1
                 const auto reactant1_index = static_cast<int>(p_pair_indices_1[i]);
                 const auto product1_index = products_np_data[2] + with_product_p_offsets[i];
+                amrex::Print() << "Creating " << pc_products[2]->getSpeciesTypeName() << " particle at index " << product1_index << " from reactant index " << reactant1_index << std::endl;
                 copy_species1[2](soa_products_data[2], soa_1, reactant1_index,
                                 static_cast<int>(product1_index), engine);
                 // Set the weight of the new particle to p_pair_reaction_weight[i]
@@ -294,6 +295,7 @@ public:
                 // create a copy of the other product species at the location of species 2
                 const auto reactant2_index = static_cast<int>(p_pair_indices_2[i]);
                 const auto product2_index = products_np_data[3] + with_product_p_offsets[i];
+                amrex::Print() << "Creating " << pc_products[3]->getSpeciesTypeName() << " particle at index " << product2_index << " from reactant index " << reactant2_index << std::endl;
                 copy_species2[3](soa_products_data[3], soa_2, reactant2_index,
                                 static_cast<int>(product2_index), engine);
                 // Set the weight of the new particle to p_pair_reaction_weight[i]
@@ -509,6 +511,7 @@ public:
         {
             const int start_index = int(products_np[i]);
             const int stop_index  = int(products_np[i] + num_added_vec[i]);
+            amrex::Print() << "Calling ParticleCreation::DefaultInitializeRuntimeAttributes for species " << pc_products[i]->getSpeciesTypeName() << std::endl;
             ParticleCreation::DefaultInitializeRuntimeAttributes(*tile_products[i],
                                        0, 0,
                                        pc_products[i]->getUserRealAttribs(), pc_products[i]->getUserIntAttribs(),

Source/Particles/ParticleCreation/DefaultInitialization.H

@@ -137,6 +137,14 @@ void DefaultInitializeRuntimeAttributes (PTile& ptile,
 {
         using namespace amrex::literals;
 
+        amrex::Print() << "Starting ParticleCreation::DefaultInitializeRuntimeAttributes" << std::endl;
+        for (auto const & comp : particle_comps) {
+            amrex::Print() << "particle_comps: " << comp << std::endl;
+        }
+        for (auto const & comp: ptile.GetStructOfArrays().GetRealNames()) {
+            amrex::Print() << "ptile.GetStructOfArrays().GetRealNames(): " << comp << std::endl;
+        }
+
         // Preparing data needed for user defined attributes
         const auto n_user_real_attribs = static_cast<int>(user_real_attribs.size());
         const auto n_user_int_attribs = static_cast<int>(user_int_attribs.size());

Source/Particles/ParticleCreation/SmartCopy.H

@@ -50,9 +50,11 @@ struct SmartCopy
     {
         // initialize the real components
         for (int j = 0; j < DstData::NAR; ++j) {
+            amrex::Print() << "Initializing real component " << j << " for dst particle " << i_dst << " with policy " << static_cast<int>(m_policy_real[j]) << std::endl;
             dst.m_rdata[j][i_dst] = initializeRealValue(m_policy_real[j], engine);
         }
         for (int j = 0; j < dst.m_num_runtime_real; ++j) {
+            amrex::Print() << "Initializing runtime real component " << j << " for dst particle " << i_dst << " with policy " << static_cast<int>(m_policy_real[j+DstData::NAR]) << std::endl;
             dst.m_runtime_rdata[j][i_dst] = initializeRealValue(m_policy_real[j+DstData::NAR], engine);
         }
 

Source/Particles/PhysicalParticleContainer.cpp

@@ -351,6 +351,7 @@ PhysicalParticleContainer::PhysicalParticleContainer (AmrCore* amr_core, int isp
     }
 
     if (m_collisions_split_position_push) {
+        amrex::Print() << "Adding average momentum components for " << species_name << std::endl;
         // Add the average momentum components to deposit the current correctly
         // if the position push is split to perform collisions
         AddRealComp("ux_avg");
@@ -428,6 +429,8 @@ PhysicalParticleContainer::DefaultInitializeRuntimeAttributes (
     int n_external_attr_real,
     int n_external_attr_int)
 {
+    amrex::Print() << "Starting PhysicalParticleContainer::DefaultInitializeRuntimeAttributes for "
+                   << species_name << std::endl;
     ParticleCreation::DefaultInitializeRuntimeAttributes(pinned_tile,
                                        n_external_attr_real, n_external_attr_int,
                                        m_user_real_attribs, m_user_int_attribs,
@@ -455,6 +458,15 @@ PhysicalParticleContainer::Evolve (ablastr::fields::MultiFabRegister& fields,
     using ablastr::fields::Direction;
     using warpx::fields::FieldType;
 
+    amrex::Print() << "Starting PhysicalParticleContainer::Evolve for " << species_name << std::endl;
+    amrex::Print() << "NumRealComps() = " << NumRealComps() << std::endl;
+    amrex::Print() << "NumRuntimeRealComps() = " << NumRuntimeRealComps() << std::endl;
+    amrex::Print() << "NStructReal = " << NStructReal << std::endl;
+    amrex::Print() << "NArrayReal = " << NArrayReal << std::endl;
+    for (auto const & comp : GetRealSoANames()) {
+        amrex::Print() << "GetRealSoANames(): " << comp << std::endl;
+    }
+
     WARPX_PROFILE("PhysicalParticleContainer::Evolve()");
     WARPX_PROFILE_VAR_NS("PhysicalParticleContainer::Evolve::GatherAndPush", blp_fg);