powiaina_num.js 0.2.0-alpha.2.7 → 0.2.0-alpha.2.8

package/README.md

@@ -16,7 +16,8 @@ The library exports a class,
 Create a PowiainaNum.js object like this:
 
 ```javascript
-import PowiainaNum from "powiaina_num.js";
+import PowiainaNum from "powiaina_num.js"; // static import 
+const { default: PowiainaNum } = await import("powiaina_num.js") // dynamic import
 
 let a = new PowiainaNum(); // create PN.js number with NaN
 let b = new PowiainaNum(3); // create PN.js number with number 3

package/dist/PowiainaNum.cjs.js

@@ -43,6 +43,7 @@ var MSI = 9007199254740991;
 var MSI_LOG10 = 15.954589770191003;
 var MSI_REC = 1.1102230246251568e-16;
 var LONG_STRING_MIN_LENGTH = 17;
+var EXP_E_REC = 1.444667861009766;
 var isPowiainaNum = /^[-\+]*(Infinity|NaN|(10(\^+|\{([1-9]\d*|!)(,([1-9]\d*|!))?(,[1-9]\d*)?\})|\(10(\^+|\{([1-9]\d*|!)(,([1-9]\d*|!))?(,[1-9]\d*)?\})\)\^[1-9]\d* )*((\d+(\.\d*)?|\d*\.\d+)?([Ee][-\+]*))*(0|\d+(\.\d*)?|\d*\.\d+))$/;
 function newOperator(r) {
   var a = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0;
@@ -162,7 +163,7 @@ function d_lambertw(z) {
   var principal = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : true;
   var w;
   var ew, wewz, wn;
-  if (z.isInfinite()) return z;
+  if (z.isInfiNaN()) return z;
   if (principal) {
     if (z.eq(PowiainaNum.ZERO)) {
       return PowiainaNum.ZERO.clone();
@@ -294,7 +295,7 @@ var PowiainaNum = /*#__PURE__*/function () {
       var y = new PowiainaNum(other);
       // inf * -inf = -inf
       if (x.eq(PowiainaNum.POSITIVE_INFINITY) && y.eq(PowiainaNum.NEGATIVE_INFINITY) || y.eq(PowiainaNum.POSITIVE_INFINITY) && x.eq(PowiainaNum.NEGATIVE_INFINITY)) return PowiainaNum.NEGATIVE_INFINITY.clone();
-      if (x.isInfinite() && y.eq(PowiainaNum.ZERO) && y.isInfinite() && x.eq(PowiainaNum.ZERO)) return PowiainaNum.NaN.clone();
+      if (x.isInfiNaN() && y.eq(PowiainaNum.ZERO) && y.isInfiNaN() && x.eq(PowiainaNum.ZERO)) return PowiainaNum.NaN.clone();
       if (x.eq(PowiainaNum.NEGATIVE_INFINITY) && y.eq(PowiainaNum.NEGATIVE_INFINITY)) return PowiainaNum.POSITIVE_INFINITY.clone();
       // inf & nan check
       if (!x.isFinite()) return x.clone();
@@ -393,6 +394,127 @@ var PowiainaNum = /*#__PURE__*/function () {
     value: function cbrt() {
       return this.root(3);
     }
+  }, {
+    key: "tetrate",
+    value:
+    // Code from ExpantaNum.js
+    function tetrate(other2) {
+      var payload = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 1;
+      var t = this.clone();
+      var other = new PowiainaNum(other2);
+      var payl = new PowiainaNum(payload);
+      if (t.isNaN() || other.isNaN() || payl.isNaN()) return PowiainaNum.NaN.clone();
+      if (payl.neq(PowiainaNum.ONE)) other = other.add(payl.slog(t));
+      var negln;
+      if (other.isInfi() && other.sign > 0) {
+        if (t.gte(EXP_E_REC)) return PowiainaNum.POSITIVE_INFINITY.clone();
+        negln = this.log().neg();
+        return negln.lambertw().div(negln);
+      }
+      if (other.lte(-2)) return PowiainaNum.NaN.clone();
+      if (t.isZero()) {
+        if (other.isZero()) return PowiainaNum.NaN.clone();
+        if (other.gte(MSI / 2) || other.toNumber() % 2 == 0) return PowiainaNum.ZERO.clone();
+        return PowiainaNum.ONE.clone();
+      }
+      if (t.eq(PowiainaNum.ONE)) {
+        if (other.eq(PowiainaNum.ONE.neg())) return PowiainaNum.NaN.clone();
+        return PowiainaNum.ONE.clone();
+      }
+      if (other.eq(PowiainaNum.ONE.neg())) return PowiainaNum.ZERO.clone();
+      if (other.eq(PowiainaNum.ZERO)) return PowiainaNum.ONE.clone();
+      if (other.eq(PowiainaNum.ONE)) return t;
+      if (other.eq(2)) return t.pow(t);
+      if (t.eq(2)) {
+        if (other.eq(3)) return PowiainaNum.fromNumber(16);
+        if (other.eq(4)) return PowiainaNum.fromNumber(65536);
+      }
+      var m = t.max(other);
+      if (m.gt(PowiainaNum.PENTATED_MSI)) return m;
+      if (m.gt(PowiainaNum.TETRATED_MSI) || other.gt(MSI)) {
+        if (this.lt(EXP_E_REC)) {
+          negln = t.ln().neg();
+          return negln.lambertw().div(negln);
+        }
+        var j = t.slog(10).add(other);
+        j.setOperator(j.getOperator(2) + 1, 2);
+        j.normalize();
+        return j;
+      }
+      var y = other.toNumber();
+      var f = Math.floor(y);
+      var r = t.pow(y - f);
+      var l = PowiainaNum.NaN;
+      var i = 0;
+      for (var w = PowiainaNum.E_MSI.clone(); f !== 0 && r.lt(w) && i < 100; ++i) {
+        if (f > 0) {
+          r = t.pow(r);
+          if (l.eq(r)) {
+            f = 0;
+            break;
+          }
+          l = r;
+          --f;
+        } else {
+          r = r.log(t);
+          if (l.eq(r)) {
+            f = 0;
+            break;
+          }
+          l = r;
+          ++f;
+        }
+      }
+      if (i == 100 || this.lt(EXP_E_REC)) f = 0;
+      r.setOperator(r.getOperator(1) + f, 1);
+      r.normalize();
+      return r;
+    }
+    // Code from ExpantaNum.js
+  }, {
+    key: "slog",
+    value: function slog() {
+      var base = arguments.length > 0 && arguments[0] !== undefined ? arguments[0] : 10;
+      var x = this.clone();
+      var b = new PowiainaNum(base);
+      if (x.isInfiNaN()) return x;
+      if (b.isNaN()) return b;
+      if (b.isInfi()) return PowiainaNum.ZERO.clone();
+      if (x.isZero()) return PowiainaNum.ONE.clone();
+      if (x.eq(PowiainaNum.ONE)) return PowiainaNum.ZERO.clone();
+      if (x.eq(b)) return PowiainaNum.ONE.clone();
+      if (b.lt(EXP_E_REC)) {
+        var a = b.tetrate(Infinity);
+        if (x.eq(a)) return PowiainaNum.POSITIVE_INFINITY.clone();
+        if (x.gt(a)) return PowiainaNum.NaN.clone();
+      }
+      if (x.max(b).gt(PowiainaNum.PENTATED_MSI)) {
+        if (x.gt(b)) return x;
+        return PowiainaNum.ZERO.clone();
+      }
+      if (x.max(b).gt(PowiainaNum.TETRATED_MSI)) {
+        if (x.gt(b)) {
+          x.setOperator(x.getOperator(2) - 1, 2);
+          x.normalize();
+          return x.sub(x.getOperator(1));
+        }
+      }
+      if (x.lt(PowiainaNum.ZERO.clone())) return b.pow(x).sub(2);
+      var r = 0;
+      var t = x.getOperator(1) - b.getOperator(1);
+      if (t > 3) {
+        var l = t - 3;
+        r += l;
+        x.setOperator(x.getOperator(1) - l, 1);
+      }
+      for (var i = 0; i < 100; ++i) {
+        if (x.lte(PowiainaNum.ONE)) return new PowiainaNum(r + x.toNumber());
+        ++r;
+        x = PowiainaNum.log(x, base);
+      }
+      console.warn("Failed to converage");
+      return PowiainaNum.NaN.clone();
+    }
   }, {
     key: "abs",
     value: function abs() {
@@ -424,6 +546,11 @@ var PowiainaNum = /*#__PURE__*/function () {
       var other = new PowiainaNum(base);
       return this.log10().div(other.log10());
     }
+  }, {
+    key: "ln",
+    value: function ln() {
+      return this.log();
+    }
   }, {
     key: "exp",
     value: function exp() {
@@ -630,6 +757,18 @@ var PowiainaNum = /*#__PURE__*/function () {
       r.sign = this.sign;
       return r;
     }
+    /**
+     * Work like `Math.trunc`,
+     *
+     * if `this > 0`, return `floor(this)`
+     *
+     * if `this < 0`, return `ceil(this)`
+     *
+     * @example
+     * new PowiainaNum(3.3).trunc() == new PowiainaNum(3)
+     * new PowiainaNum(-1.114514).trunc() == new PowiainaNum(-1)
+     * @returns
+     */
   }, {
     key: "trunc",
     value: function trunc() {
@@ -736,9 +875,14 @@ var PowiainaNum = /*#__PURE__*/function () {
       return Boolean(this.small || isFinite(this.getOperator(0))) && !this.isNaN();
     })
   }, {
-    key: "isInfinite",
-    value: function isInfinite() {
-      return Boolean(!this.small && !isFinite(this.getOperator(0))) || this.isNaN();
+    key: "isInfi",
+    value: function isInfi() {
+      return this.rec().isZero();
+    }
+  }, {
+    key: "isInfiNaN",
+    value: function isInfiNaN() {
+      return this.isInfi() || this.isNaN();
     }
   }, {
     key: "isInt",
@@ -770,6 +914,7 @@ var PowiainaNum = /*#__PURE__*/function () {
       var renormalize = true;
       var x = this;
       for (var _i = 0; _i < this.array.length; _i++) {
+        // Check what is infinity
         if (this.array[_i].repeat == Infinity) {
           this.array = [{
             arrow: 0,
@@ -805,36 +950,42 @@ var PowiainaNum = /*#__PURE__*/function () {
         if (!Number.isInteger(e.expans)) e.expans = Math.floor(e.expans);
         if (!Number.isInteger(e.megota)) e.megota = Math.floor(e.megota);
       }
+      if (!x.array.length) {
+        x.small = !x.small;
+        x.array = [newOperator(Infinity)]; // if no array set zero
+      }
       do {
         renormalize = false;
+        // Sort arrays.
         this.array.sort(function (a, b) {
           return compareTuples([a.megota, a.expans, a.arrow], [b.megota, b.expans, b.arrow]);
         });
+        for (i = 1; i < x.array.length - 1; ++i) {
+          if (x.array[i].arrow == x.array[i + 1].arrow && x.array[i].expans == x.array[i + 1].expans && x.array[i].megota == x.array[i + 1].megota) {
+            // same array's merge
+            x.array[i].repeat += x.array[i + 1].repeat;
+            x.array.splice(i + 1, 1);
+            --i;
+            renormalize = true;
+          }
+        }
         for (i = 1; i < x.array.length; ++i) {
-          // check 0 repeat count
+          // If there is a 0 repeat operator, remove it.
           if (x.array[i].arrow !== 0 && (x.array[i].repeat === 0 || x.array[i].repeat === null || x.array[i].repeat === undefined)) {
             x.array.splice(i, 1);
             --i;
             continue;
           }
-          // check arrow 0 and brace count >=2
+          // If there is a operator which arrow 0 and brace count >=2
+          // replace it as arrow replacement operaotr
           if (x.array[i].arrow == 0 && x.array[i].expans >= 2) {
             x.array[i].arrow = Infinity;
             x.array[i].valuereplaced = 0;
             x.array[i].expans = x.array[i].expans - 1;
           }
         }
-        for (var _i2 = 1; _i2 < this.array.length; _i2++) {
-          if (this.array[_i2].repeat == 0) {
-            this.array.splice(_i2, 1);
-            _i2--;
-          }
-        }
-        if (!x.array.length) {
-          x.small = !x.small;
-          x.array = [newOperator(Infinity)]; // if no array set zero
-        }
         if (x.array.length > PowiainaNum.maxOps) x.array.splice(1, x.array.length - PowiainaNum.maxOps); // max operators check
+        // for any 10^a but a >log10(MSI), replace to regular 10^a
         if (this.getOperator(1) >= 1 && this.getOperator(0) < MSI_LOG10) {
           this.setOperator(this.getOperator(1) - 1, 1);
           this.setOperator(Math.pow(10, this.getOperator(0)), 0);
@@ -845,13 +996,15 @@ var PowiainaNum = /*#__PURE__*/function () {
           this.setOperator(Math.log10(this.getOperator(0)), 0);
           renormalize = true;
         }
+        // for a<1, turn into reciprocate
         if (this.array.length == 1 && this.array[0].repeat < 1) {
           this.array[0].repeat = 1 / this.array[0].repeat;
           this.small = !this.small;
           renormalize = true;
         }
         while (x.array.length >= 2 && x.array[0].repeat == 1 && x.array[1].repeat) {
-          // [1, [R=1, A=sth, E=sth, M=sth]]
+          // for any 10{X}10{X} 1, turn into 10{x}10  
+          // [1, [R=sth, A=sth, E=sth, M=sth]]
           if (x.array[1].repeat > 1) {
             x.array[1].repeat--;
           } else {
@@ -861,19 +1014,18 @@ var PowiainaNum = /*#__PURE__*/function () {
           renormalize = true;
         }
         if (x.array.length >= 2 && x.array[0].repeat < MSI && x.array[1].arrow >= 2 && x.array[1].repeat == 1) {
+          // for any 10{A sample=2}1e9, turn into (10{A-1})^1e9-1 10
           // [1e9, [R=1, A=2, sth, sth]]
           x.array.splice(1, 1, newOperator(x.array[0].repeat, x.array[1].arrow - 1, x.array[1].expans, x.array[1].megota));
           x.array[0].repeat = 10;
           renormalize = true;
         }
-        for (i = 1; i < x.array.length - 1; ++i) {
-          if (x.array[i].arrow == x.array[i + 1].arrow && x.array[i].expans == x.array[i + 1].expans && x.array[i].megota == x.array[i + 1].megota) {
-            // same array's merge
-            x.array[i].repeat += x.array[i + 1].repeat;
-            x.array.splice(i + 1, 1);
-            --i;
-            renormalize = true;
-          }
+        // for any (10{A=2})^1e16 10, turn into (10{A+1}) 1e16
+        if (x.array.length >= 2 && x.array[1].repeat > MSI) {
+          x.array[1].arrow++;
+          x.array[0].repeat = x.array[1].repeat;
+          x.array[1].repeat = 1;
+          renormalize = true;
         }
       } while (renormalize);
       return this;
@@ -981,7 +1133,7 @@ var PowiainaNum = /*#__PURE__*/function () {
         if (this.eq(PowiainaNum.ZERO)) return "0";
         return "/".concat(this.rec().toString());
       }
-      if (this.isInfinite()) return "Infinity";
+      if (this.isInfi()) return "Infinity";
       // O^a (10{arrow,expans,megota})^repeation base
       var res = "";
       for (var i = this.array.length - 1; i >= 0; i--) {
@@ -1064,6 +1216,12 @@ var PowiainaNum = /*#__PURE__*/function () {
     value: function cbrt(t) {
       return new PowiainaNum(t).cbrt();
     }
+  }, {
+    key: "tetrate",
+    value: function tetrate(t, other2) {
+      var payload = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : 1;
+      return new PowiainaNum(t).tetrate(other2, payload);
+    }
   }, {
     key: "log10",
     value: function log10(t) {
@@ -1100,26 +1258,37 @@ var PowiainaNum = /*#__PURE__*/function () {
   }, {
     key: "tetrate_10",
     value: function tetrate_10(other2) {
-      var other = new PowiainaNum(other2);
-      var height_int = other.trunc().toNumber();
-      var height_frac = other.sub(height_int).toNumber();
-      if (other.gt(PowiainaNum.PENTATED_MSI)) return other.clone();else if (other.gt(PowiainaNum.MSI)) {
-        other.setOperator(other.getOperator(2) + 1, 2);
-      } else if (other.lt(-2)) return PowiainaNum.NaN.clone();else if (other.lt(-1)) {
+      return PowiainaNum.fromNumber(10).tetrate(other2);
+      /*
+      const other = new PowiainaNum(other2);
+             const height_int = other.trunc().toNumber();
+      const height_frac = other.sub(height_int).toNumber();
+           if (other.gt(PowiainaNum.PENTATED_MSI))
+        return other.clone();
+      else if (other.gt(PowiainaNum.MSI)) {
+        other.setOperator(other.getOperator(2)+1, 2);
+      } else if (other.lt(-2))
+        return PowiainaNum.NaN.clone();
+      else if (other.lt(-1)) {
         return other.add(2).log10();
-      } else if (other.lt(0)) {
+      }
+      else if (other.lt(0)) {
         return other.add(1);
-      } else if (other.lt(1)) {
-        return other.pow10(); // 10^x
-      } else if (height_int == 1) return PowiainaNum.pow(10, PowiainaNum.pow(10, height_frac));else if (height_int == 2) return PowiainaNum.pow(10, PowiainaNum.pow(10, PowiainaNum.pow(10, height_frac)));else {
-        var remain = height_int - 2;
-        var a = PowiainaNum.pow(10, PowiainaNum.pow(10, PowiainaNum.pow(10, height_frac)));
-        a.setOperator(a.getOperator(1) + remain, 1);
+      }else if (other.lt(1)) {
+        return other.pow10() // 10^x
+      }else if (height_int==1)
+        return PowiainaNum.pow(10,PowiainaNum.pow(10,height_frac));
+      else if (height_int==2)
+        return PowiainaNum.pow(10,PowiainaNum.pow(10,PowiainaNum.pow(10,height_frac)));
+      else {
+        const remain = height_int-2;
+        let a = PowiainaNum.pow(10,PowiainaNum.pow(10,PowiainaNum.pow(10,height_frac)));
+        a.setOperator(a.getOperator(1)+remain, 1);
         return a;
       }
       // 1--2, 10-<1e10, 10^10^0->1
       // 2--3, 1e10-<e1e10, 10^10^10^0->1
-      return PowiainaNum.NaN.clone();
+      return PowiainaNum.NaN.clone();*/
     }
   }, {
     key: "isNaN",
@@ -1259,7 +1428,7 @@ var PowiainaNum = /*#__PURE__*/function () {
         a = input.split(/[Ee]/);
         b = [x.array[0].repeat, 0];
         c = 1;
-        for (var _i3 = a.length - 1; _i3 >= 0; --_i3) {
+        for (var _i2 = a.length - 1; _i2 >= 0; --_i2) {
           //The things that are already there
           if (b[0] < MSI_LOG10 && b[1] === 0) {
             b[0] = Math.pow(10, c * b[0]);
@@ -1276,12 +1445,12 @@ var PowiainaNum = /*#__PURE__*/function () {
             b[1]++;
           }
           //Multiplying coefficient
-          var decimalPointPos = a[_i3].indexOf(".");
-          var intPartLen = decimalPointPos == -1 ? a[_i3].length : decimalPointPos;
+          var decimalPointPos = a[_i2].indexOf(".");
+          var intPartLen = decimalPointPos == -1 ? a[_i2].length : decimalPointPos;
           if (b[1] === 0) {
-            if (intPartLen >= LONG_STRING_MIN_LENGTH) b[0] = Math.log10(b[0]) + log10LongString(a[_i3].substring(0, intPartLen)), b[1] = 1;else if (a[_i3]) b[0] *= Number(a[_i3]);
+            if (intPartLen >= LONG_STRING_MIN_LENGTH) b[0] = Math.log10(b[0]) + log10LongString(a[_i2].substring(0, intPartLen)), b[1] = 1;else if (a[_i2]) b[0] *= Number(a[_i2]);
           } else {
-            d = intPartLen >= LONG_STRING_MIN_LENGTH ? log10LongString(a[_i3].substring(0, intPartLen)) : a[_i3] ? Math.log10(Number(a[_i3])) : 0;
+            d = intPartLen >= LONG_STRING_MIN_LENGTH ? log10LongString(a[_i2].substring(0, intPartLen)) : a[_i2] ? Math.log10(Number(a[_i2])) : 0;
             if (b[1] == 1) {
               b[0] += d;
             } else if (b[1] == 2 && b[0] < MSI_LOG10 + Math.log10(d)) {