mathjs/lib/function/trigonometry/acsch.js

'use strict';

module.exports = function (math) {
  var util = require('../../util/index'),

      BigNumber = math.type.BigNumber,
      Complex = require('../../type/Complex'),
      Unit = require('../../type/Unit'),
      collection = require('../../type/collection'),

      isNumber = util.number.isNumber,
      isBoolean = util['boolean'].isBoolean,
      isComplex = Complex.isComplex,
      isUnit = Unit.isUnit,
      isCollection = collection.isCollection,

      bigAcsch = util.bignumber.acosh_asinh_asech_acsch;

  /**
   * Calculate the hyperbolic arccosecant of a value,
   * defined as `acsch(x) = ln(1/x + sqrt(1/x^2 + 1))`.
   *
   * For matrices, the function is evaluated element wise.
   *
   * Syntax:
   *
   *    math.acsch(x)
   *
   * Examples:
   *
   *    math.acsch(0.5);       // returns 1.4436354751788103
   *
   * See also:
   *
   *    asech, acoth
   *
   * @param {Number | Boolean | Complex | Unit | Array | Matrix | null} x  Function input
   * @return {Number | Complex | Array | Matrix} Hyperbolic arccosecant of x
   */
  math.acsch = function acsch(x) {
    if (arguments.length != 1) {
      throw new math.error.ArgumentsError('acsch', arguments.length, 1);
    }

    if (isNumber(x)) {
      x = 1 / x;
      return Math.log(x + Math.sqrt(x*x + 1));
    }

    if (isComplex(x)) {
      // acsch(z) = -i*asinh(1/z)
      var den = x.re*x.re + x.im*x.im;
      if (den != 0) {
        x.re =  x.re / den;
        x.im = -x.im / den;
      } else {
        x.re = (x.re != 0) ?  (x.re / 0) : 0;
        x.im = (x.im != 0) ? -(x.im / 0) : 0;
      }

      if (x.im) {
        return math.asinh(x);
      }

      x = 1 / x.re;
      return new Complex(Math.log(x + Math.sqrt(x*x + 1)), 0);
    }

    if (isCollection(x)) {
      return collection.deepMap(x, acsch);
    }

    if (isBoolean(x) || x === null) {
      return (x) ? 0.881373587019543 : Infinity;
    }

    if (x instanceof BigNumber) {
      return bigAcsch(x, true, true);
    }

    throw new math.error.UnsupportedTypeError('acsch', math['typeof'](x));
  };
};