diff --git a/src/io.c b/src/io.c index 3287233..4f11dea 100644 --- a/src/io.c +++ b/src/io.c @@ -71,7 +71,7 @@ int read_vec(_Complex double* out, int K1, int K2, FILE* file){ } // allocate line buffer - line=calloc(sizeof(char), len); + line=calloc(len,sizeof(char)); while(1){ c=fgetc(file); @@ -126,7 +126,7 @@ int read_vec(_Complex double* out, int K1, int K2, FILE* file){ // check that there is room in buffer if(pos==len){ // too short: reallocate - line_realloc=calloc(sizeof(char), 2*len); + line_realloc=calloc(2*len,sizeof(char)); for(pos=0;posK1*(2*parameters->K2+1)+parameters->K2); + _Complex double* u0=calloc(parameters->K1*(2*parameters->K2+1)+parameters->K2,sizeof(_Complex double)); switch(parameters->init){ case INIT_RANDOM: diff --git a/src/navier-stokes.c b/src/navier-stokes.c index 3e87f5e..78c0ded 100644 --- a/src/navier-stokes.c +++ b/src/navier-stokes.c @@ -346,30 +346,30 @@ int ns_init_tmps( unsigned int algorithm ){ // velocity field - *u=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); + *u=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); // allocate tmp vectors for computation if(algorithm==ALGORITHM_RK2){ - *tmp1=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp2=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); + *tmp1=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp2=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); } else if (algorithm==ALGORITHM_RK4){ - *tmp1=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp2=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp3=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); + *tmp1=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp2=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp3=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); } else if (algorithm==ALGORITHM_RKF45 || algorithm==ALGORITHM_RKDP54){ - *tmp1=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp2=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp3=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp4=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp5=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp6=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp7=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); + *tmp1=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp2=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp3=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp4=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp5=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp6=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp7=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); } else if (algorithm==ALGORITHM_RKBS32){ - *tmp1=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp2=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp3=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp4=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); - *tmp5=calloc(sizeof(_Complex double),K1*(2*K2+1)+K2); + *tmp1=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp2=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp3=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp4=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); + *tmp5=calloc(K1*(2*K2+1)+K2,sizeof(_Complex double)); } else { fprintf(stderr,"bug: unknown algorithm: %u, contact ian.jauslin@rutgers,edu\n",algorithm); };