fusl/src/math/__expo2f.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/__expo2f.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(FLT_MIN) /

	4 static const int k = 235;

	5 static const float kln2 = 0x1.45c778p+7f;

	6

	7 /* expf(x)/2 for x >= log(FLT_MAX), slightly better than 0.5fexpf(x/2)expf(x/2 ) */

	8 float __expo2f(float x)

	9 {

	10 float scale;

	11

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

	13 SET_FLOAT_WORD(scale, (uint32_t)(0x7f + k/2) << 23);

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

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

	16 }

OLD	NEW

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