[fpc-pascal] Re: fast sqrt routine

2006-03-22 Thread Thomas Schatzl
>Andrew Bennett >Sat, 11 Mar 2006 17:24:51 -0800 > >On Sat, 11 Mar 2006 13:44:52 +0100, Pianoman wrote: > >>Hi, I need to perform very fast squareroot (sqrt) operations on double >> type nukmbers. >>I tried this >>function mysqrt(a:double):double; >> [...] >>yn:=(y*y+a)/(2*y); >>until abs(yn-y) <

[fpc-pascal] Re: fast sqrt routine

2006-03-11 Thread Andrew Bennett
On Sat, 11 Mar 2006 13:44:52 +0100, Pianoman wrote: >Hi, I need to perform very fast squareroot (sqrt) operations on double type >nukmbers. >I tried this >function mysqrt(a:double):double; >var y,yn:double; >begin >yn:=a; >repeat >y:=yn; >yn:=(y*y+a)/(2*y); >until abs(yn-y) < 10e-16; >mysqrt:=yn;