for the record all the computation in the C code is done with float32 so it
can explain numerical errors...