diff options
Diffstat (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJLinearInTime.cpp')
| -rw-r--r-- | Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmJLinearInTime.cpp | 30 |
1 files changed, 15 insertions, 15 deletions
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 |