powiaina_num.js 0.2.0-alpha.2.2 → 0.2.0-alpha.2.4

package/dist/PowiainaNum.cjs.js

@@ -38,10 +38,12 @@ function _typeof(o) {
   }, _typeof(o);
 }
 
+var powiainaNumError = "[PowiainaNum 0.2 error]";
 var MSI = 9007199254740991;
 var MSI_LOG10 = 15.954589770191003;
 var MSI_REC = 1.1102230246251568e-16;
 var LONG_STRING_MIN_LENGTH = 17;
+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;
   var e = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : 1;
@@ -69,6 +71,125 @@ function compareTuples() {
 function log10LongString(str) {
   return Math.log10(Number(str.substring(0, LONG_STRING_MIN_LENGTH))) + (str.length - LONG_STRING_MIN_LENGTH);
 }
+// Code from break_eternity.js
+function f_gamma(n) {
+  if (!isFinite(n)) {
+    return n;
+  }
+  if (n < -50) {
+    if (n === Math.trunc(n)) {
+      return Number.NEGATIVE_INFINITY;
+    }
+    return 0;
+  }
+  var scal1 = 1;
+  while (n < 10) {
+    scal1 = scal1 * n;
+    ++n;
+  }
+  n -= 1;
+  var l = 0.9189385332046727; //0.5*Math.log(2*Math.PI)
+  l = l + (n + 0.5) * Math.log(n);
+  l = l - n;
+  var n2 = n * n;
+  var np = n;
+  l = l + 1 / (12 * np);
+  np = np * n2;
+  l = l - 1 / (360 * np);
+  np = np * n2;
+  l = l + 1 / (1260 * np);
+  np = np * n2;
+  l = l - 1 / (1680 * np);
+  np = np * n2;
+  l = l + 1 / (1188 * np);
+  np = np * n2;
+  l = l - 691 / (360360 * np);
+  np = np * n2;
+  l = l + 7 / (1092 * np);
+  np = np * n2;
+  l = l - 3617 / (122400 * np);
+  return Math.exp(l) / scal1;
+}
+var OMEGA = 0.56714329040978387299997; // W(1, 0)
+//from https://math.stackexchange.com/a/465183
+// The evaluation can become inaccurate very close to the branch point
+// Evaluates W(x, 0) if principal is true, W(x, -1) if principal is false
+function f_lambertw(z) {
+  var tol = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 1e-10;
+  var principal = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : true;
+  var w;
+  var wn;
+  if (!Number.isFinite(z)) {
+    return z;
+  }
+  if (principal) {
+    if (z === 0) {
+      return z;
+    }
+    if (z === 1) {
+      return OMEGA;
+    }
+    if (z < 10) {
+      w = 0;
+    } else {
+      w = Math.log(z) - Math.log(Math.log(z));
+    }
+  } else {
+    if (z === 0) return -Infinity;
+    if (z <= -0.1) {
+      w = -2;
+    } else {
+      w = Math.log(-z) - Math.log(-Math.log(-z));
+    }
+  }
+  for (var i = 0; i < 100; ++i) {
+    wn = (z * Math.exp(-w) + w * w) / (w + 1);
+    if (Math.abs(wn - w) < tol * Math.abs(wn)) {
+      return wn;
+    } else {
+      w = wn;
+    }
+  }
+  throw Error("Iteration failed to converge: ".concat(z.toString())); //return Number.NaN;
+}
+//from https://github.com/scipy/scipy/blob/8dba340293fe20e62e173bdf2c10ae208286692f/scipy/special/lambertw.pxd
+// The evaluation can become inaccurate very close to the branch point
+// at ``-1/e``. In some corner cases, `lambertw` might currently
+// fail to converge, or can end up on the wrong branch.
+// Evaluates W(x, 0) if principal is true, W(x, -1) if principal is false
+function d_lambertw(z) {
+  var tol = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 1e-10;
+  var principal = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : true;
+  var w;
+  var ew, wewz, wn;
+  if (z.isInfinite()) return z;
+  if (principal) {
+    if (z.eq(PowiainaNum.ZERO)) {
+      return PowiainaNum.ZERO.clone();
+    }
+    if (z.eq(PowiainaNum.ONE)) {
+      //Split out this case because the asymptotic series blows up
+      return PowiainaNum.fromNumber(OMEGA);
+    } //Get an initial guess for Halley's method
+    w = z.log();
+  } else {
+    if (z.eq(PowiainaNum.ZERO)) {
+      return PowiainaNum.NEGATIVE_INFINITY.clone();
+    } //Get an initial guess for Halley's method
+    w = z.neg().log();
+  } //Halley's method; see 5.9 in [1]
+  for (var i = 0; i < 100; ++i) {
+    ew = w.neg().exp();
+    wewz = w.sub(z.mul(ew));
+    wn = w.sub(wewz.div(w.add(1).sub(w.add(2).mul(wewz).div(w.mul(2).add(2)))));
+    if (wn.sub(w).abs().lt(wn.abs().mul(tol))) {
+      return wn;
+    } else {
+      w = wn;
+    }
+  }
+  throw Error("Iteration failed to converge: ".concat(z.toString())); //return Decimal.dNaN;
+}
 var PowiainaNum = /*#__PURE__*/function () {
   function PowiainaNum(arg1) {
     _classCallCheck(this, PowiainaNum);
@@ -99,6 +220,7 @@ var PowiainaNum = /*#__PURE__*/function () {
   return _createClass(PowiainaNum, [{
     key: "add",
     value: function add(other) {
+      var _a, _b, _c, _d;
       var x = this.clone();
       var y = new PowiainaNum(other);
       // inf + -inf = nan
@@ -132,11 +254,11 @@ var PowiainaNum = /*#__PURE__*/function () {
         a = y;
       }
       var mult = 1;
-      if (!a.small && !b.small && !a.getOperator(1) && !b.getOperator(1) && a.sign == b.sign) {
-        return new PowiainaNum((a.getOperator(0) + b.getOperator(0)) * a.sign);
+      if (!a.small && !b.small && !((_a = a.array[1]) === null || _a === void 0 ? void 0 : _a.repeat) && !((_b = b.array[1]) === null || _b === void 0 ? void 0 : _b.repeat) && a.sign == b.sign) {
+        return new PowiainaNum((a.array[0].repeat + b.array[0].repeat) * a.sign);
       }
-      var alog10 = (a.small ? -1 : 1) * (a.getOperator(1) ? a.getOperator(0) : Math.log10(a.getOperator(0)));
-      var blog10 = (b.small ? -1 : 1) * (b.getOperator(1) ? b.getOperator(0) : Math.log10(b.getOperator(0)));
+      var alog10 = (a.small ? -1 : 1) * (((_c = a.array[1]) === null || _c === void 0 ? void 0 : _c.repeat) ? a.array[0].repeat : Math.log10(a.array[0].repeat));
+      var blog10 = (b.small ? -1 : 1) * (((_d = b.array[1]) === null || _d === void 0 ? void 0 : _d.repeat) ? b.array[0].repeat : Math.log10(b.array[0].repeat));
       if (alog10 - blog10 > MSI_LOG10) return a;
       var offset = -Math.floor(alog10); //a number can make a+off in [0,1)
       var r,
@@ -152,10 +274,9 @@ var PowiainaNum = /*#__PURE__*/function () {
       r = new PowiainaNum();
       r.sign = 1;
       if (l > MSI_LOG10 || l < -MSI_LOG10) {
-        r.setOperator(1, 1);
-        r.setOperator(Math.log10(Math.abs(l)), 1);
+        r.array = [newOperator(l, 0), newOperator(1, 1)];
       } else {
-        r.setOperator(Math.pow(10, Math.abs(l)), 0);
+        r.array = [newOperator(Math.pow(10, l), 0)];
       }
       r.small = l < 0 ? true : false;
       r.sign *= mult;
@@ -286,11 +407,17 @@ var PowiainaNum = /*#__PURE__*/function () {
     }
   }, {
     key: "log",
-    value: function log(base) {
+    value: function log() {
+      var base = arguments.length > 0 && arguments[0] !== undefined ? arguments[0] : Math.E;
       // log_a b = log_x b / log_x a;
       var other = new PowiainaNum(base);
       return this.log10().div(other.log10());
     }
+  }, {
+    key: "exp",
+    value: function exp() {
+      return this.pow_base(Math.E);
+    }
   }, {
     key: "mod",
     value: function mod(x) {
@@ -298,6 +425,116 @@ var PowiainaNum = /*#__PURE__*/function () {
       var division = this.div(other);
       return division.sub(division.floor()).mul(other);
     }
+    /**
+     * For positive integers, X factorial (written as X!) equals X * (X - 1) * (X - 2) *... * 3 * 2 * 1. 0! equals 1.
+     * This can be extended to real numbers (except for negative integers) via the gamma function, which is what this function does.
+     */
+    //[Code from break_eternity.js]
+  }, {
+    key: "factorial",
+    value: function factorial() {
+      if (this.abs().lt(MSI)) {
+        return this.add(1).gamma();
+      } else if (this.abs().lt(PowiainaNum.E_MSI)) {
+        return PowiainaNum.exp(this.mul(this.log10().sub(1)));
+      } else {
+        return PowiainaNum.exp(this);
+      }
+    }
+  }, {
+    key: "gamma",
+    value:
+    /**
+     * The gamma function extends the idea of factorials to non-whole numbers using some calculus.
+     * Gamma(x) is defined as the integral of t^(x-1) * e^-t dt from t = 0 to t = infinity,
+     * and gamma(x) = (x - 1)! for nonnegative integer x, so the factorial for non-whole numbers is defined using the gamma function.
+     */
+    //[Code from break_eternity.js]
+    //from HyperCalc source code  
+    function gamma() {
+      if (this.small) {
+        return this.rec();
+      } else if (this.lte(MSI)) {
+        if (this.lt(24)) {
+          return PowiainaNum.fromNumber(f_gamma(this.sign * this.getOperator(0)));
+        }
+        var t = this.getOperator(0) - 1;
+        var l = 0.9189385332046727; //0.5*Math.log(2*Math.PI)
+        l = l + (t + 0.5) * Math.log(t);
+        l = l - t;
+        var n2 = t * t;
+        var np = t;
+        var lm = 12 * np;
+        var adj = 1 / lm;
+        var l2 = l + adj;
+        if (l2 === l) {
+          return PowiainaNum.exp(l);
+        }
+        l = l2;
+        np = np * n2;
+        lm = 360 * np;
+        adj = 1 / lm;
+        l2 = l - adj;
+        if (l2 === l) {
+          return PowiainaNum.exp(l);
+        }
+        l = l2;
+        np = np * n2;
+        lm = 1260 * np;
+        var lt = 1 / lm;
+        l = l + lt;
+        np = np * n2;
+        lm = 1680 * np;
+        lt = 1 / lm;
+        l = l - lt;
+        return PowiainaNum.exp(l);
+      } else if (this.gt(MSI)) {
+        return PowiainaNum.exp(this.mul(this.log().sub(1)));
+      } else {
+        return PowiainaNum.exp(this);
+      }
+    }
+  }, {
+    key: "lambertw",
+    value:
+    /**
+       * The Lambert W function, also called the omega function or product logarithm, is the solution W(x) === x*e^x.
+       * https://en.wikipedia.org/wiki/Lambert_W_function
+       *
+       * This is a multi-valued function in the complex plane, but only two branches matter for real numbers: the "principal branch" W0, and the "non-principal branch" W_-1.
+       * W_0 works for any number >= -1/e, but W_-1 only works for nonpositive numbers >= -1/e.
+       * The "principal" parameter, which is true by default, decides which branch we're looking for: W_0 is used if principal is true, W_-1 is used if principal is false.
+       */
+    //Code from break_eternity.js
+    //Some special values, for testing: https://en.wikipedia.org/wiki/Lambert_W_function#Special_values
+    function lambertw() {
+      var principal = arguments.length > 0 && arguments[0] !== undefined ? arguments[0] : true;
+      if (this.lt(-0.3678794411710499)) {
+        return PowiainaNum.NaN.clone(); //complex
+      } else if (principal) {
+        if (this.abs().lt("1e-300")) return new PowiainaNum(this);else if (this.small) {
+          return PowiainaNum.fromNumber(f_lambertw(this.toNumber()));
+        } else if (this.layer === 0) {
+          return PowiainaNum.fromNumber(f_lambertw(this.sign * this.getOperator(0)));
+        } else if (this.lt("eee15")) {
+          return d_lambertw(this);
+        } else {
+          // Numbers this large would sometimes fail to converge using d_lambertw, and at this size this.ln() is close enough
+          return this.log();
+        }
+      } else {
+        if (this.sign === 1) {
+          return PowiainaNum.NaN.clone(); //complex
+        }
+        if (this.layer === 0) {
+          return PowiainaNum.fromNumber(f_lambertw(this.sign * this.getOperator(0), 1e-10, false));
+        } else if (this.layer == 1) {
+          return d_lambertw(this, 1e-10, false);
+        } else {
+          return this.neg().rec().lambertw().neg();
+        }
+      }
+    }
   }, {
     key: "max",
     value: function max(x) {
@@ -494,13 +731,14 @@ var PowiainaNum = /*#__PURE__*/function () {
       if (this.abs().gte(Math.pow(2, 52))) return true;
       return false;
     }
+  }, {
+    key: "normalize",
+    value:
     /**
      * Normalize functions will make this number convert into standard format.(it also change `this`, like [].sort)
      * @returns normalized number
      */
-  }, {
-    key: "normalize",
-    value: function normalize() {
+    function normalize() {
       //TODO: normalize
       var renormalize = true;
       var x = this;
@@ -731,7 +969,109 @@ var PowiainaNum = /*#__PURE__*/function () {
       }
       return res;
     }
+  }, {
+    key: "arr01",
+    get:
+    /**
+     * A property arary value for version 0.1.x PowiainaNum.
+     */
+    function get() {
+      var res = [0];
+      for (var i = 0; i < this.array.length; i++) {
+        if (i == 0) res[0] = this.array[i].repeat;else {
+          // @ts-ignore
+          res[i] = [0, 0, 0, 0];
+          // @ts-ignore
+          res[i][0] = this.array[i].arrow == Infinity ? "x" : this.array[i].arrow;
+          // @ts-ignore
+          res[i][1] = this.array[i].repeat;
+          // @ts-ignore
+          res[i][2] = this.array[i].expans == Infinity ? "x" : this.array[i].expans;
+          // @ts-ignore
+          res[i][3] = this.array[i].megota;
+        }
+      }
+      return res;
+    }
   }], [{
+    key: "add",
+    value: function add(t, other) {
+      return new PowiainaNum(t).add(other);
+    }
+  }, {
+    key: "sub",
+    value: function sub(t, other) {
+      return new PowiainaNum(t).sub(other);
+    }
+  }, {
+    key: "mul",
+    value: function mul(t, other) {
+      return new PowiainaNum(t).mul(other);
+    }
+  }, {
+    key: "div",
+    value: function div(t, other) {
+      return new PowiainaNum(t).div(other);
+    }
+  }, {
+    key: "pow",
+    value: function pow(t, other) {
+      return new PowiainaNum(t).pow(other);
+    }
+  }, {
+    key: "root",
+    value: function root(t, other) {
+      return new PowiainaNum(t).root(other);
+    }
+  }, {
+    key: "sqrt",
+    value: function sqrt(t) {
+      return new PowiainaNum(t).sqrt();
+    }
+  }, {
+    key: "cbrt",
+    value: function cbrt(t) {
+      return new PowiainaNum(t).cbrt();
+    }
+  }, {
+    key: "log10",
+    value: function log10(t) {
+      return new PowiainaNum(t).log10();
+    }
+  }, {
+    key: "log",
+    value: function log(t) {
+      var base = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : Math.E;
+      return new PowiainaNum(t).log(base);
+    }
+  }, {
+    key: "exp",
+    value: function exp(x) {
+      var y = new PowiainaNum(x);
+      return y.pow_base(Math.E);
+    }
+  }, {
+    key: "factorial",
+    value: function factorial(x) {
+      return new PowiainaNum(x).factorial();
+    }
+  }, {
+    key: "gamma",
+    value: function gamma(x) {
+      return new PowiainaNum(x).gamma();
+    }
+  }, {
+    key: "lambertw",
+    value: function lambertw(x) {
+      var principal = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : true;
+      return new PowiainaNum(x).lambertw(principal);
+    }
+  }, {
+    key: "isNaN",
+    value: function isNaN(x) {
+      return new PowiainaNum(x).isNaN();
+    }
+  }, {
     key: "fromNumber",
     value: function fromNumber(x) {
       var obj = new PowiainaNum(); // NaN
@@ -740,22 +1080,35 @@ var PowiainaNum = /*#__PURE__*/function () {
       var y = Math.abs(x);
       if (y >= MSI_REC && y < 1) {
         obj.small = true;
-        obj.setOperator(1 / y, 0);
+        obj.array = [newOperator(1 / y, 0)];
       } else if (y < MSI_REC) {
         obj.small = true;
-        obj.setOperator(-Math.log10(y), 0);
-        obj.setOperator(1, 1);
+        obj.array = [newOperator(-Math.log10(y), 0), newOperator(1, 1)];
+      } else if (y <= MSI) {
+        obj.array = [newOperator(y, 0)];
       } else {
-        obj.setOperator(y, 0);
+        obj.setOperator(Math.log10(y), 0);
+        obj.array = [newOperator(Math.log10(y), 0), newOperator(1, 1)];
       }
-      obj.normalize();
       return obj;
     }
   }, {
     key: "fromString",
     value: function fromString(input) {
       var _a, _b, _c, _d, _e, _f;
+      if (!isPowiainaNum.test(input)) {
+        throw powiainaNumError + "malformed input: " + input;
+      }
       var x = new PowiainaNum();
+      // Judge the string was a number
+      // @ts-ignore
+      if (!isNaN(Number(input))) {
+        // @ts-ignore
+        if (isFinite(Number(input))) {
+          x.resetFromObject(PowiainaNum.fromNumber(Number(input)));
+          return x;
+        }
+      }
       var negateIt = false;
       var recipIt = false;
       if (input[0] == "-" || input[0] == "+") {
@@ -976,6 +1329,13 @@ PowiainaNum.E_MSI_REC = new PowiainaNum({
   layer: 0,
   sign: 1
 });
+PowiainaNum.TRITRI = new PowiainaNum({
+  small: false,
+  layer: 0,
+  sign: 1,
+  array: [newOperator(3638334640023.7783, 0, 1, 1), newOperator(7625587484984, 1, 1, 1)]
+});
+PowiainaNum.GRAHAMS_NUMBER = new PowiainaNum("(10{!})^63 10^^^(10^)^7625597484984 3638334640023.7783");
 PowiainaNum.POSITIVE_INFINITY = new PowiainaNum(Infinity);
 PowiainaNum.NEGATIVE_INFINITY = new PowiainaNum(-Infinity);
 PowiainaNum.NaN = new PowiainaNum({