From 731244a83d732be84d96763e4747bfe0a8f8a5be Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Fri, 21 Apr 2023 17:36:17 -0400 Subject: Only use T term in alpha for ns_step computation --- src/navier-stokes.h | 2 ++ 1 file changed, 2 insertions(+) (limited to 'src/navier-stokes.h') diff --git a/src/navier-stokes.h b/src/navier-stokes.h index da3894e..ff292cc 100644 --- a/src/navier-stokes.h +++ b/src/navier-stokes.h @@ -49,6 +49,8 @@ int ns_T_nofft( _Complex double* out, _Complex double* u, int K1, int K2, int N1 // 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); // compute energy double compute_energy( _Complex double* u, int K1, int K2); // compute enstrophy -- cgit v1.2.3-54-g00ecf