diff options
Diffstat (limited to 'Source')
| -rw-r--r-- | Source/Diagnostics/BoostedFrameDiagnostic.cpp | 17 | ||||
| -rw-r--r-- | Source/Evolve/WarpXEvolveEM.cpp | 4 | ||||
| -rw-r--r-- | Source/Particles/MultiParticleContainer.cpp | 1 | ||||
| -rw-r--r-- | Source/Particles/PhysicalParticleContainer.cpp | 2 |
4 files changed, 19 insertions, 5 deletions
diff --git a/Source/Diagnostics/BoostedFrameDiagnostic.cpp b/Source/Diagnostics/BoostedFrameDiagnostic.cpp index ad33b68d8..934610b9e 100644 --- a/Source/Diagnostics/BoostedFrameDiagnostic.cpp +++ b/Source/Diagnostics/BoostedFrameDiagnostic.cpp @@ -498,7 +498,7 @@ void BoostedFrameDiagnostic::Flush(const Geometry& geom) { BL_PROFILE("BoostedFrameDiagnostic::Flush"); - std::cout<<"in Flush"\n; + std::cout<<"in Flush\n"; VisMF::Header::Version current_version = VisMF::GetHeaderVersion(); VisMF::SetHeaderVersion(amrex::VisMF::Header::NoFabHeader_v1); @@ -546,6 +546,7 @@ void BoostedFrameDiagnostic::Flush(const Geometry& geom) int js = mypc.map_species_lab_diags[j]; std::string species_name = species_names[js]; #ifdef WARPX_USE_HDF5 + std::cout<<"particles_buffer_ j "<<j<<'\n'; writeParticleDataHDF5(particles_buffer_[i][j], snapshots_[i].file_name, species_name); @@ -585,17 +586,21 @@ writeLabFrameData(const MultiFab* cell_centered_data, for (int i = 0; i < N_snapshots_; ++i) { std::cout<<"in BoostedFrameDiagnostic::writeLabFrameData 2\n"; + std::cout<<"i "<<i<<std::endl; const Real old_z_boost = snapshots_[i].current_z_boost; + std::cout<<"toto 1\n"; snapshots_[i].updateCurrentZPositions(t_boost, inv_gamma_boost_, inv_beta_boost_); + std::cout<<"toto 2\n"; if ( (snapshots_[i].current_z_boost < zlo_boost) or (snapshots_[i].current_z_boost > zhi_boost) or (snapshots_[i].current_z_lab < snapshots_[i].zmin_lab) or (snapshots_[i].current_z_lab > snapshots_[i].zmax_lab) ) continue; + std::cout<<"toto 3\n"; int i_lab = (snapshots_[i].current_z_lab - snapshots_[i].zmin_lab) / dz_lab_; std::cout<<"in BoostedFrameDiagnostic::writeLabFrameData 3\n"; @@ -656,9 +661,9 @@ writeLabFrameData(const MultiFab* cell_centered_data, if (WarpX::do_boosted_frame_particles) { std::cout<<"in BoostedFrameDiagnostic::writeLabFrameData 6\n"; - mypc.GetLabFrameData(snapshots_[i].file_name, i_lab, boost_direction_, - old_z_boost, snapshots_[i].current_z_boost, - t_boost, snapshots_[i].t_lab, dt, particles_buffer_[i]); + //mypc.GetLabFrameData(snapshots_[i].file_name, i_lab, boost_direction_, + // old_z_boost, snapshots_[i].current_z_boost, + // t_boost, snapshots_[i].t_lab, dt, particles_buffer_[i]); std::cout<<"in BoostedFrameDiagnostic::writeLabFrameData 7\n"; } @@ -688,6 +693,7 @@ writeLabFrameData(const MultiFab* cell_centered_data, const std::string species_name = species_names[mypc.map_species_lab_diags[j]]; #ifdef WARPX_USE_HDF5 + std::cout<<"particles_buffer_ j "<<j<<'\n'; writeParticleDataHDF5(particles_buffer_[i][j], snapshots_[i].file_name, species_name); @@ -697,8 +703,9 @@ writeLabFrameData(const MultiFab* cell_centered_data, part_ss << snapshots_[i].file_name + "/" + species_name + "/"; std::cout<<"in BoostedFrameDiagnostic::writeLabFrameData 11\n"; + std::cout<<"particles_buffer_ j "<<j<<'\n'; - writeParticleData(particles_buffer_[i][j], part_ss.str(), i_lab); +// writeParticleData(particles_buffer_[i][j], part_ss.str(), i_lab); #endif std::cout<<"in BoostedFrameDiagnostic::writeLabFrameData 12\n"; } diff --git a/Source/Evolve/WarpXEvolveEM.cpp b/Source/Evolve/WarpXEvolveEM.cpp index dab58f95b..5f6760f22 100644 --- a/Source/Evolve/WarpXEvolveEM.cpp +++ b/Source/Evolve/WarpXEvolveEM.cpp @@ -173,7 +173,9 @@ WarpX::EvolveEM (int numsteps) if (WarpX::do_boosted_frame_fields) { cell_centered_data = GetCellCenteredData(); } + std::cout<<"before myBFD->writeLabFrameData"<<std::endl; myBFD->writeLabFrameData(cell_centered_data.get(), *mypc, geom[0], cur_time, dt[0]); + std::cout<<"after myBFD->writeLabFrameData"<<std::endl; } if (to_make_plot || do_insitu) @@ -247,7 +249,9 @@ WarpX::EvolveEM (int numsteps) } if (do_boosted_frame_diagnostic) { + std::cout<<"before myBFD->Flush"<<std::endl; myBFD->Flush(geom[0]); + std::cout<<"after myBFD->Flush"<<std::endl; } #ifdef BL_USE_SENSEI_INSITU diff --git a/Source/Particles/MultiParticleContainer.cpp b/Source/Particles/MultiParticleContainer.cpp index 4726cb4cf..841c41835 100644 --- a/Source/Particles/MultiParticleContainer.cpp +++ b/Source/Particles/MultiParticleContainer.cpp @@ -396,6 +396,7 @@ MultiParticleContainer int isp = map_species_lab_diags[i]; std::cout<<"GetLabFrameData 2\n"; WarpXParticleContainer* pc = allcontainers[isp].get(); + std::cout<<"getparticleslice "<< isp<<std::endl; WarpXParticleContainer::DiagnosticParticles diagnostic_particles; pc->GetParticleSlice(direction, z_old, z_new, t_boost, t_lab, dt, diagnostic_particles); std::cout<<"GetLabFrameData 3\n"; diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index d99cc9c66..c4d4abfd6 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -1863,6 +1863,8 @@ void PhysicalParticleContainer::GetParticleSlice(const int direction, const Real // Note the the slice should always move in the negative boost direction. AMREX_ALWAYS_ASSERT(z_new < z_old); + AMREX_ALWAYS_ASSERT(do_boosted_frame_diags == 1); + const int nlevs = std::max(0, finestLevel()+1); // we figure out a box for coarse-grained rejection. If the RealBox corresponding to a |