software/Nstrophy
1
Fork 0
Code Issues Pull Requests Releases Activity

Remove avg from energy

Browse Source
This commit is contained in:
Ian Jauslin
2022-05-27 00:38:05 -04:00
parent c62ce07579
commit f9d9d48154
1 changed files with 3 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
src/navier-stokes.c
View File

@@ -105,17 +105,12 @@ int energy(
// iterate // iterate
for(t=0;t<nsteps;t++){ for(t=0;t<nsteps;t++){
ins_step(u, K1, K2, N1, N2, nu, delta, L, g, fft1, fft2, ifft, tmp1, tmp2, tmp3); ins_step(u, K1, K2, N1, N2, nu, delta, L, g, fft1, fft2, ifft, tmp1, tmp2, tmp3);
if(t%print_freq==0){ if(t%print_freq==0){
energy=compute_energy(u, K1, K2); energy=compute_energy(u, K1, K2);
// running average fprintf(stderr,"%d % .8e % .8e\n",t,t*delta, energy);
if(t>0){ printf("%8d % .15e % .15e\n",t,t*delta,energy);
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);
} }
} }