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

Bug with computing averages

Browse Source
This commit is contained in:
Ian Jauslin 2023-06-20 11:07:35 -04:00
parent 870196ec07
commit 9c3ad2c400
1 changed files with 6 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/navier-stokes.c
View File

@ -243,19 +243,17 @@ int enstrophy(
fprintf(stderr,"% .8e % .8e % .8e % .8e % .8e % .8e % .8e\n",time, avg_a, avg_en, avg_en_x_a, alpha, enstrophy, alpha*enstrophy); fprintf(stderr,"% .8e % .8e % .8e % .8e % .8e % .8e % .8e\n",time, avg_a, avg_en, avg_en_x_a, alpha, enstrophy, alpha*enstrophy);
printf("% .15e % .15e % .15e % .15e % .15e % .15e % .15e\n",time, avg_a, avg_en, avg_en_x_a, alpha, enstrophy, alpha*enstrophy); printf("% .15e % .15e % .15e % .15e % .15e % .15e % .15e\n",time, avg_a, avg_en, avg_en_x_a, alpha, enstrophy, alpha*enstrophy);
} }
// print to stdout
}
// reset averages // reset averages
avg_a=0; avg_a=0;
avg_en=0; avg_en=0;
avg_en_x_a=0; avg_en_x_a=0;
prevtime=time; prevtime=time;
}
// get ready for next step // get ready for next step
step=next_step; step=next_step;
// catch abort signal // catch abort signal
if (g_abort){ if (g_abort){
// print u to stderr if no savefile // print u to stderr if no savefile