diff options
| author | Ian Jauslin <ian.jauslin@rutgers.edu> | 2023-05-08 15:04:00 -0400 |
|---|---|---|
| committer | Ian Jauslin <ian.jauslin@rutgers.edu> | 2023-05-08 15:04:00 -0400 |
| commit | e5a126f45a072cc1e1b533cfacc870863db763a5 (patch) | |
| tree | ccf67d7494c23f182cf9d3d5ca93cfc739e9a5ac | |
| parent | 0b6fe9bd6bef5a3cba2cb61920de3ad3dc6b57c2 (diff) | |
Replace energy with energy*alpha
| -rw-r--r-- | src/navier-stokes.c | 13 |
1 files changed, 6 insertions, 7 deletions
diff --git a/src/navier-stokes.c b/src/navier-stokes.c index afcb1ee..39f0fad 100644 --- a/src/navier-stokes.c +++ b/src/navier-stokes.c @@ -111,8 +111,8 @@ int eea( _Complex double* tmp1; _Complex double* tmp2; _Complex double* tmp3; - double alpha, energy, enstrophy; - double avg_e,avg_a,avg_en; + double alpha, enstrophy; + double avg_a,avg_en,avg_en_x_a; // index uint64_t t; fft_vect fft1; @@ -125,9 +125,9 @@ int eea( // init running average - avg_e=0; avg_a=0; avg_en=0; + avg_en_x_a=0; // special first case when starting_time is not a multiple of print_freq uint64_t first_box = print_freq - (starting_time % print_freq); @@ -140,17 +140,16 @@ int eea( ns_step_rk4(u, K1, K2, N1, N2, nu, delta, L, g, fft1, fft2, ifft, tmp1, tmp2, tmp3, irreversible); } - energy=compute_energy(u, K1, K2); alpha=compute_alpha(u, K1, K2, g, L); enstrophy=compute_enstrophy(u, K1, K2, L); - avg_e=average_step(energy, avg_e, t, starting_time, print_freq, first_box); avg_a=average_step(alpha, avg_a, t, starting_time, print_freq, first_box); avg_en=average_step(enstrophy, avg_en, t, starting_time, print_freq, first_box); + avg_en_x_a=average_step(enstrophy*alpha, avg_en_x_a, t, starting_time, print_freq, first_box); if(t>starting_time && t%print_freq==0){ - fprintf(stderr,"%lu % .8e % .8e % .8e % .8e % .8e % .8e % .8e\n",t,t*delta, avg_a, avg_e, avg_en, alpha, energy, enstrophy); - printf("%8lu % .15e % .15e % .15e % .15e % .15e % .15e % .15e\n",t,t*delta, avg_a, avg_e, avg_en, alpha, energy, enstrophy); + fprintf(stderr,"%lu % .8e % .8e % .8e % .8e % .8e % .8e % .8e\n",t,t*delta, avg_a, avg_en_x_a, avg_en, alpha, alpha*enstrophy, enstrophy); + printf("%8lu % .15e % .15e % .15e % .15e % .15e % .15e % .15e\n",t,t*delta, avg_a, avg_en_x_a, avg_en, alpha, alpha*enstrophy, enstrophy); } // catch abort signal |