Square root by hand: Difference between revisions
→{{header|Phix}}: stress test?
m (→{{header|Phix}}: extended tests) |
(→{{header|Phix}}: stress test?) |
||
Line 509:
8968725339654633180882964062061525835239505474575028775996172983557522033753185701135437460340849884
71603868999706990048150305440277903164542478230684929369186215805784631115966687130130156185689872372
</pre>
stress test?:
<lang Phix>include mpfr.e
mpfr_set_default_prec(-100) -- 100 d.p precision
mpfr pi = mpfr_init()
mpfr_const_pi(pi)
string ps = mpfr_sprintf("%.100Rf", pi),
rs = spigot_sqrt(ps,102) -- (<=101 is not enough)
mpfr_set_str(pi,rs)
mpfr_mul(pi,pi,pi)
rs = mpfr_sprintf("%.100Rf", pi)
printf(1,"Pi (builtin) vs spigot_sqrt(pi) squared:\n %s\n %s\n",{ps,rs})</lang>
{{out}}
<pre>
Pi (builtin) vs spigot_sqrt(pi) squared:
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170680
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170680
</pre>
| |||