diff options
| author |  Remi Lehe <remi.lehe@normalesup.org>
| 2019-04-30 09:41:21 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2019-04-30 09:41:21 -0700 |
| commit | bb1289f3879155244372f791d6854107150f4266 (patch) | |
| tree | 5c2697eeb93fa75b6f82b5b3603c393b6c57e39a /Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp | |
| parent | 9a7a87fcd64d7faab8c57cbc769840f8ee74454a (diff) | |
| parent | fcab5379ddadb02f81d14c9583000c0aeb5b7db6 (diff) | |
| download | WarpX-bb1289f3879155244372f791d6854107150f4266.tar.gz WarpX-bb1289f3879155244372f791d6854107150f4266.tar.zst WarpX-bb1289f3879155244372f791d6854107150f4266.zip | |
Merge pull request #104 from ECP-WarpX/fix_spectral
Fix bugs in the spectral solver
Diffstat (limited to 'Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp')
| -rw-r--r-- | Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp b/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp index d91891a30..ddb2020d8 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp @@ -161,12 +161,13 @@ SpectralKSpace::getModifiedKComponent( const DistributionMapping& dm, // Fill the modified k vector for (int i=0; i<k.size(); i++ ){ + modified_k[i] = 0; for (int n=1; n<stencil_coef.size(); n++){ if (nodal){ - modified_k[i] = stencil_coef[n]* \ + modified_k[i] += stencil_coef[n]* \ std::sin( k[i]*n*delta_x )/( n*delta_x ); } else { - modified_k[i] = stencil_coef[n]* \ + modified_k[i] += stencil_coef[n]* \ std::sin( k[i]*(n-0.5)*delta_x )/( (n-0.5)*delta_x ); } } |