Ian Jauslin
summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorIan Jauslin <ian@jauslin.org>2022-05-27 00:38:05 -0400
committerIan Jauslin <ian@jauslin.org>2022-05-27 00:38:05 -0400
commitf9d9d481545d4769c3cc73d3c7b516725213f85a (patch)
treea3242dee402e3bea971c9c83b8c7d030ebcf8dc2 /src/navier-stokes.c
parentc62ce075798fd8db209ca3f6a6749e8bd740dc67 (diff)
Remove avg from energy
Diffstat (limited to 'src/navier-stokes.c')
-rw-r--r--src/navier-stokes.c11
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);
}
}