From 091c8d6571d359fc88ce3086ebf4196873797e14 Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Wed, 18 Sep 2024 06:03:33 -0700 Subject: [PATCH] Copy rho from guard cells in IGF solver (#5284) * Copy rho from guard cells in IGF solver * Reset checksum --- .../test_3d_open_bc_poisson_solver.json | 24 +++++++++---------- .../fields/IntegratedGreenFunctionSolver.cpp | 2 +- 2 files changed, 13 insertions(+), 13 deletions(-) diff --git a/Regression/Checksum/benchmarks_json/test_3d_open_bc_poisson_solver.json b/Regression/Checksum/benchmarks_json/test_3d_open_bc_poisson_solver.json index e4ff1fc68a8..0453481ec60 100644 --- a/Regression/Checksum/benchmarks_json/test_3d_open_bc_poisson_solver.json +++ b/Regression/Checksum/benchmarks_json/test_3d_open_bc_poisson_solver.json @@ -1,20 +1,20 @@ { "lev=0": { - "Bx": 100915933.44993827, - "By": 157610622.1855512, - "Bz": 9.717358898362187e-14, - "Ex": 4.7250652706211096e+16, - "Ey": 3.0253948990559976e+16, - "Ez": 3276573.9514776524, - "rho": 10994013582437.193 + "Bx": 100915975.15792876, + "By": 157610677.31483692, + "Bz": 2.404060922276648e-13, + "Ex": 4.725066923361703e+16, + "Ey": 3.0253961494347724e+16, + "Ez": 3276584.4383433666, + "rho": 10994013582437.197 }, "electron": { - "particle_momentum_x": 5.701277606050295e-19, - "particle_momentum_y": 3.6504516641520437e-19, + "particle_momentum_x": 5.701279599504008e-19, + "particle_momentum_y": 3.650453172860547e-19, "particle_momentum_z": 1.145432768297242e-10, - "particle_position_x": 17.314086912497864, - "particle_position_y": 0.2583691267187796, + "particle_position_x": 17.31408691249785, + "particle_position_y": 0.2583691267187801, "particle_position_z": 10066.329600000008, "particle_weight": 19969036501.910976 } -} \ No newline at end of file +} diff --git a/Source/ablastr/fields/IntegratedGreenFunctionSolver.cpp b/Source/ablastr/fields/IntegratedGreenFunctionSolver.cpp index 0767ecfb2f3..1aeee9d81d2 100644 --- a/Source/ablastr/fields/IntegratedGreenFunctionSolver.cpp +++ b/Source/ablastr/fields/IntegratedGreenFunctionSolver.cpp @@ -75,7 +75,7 @@ computePhiIGF ( amrex::MultiFab const & rho, SpectralField tmp_G_fft = SpectralField( spectralspace_ba, dm_global_fft, 1, 0 ); // Copy from rho to tmp_rho - tmp_rho.ParallelCopy( rho, 0, 0, 1, amrex::IntVect::TheZeroVector(), amrex::IntVect::TheZeroVector() ); + tmp_rho.ParallelCopy( rho, 0, 0, 1, rho.nGrowVect(), amrex::IntVect::TheZeroVector() ); // Compute the integrated Green function {