diff options
author | Ian Jauslin <ian@jauslin.org> | 2022-05-27 00:38:05 -0400 |
---|---|---|
committer | Ian Jauslin <ian@jauslin.org> | 2022-05-27 00:38:05 -0400 |
commit | f9d9d481545d4769c3cc73d3c7b516725213f85a (patch) | |
tree | a3242dee402e3bea971c9c83b8c7d030ebcf8dc2 /src | |
parent | c62ce075798fd8db209ca3f6a6749e8bd740dc67 (diff) |
Remove avg from energy
Diffstat (limited to 'src')
-rw-r--r-- | src/navier-stokes.c | 11 |
1 files changed, 3 insertions, 8 deletions
diff --git a/src/navier-stokes.c b/src/navier-stokes.c index 588dcf4..5082989 100644 --- a/src/navier-stokes.c +++ b/src/navier-stokes.c @@ -105,17 +105,12 @@ int energy( // iterate for(t=0;t<nsteps;t++){ ins_step(u, K1, K2, N1, N2, nu, delta, L, g, fft1, fft2, ifft, tmp1, tmp2, tmp3); - + if(t%print_freq==0){ energy=compute_energy(u, K1, K2); - // running average - if(t>0){ - avg=avg-(avg-energy)/t; - } - - fprintf(stderr,"%d % .8e % .8e % .8e\n",t,t*delta, energy, avg); - printf("%8d % .15e % .15e % .15e\n",t,t*delta,energy,avg); + fprintf(stderr,"%d % .8e % .8e\n",t,t*delta, energy); + printf("%8d % .15e % .15e\n",t,t*delta,energy); } } |