From 80787133443c9adebdf6a9c6cdc6538bb2bcd2df Mon Sep 17 00:00:00 2001
From: Remi Lehe <remi.lehe@normalesup.org>
Date: Sat, 20 Apr 2019 22:53:00 -0700
Subject: Added shift function

---
 Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp | 8 ++++++++
 1 file changed, 8 insertions(+)

(limited to 'Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp')

diff --git a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
index 6b39e9fa4..78b7d7c67 100644
--- a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
+++ b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
@@ -25,6 +25,14 @@ SpectralFieldData::SpectralFieldData( const BoxArray& realspace_ba,
     tmpRealField = SpectralField(realspace_ba, dm, 1, 0);
     tmpSpectralField = SpectralField(spectralspace_ba, dm, 1, 0);
 
+    // Allocate the vectors that allow to shift between nodal and cell-centered
+    for (int i_dim=0; i_dim<AMREX_SPACEDIM; i_dim++) {
+        shift_C2N[i_dim] = k_space.AllocateAndFillSpectralShiftFactor(
+                            dm, i_dim, ShiftType::CenteredToNodal );
+        shift_N2C[i_dim] = k_space.AllocateAndFillSpectralShiftFactor(
+                            dm, i_dim, ShiftType::NodalToCentered );        
+    }
+
     // Allocate and initialize the FFT plans
     forward_plan = FFTplans(spectralspace_ba, dm);
     backward_plan = FFTplans(spectralspace_ba, dm);
-- 
cgit v1.2.3