diff options
author | Ian Jauslin <ian.jauslin@rutgers.edu> | 2023-04-21 17:36:17 -0400 |
---|---|---|
committer | Ian Jauslin <ian.jauslin@rutgers.edu> | 2023-04-21 17:36:17 -0400 |
commit | 731244a83d732be84d96763e4747bfe0a8f8a5be (patch) | |
tree | d9c10b9d7b63865275c8f522b4200f5789d266cf /src/navier-stokes.h | |
parent | a53aec9b225ce5f56b18d51552ee5d1df4403008 (diff) |
Only use T term in alpha for ns_step computation
Diffstat (limited to 'src/navier-stokes.h')
-rw-r--r-- | src/navier-stokes.h | 2 |
1 files changed, 2 insertions, 0 deletions
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 |