diff options
Diffstat (limited to 'src/navier-stokes.h')
| -rw-r--r-- | src/navier-stokes.h | 7 |
1 files changed, 7 insertions, 0 deletions
diff --git a/src/navier-stokes.h b/src/navier-stokes.h index 2ff4e3f..d144fce 100644 --- a/src/navier-stokes.h +++ b/src/navier-stokes.h @@ -17,6 +17,10 @@ int uk( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double d // compute enstrophy int enstrophy( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double delta, _Complex double (*g)(int,int), unsigned int print_freq); +// compute solution as a function of time, but do not print anything (useful for debugging) +int quiet( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double delta, _Complex double (*g)(int,int)); + + // initialize vectors for computation int ns_init_tmps( _Complex double **u, _Complex double ** tmp1, _Complex double **tmp2, _Complex double **tmp3, fft_vect* fft1, fft_vect *fft2, fft_vect *ifft, int K1, int K2, int N1, int N2); // release vectors @@ -34,6 +38,9 @@ int ins_rhs( _Complex double* out, _Complex double* u, int K1, int K2, int N1, i // convolution term in right side of equation int ns_T( _Complex double* u, int K1, int K2, int N1, int N2, fft_vect fft1, fft_vect fft2, fft_vect ifft); +// convolution term in right side of equation +int ns_T_nofft( _Complex double* out, _Complex double* u, int K1, int K2, int N1, int N2); + // compute alpha _Complex double compute_alpha( _Complex double* u, int K1, int K2, _Complex double (*g)(int,int)); |