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
JavaScript
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
})