Further... If you need double precision you might have to do it yourself. In absence of comments in the code I posted, the algorithm I used for single precision is from here: http://www.utdallas.edu/~cantrell/ee6481/lectures/float_comp.pdf Eventually I'll get around to double precision, but not soon. Best wishes Stephen