Fix order of arguments in calling get_TrueFIS. While here, use stat['exp_size'] rather than exp_size when calling pvalue as the latter does not exist.
Use the lowered delta when computing the epsilon, and fix calls to bc by including -l (although I'm not sure what it does)
Use right order of the arguments when calling get_trueFIs in main(). Also only do the optimization part if we actually have candidates.