@@ -24,7 +24,7 @@ void ThetaImplicitEM::Define ( WarpX* const a_WarpX )
2424 m_E.Define ( m_WarpX, " Efield_fp" );
2525 m_Eold.Define ( m_E );
2626
27- // Define Bold MultiFabs
27+ // Define B_old MultiFabs
2828 using ablastr::fields::Direction;
2929 const int num_levels = 1 ;
3030 for (int lev = 0 ; lev < num_levels; ++lev) {
@@ -33,9 +33,9 @@ void ThetaImplicitEM::Define ( WarpX* const a_WarpX )
3333 const auto & ba_Bz = m_WarpX->m_fields .get (FieldType::Bfield_fp, Direction{2 }, lev)->boxArray ();
3434 const auto & dm = m_WarpX->m_fields .get (FieldType::Bfield_fp, Direction{0 }, lev)->DistributionMap ();
3535 const amrex::IntVect ngb = m_WarpX->m_fields .get (FieldType::Bfield_fp, Direction{0 }, lev)->nGrowVect ();
36- m_WarpX->m_fields .alloc_init (FieldType::Bold , Direction{0 }, lev, ba_Bx, dm, 1 , ngb, 0 .0_rt);
37- m_WarpX->m_fields .alloc_init (FieldType::Bold , Direction{1 }, lev, ba_By, dm, 1 , ngb, 0 .0_rt);
38- m_WarpX->m_fields .alloc_init (FieldType::Bold , Direction{2 }, lev, ba_Bz, dm, 1 , ngb, 0 .0_rt);
36+ m_WarpX->m_fields .alloc_init (FieldType::B_old , Direction{0 }, lev, ba_Bx, dm, 1 , ngb, 0 .0_rt);
37+ m_WarpX->m_fields .alloc_init (FieldType::B_old , Direction{1 }, lev, ba_By, dm, 1 , ngb, 0 .0_rt);
38+ m_WarpX->m_fields .alloc_init (FieldType::B_old , Direction{2 }, lev, ba_Bz, dm, 1 , ngb, 0 .0_rt);
3939 }
4040
4141 // Parse theta-implicit solver specific parameters
@@ -92,10 +92,10 @@ void ThetaImplicitEM::OneStep ( const amrex::Real a_time,
9292 const int num_levels = 1 ;
9393 for (int lev = 0 ; lev < num_levels; ++lev) {
9494 const ablastr::fields::VectorField Bfp = m_WarpX->m_fields .get_alldirs (FieldType::Bfield_fp, lev);
95- ablastr::fields::VectorField Bold = m_WarpX->m_fields .get_alldirs (FieldType::Bold , lev);
95+ ablastr::fields::VectorField B_old = m_WarpX->m_fields .get_alldirs (FieldType::B_old , lev);
9696 for (int n = 0 ; n < 3 ; ++n) {
97- amrex::MultiFab::Copy ( *Bold [n], *Bfp[n], 0 , 0 , Bold [n]->nComp (),
98- Bold [n]->nGrowVect () );
97+ amrex::MultiFab::Copy (*B_old [n], *Bfp[n], 0 , 0 , B_old [n]->nComp (),
98+ B_old [n]->nGrowVect () );
9999 }
100100 }
101101
@@ -147,8 +147,8 @@ void ThetaImplicitEM::UpdateWarpXFields ( const WarpXSolverVec& a_E,
147147 m_WarpX->SetElectricFieldAndApplyBCs ( a_E );
148148
149149 // Update Bfield_fp owned by WarpX
150- ablastr::fields::MultiLevelVectorField const & Bold = m_WarpX->m_fields .get_mr_levels_alldirs (FieldType::Bold , 0 );
151- m_WarpX->UpdateMagneticFieldAndApplyBCs ( Bold , m_theta* a_dt );
150+ ablastr::fields::MultiLevelVectorField const & B_old = m_WarpX->m_fields .get_mr_levels_alldirs (FieldType::B_old , 0 );
151+ m_WarpX->UpdateMagneticFieldAndApplyBCs (B_old , m_theta * a_dt );
152152
153153}
154154
@@ -163,7 +163,7 @@ void ThetaImplicitEM::FinishFieldUpdate ( amrex::Real a_new_time )
163163 const amrex::Real c1 = 1 ._rt - c0;
164164 m_E.linComb ( c0, m_E, c1, m_Eold );
165165 m_WarpX->SetElectricFieldAndApplyBCs ( m_E );
166- ablastr::fields::MultiLevelVectorField const & Bold = m_WarpX->m_fields .get_mr_levels_alldirs (FieldType::Bold , 0 );
167- m_WarpX->FinishMagneticFieldAndApplyBCs ( Bold , m_theta );
166+ ablastr::fields::MultiLevelVectorField const & B_old = m_WarpX->m_fields .get_mr_levels_alldirs (FieldType::B_old , 0 );
167+ m_WarpX->FinishMagneticFieldAndApplyBCs (B_old , m_theta );
168168
169169}
0 commit comments