From d16c42d9f5a40b94406a859fa510bba96480d5e8 Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Tue, 11 Apr 2023 18:45:45 -0400 Subject: Only store u[kx,ky] with kx>=0 --- src/navier-stokes.h | 4 ++++ 1 file changed, 4 insertions(+) (limited to 'src/navier-stokes.h') diff --git a/src/navier-stokes.h b/src/navier-stokes.h index 582212a..d100f44 100644 --- a/src/navier-stokes.h +++ b/src/navier-stokes.h @@ -52,6 +52,10 @@ double compute_enstrophy( _Complex double* u, int K1, int K2, double L); // get index for kx,ky in array of size S int klookup( int kx, int ky, int S1, int S2); +// get index for kx,ky in array of size K1,K2 in which only the terms with kx>=0 are stored +int klookup_sym( int kx, int ky, int K2); +// get u_{kx,ky} from a vector u in which only the values for kx>=0 are stored, assuming u_{-k}=u_k^* +_Complex double getval_sym( _Complex double* u, int kx, int ky, int K2); #endif -- cgit v1.2.3-54-g00ecf