I have a newton raphson loop where I have the following:

lambda -= delta;

typically, I have to do an operation such as:

lambda = 52347.592

delta = 0.001542

and then I end up with terribad floating point rounding errors since floating point operations can only support approximately 7 significant digits.

I’d very much prefer to stay away from doubles if possible due to its terrible floating point performance.

edit: after some more investigation - seems like I need to change a lot of floats to doubles - give me a 50% penalty in running time… sigh.

Apparently doubles run at 1/10th the speed of single precision floats