YourWishes/jerryscript
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Implement missing Math Inverse Hyperbolic functions (#3675)

Browse Source 
Math.acosh, Math.asinh, Math.acosh

C implementation ported from fdlibm

Part of Issue #3568

JerryScript-DCO-1.0-Signed-off-by: Rafal Walczyna r.walczyna@samsung.com
This commit is contained in:
Rafal Walczyna
2020-04-20 15:41:06 +02:00
committed by GitHub
parent 453da11398
commit 9c7a699d10
11 changed files with 554 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files
jerry-libm/acosh.c
+92
View File
@@ -0,0 +1,92 @@
/* Copyright JS Foundation and other contributors, http://js.foundation
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
* This file is based on work under the following copyright and permission
* notice:
*
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
*
* Developed at SunSoft, a Sun Microsystems, Inc. business.
* Permission to use, copy, modify, and distribute this
* software is freely granted, provided that this notice
* is preserved.
*
* @(#)e_acosh.c 1.3 95/01/18
*/
#include "jerry-libm-internal.h"
/* acosh(x)
* Method :
* Based on
* acosh(x) = log [ x + sqrt(x * x - 1) ]
* we have
* acosh(x) := log(x) + ln2, if x is large; else
* acosh(x) := log(2x - 1 / (sqrt(x * x - 1) + x)), if x > 2; else
* acosh(x) := log1p(t + sqrt(2.0 * t + t * t)); where t = x - 1.
*
* Special cases:
* acosh(x) is NaN with signal if x < 1.
* acosh(NaN) is NaN without signal.
*/
#define one 1.0
#define ln2 6.93147180559945286227e-01 /* 0x3FE62E42, 0xFEFA39EF */
double
acosh (double x)
{
double t;
int hx;
hx = __HI (x);
if (hx < 0x3ff00000)
{
/* x < 1 */
return NAN;
}
else if (hx >= 0x41b00000)
{
/* x > 2**28 */
if (hx >= 0x7ff00000)
{
/* x is inf of NaN */
return x + x;
}
else
{
/* acosh(huge) = log(2x) */
return log (x) + ln2;
}
}
else if (((hx - 0x3ff00000) | __LO (x)) == 0)
{
/* acosh(1) = 0 */
return 0.0;
}
else if (hx > 0x40000000)
{
/* 2**28 > x > 2 */
t = x * x;
return log (2.0 * x - one / (x + sqrt (t - one)));
}
else
{
/* 1 < x < 2 */
t = x - one;
return log1p (t + sqrt (2.0 * t + t * t));
}
} /* acosh */
#undef one
#undef ln2
jerry-libm/asinh.c
+95
View File
@@ -0,0 +1,95 @@
/* Copyright JS Foundation and other contributors, http://js.foundation
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
* This file is based on work under the following copyright and permission
* notice:
*
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
*
* Developed at SunSoft, a Sun Microsystems, Inc. business.
* Permission to use, copy, modify, and distribute this
* software is freely granted, provided that this notice
* is preserved.
*
* @(#)s_asinh.c 1.3 95/01/18
*/
#include "jerry-libm-internal.h"
/* asinh(x)
* Method :
* Based on
* asinh(x) = sign(x) * log [ |x| + sqrt(x*x+1) ]
* we have
* asinh(x) := x if 1 + x * x = 1,
* := sign(x) * (log(x)+ln2)) for large |x|, else
* := sign(x) * log(2|x| + 1 / (|x| + sqrt(x * x + 1))) if|x| > 2, else
* := sign(x) * log1p(|x| + x^2 / (1 + sqrt(1 + x^2)))
*/
#define one 1.0
#define ln2 6.93147180559945286227e-01 /* 0x3FE62E42, 0xFEFA39EF */
#define huge 1.0e+300
double
asinh (double x)
{
double t, w;
int hx, ix;
hx = __HI (x);
ix = hx & 0x7fffffff;
if (ix >= 0x7ff00000)
{
/* x is inf or NaN */
return x + x;
}
if (ix < 0x3e300000)
{
/* |x| < 2**-28 */
if (huge + x > one)
{
/* return x inexact except 0 */
return x;
}
}
if (ix > 0x41b00000)
{
/* |x| > 2**28 */
w = log (fabs (x)) + ln2;
}
else if (ix > 0x40000000)
{
/* 2**28 > |x| > 2.0 */
t = fabs (x);
w = log (2.0 * t + one / (sqrt (x * x + one) + t));
}
else
{
/* 2.0 > |x| > 2**-28 */
t = x * x;
w = log1p (fabs (x) + t / (one + sqrt (one + t)));
}
if (hx > 0)
{
return w;
}
else
{
return -w;
}
} /* asinh */
#undef one
#undef ln2
#undef huge
jerry-libm/atanh.c
+100
View File
@@ -0,0 +1,100 @@
/* Copyright JS Foundation and other contributors, http://js.foundation
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
* This file is based on work under the following copyright and permission
* notice:
*
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
*
* Developed at SunSoft, a Sun Microsystems, Inc. business.
* Permission to use, copy, modify, and distribute this
* software is freely granted, provided that this notice
* is preserved.
*
* @(#)e_atanh.c 1.3 95/01/18
*/
#include "jerry-libm-internal.h"
/* atanh(x)
* Method :
* 1.Reduced x to positive by atanh(-x) = -atanh(x)
* 2.For x >= 0.5
* 1 2x x
* atanh(x) = --- * log(1 + -------) = 0.5 * log1p(2 * --------)
* 2 1 - x 1 - x
*
* For x < 0.5
* atanh(x) = 0.5 * log1p(2x + 2x * x / (1 - x))
*
* Special cases:
* atanh(x) is NaN if |x| > 1 with signal;
* atanh(NaN) is that NaN with no signal;
* atanh(+-1) is +-INF with signal.
*
*/
#define zero 0.0
#define one 1.0
#define huge 1.0e+300
double
atanh (double x)
{
double t;
int hx, ix;
double_accessor temp;
temp.dbl = x;
hx = temp.as_int.hi;
ix = hx & 0x7fffffff;
/* |x| > 1 */
if ((ix | ((unsigned int)(temp.as_int.lo | (-temp.as_int.lo)) >> 31)) > 0x3ff00000)
{
return NAN;
}
if (ix == 0x3ff00000)
{
return x / zero;
}
if (ix < 0x3e300000 && (huge + x) > zero)
{
return x; /* x<2**-28 */
}
/* x <- |x| */
temp.as_int.hi = ix;
if (ix < 0x3fe00000)
{
/* x < 0.5 */
t = temp.dbl + temp.dbl;
t = 0.5 * log1p (t + t * temp.dbl / (one - temp.dbl));
}
else
{
t = 0.5 * log1p ((temp.dbl + temp.dbl) / (one - temp.dbl));
}
if (hx >= 0)
{
return t;
}
else
{
return -t;
}
} /* atanh */
#undef zero
#undef one
#undef huge
jerry-libm/include/math.h
+5
View File
@@ -56,6 +56,11 @@ double asin (double);
double atan (double);
double atan2 (double, double);
/* Inverse hyperbolic functions */
double acosh (double);
double asinh (double);
double atanh (double);
/* Exponential and logarithmic functions. */
double exp (double);
double expm1 (double);
jerry-libm/jerry-libm-internal.h
+4
View File
@@ -88,6 +88,10 @@ double cos (double x);
double sin (double x);
double tan (double x);
double acosh (double x);
double asinh (double x);
double atanh (double x);
double exp (double x);
double expm1 (double x);
double log (double x);