From 60fd6699d71212a71b5c2373ac2d13c21296ae1b Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Thu, 18 Apr 2019 16:31:30 -0700 Subject: Added modified k vector --- Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp | 16 ++++++++++++++++ 1 file changed, 16 insertions(+) (limited to 'Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp') diff --git a/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp b/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp index ab684444d..03d62892f 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralKSpace.cpp @@ -60,3 +60,19 @@ AllocateAndFillKvector( ManagedVector& k, const Box& bx, const Real* dx, c // TODO: For real-to-complex, } + +void +ComputeModifiedKVector( ManagedVector& modified_k, + const ManagedVector& k, + const Box& bx, const Real dx, const int norder ) +{ + // Allocate modified_k to the right size + int N = k.size(); + modified_k.resize( N ); + + // For now, this simply copies the infinite order k + for (int i=0; i