redshift_list is a double but has integer comparison
jacobseiler opened this issue · 1 comments
jacobseiler commented
In main.c we initialize double *redshift_list = NULL; and then read values in input_redshifts.c.
However then in cifog.c we attempt to make an integer comparison if((redshift_list[2*cycle+1] == 1) || (cycle == 0)) snap++;.
I saved my redshift file using np.savetxt() which saved as 1.000000 instead of 1 resulting in snap not being incremented properly.
jacobseiler commented
Oops it turns out I wasn't reading in the redshift file properly....