Mercurial > mxe-octave
view src/build-msvctools/math/asinhf.c @ 6161:8fcac4d6d983 release
of-queueing: Fix syntax error with Octave 7 (bug #62314).
* src/of-queueing-1-octave7.patch: Add new file.
* dist-files.mk: Include new patch.
author | Markus Mützel <markus.muetzel@gmx.de> |
---|---|
date | Thu, 14 Apr 2022 19:27:13 +0200 |
parents | f8299bb6c872 |
children |
line wrap: on
line source
#include <math.h> #include <errno.h> #include "fastmath.h" /* asinh(x) = copysign(log(fabs(x) + sqrt(x * x + 1.0)), x) */ float asinhf(float x) { float z; if (!isfinite (x)) return x; z = fabsf (x); /* Avoid setting FPU underflow exception flag in x * x. */ #if 0 if ( z < 0x1p-32) return x; #endif /* Use log1p to avoid cancellation with small x. Put x * x in denom, so overflow is harmless. asinh(x) = log1p (x + sqrt (x * x + 1.0) - 1.0) = log1p (x + x * x / (sqrt (x * x + 1.0) + 1.0)) */ z = __fast_log1p (z + z * z / (__fast_sqrt (z * z + 1.0) + 1.0)); return ( x > 0.0 ? z : -z); }