diff options
Diffstat (limited to 'src/navier-stokes.c')
-rw-r--r-- | src/navier-stokes.c | 12 |
1 files changed, 12 insertions, 0 deletions
diff --git a/src/navier-stokes.c b/src/navier-stokes.c index 7f22c31..18c2fc5 100644 --- a/src/navier-stokes.c +++ b/src/navier-stokes.c @@ -30,6 +30,18 @@ int uk( ns_init_tmps(&u, &tmp1, &tmp2, &tmp3, &fft1, &fft2, &ifft, K1, K2, N1, N2, nthreads); ns_init_u(u, K1, K2); + // print column headers + printf("# 1:i 2:t "); + t=3; + for(kx=-K1;kx<=K1;kx++){ + for (ky=-K2;ky<=K2;ky++){ + printf(" %6d:(%4d,%4d)r ",t,kx,ky); + t++; + printf(" %6d:(%4d,%4d)i ",t,kx,ky); + t++; + } + } + // iterate for(t=0;t<nsteps;t++){ ins_step(u, K1, K2, N1, N2, nu, delta, g, fft1, fft2, ifft, tmp1, tmp2, tmp3); |