diff options
Diffstat (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms')
6 files changed, 65 insertions, 65 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmComoving.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmComoving.cpp index eb7e1578e..37bbc2ca3 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmComoving.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmComoving.cpp @@ -75,15 +75,15 @@ PsatdAlgorithmComoving::pushSpectralFields (SpectralFieldData& f) const const amrex::Box& bx = f.fields[mfi].box(); // Extract arrays for the fields to be updated - amrex::Array4<Complex> fields = f.fields[mfi].array(); + const amrex::Array4<Complex> fields = f.fields[mfi].array(); // Extract arrays for the coefficients - amrex::Array4<const amrex::Real> C_arr = C_coef [mfi].array(); - amrex::Array4<const amrex::Real> S_ck_arr = S_ck_coef[mfi].array(); - amrex::Array4<const Complex> X1_arr = X1_coef [mfi].array(); - amrex::Array4<const Complex> X2_arr = X2_coef [mfi].array(); - amrex::Array4<const Complex> X3_arr = X3_coef [mfi].array(); - amrex::Array4<const Complex> X4_arr = X4_coef [mfi].array(); + const amrex::Array4<const amrex::Real> C_arr = C_coef [mfi].array(); + const amrex::Array4<const amrex::Real> S_ck_arr = S_ck_coef[mfi].array(); + const amrex::Array4<const Complex> X1_arr = X1_coef [mfi].array(); + const amrex::Array4<const Complex> X2_arr = X2_coef [mfi].array(); + const amrex::Array4<const Complex> X3_arr = X3_coef [mfi].array(); + const amrex::Array4<const Complex> X4_arr = X4_coef [mfi].array(); // Extract pointers for the k vectors const amrex::Real* modified_kx_arr = modified_kx_vec[mfi].dataPtr(); @@ -178,13 +178,13 @@ void PsatdAlgorithmComoving::InitializeSpectralCoefficients (const SpectralKSpac const amrex::Real* kz = kz_vec[mfi].dataPtr(); // Extract arrays for the coefficients - amrex::Array4<amrex::Real> C = C_coef [mfi].array(); - amrex::Array4<amrex::Real> S_ck = S_ck_coef [mfi].array(); - amrex::Array4<Complex> X1 = X1_coef [mfi].array(); - amrex::Array4<Complex> X2 = X2_coef [mfi].array(); - amrex::Array4<Complex> X3 = X3_coef [mfi].array(); - amrex::Array4<Complex> X4 = X4_coef [mfi].array(); - amrex::Array4<Complex> T2 = Theta2_coef[mfi].array(); + const amrex::Array4<amrex::Real> C = C_coef [mfi].array(); + const amrex::Array4<amrex::Real> S_ck = S_ck_coef [mfi].array(); + const amrex::Array4<Complex> X1 = X1_coef [mfi].array(); + const amrex::Array4<Complex> X2 = X2_coef [mfi].array(); + const amrex::Array4<Complex> X3 = X3_coef [mfi].array(); + const amrex::Array4<Complex> X4 = X4_coef [mfi].array(); + const amrex::Array4<Complex> T2 = Theta2_coef[mfi].array(); // Store comoving velocity const amrex::Real vx = m_v_comoving[0]; @@ -254,7 +254,7 @@ void PsatdAlgorithmComoving::InitializeSpectralCoefficients (const SpectralKSpac if ( (nu != om_mod/om) && (nu != -om_mod/om) && (nu != 0.) ) { - Complex x1 = om2 / (om2_mod - nu * nu * om2) + const Complex x1 = om2 / (om2_mod - nu * nu * om2) * (theta_star - theta * C(i,j,k) + I * nu * om * theta * S_ck(i,j,k)); // X1 multiplies i*(k \times J) in the update equation for B @@ -424,7 +424,7 @@ void PsatdAlgorithmComoving::CurrentCorrection (SpectralFieldData& field_data) const amrex::Box& bx = field_data.fields[mfi].box(); // Extract arrays for the fields to be updated - amrex::Array4<Complex> fields = field_data.fields[mfi].array(); + const amrex::Array4<Complex> fields = field_data.fields[mfi].array(); // Extract pointers for the k vectors const amrex::Real* const modified_kx_arr = modified_kx_vec[mfi].dataPtr(); diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmFirstOrder.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmFirstOrder.cpp index 113f88fca..946d938c7 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmFirstOrder.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmFirstOrder.cpp @@ -69,7 +69,7 @@ PsatdAlgorithmFirstOrder::pushSpectralFields (SpectralFieldData& f) const const amrex::Box& bx = f.fields[mfi].box(); // Extract arrays for the fields to be updated - amrex::Array4<Complex> fields = f.fields[mfi].array(); + const amrex::Array4<Complex> fields = f.fields[mfi].array(); // Extract pointers for the k vectors const amrex::Real* modified_kx_arr = modified_kx_vec[mfi].dataPtr(); diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJConstantInTime.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJConstantInTime.cpp index 6aad1831e..72c2b7a28 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJConstantInTime.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJConstantInTime.cpp @@ -129,14 +129,14 @@ PsatdAlgorithmJConstantInTime::pushSpectralFields (SpectralFieldData& f) const const amrex::Box& bx = f.fields[mfi].box(); // Extract arrays for the fields to be updated - amrex::Array4<Complex> fields = f.fields[mfi].array(); + const amrex::Array4<Complex> fields = f.fields[mfi].array(); // These coefficients are always allocated - amrex::Array4<const amrex::Real> C_arr = C_coef[mfi].array(); - amrex::Array4<const amrex::Real> S_ck_arr = S_ck_coef[mfi].array(); - amrex::Array4<const Complex> X1_arr = X1_coef[mfi].array(); - amrex::Array4<const Complex> X2_arr = X2_coef[mfi].array(); - amrex::Array4<const Complex> X3_arr = X3_coef[mfi].array(); + const amrex::Array4<const amrex::Real> C_arr = C_coef[mfi].array(); + const amrex::Array4<const amrex::Real> S_ck_arr = S_ck_coef[mfi].array(); + const amrex::Array4<const Complex> X1_arr = X1_coef[mfi].array(); + const amrex::Array4<const Complex> X2_arr = X2_coef[mfi].array(); + const amrex::Array4<const Complex> X3_arr = X3_coef[mfi].array(); amrex::Array4<const Complex> X4_arr; amrex::Array4<const Complex> T2_arr; @@ -379,11 +379,11 @@ void PsatdAlgorithmJConstantInTime::InitializeSpectralCoefficients ( const amrex::Real* kz_c = modified_kz_vec_centered[mfi].dataPtr(); // Coefficients always allocated - amrex::Array4<amrex::Real> C = C_coef[mfi].array(); - amrex::Array4<amrex::Real> S_ck = S_ck_coef[mfi].array(); - amrex::Array4<Complex> X1 = X1_coef[mfi].array(); - amrex::Array4<Complex> X2 = X2_coef[mfi].array(); - amrex::Array4<Complex> X3 = X3_coef[mfi].array(); + const amrex::Array4<amrex::Real> C = C_coef[mfi].array(); + const amrex::Array4<amrex::Real> S_ck = S_ck_coef[mfi].array(); + const amrex::Array4<Complex> X1 = X1_coef[mfi].array(); + const amrex::Array4<Complex> X2 = X2_coef[mfi].array(); + const amrex::Array4<Complex> X3 = X3_coef[mfi].array(); amrex::Array4<Complex> X4; amrex::Array4<Complex> T2; @@ -394,11 +394,11 @@ void PsatdAlgorithmJConstantInTime::InitializeSpectralCoefficients ( } // Extract Galilean velocity - amrex::Real vg_x = m_v_galilean[0]; + const amrex::Real vg_x = m_v_galilean[0]; #if defined(WARPX_DIM_3D) - amrex::Real vg_y = m_v_galilean[1]; + const amrex::Real vg_y = m_v_galilean[1]; #endif - amrex::Real vg_z = m_v_galilean[2]; + const amrex::Real vg_z = m_v_galilean[2]; // Loop over indices within one box ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept @@ -547,19 +547,19 @@ void PsatdAlgorithmJConstantInTime::InitializeSpectralCoefficientsAveraging ( const amrex::Real* kz_c = modified_kz_vec_centered[mfi].dataPtr(); // Coefficients allocated only with averaged Galilean PSATD - amrex::Array4<Complex> Psi1 = Psi1_coef[mfi].array(); - amrex::Array4<Complex> Psi2 = Psi2_coef[mfi].array(); - amrex::Array4<Complex> Y1 = Y1_coef[mfi].array(); - amrex::Array4<Complex> Y3 = Y3_coef[mfi].array(); - amrex::Array4<Complex> Y2 = Y2_coef[mfi].array(); - amrex::Array4<Complex> Y4 = Y4_coef[mfi].array(); + const amrex::Array4<Complex> Psi1 = Psi1_coef[mfi].array(); + const amrex::Array4<Complex> Psi2 = Psi2_coef[mfi].array(); + const amrex::Array4<Complex> Y1 = Y1_coef[mfi].array(); + const amrex::Array4<Complex> Y3 = Y3_coef[mfi].array(); + const amrex::Array4<Complex> Y2 = Y2_coef[mfi].array(); + const amrex::Array4<Complex> Y4 = Y4_coef[mfi].array(); // Extract Galilean velocity - amrex::Real vg_x = m_v_galilean[0]; + const amrex::Real vg_x = m_v_galilean[0]; #if defined(WARPX_DIM_3D) - amrex::Real vg_y = m_v_galilean[1]; + const amrex::Real vg_y = m_v_galilean[1]; #endif - amrex::Real vg_z = m_v_galilean[2]; + const amrex::Real vg_z = m_v_galilean[2]; // Loop over indices within one box ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept @@ -728,7 +728,7 @@ void PsatdAlgorithmJConstantInTime::CurrentCorrection (SpectralFieldData& field_ const amrex::Box& bx = field_data.fields[mfi].box(); // Extract arrays for the fields to be updated - amrex::Array4<Complex> fields = field_data.fields[mfi].array(); + const amrex::Array4<Complex> fields = field_data.fields[mfi].array(); // Extract pointers for the k vectors const amrex::Real* const modified_kx_arr = modified_kx_vec[mfi].dataPtr(); @@ -828,7 +828,7 @@ PsatdAlgorithmJConstantInTime::VayDeposition (SpectralFieldData& field_data) const amrex::Box& bx = field_data.fields[mfi].box(); // Extract arrays for the fields to be updated - amrex::Array4<Complex> fields = field_data.fields[mfi].array(); + const amrex::Array4<Complex> fields = field_data.fields[mfi].array(); // Extract pointers for the modified k vectors const amrex::Real* const modified_kx_arr = modified_kx_vec[mfi].dataPtr(); diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJLinearInTime.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJLinearInTime.cpp index c03ebcf69..30243e630 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJLinearInTime.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJLinearInTime.cpp @@ -84,14 +84,14 @@ PsatdAlgorithmJLinearInTime::pushSpectralFields (SpectralFieldData& f) const const amrex::Box& bx = f.fields[mfi].box(); // Extract arrays for the fields to be updated - amrex::Array4<Complex> fields = f.fields[mfi].array(); + const amrex::Array4<Complex> fields = f.fields[mfi].array(); // These coefficients are always allocated - amrex::Array4<const amrex::Real> C_arr = C_coef[mfi].array(); - amrex::Array4<const amrex::Real> S_ck_arr = S_ck_coef[mfi].array(); - amrex::Array4<const amrex::Real> X1_arr = X1_coef[mfi].array(); - amrex::Array4<const amrex::Real> X2_arr = X2_coef[mfi].array(); - amrex::Array4<const amrex::Real> X3_arr = X3_coef[mfi].array(); + const amrex::Array4<const amrex::Real> C_arr = C_coef[mfi].array(); + const amrex::Array4<const amrex::Real> S_ck_arr = S_ck_coef[mfi].array(); + const amrex::Array4<const amrex::Real> X1_arr = X1_coef[mfi].array(); + const amrex::Array4<const amrex::Real> X2_arr = X2_coef[mfi].array(); + const amrex::Array4<const amrex::Real> X3_arr = X3_coef[mfi].array(); amrex::Array4<const amrex::Real> X5_arr; amrex::Array4<const amrex::Real> X6_arr; @@ -280,11 +280,11 @@ void PsatdAlgorithmJLinearInTime::InitializeSpectralCoefficients ( const amrex::Real* kz_s = modified_kz_vec[mfi].dataPtr(); // Coefficients always allocated - amrex::Array4<amrex::Real> C = C_coef[mfi].array(); - amrex::Array4<amrex::Real> S_ck = S_ck_coef[mfi].array(); - amrex::Array4<amrex::Real> X1 = X1_coef[mfi].array(); - amrex::Array4<amrex::Real> X2 = X2_coef[mfi].array(); - amrex::Array4<amrex::Real> X3 = X3_coef[mfi].array(); + const amrex::Array4<amrex::Real> C = C_coef[mfi].array(); + const amrex::Array4<amrex::Real> S_ck = S_ck_coef[mfi].array(); + const amrex::Array4<amrex::Real> X1 = X1_coef[mfi].array(); + const amrex::Array4<amrex::Real> X2 = X2_coef[mfi].array(); + const amrex::Array4<amrex::Real> X3 = X3_coef[mfi].array(); // Loop over indices within one box amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept @@ -372,11 +372,11 @@ void PsatdAlgorithmJLinearInTime::InitializeSpectralCoefficientsAveraging ( #endif const amrex::Real* kz_s = modified_kz_vec[mfi].dataPtr(); - amrex::Array4<amrex::Real const> C = C_coef[mfi].array(); - amrex::Array4<amrex::Real const> S_ck = S_ck_coef[mfi].array(); + const amrex::Array4<amrex::Real const> C = C_coef[mfi].array(); + const amrex::Array4<amrex::Real const> S_ck = S_ck_coef[mfi].array(); - amrex::Array4<amrex::Real> X5 = X5_coef[mfi].array(); - amrex::Array4<amrex::Real> X6 = X6_coef[mfi].array(); + const amrex::Array4<amrex::Real> X5 = X5_coef[mfi].array(); + const amrex::Array4<amrex::Real> X6 = X6_coef[mfi].array(); // Loop over indices within one box amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmPml.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmPml.cpp index 2b48bb99e..1a4039915 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmPml.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmPml.cpp @@ -69,12 +69,12 @@ PsatdAlgorithmPml::pushSpectralFields(SpectralFieldData& f) const { const Box& bx = f.fields[mfi].box(); // Extract arrays for the fields to be updated - Array4<Complex> fields = f.fields[mfi].array(); + const Array4<Complex> fields = f.fields[mfi].array(); // Extract arrays for the coefficients - Array4<const Real> C_arr = C_coef[mfi].array(); - Array4<const Real> S_ck_arr = S_ck_coef[mfi].array(); - Array4<const Real> inv_k2_arr = inv_k2_coef[mfi].array(); + const Array4<const Real> C_arr = C_coef[mfi].array(); + const Array4<const Real> S_ck_arr = S_ck_coef[mfi].array(); + const Array4<const Real> inv_k2_arr = inv_k2_coef[mfi].array(); // Extract pointers for the k vectors const Real* modified_kx_arr = modified_kx_vec[mfi].dataPtr(); @@ -368,9 +368,9 @@ void PsatdAlgorithmPml::InitializeSpectralCoefficients ( const Real* modified_kz = modified_kz_vec[mfi].dataPtr(); // Extract arrays for the coefficients - Array4<Real> C = C_coef[mfi].array(); - Array4<Real> S_ck = S_ck_coef[mfi].array(); - Array4<Real> inv_k2 = inv_k2_coef[mfi].array(); + const Array4<Real> C = C_coef[mfi].array(); + const Array4<Real> S_ck = S_ck_coef[mfi].array(); + const Array4<Real> inv_k2 = inv_k2_coef[mfi].array(); // Loop over indices within one box ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/SpectralBaseAlgorithm.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/SpectralBaseAlgorithm.cpp index 099f6edf3..172715b98 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/SpectralBaseAlgorithm.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/SpectralBaseAlgorithm.cpp @@ -70,7 +70,7 @@ SpectralBaseAlgorithm::ComputeSpectralDivE ( const Box& bx = field_data.fields[mfi].box(); // Extract arrays for the fields to be updated - Array4<Complex> fields = field_data.fields[mfi].array(); + const Array4<Complex> fields = field_data.fields[mfi].array(); // Extract pointers for the k vectors const Real* modified_kx_arr = modified_kx_vec[mfi].dataPtr(); #if defined(WARPX_DIM_3D) |