diff options
| author | Ian Jauslin <ian@jauslin.org> | 2022-05-19 17:51:45 +0200 |
|---|---|---|
| committer | Ian Jauslin <ian@jauslin.org> | 2022-05-19 17:51:45 +0200 |
| commit | 146903265abcfb0dce8ee6075b1cd355cab3ab94 (patch) | |
| tree | b91c91c57e484192500c53d26fa4e4d18a7ed880 /src/navier-stokes.c | |
| parent | 8bce8632c520acea86b4a6b2b7cf3eafe7964124 (diff) | |
no driving option
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); |