It sounds like my method will not help you in this case. I have used it in solving triangles. I want $b=sqrt($c**2-$a**2). The mathematically equivalent statement $b=sqrt(($c+$a)*($c-$a)) allocates almost twice as many bits for the subtraction. Not knowing this simple trick added weeks of effort to my first real FORTRAN program.