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

Bug fix

Browse Source
This commit is contained in:
Ian Jauslin 2025-01-31 15:04:41 -05:00
parent 0244f03d02
commit 3fa3a86066
1 changed files with 2 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/lyapunov.c
View File

@ -221,8 +221,8 @@ int lyapunov_D_step(
// compute T // compute T
lyapunov_T(N1,N2,fftu1,fftu2,fftu3,fftu4,ifft); lyapunov_T(N1,N2,fftu1,fftu2,fftu3,fftu4,ifft);
// save to vect // save to vect
for(lx=-K1;lx<=K1;lx++){ for(lx=0;lx<=K1;lx++){
for(ly=-K2;ly<=K2;ly++){ for(ly=(lx>0 ? -K2 : 1);ly<=K2;ly++){
tmp4[klookup_sym(lx,ly,K2)]=ifft.fft[klookup(lx,ly,N1,N2)]; tmp4[klookup_sym(lx,ly,K2)]=ifft.fft[klookup(lx,ly,N1,N2)];
} }
} }
@ -664,9 +664,6 @@ int lyapunov_init_tmps(
unsigned int algorithm, unsigned int algorithm,
unsigned int algorithm_lyapunov unsigned int algorithm_lyapunov
){ ){
// velocity field
*u=calloc(USIZE,sizeof(_Complex double));
// allocate tmp vectors for computation // allocate tmp vectors for computation
if(algorithm_lyapunov==ALGORITHM_RK2){ if(algorithm_lyapunov==ALGORITHM_RK2){
*tmp1=calloc(MATSIZE,sizeof(double)); *tmp1=calloc(MATSIZE,sizeof(double));