diff options
| -rw-r--r-- | src/navier-stokes.c | 9 |
1 files changed, 9 insertions, 0 deletions
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); } |