UNPKG

mathjs

Version:

Math.js is an extensive math library for JavaScript and Node.js. It features a flexible expression parser with support for symbolic computation, comes with a large set of built-in functions and constants, and offers an integrated solution to work with dif

199 lines (171 loc) 6.65 kB
"use strict"; Object.defineProperty(exports, "__esModule", { value: true }); exports.createSimplifyCore = void 0; var _is = require("../../../utils/is.js"); var _factory = require("../../../utils/factory.js"); var name = 'simplifyCore'; var dependencies = ['equal', 'isZero', 'add', 'subtract', 'multiply', 'divide', 'pow', 'ConstantNode', 'OperatorNode', 'FunctionNode', 'ParenthesisNode']; var createSimplifyCore = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) { var equal = _ref.equal, isZero = _ref.isZero, add = _ref.add, subtract = _ref.subtract, multiply = _ref.multiply, divide = _ref.divide, pow = _ref.pow, ConstantNode = _ref.ConstantNode, OperatorNode = _ref.OperatorNode, FunctionNode = _ref.FunctionNode, ParenthesisNode = _ref.ParenthesisNode; var node0 = new ConstantNode(0); var node1 = new ConstantNode(1); /** * simplifyCore() performs single pass simplification suitable for * applications requiring ultimate performance. In contrast, simplify() * extends simplifyCore() with additional passes to provide deeper * simplification. * * Syntax: * * simplify.simplifyCore(expr) * * Examples: * * const f = math.parse('2 * 1 * x ^ (2 - 1)') * math.simplify.simpifyCore(f) // Node {2 * x} * math.simplify('2 * 1 * x ^ (2 - 1)', [math.simplify.simpifyCore]) // Node {2 * x} * * See also: * * derivative * * @param {Node} node * The expression to be simplified */ function simplifyCore(node) { if ((0, _is.isOperatorNode)(node) && node.isUnary()) { var a0 = simplifyCore(node.args[0]); if (node.op === '+') { // unary plus return a0; } if (node.op === '-') { // unary minus if ((0, _is.isOperatorNode)(a0)) { if (a0.isUnary() && a0.op === '-') { return a0.args[0]; } else if (a0.isBinary() && a0.fn === 'subtract') { return new OperatorNode('-', 'subtract', [a0.args[1], a0.args[0]]); } } return new OperatorNode(node.op, node.fn, [a0]); } } else if ((0, _is.isOperatorNode)(node) && node.isBinary()) { var _a = simplifyCore(node.args[0]); var a1 = simplifyCore(node.args[1]); if (node.op === '+') { if ((0, _is.isConstantNode)(_a)) { if (isZero(_a.value)) { return a1; } else if ((0, _is.isConstantNode)(a1)) { return new ConstantNode(add(_a.value, a1.value)); } } if ((0, _is.isConstantNode)(a1) && isZero(a1.value)) { return _a; } if ((0, _is.isOperatorNode)(a1) && a1.isUnary() && a1.op === '-') { return new OperatorNode('-', 'subtract', [_a, a1.args[0]]); } return new OperatorNode(node.op, node.fn, a1 ? [_a, a1] : [_a]); } else if (node.op === '-') { if ((0, _is.isConstantNode)(_a) && a1) { if ((0, _is.isConstantNode)(a1)) { return new ConstantNode(subtract(_a.value, a1.value)); } else if (isZero(_a.value)) { return new OperatorNode('-', 'unaryMinus', [a1]); } } // if (node.fn === "subtract" && node.args.length === 2) { if (node.fn === 'subtract') { if ((0, _is.isConstantNode)(a1) && isZero(a1.value)) { return _a; } if ((0, _is.isOperatorNode)(a1) && a1.isUnary() && a1.op === '-') { return simplifyCore(new OperatorNode('+', 'add', [_a, a1.args[0]])); } return new OperatorNode(node.op, node.fn, [_a, a1]); } } else if (node.op === '*') { if ((0, _is.isConstantNode)(_a)) { if (isZero(_a.value)) { return node0; } else if (equal(_a.value, 1)) { return a1; } else if ((0, _is.isConstantNode)(a1)) { return new ConstantNode(multiply(_a.value, a1.value)); } } if ((0, _is.isConstantNode)(a1)) { if (isZero(a1.value)) { return node0; } else if (equal(a1.value, 1)) { return _a; } else if ((0, _is.isOperatorNode)(_a) && _a.isBinary() && _a.op === node.op) { var a00 = _a.args[0]; if ((0, _is.isConstantNode)(a00)) { var a00a1 = new ConstantNode(multiply(a00.value, a1.value)); return new OperatorNode(node.op, node.fn, [a00a1, _a.args[1]], node.implicit); // constants on left } } return new OperatorNode(node.op, node.fn, [a1, _a], node.implicit); // constants on left } return new OperatorNode(node.op, node.fn, [_a, a1], node.implicit); } else if (node.op === '/') { if ((0, _is.isConstantNode)(_a)) { if (isZero(_a.value)) { return node0; } else if ((0, _is.isConstantNode)(a1) && (equal(a1.value, 1) || equal(a1.value, 2) || equal(a1.value, 4))) { return new ConstantNode(divide(_a.value, a1.value)); } } return new OperatorNode(node.op, node.fn, [_a, a1]); } else if (node.op === '^') { if ((0, _is.isConstantNode)(a1)) { if (isZero(a1.value)) { return node1; } else if (equal(a1.value, 1)) { return _a; } else { if ((0, _is.isConstantNode)(_a)) { // fold constant return new ConstantNode(pow(_a.value, a1.value)); } else if ((0, _is.isOperatorNode)(_a) && _a.isBinary() && _a.op === '^') { var a01 = _a.args[1]; if ((0, _is.isConstantNode)(a01)) { return new OperatorNode(node.op, node.fn, [_a.args[0], new ConstantNode(multiply(a01.value, a1.value))]); } } } } return new OperatorNode(node.op, node.fn, [_a, a1]); } } else if ((0, _is.isParenthesisNode)(node)) { var c = simplifyCore(node.content); if ((0, _is.isParenthesisNode)(c) || (0, _is.isSymbolNode)(c) || (0, _is.isConstantNode)(c)) { return c; } return new ParenthesisNode(c); } else if ((0, _is.isFunctionNode)(node)) { var args = node.args.map(simplifyCore).map(function (arg) { return (0, _is.isParenthesisNode)(arg) ? arg.content : arg; }); return new FunctionNode(simplifyCore(node.fn), args); } else {// cannot simplify } return node; } return simplifyCore; }); exports.createSimplifyCore = createSimplifyCore;