mathjs 12.4.2 → 13.0.0

package/lib/cjs/function/complex/im.js

@@ -6,10 +6,12 @@ Object.defineProperty(exports, "__esModule", {
 exports.createIm = void 0;
 var _factory = require("../../utils/factory.js");
 var _collection = require("../../utils/collection.js");
-var name = 'im';
-var dependencies = ['typed'];
-var createIm = exports.createIm = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
-  var typed = _ref.typed;
+const name = 'im';
+const dependencies = ['typed'];
+const createIm = exports.createIm = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+  let {
+    typed
+  } = _ref;
   /**
    * Get the imaginary part of a complex number.
    * For a complex number `a + bi`, the function returns `b`.
@@ -38,19 +40,9 @@ var createIm = exports.createIm = /* #__PURE__ */(0, _factory.factory)(name, dep
    * @return {number | BigNumber | Array | Matrix} The imaginary part of x
    */
   return typed(name, {
-    number: function number() {
-      return 0;
-    },
-    'BigNumber | Fraction': function BigNumberFraction(x) {
-      return x.mul(0);
-    },
-    Complex: function Complex(x) {
-      return x.im;
-    },
-    'Array | Matrix': typed.referToSelf(function (self) {
-      return function (x) {
-        return (0, _collection.deepMap)(x, self);
-      };
-    })
+    number: () => 0,
+    'BigNumber | Fraction': x => x.mul(0),
+    Complex: x => x.im,
+    'Array | Matrix': typed.referToSelf(self => x => (0, _collection.deepMap)(x, self))
   });
 });

package/lib/cjs/function/complex/re.js

@@ -6,10 +6,12 @@ Object.defineProperty(exports, "__esModule", {
 exports.createRe = void 0;
 var _factory = require("../../utils/factory.js");
 var _collection = require("../../utils/collection.js");
-var name = 're';
-var dependencies = ['typed'];
-var createRe = exports.createRe = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
-  var typed = _ref.typed;
+const name = 're';
+const dependencies = ['typed'];
+const createRe = exports.createRe = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+  let {
+    typed
+  } = _ref;
   /**
    * Get the real part of a complex number.
    * For a complex number `a + bi`, the function returns `a`.
@@ -38,16 +40,8 @@ var createRe = exports.createRe = /* #__PURE__ */(0, _factory.factory)(name, dep
    * @return {number | BigNumber | Array | Matrix} The real part of x
    */
   return typed(name, {
-    'number | BigNumber | Fraction': function numberBigNumberFraction(x) {
-      return x;
-    },
-    Complex: function Complex(x) {
-      return x.re;
-    },
-    'Array | Matrix': typed.referToSelf(function (self) {
-      return function (x) {
-        return (0, _collection.deepMap)(x, self);
-      };
-    })
+    'number | BigNumber | Fraction': x => x,
+    Complex: x => x.re,
+    'Array | Matrix': typed.referToSelf(self => x => (0, _collection.deepMap)(x, self))
   });
 });

package/lib/cjs/function/geometry/distance.js

@@ -6,17 +6,19 @@ Object.defineProperty(exports, "__esModule", {
 exports.createDistance = void 0;
 var _is = require("../../utils/is.js");
 var _factory = require("../../utils/factory.js");
-var name = 'distance';
-var dependencies = ['typed', 'addScalar', 'subtractScalar', 'divideScalar', 'multiplyScalar', 'deepEqual', 'sqrt', 'abs'];
-var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
-  var typed = _ref.typed,
-    addScalar = _ref.addScalar,
-    subtractScalar = _ref.subtractScalar,
-    multiplyScalar = _ref.multiplyScalar,
-    divideScalar = _ref.divideScalar,
-    deepEqual = _ref.deepEqual,
-    sqrt = _ref.sqrt,
-    abs = _ref.abs;
+const name = 'distance';
+const dependencies = ['typed', 'addScalar', 'subtractScalar', 'divideScalar', 'multiplyScalar', 'deepEqual', 'sqrt', 'abs'];
+const createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+  let {
+    typed,
+    addScalar,
+    subtractScalar,
+    multiplyScalar,
+    divideScalar,
+    deepEqual,
+    sqrt,
+    abs
+  } = _ref;
   /**
     * Calculates:
     *    The eucledian distance between two points in N-dimensional spaces.
@@ -69,7 +71,7 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
     * @return {Number | BigNumber} Returns the distance from two/three points
   */
   return typed(name, {
-    'Array, Array, Array': function ArrayArrayArray(x, y, z) {
+    'Array, Array, Array': function (x, y, z) {
       // Point to Line 2D (x=Point, y=LinePoint1, z=LinePoint2)
       if (x.length === 2 && y.length === 2 && z.length === 2) {
         if (!_2d(x)) {
@@ -84,15 +86,15 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
         if (deepEqual(y, z)) {
           throw new TypeError('LinePoint1 should not be same with LinePoint2');
         }
-        var xCoeff = subtractScalar(z[1], y[1]);
-        var yCoeff = subtractScalar(y[0], z[0]);
-        var constant = subtractScalar(multiplyScalar(z[0], y[1]), multiplyScalar(y[0], z[1]));
+        const xCoeff = subtractScalar(z[1], y[1]);
+        const yCoeff = subtractScalar(y[0], z[0]);
+        const constant = subtractScalar(multiplyScalar(z[0], y[1]), multiplyScalar(y[0], z[1]));
         return _distancePointLine2D(x[0], x[1], xCoeff, yCoeff, constant);
       } else {
         throw new TypeError('Invalid Arguments: Try again');
       }
     },
-    'Object, Object, Object': function ObjectObjectObject(x, y, z) {
+    'Object, Object, Object': function (x, y, z) {
       if (Object.keys(x).length === 2 && Object.keys(y).length === 2 && Object.keys(z).length === 2) {
         if (!_2d(x)) {
           throw new TypeError('Values of pointX and pointY should be numbers or BigNumbers');
@@ -107,9 +109,9 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
           throw new TypeError('LinePoint1 should not be same with LinePoint2');
         }
         if ('pointX' in x && 'pointY' in x && 'lineOnePtX' in y && 'lineOnePtY' in y && 'lineTwoPtX' in z && 'lineTwoPtY' in z) {
-          var xCoeff = subtractScalar(z.lineTwoPtY, y.lineOnePtY);
-          var yCoeff = subtractScalar(y.lineOnePtX, z.lineTwoPtX);
-          var constant = subtractScalar(multiplyScalar(z.lineTwoPtX, y.lineOnePtY), multiplyScalar(y.lineOnePtX, z.lineTwoPtY));
+          const xCoeff = subtractScalar(z.lineTwoPtY, y.lineOnePtY);
+          const yCoeff = subtractScalar(y.lineOnePtX, z.lineTwoPtX);
+          const constant = subtractScalar(multiplyScalar(z.lineTwoPtX, y.lineOnePtY), multiplyScalar(y.lineOnePtX, z.lineTwoPtY));
           return _distancePointLine2D(x.pointX, x.pointY, xCoeff, yCoeff, constant);
         } else {
           throw new TypeError('Key names do not match');
@@ -118,7 +120,7 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
         throw new TypeError('Invalid Arguments: Try again');
       }
     },
-    'Array, Array': function ArrayArray(x, y) {
+    'Array, Array': function (x, y) {
       // Point to Line 2D (x=[pointX, pointY], y=[x-coeff, y-coeff, const])
       if (x.length === 2 && y.length === 3) {
         if (!_2d(x)) {
@@ -150,7 +152,7 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
         throw new TypeError('Invalid Arguments: Try again');
       }
     },
-    'Object, Object': function ObjectObject(x, y) {
+    'Object, Object': function (x, y) {
       if (Object.keys(x).length === 2 && Object.keys(y).length === 3) {
         if (!_2d(x)) {
           throw new TypeError('Values of pointX and pointY should be numbers or BigNumbers');
@@ -206,7 +208,7 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
         throw new TypeError('Invalid Arguments: Try again');
       }
     },
-    Array: function Array(arr) {
+    Array: function (arr) {
       if (!_pairwise(arr)) {
         throw new TypeError('Incorrect array format entered for pairwise distance calculation');
       }
@@ -245,9 +247,9 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
     return _isNumber(a[0]) && _isNumber(a[1]) && _isNumber(a[2]) && _isNumber(a[3]) && _isNumber(a[4]) && _isNumber(a[5]);
   }
   function _objectToArray(o) {
-    var keys = Object.keys(o);
-    var a = [];
-    for (var i = 0; i < keys.length; i++) {
+    const keys = Object.keys(o);
+    const a = [];
+    for (let i = 0; i < keys.length; i++) {
       a.push(o[keys[i]]);
     }
     return a;
@@ -255,15 +257,11 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
   function _pairwise(a) {
     // checks for valid arguments passed to _distancePairwise(Array)
     if (a[0].length === 2 && _isNumber(a[0][0]) && _isNumber(a[0][1])) {
-      if (a.some(function (aI) {
-        return aI.length !== 2 || !_isNumber(aI[0]) || !_isNumber(aI[1]);
-      })) {
+      if (a.some(aI => aI.length !== 2 || !_isNumber(aI[0]) || !_isNumber(aI[1]))) {
         return false;
       }
     } else if (a[0].length === 3 && _isNumber(a[0][0]) && _isNumber(a[0][1]) && _isNumber(a[0][2])) {
-      if (a.some(function (aI) {
-        return aI.length !== 3 || !_isNumber(aI[0]) || !_isNumber(aI[1]) || !_isNumber(aI[2]);
-      })) {
+      if (a.some(aI => aI.length !== 3 || !_isNumber(aI[0]) || !_isNumber(aI[1]) || !_isNumber(aI[2]))) {
         return false;
       }
     } else {
@@ -272,32 +270,32 @@ var createDistance = exports.createDistance = /* #__PURE__ */(0, _factory.factor
     return true;
   }
   function _distancePointLine2D(x, y, a, b, c) {
-    var num = abs(addScalar(addScalar(multiplyScalar(a, x), multiplyScalar(b, y)), c));
-    var den = sqrt(addScalar(multiplyScalar(a, a), multiplyScalar(b, b)));
+    const num = abs(addScalar(addScalar(multiplyScalar(a, x), multiplyScalar(b, y)), c));
+    const den = sqrt(addScalar(multiplyScalar(a, a), multiplyScalar(b, b)));
     return divideScalar(num, den);
   }
   function _distancePointLine3D(x, y, z, x0, y0, z0, a, b, c) {
-    var num = [subtractScalar(multiplyScalar(subtractScalar(y0, y), c), multiplyScalar(subtractScalar(z0, z), b)), subtractScalar(multiplyScalar(subtractScalar(z0, z), a), multiplyScalar(subtractScalar(x0, x), c)), subtractScalar(multiplyScalar(subtractScalar(x0, x), b), multiplyScalar(subtractScalar(y0, y), a))];
+    let num = [subtractScalar(multiplyScalar(subtractScalar(y0, y), c), multiplyScalar(subtractScalar(z0, z), b)), subtractScalar(multiplyScalar(subtractScalar(z0, z), a), multiplyScalar(subtractScalar(x0, x), c)), subtractScalar(multiplyScalar(subtractScalar(x0, x), b), multiplyScalar(subtractScalar(y0, y), a))];
     num = sqrt(addScalar(addScalar(multiplyScalar(num[0], num[0]), multiplyScalar(num[1], num[1])), multiplyScalar(num[2], num[2])));
-    var den = sqrt(addScalar(addScalar(multiplyScalar(a, a), multiplyScalar(b, b)), multiplyScalar(c, c)));
+    const den = sqrt(addScalar(addScalar(multiplyScalar(a, a), multiplyScalar(b, b)), multiplyScalar(c, c)));
     return divideScalar(num, den);
   }
   function _euclideanDistance(x, y) {
-    var vectorSize = x.length;
-    var result = 0;
-    var diff = 0;
-    for (var i = 0; i < vectorSize; i++) {
+    const vectorSize = x.length;
+    let result = 0;
+    let diff = 0;
+    for (let i = 0; i < vectorSize; i++) {
       diff = subtractScalar(x[i], y[i]);
       result = addScalar(multiplyScalar(diff, diff), result);
     }
     return sqrt(result);
   }
   function _distancePairwise(a) {
-    var result = [];
-    var pointA = [];
-    var pointB = [];
-    for (var i = 0; i < a.length - 1; i++) {
-      for (var j = i + 1; j < a.length; j++) {
+    const result = [];
+    let pointA = [];
+    let pointB = [];
+    for (let i = 0; i < a.length - 1; i++) {
+      for (let j = i + 1; j < a.length; j++) {
         if (a[0].length === 2) {
           pointA = [a[i][0], a[i][1]];
           pointB = [a[j][0], a[j][1]];

package/lib/cjs/function/geometry/intersect.js

@@ -5,24 +5,26 @@ Object.defineProperty(exports, "__esModule", {
 });
 exports.createIntersect = void 0;
 var _factory = require("../../utils/factory.js");
-var name = 'intersect';
-var dependencies = ['typed', 'config', 'abs', 'add', 'addScalar', 'matrix', 'multiply', 'multiplyScalar', 'divideScalar', 'subtract', 'smaller', 'equalScalar', 'flatten', 'isZero', 'isNumeric'];
-var createIntersect = exports.createIntersect = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
-  var typed = _ref.typed,
-    config = _ref.config,
-    abs = _ref.abs,
-    add = _ref.add,
-    addScalar = _ref.addScalar,
-    matrix = _ref.matrix,
-    multiply = _ref.multiply,
-    multiplyScalar = _ref.multiplyScalar,
-    divideScalar = _ref.divideScalar,
-    subtract = _ref.subtract,
-    smaller = _ref.smaller,
-    equalScalar = _ref.equalScalar,
-    flatten = _ref.flatten,
-    isZero = _ref.isZero,
-    isNumeric = _ref.isNumeric;
+const name = 'intersect';
+const dependencies = ['typed', 'config', 'abs', 'add', 'addScalar', 'matrix', 'multiply', 'multiplyScalar', 'divideScalar', 'subtract', 'smaller', 'equalScalar', 'flatten', 'isZero', 'isNumeric'];
+const createIntersect = exports.createIntersect = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+  let {
+    typed,
+    config,
+    abs,
+    add,
+    addScalar,
+    matrix,
+    multiply,
+    multiplyScalar,
+    divideScalar,
+    subtract,
+    smaller,
+    equalScalar,
+    flatten,
+    isZero,
+    isNumeric
+  } = _ref;
   /**
    * Calculates the point of intersection of two lines in two or three dimensions
    * and of a line and a plane in three dimensions. The inputs are in the form of
@@ -53,13 +55,13 @@ var createIntersect = exports.createIntersect = /* #__PURE__ */(0, _factory.fact
   return typed('intersect', {
     'Array, Array, Array': _AAA,
     'Array, Array, Array, Array': _AAAA,
-    'Matrix, Matrix, Matrix': function MatrixMatrixMatrix(x, y, plane) {
-      var arr = _AAA(x.valueOf(), y.valueOf(), plane.valueOf());
+    'Matrix, Matrix, Matrix': function (x, y, plane) {
+      const arr = _AAA(x.valueOf(), y.valueOf(), plane.valueOf());
       return arr === null ? null : matrix(arr);
     },
-    'Matrix, Matrix, Matrix, Matrix': function MatrixMatrixMatrixMatrix(w, x, y, z) {
+    'Matrix, Matrix, Matrix, Matrix': function (w, x, y, z) {
       // TODO: output matrix type should match input matrix type
-      var arr = _AAAA(w.valueOf(), x.valueOf(), y.valueOf(), z.valueOf());
+      const arr = _AAAA(w.valueOf(), x.valueOf(), y.valueOf(), z.valueOf());
       return arr === null ? null : matrix(arr);
     }
   });
@@ -123,9 +125,7 @@ var createIntersect = exports.createIntersect = /* #__PURE__ */(0, _factory.fact
 
     // column matrix
     if (arr.length > 1 && Array.isArray(arr[0])) {
-      if (arr.every(function (el) {
-        return Array.isArray(el) && el.length === 1;
-      })) return flatten(arr);
+      if (arr.every(el => Array.isArray(el) && el.length === 1)) return flatten(arr);
     }
     return arr;
   }
@@ -139,46 +139,46 @@ var createIntersect = exports.createIntersect = /* #__PURE__ */(0, _factory.fact
     return x.length === 4 && isNumeric(x[0]) && isNumeric(x[1]) && isNumeric(x[2]) && isNumeric(x[3]);
   }
   function _intersect2d(p1a, p1b, p2a, p2b) {
-    var o1 = p1a;
-    var o2 = p2a;
-    var d1 = subtract(o1, p1b);
-    var d2 = subtract(o2, p2b);
-    var det = subtract(multiplyScalar(d1[0], d2[1]), multiplyScalar(d2[0], d1[1]));
+    const o1 = p1a;
+    const o2 = p2a;
+    const d1 = subtract(o1, p1b);
+    const d2 = subtract(o2, p2b);
+    const det = subtract(multiplyScalar(d1[0], d2[1]), multiplyScalar(d2[0], d1[1]));
     if (isZero(det)) return null;
-    if (smaller(abs(det), config.epsilon)) {
+    if (smaller(abs(det), config.relTol)) {
       return null;
     }
-    var d20o11 = multiplyScalar(d2[0], o1[1]);
-    var d21o10 = multiplyScalar(d2[1], o1[0]);
-    var d20o21 = multiplyScalar(d2[0], o2[1]);
-    var d21o20 = multiplyScalar(d2[1], o2[0]);
-    var t = divideScalar(addScalar(subtract(subtract(d20o11, d21o10), d20o21), d21o20), det);
+    const d20o11 = multiplyScalar(d2[0], o1[1]);
+    const d21o10 = multiplyScalar(d2[1], o1[0]);
+    const d20o21 = multiplyScalar(d2[0], o2[1]);
+    const d21o20 = multiplyScalar(d2[1], o2[0]);
+    const t = divideScalar(addScalar(subtract(subtract(d20o11, d21o10), d20o21), d21o20), det);
     return add(multiply(d1, t), o1);
   }
   function _intersect3dHelper(a, b, c, d, e, f, g, h, i, j, k, l) {
     // (a - b)*(c - d) + (e - f)*(g - h) + (i - j)*(k - l)
-    var add1 = multiplyScalar(subtract(a, b), subtract(c, d));
-    var add2 = multiplyScalar(subtract(e, f), subtract(g, h));
-    var add3 = multiplyScalar(subtract(i, j), subtract(k, l));
+    const add1 = multiplyScalar(subtract(a, b), subtract(c, d));
+    const add2 = multiplyScalar(subtract(e, f), subtract(g, h));
+    const add3 = multiplyScalar(subtract(i, j), subtract(k, l));
     return addScalar(addScalar(add1, add2), add3);
   }
   function _intersect3d(x1, y1, z1, x2, y2, z2, x3, y3, z3, x4, y4, z4) {
-    var d1343 = _intersect3dHelper(x1, x3, x4, x3, y1, y3, y4, y3, z1, z3, z4, z3);
-    var d4321 = _intersect3dHelper(x4, x3, x2, x1, y4, y3, y2, y1, z4, z3, z2, z1);
-    var d1321 = _intersect3dHelper(x1, x3, x2, x1, y1, y3, y2, y1, z1, z3, z2, z1);
-    var d4343 = _intersect3dHelper(x4, x3, x4, x3, y4, y3, y4, y3, z4, z3, z4, z3);
-    var d2121 = _intersect3dHelper(x2, x1, x2, x1, y2, y1, y2, y1, z2, z1, z2, z1);
-    var numerator = subtract(multiplyScalar(d1343, d4321), multiplyScalar(d1321, d4343));
-    var denominator = subtract(multiplyScalar(d2121, d4343), multiplyScalar(d4321, d4321));
+    const d1343 = _intersect3dHelper(x1, x3, x4, x3, y1, y3, y4, y3, z1, z3, z4, z3);
+    const d4321 = _intersect3dHelper(x4, x3, x2, x1, y4, y3, y2, y1, z4, z3, z2, z1);
+    const d1321 = _intersect3dHelper(x1, x3, x2, x1, y1, y3, y2, y1, z1, z3, z2, z1);
+    const d4343 = _intersect3dHelper(x4, x3, x4, x3, y4, y3, y4, y3, z4, z3, z4, z3);
+    const d2121 = _intersect3dHelper(x2, x1, x2, x1, y2, y1, y2, y1, z2, z1, z2, z1);
+    const numerator = subtract(multiplyScalar(d1343, d4321), multiplyScalar(d1321, d4343));
+    const denominator = subtract(multiplyScalar(d2121, d4343), multiplyScalar(d4321, d4321));
     if (isZero(denominator)) return null;
-    var ta = divideScalar(numerator, denominator);
-    var tb = divideScalar(addScalar(d1343, multiplyScalar(ta, d4321)), d4343);
-    var pax = addScalar(x1, multiplyScalar(ta, subtract(x2, x1)));
-    var pay = addScalar(y1, multiplyScalar(ta, subtract(y2, y1)));
-    var paz = addScalar(z1, multiplyScalar(ta, subtract(z2, z1)));
-    var pbx = addScalar(x3, multiplyScalar(tb, subtract(x4, x3)));
-    var pby = addScalar(y3, multiplyScalar(tb, subtract(y4, y3)));
-    var pbz = addScalar(z3, multiplyScalar(tb, subtract(z4, z3)));
+    const ta = divideScalar(numerator, denominator);
+    const tb = divideScalar(addScalar(d1343, multiplyScalar(ta, d4321)), d4343);
+    const pax = addScalar(x1, multiplyScalar(ta, subtract(x2, x1)));
+    const pay = addScalar(y1, multiplyScalar(ta, subtract(y2, y1)));
+    const paz = addScalar(z1, multiplyScalar(ta, subtract(z2, z1)));
+    const pbx = addScalar(x3, multiplyScalar(tb, subtract(x4, x3)));
+    const pby = addScalar(y3, multiplyScalar(tb, subtract(y4, y3)));
+    const pbz = addScalar(z3, multiplyScalar(tb, subtract(z4, z3)));
     if (equalScalar(pax, pbx) && equalScalar(pay, pby) && equalScalar(paz, pbz)) {
       return [pax, pay, paz];
     } else {
@@ -186,18 +186,18 @@ var createIntersect = exports.createIntersect = /* #__PURE__ */(0, _factory.fact
     }
   }
   function _intersectLinePlane(x1, y1, z1, x2, y2, z2, x, y, z, c) {
-    var x1x = multiplyScalar(x1, x);
-    var x2x = multiplyScalar(x2, x);
-    var y1y = multiplyScalar(y1, y);
-    var y2y = multiplyScalar(y2, y);
-    var z1z = multiplyScalar(z1, z);
-    var z2z = multiplyScalar(z2, z);
-    var numerator = subtract(subtract(subtract(c, x1x), y1y), z1z);
-    var denominator = subtract(subtract(subtract(addScalar(addScalar(x2x, y2y), z2z), x1x), y1y), z1z);
-    var t = divideScalar(numerator, denominator);
-    var px = addScalar(x1, multiplyScalar(t, subtract(x2, x1)));
-    var py = addScalar(y1, multiplyScalar(t, subtract(y2, y1)));
-    var pz = addScalar(z1, multiplyScalar(t, subtract(z2, z1)));
+    const x1x = multiplyScalar(x1, x);
+    const x2x = multiplyScalar(x2, x);
+    const y1y = multiplyScalar(y1, y);
+    const y2y = multiplyScalar(y2, y);
+    const z1z = multiplyScalar(z1, z);
+    const z2z = multiplyScalar(z2, z);
+    const numerator = subtract(subtract(subtract(c, x1x), y1y), z1z);
+    const denominator = subtract(subtract(subtract(addScalar(addScalar(x2x, y2y), z2z), x1x), y1y), z1z);
+    const t = divideScalar(numerator, denominator);
+    const px = addScalar(x1, multiplyScalar(t, subtract(x2, x1)));
+    const py = addScalar(y1, multiplyScalar(t, subtract(y2, y1)));
+    const pz = addScalar(z1, multiplyScalar(t, subtract(z2, z1)));
     return [px, py, pz];
     // TODO: Add cases when line is parallel to the plane:
     //       (a) no intersection,

package/lib/cjs/function/logical/and.js

@@ -11,34 +11,36 @@ var _matAlgo06xS0S = require("../../type/matrix/utils/matAlgo06xS0S0.js");
 var _factory = require("../../utils/factory.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var _index = require("../../plain/number/index.js");
-var name = 'and';
-var dependencies = ['typed', 'matrix', 'equalScalar', 'zeros', 'not', 'concat'];
-var createAnd = exports.createAnd = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
-  var typed = _ref.typed,
-    matrix = _ref.matrix,
-    equalScalar = _ref.equalScalar,
-    zeros = _ref.zeros,
-    not = _ref.not,
-    concat = _ref.concat;
-  var matAlgo02xDS0 = (0, _matAlgo02xDS.createMatAlgo02xDS0)({
-    typed: typed,
-    equalScalar: equalScalar
+const name = 'and';
+const dependencies = ['typed', 'matrix', 'equalScalar', 'zeros', 'not', 'concat'];
+const createAnd = exports.createAnd = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+  let {
+    typed,
+    matrix,
+    equalScalar,
+    zeros,
+    not,
+    concat
+  } = _ref;
+  const matAlgo02xDS0 = (0, _matAlgo02xDS.createMatAlgo02xDS0)({
+    typed,
+    equalScalar
   });
-  var matAlgo06xS0S0 = (0, _matAlgo06xS0S.createMatAlgo06xS0S0)({
-    typed: typed,
-    equalScalar: equalScalar
+  const matAlgo06xS0S0 = (0, _matAlgo06xS0S.createMatAlgo06xS0S0)({
+    typed,
+    equalScalar
   });
-  var matAlgo11xS0s = (0, _matAlgo11xS0s.createMatAlgo11xS0s)({
-    typed: typed,
-    equalScalar: equalScalar
+  const matAlgo11xS0s = (0, _matAlgo11xS0s.createMatAlgo11xS0s)({
+    typed,
+    equalScalar
   });
-  var matAlgo14xDs = (0, _matAlgo14xDs.createMatAlgo14xDs)({
-    typed: typed
+  const matAlgo14xDs = (0, _matAlgo14xDs.createMatAlgo14xDs)({
+    typed
   });
-  var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
-    typed: typed,
-    matrix: matrix,
-    concat: concat
+  const matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
+    typed,
+    matrix,
+    concat
   });
 
   /**
@@ -64,75 +66,60 @@ var createAnd = exports.createAnd = /* #__PURE__ */(0, _factory.factory)(name, d
    *
    *    not, or, xor
    *
-   * @param  {number | BigNumber | Complex | Unit | Array | Matrix} x First value to check
-   * @param  {number | BigNumber | Complex | Unit | Array | Matrix} y Second value to check
+   * @param  {number | BigNumber | bigint | Complex | Unit | Array | Matrix} x First value to check
+   * @param  {number | BigNumber | bigint | Complex | Unit | Array | Matrix} y Second value to check
    * @return {boolean | Array | Matrix}
    *            Returns true when both inputs are defined with a nonzero/nonempty value.
    */
   return typed(name, {
     'number, number': _index.andNumber,
-    'Complex, Complex': function ComplexComplex(x, y) {
+    'Complex, Complex': function (x, y) {
       return (x.re !== 0 || x.im !== 0) && (y.re !== 0 || y.im !== 0);
     },
-    'BigNumber, BigNumber': function BigNumberBigNumber(x, y) {
+    'BigNumber, BigNumber': function (x, y) {
       return !x.isZero() && !y.isZero() && !x.isNaN() && !y.isNaN();
     },
-    'Unit, Unit': typed.referToSelf(function (self) {
-      return function (x, y) {
-        return self(x.value || 0, y.value || 0);
-      };
+    'bigint, bigint': _index.andNumber,
+    'Unit, Unit': typed.referToSelf(self => (x, y) => self(x.value || 0, y.value || 0)),
+    'SparseMatrix, any': typed.referToSelf(self => (x, y) => {
+      // check scalar
+      if (not(y)) {
+        // return zero matrix
+        return zeros(x.size(), x.storage());
+      }
+      return matAlgo11xS0s(x, y, self, false);
     }),
-    'SparseMatrix, any': typed.referToSelf(function (self) {
-      return function (x, y) {
-        // check scalar
-        if (not(y)) {
-          // return zero matrix
-          return zeros(x.size(), x.storage());
-        }
-        return matAlgo11xS0s(x, y, self, false);
-      };
+    'DenseMatrix, any': typed.referToSelf(self => (x, y) => {
+      // check scalar
+      if (not(y)) {
+        // return zero matrix
+        return zeros(x.size(), x.storage());
+      }
+      return matAlgo14xDs(x, y, self, false);
     }),
-    'DenseMatrix, any': typed.referToSelf(function (self) {
-      return function (x, y) {
-        // check scalar
-        if (not(y)) {
-          // return zero matrix
-          return zeros(x.size(), x.storage());
-        }
-        return matAlgo14xDs(x, y, self, false);
-      };
+    'any, SparseMatrix': typed.referToSelf(self => (x, y) => {
+      // check scalar
+      if (not(x)) {
+        // return zero matrix
+        return zeros(x.size(), x.storage());
+      }
+      return matAlgo11xS0s(y, x, self, true);
     }),
-    'any, SparseMatrix': typed.referToSelf(function (self) {
-      return function (x, y) {
-        // check scalar
-        if (not(x)) {
-          // return zero matrix
-          return zeros(x.size(), x.storage());
-        }
-        return matAlgo11xS0s(y, x, self, true);
-      };
+    'any, DenseMatrix': typed.referToSelf(self => (x, y) => {
+      // check scalar
+      if (not(x)) {
+        // return zero matrix
+        return zeros(x.size(), x.storage());
+      }
+      return matAlgo14xDs(y, x, self, true);
     }),
-    'any, DenseMatrix': typed.referToSelf(function (self) {
-      return function (x, y) {
-        // check scalar
-        if (not(x)) {
-          // return zero matrix
-          return zeros(x.size(), x.storage());
-        }
-        return matAlgo14xDs(y, x, self, true);
-      };
+    'Array, any': typed.referToSelf(self => (x, y) => {
+      // use matrix implementation
+      return self(matrix(x), y).valueOf();
     }),
-    'Array, any': typed.referToSelf(function (self) {
-      return function (x, y) {
-        // use matrix implementation
-        return self(matrix(x), y).valueOf();
-      };
-    }),
-    'any, Array': typed.referToSelf(function (self) {
-      return function (x, y) {
-        // use matrix implementation
-        return self(x, matrix(y)).valueOf();
-      };
+    'any, Array': typed.referToSelf(self => (x, y) => {
+      // use matrix implementation
+      return self(x, matrix(y)).valueOf();
     })
   }, matrixAlgorithmSuite({
     SS: matAlgo06xS0S0,