From 184b46f75614308478cba83a6de96c5b2f18f372 Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Thu, 6 Apr 2023 11:28:18 -0400 Subject: Enforce reality directly in RK --- src/navier-stokes.c | 9 +++++++++ 1 file changed, 9 insertions(+) (limited to 'src') diff --git a/src/navier-stokes.c b/src/navier-stokes.c index 07eb343..0fc541e 100644 --- a/src/navier-stokes.c +++ b/src/navier-stokes.c @@ -348,6 +348,13 @@ int ns_step( } } + // enforce u(-k)=u(k)^* + for(kx=-K1;kx<=K1;kx++){ + for(ky=-K2;ky<=K2;ky++){ + u[klookup(kx,ky,2*K1+1,2*K2+1)]=(u[klookup(kx,ky,2*K1+1,2*K2+1)]+conj(u[klookup(-kx,-ky,2*K1+1,2*K2+1)]))/2; + } + } + return(0); } @@ -470,12 +477,14 @@ int ns_T( // inverse fft fftw_execute(ifft.fft_plan); + /* // enforce T(u,-k)=T(u,k)^* for(kx=-K1;kx<=K1;kx++){ for(ky=-K2;ky<=K2;ky++){ ifft.fft[klookup(kx,ky,N1,N2)]=(ifft.fft[klookup(kx,ky,N1,N2)]+conj(ifft.fft[klookup(-kx,-ky,N1,N2)]))/2; } } + */ return(0); } -- cgit v1.2.3-54-g00ecf