I don’t know what you think but for me it is not possible to get negative values …
The only way to do it is to get a variable overflow. No?
irej[channel][0] / nb_tot * 100.0
while irej is unsigned, and double nb_tot = (double)(naxes[0] * naxes[1] * nb_frames);
That makes me crazy.