aboutsummaryrefslogtreecommitdiff
path: root/Source/Diagnostics/BoostedFrameDiagnostic.cpp
diff options
context:
space:
mode:
authorGravatar MaxThevenet <mthevenet@lbl.gov> 2019-05-12 06:54:06 -0700
committerGravatar MaxThevenet <mthevenet@lbl.gov> 2019-05-12 06:54:06 -0700
commit71cdf5ca0bfe5a239382cdfb56df1d2d9e68c65f (patch)
tree68d1c8d8572e826616077011fb9f5f1fad4d4661 /Source/Diagnostics/BoostedFrameDiagnostic.cpp
parent5ebe5c35e5b6cac0428173291ece3e431954037f (diff)
downloadWarpX-71cdf5ca0bfe5a239382cdfb56df1d2d9e68c65f.tar.gz
WarpX-71cdf5ca0bfe5a239382cdfb56df1d2d9e68c65f.tar.zst
WarpX-71cdf5ca0bfe5a239382cdfb56df1d2d9e68c65f.zip
debugging
Diffstat (limited to 'Source/Diagnostics/BoostedFrameDiagnostic.cpp')
-rw-r--r--Source/Diagnostics/BoostedFrameDiagnostic.cpp17
1 files changed, 12 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";
}