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

Print Lyapunov exponents

Browse Source
This commit is contained in:
Ian Jauslin 2024-02-19 19:09:26 -05:00
parent 15bcb8b467
commit 1c3f4b9f2f
1 changed files with 6 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/lyapunov.c
View File

@ -132,6 +132,12 @@ int lyapunov(
lyapunov_avg[i]=lyapunov_avg[i]*prevtime/time+lyapunov[i]*(time-prevtime)/time; lyapunov_avg[i]=lyapunov_avg[i]*prevtime/time+lyapunov[i]*(time-prevtime)/time;
} }
// print
for(i=0; i<MATSIZE; i++){
printf("% .15e % .15e % .15e\n",time, lyapunov[i], lyapunov_avg[i]);
}
printf("\n");
// set Du_prod to Q: // set Du_prod to Q:
LAPACKE_zungrq(LAPACK_COL_MAJOR, MATSIZE, MATSIZE, MATSIZE, Du_prod, MATSIZE, tmp2); LAPACKE_zungrq(LAPACK_COL_MAJOR, MATSIZE, MATSIZE, MATSIZE, Du_prod, MATSIZE, tmp2);