fusl/src/math/__expo2.c - Issue 1573973002: Add a "fork" of musl as //fusl.

Side by Side Diff

Use n/p to move between diff chunks; N/P to move between comments. Draft comments are only viewable by you.

Keyboard Shortcuts

	File
u :	up to issue
j / k :	jump to file after / before current file
J / K :	jump to next file with a comment after / before current file
	Side-by-side diff
i :	toggle intra-line diffs
e :	expand all comments
c :	collapse all comments
s :	toggle showing all comments
n / p :	next / previous diff chunk or comment
N / P :	next / previous comment
<Up> / <Down> :	next / previous line

	Issue
u :	up to list of issues
j / k :	jump to patch after / before current patch
o / <Enter> :	open current patch in side-by-side view
i :	open current patch in unified diff view

	Issue List
j / k :	jump to issue after / before current issue
o / <Enter> :	open current issue

Side by Side Diff: fusl/src/math/__expo2.c

Issue 1573973002: Add a "fork" of musl as //fusl. (Closed) Base URL: https://github.com/domokit/mojo.git@master

Patch Set: Created 4 years, 11 months ago

Use n/p to move between diff chunks; N/P to move between comments. Draft comments are only viewable by you.

Jump to:

OLD	NEW
(Empty)
	1 #include "libm.h"

	2

	3 /* k is such that kln2 has minimal relative error and x - kln2 > log(DBL_MIN) /

	4 static const int k = 2043;

	5 static const double kln2 = 0x1.62066151add8bp+10;

	6

	7 /* exp(x)/2 for x >= log(DBL_MAX), slightly better than 0.5exp(x/2)exp(x/2) */

	8 double __expo2(double x)

	9 {

	10 double scale;

	11

	12 /* note that k is odd and scalescale overflows /

	13 INSERT_WORDS(scale, (uint32_t)(0x3ff + k/2) << 20, 0);

	14 /* exp(x - k ln2) * 2*(k-1) /

	15 return exp(x - kln2) * scale * scale;

	16 }

OLD	NEW

« no previous file with comments | « fusl/src/math/__cosl.c ('k') | fusl/src/math/__expo2f.c » ('j') | no next file with comments »