diff options
Diffstat (limited to 'Source/Evolve/WarpXEvolveEM.cpp')
| -rw-r--r-- | Source/Evolve/WarpXEvolveEM.cpp | 27 |
1 files changed, 20 insertions, 7 deletions
diff --git a/Source/Evolve/WarpXEvolveEM.cpp b/Source/Evolve/WarpXEvolveEM.cpp index a517a91f5..7a3262703 100644 --- a/Source/Evolve/WarpXEvolveEM.cpp +++ b/Source/Evolve/WarpXEvolveEM.cpp @@ -562,6 +562,18 @@ WarpX::ComputeDt () if (do_electrostatic) { dt[0] = const_dt; } + + for (int lev=0; lev <= max_level; lev++) { + const Real* dx_lev = geom[lev].CellSize(); + Print()<<"Level "<<lev<<": dt = "<<dt[lev] + <<" ; dx = "<<dx_lev[0] +#if (defined WARPX_DIM_XZ) || (defined WARPX_DIM_RZ) + <<" ; dz = "<<dx_lev[1]<<'\n'; +#elif (defined WARPX_DIM_3D) + <<" ; dy = "<<dx_lev[1] + <<" ; dz = "<<dx_lev[2]<<'\n'; +#endif + } } /* \brief computes max_step for wakefield simulation in boosted frame. @@ -578,13 +590,14 @@ WarpX::computeMaxStepBoostAccelerator(amrex::Geometry a_geom){ WarpX::moving_window_dir == AMREX_SPACEDIM-1, "Can use zmax_plasma_to_compute_max_step only if " + "moving window along z. TODO: all directions."); - - AMREX_ALWAYS_ASSERT_WITH_MESSAGE( - (WarpX::boost_direction[0]-0)*(WarpX::boost_direction[0]-0) + - (WarpX::boost_direction[1]-0)*(WarpX::boost_direction[1]-0) + - (WarpX::boost_direction[2]-1)*(WarpX::boost_direction[2]-1) < 1.e-12, - "Can use zmax_plasma_to_compute_max_step only if " + - "warpx.boost_direction = z. TODO: all directions."); + if (gamma_boost > 1){ + AMREX_ALWAYS_ASSERT_WITH_MESSAGE( + (WarpX::boost_direction[0]-0)*(WarpX::boost_direction[0]-0) + + (WarpX::boost_direction[1]-0)*(WarpX::boost_direction[1]-0) + + (WarpX::boost_direction[2]-1)*(WarpX::boost_direction[2]-1) < 1.e-12, + "Can use zmax_plasma_to_compute_max_step in boosted frame only if " + + "warpx.boost_direction = z. TODO: all directions."); + } // Lower end of the simulation domain. All quantities are given in boosted // frame except zmax_plasma_to_compute_max_step. |