From cbcb39ff19c441a770671e151218330d1ac7bc91 Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Sat, 22 Apr 2023 15:04:02 -0400 Subject: Remove term that sums up to 0 from alpha --- src/navier-stokes.h | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) (limited to 'src/navier-stokes.h') diff --git a/src/navier-stokes.h b/src/navier-stokes.h index 6f14bc8..d12f908 100644 --- a/src/navier-stokes.h +++ b/src/navier-stokes.h @@ -48,9 +48,7 @@ int ns_T( _Complex double* u, int K1, int K2, int N1, int N2, fft_vect fft1, fft int ns_T_nofft( _Complex double* out, _Complex double* u, int K1, int K2, int N1, int N2); // compute alpha -double compute_alpha( _Complex double* u, int K1, int K2, int N1, int N2, _Complex double* g, double L, _Complex double* T); -// do not include the term involving T, which, in the limit K->infty, vanishes -double compute_alpha_fast( _Complex double* u, int K1, int K2, _Complex double* g, double L); +double compute_alpha( _Complex double* u, int K1, int K2, _Complex double* g, double L); // compute energy double compute_energy( _Complex double* u, int K1, int K2); // compute enstrophy -- cgit v1.2.3-54-g00ecf