Correct sqrt initial estimate

MikeMcl · MikeMcl · commit 9a33edc8df72 · 2020-09-30T21:02:31.000+01:00
diff --git a/big.js b/big.js
@@ -782,7 +782,7 @@
       if (!(c.length + e & 1)) c += '0';
       s = Math.sqrt(c);
       e = ((e + 1) / 2 | 0) - (e < 0 || e & 1);
-      r = new Big((s == 1 / 0 ? '1e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);
+      r = new Big((s == 1 / 0 ? '5e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);
     } else {
       r = new Big(s);
     }
diff --git a/big.mjs b/big.mjs
@@ -779,7 +779,7 @@ P.sqrt = function () {
     if (!(c.length + e & 1)) c += '0';
     s = Math.sqrt(c);
     e = ((e + 1) / 2 | 0) - (e < 0 || e & 1);
-    r = new Big((s == 1 / 0 ? '1e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);
+    r = new Big((s == 1 / 0 ? '5e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);
   } else {
     r = new Big(s);
   }

Original file line number	Diff line number	Diff line change
`@@ -782,7 +782,7 @@`
`782`	`782`	`if (!(c.length + e & 1)) c += '0';`
`783`	`783`	`s = Math.sqrt(c);`
`784`	`784`	`e = ((e + 1) / 2 \| 0) - (e < 0 \|\| e & 1);`
`785`		`- r = new Big((s == 1 / 0 ? '1e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);`
	`785`	`+ r = new Big((s == 1 / 0 ? '5e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);`
`786`	`786`	`} else {`
`787`	`787`	`r = new Big(s);`
`788`	`788`	`}`
Original file line number	Diff line number	Diff line change
`@@ -779,7 +779,7 @@ P.sqrt = function () {`
`779`	`779`	`if (!(c.length + e & 1)) c += '0';`
`780`	`780`	`s = Math.sqrt(c);`
`781`	`781`	`e = ((e + 1) / 2 \| 0) - (e < 0 \|\| e & 1);`
`782`		`- r = new Big((s == 1 / 0 ? '1e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);`
	`782`	`+ r = new Big((s == 1 / 0 ? '5e' : (s = s.toExponential()).slice(0, s.indexOf('e') + 1)) + e);`
`783`	`783`	`} else {`
`784`	`784`	`r = new Big(s);`
`785`	`785`	`}`