diff options
author | Ian Jauslin <ian.jauslin@rutgers.edu> | 2023-04-22 15:04:02 -0400 |
---|---|---|
committer | Ian Jauslin <ian.jauslin@rutgers.edu> | 2023-04-22 15:04:02 -0400 |
commit | cbcb39ff19c441a770671e151218330d1ac7bc91 (patch) | |
tree | 863ef320729432ab6d6a44aa17408dd844dbc5fe /src/navier-stokes.h | |
parent | 5d0a1bcc6f952e31be072b2aa24487486c483509 (diff) |
Remove term that sums up to 0 from alpha
Diffstat (limited to 'src/navier-stokes.h')
-rw-r--r-- | src/navier-stokes.h | 4 |
1 files changed, 1 insertions, 3 deletions
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 |