orastron/brickworks
2
1
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases 9 Wiki Activity

bw_math: more improvements, more tests (prev commit: also fixes)

Browse Source
This commit is contained in:
Stefano D'Angelo 2023-08-09 19:17:54 +02:00
parent c25bd5bc78
commit 88fad66dd6
2 changed files with 624 additions and 482 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
include/bw_math.h
View File

@ -265,6 +265,8 @@ static inline float bw_sin2pif_3(float x);
* *
* Absolute error < 0.011, relative error < 1.7%. * Absolute error < 0.011, relative error < 1.7%.
* *
* `x` must be finite.
*
* #### bw_sinf_3() * #### bw_sinf_3()
* ```>>> */ * ```>>> */
static inline float bw_sinf_3(float x); static inline float bw_sinf_3(float x);
@ -274,6 +276,8 @@ static inline float bw_sinf_3(float x);
* *
* Absolute error < 0.011, relative error < 1.7%. * Absolute error < 0.011, relative error < 1.7%.
* *
* `x` must be finite.
*
* #### bw_cos2pif_3() * #### bw_cos2pif_3()
* ```>>> */ * ```>>> */
static inline float bw_cos2pif_3(float x); static inline float bw_cos2pif_3(float x);
@ -283,6 +287,8 @@ static inline float bw_cos2pif_3(float x);
* *
* Absolute error < 0.011, relative error < 1.7%. * Absolute error < 0.011, relative error < 1.7%.
* *
* `x` must be finite.
*
* #### bw_cosf_3() * #### bw_cosf_3()
* ```>>> */ * ```>>> */
static inline float bw_cosf_3(float x); static inline float bw_cosf_3(float x);
@ -292,6 +298,8 @@ static inline float bw_cosf_3(float x);
* *
* Absolute error < 0.011, relative error < 1.7%. * Absolute error < 0.011, relative error < 1.7%.
* *
* `x` must be finite.
*
* #### bw_tan2pif_3() * #### bw_tan2pif_3()
* ```>>> */ * ```>>> */
static inline float bw_tan2pif_3(float x); static inline float bw_tan2pif_3(float x);
@ -305,6 +313,8 @@ static inline float bw_tan2pif_3(float x);
* *
* Absolute error < 0.06, relative error < 0.8%. * Absolute error < 0.06, relative error < 0.8%.
* *
* `x` must be finite.
*
* #### bw_tanf_3() * #### bw_tanf_3()
* ```>>> */ * ```>>> */
static inline float bw_tanf_3(float x); static inline float bw_tanf_3(float x);
@ -318,6 +328,8 @@ static inline float bw_tanf_3(float x);
* *
* Absolute error < 0.06, relative error < 0.8%. * Absolute error < 0.06, relative error < 0.8%.
* *
* `x` must be finite.
*
* #### bw_log2f_3() * #### bw_log2f_3()
* ```>>> */ * ```>>> */
static inline float bw_log2f_3(float x); static inline float bw_log2f_3(float x);
@ -611,11 +623,14 @@ static inline void bw_intfracf(float x, float *i, float *f) {
} }
static inline float bw_rcpf_2(float x) { static inline float bw_rcpf_2(float x) {
BW_ASSERT(bw_is_finite(x));
BW_ASSERT((x >= 8.077935669e-28f && x <= 1.237940039e27f) || (x <= -8.077935669e-28f && x >= -1.237940039e27f));
union { float f; int32_t i; } v; union { float f; int32_t i; } v;
v.f = x; v.f = x;
v.i = 0x7ef0e840 - v.i; v.i = 0x7ef0e840 - v.i;
v.f = v.f + v.f - x * v.f * v.f; v.f = v.f + v.f - x * v.f * v.f;
v.f = v.f + v.f - x * v.f * v.f; v.f = v.f + v.f - x * v.f * v.f;
BW_ASSERT(bw_is_finite(v.f));
return v.f; return v.f;
} }

1091
test/bw_math.c
View File

File diff suppressed because it is too large Load Diff