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

88 lines (79 loc) 2.57 kB
import { factory } from '../../utils/factory' import { flatten } from '../../utils/array' const name = 'hypot' const dependencies = [ 'typed', 'abs', 'addScalar', 'divideScalar', 'multiplyScalar', 'sqrt', 'smaller', 'isPositive' ] export const createHypot = /* #__PURE__ */ factory(name, dependencies, ({ typed, abs, addScalar, divideScalar, multiplyScalar, sqrt, smaller, isPositive }) => { /** * Calculate the hypotenusa of a list with values. The hypotenusa is defined as: * * hypot(a, b, c, ...) = sqrt(a^2 + b^2 + c^2 + ...) * * For matrix input, the hypotenusa is calculated for all values in the matrix. * * Syntax: * * math.hypot(a, b, ...) * math.hypot([a, b, c, ...]) * * Examples: * * math.hypot(3, 4) // 5 * math.hypot(3, 4, 5) // 7.0710678118654755 * math.hypot([3, 4, 5]) // 7.0710678118654755 * math.hypot(-2) // 2 * * See also: * * abs, norm * * @param {... number | BigNumber | Array | Matrix} args A list with numeric values or an Array or Matrix. * Matrix and Array input is flattened and returns a * single number for the whole matrix. * @return {number | BigNumber} Returns the hypothenusa of the input values. */ const hypot = typed(name, { '... number | BigNumber': _hypot, Array: function (x) { return hypot.apply(hypot, flatten(x)) }, Matrix: function (x) { return hypot.apply(hypot, flatten(x.toArray())) } }) /** * Calculate the hypotenusa for an Array with values * @param {Array.<number | BigNumber>} args * @return {number | BigNumber} Returns the result * @private */ function _hypot (args) { // code based on `hypot` from es6-shim: // https://github.com/paulmillr/es6-shim/blob/master/es6-shim.js#L1619-L1633 let result = 0 let largest = 0 for (let i = 0; i < args.length; i++) { const value = abs(args[i]) if (smaller(largest, value)) { result = multiplyScalar(result, multiplyScalar(divideScalar(largest, value), divideScalar(largest, value))) result = addScalar(result, 1) largest = value } else { result = addScalar(result, isPositive(value) ? multiplyScalar(divideScalar(value, largest), divideScalar(value, largest)) : value) } } return multiplyScalar(largest, sqrt(result)) } return hypot })