jepler/micropython-ulab
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

numpy/vector.c: remove usage of fpclassify (#636)

Browse source 
Fixes https://github.com/v923z/micropython-ulab/issues/635

Verified by re-compiling circuitpython with this change.
This commit is contained in:
KB Sriram 2023-07-17 13:03:30 -07:00 committed by GitHub
parent d025aa3cf6
commit 84f99f17fc
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
3 changed files with 9 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
code/numpy/vector.c
View file

@ -827,7 +827,7 @@ MP_DEFINE_CONST_FUN_OBJ_1(vector_sin_obj, vector_sin);
//| //|
static mp_float_t vector_sinc1(mp_float_t x) { static mp_float_t vector_sinc1(mp_float_t x) {
if (fpclassify(x) == FP_ZERO) { if (x == MICROPY_FLOAT_CONST(0.)) {
return MICROPY_FLOAT_CONST(1.); return MICROPY_FLOAT_CONST(1.);
} }
x *= MP_PI; x *= MP_PI;

7
tests/1d/numpy/universal_functions.py
View file

@ -112,6 +112,13 @@ for i in range(len(x)):
cmp_result.append(math.isclose(result[i], ref_result[i], rel_tol=1E-9, abs_tol=1E-9)) cmp_result.append(math.isclose(result[i], ref_result[i], rel_tol=1E-9, abs_tol=1E-9))
print(cmp_result) print(cmp_result)
result = np.sinc(x)
ref_result = np.array([0.03935584386392389, -0.04359862862918773, 1.0, -0.04359862862918773, 0.03935584386392389])
cmp_result = []
for i in range(len(x)):
cmp_result.append(math.isclose(result[i], ref_result[i], rel_tol=1E-9, abs_tol=1E-9))
print(cmp_result)
result = (spy.special.erf(np.linspace(-3, 3, num=5))) result = (spy.special.erf(np.linspace(-3, 3, num=5)))
ref_result = np.array([-0.9999779095030014, -0.9661051464753108, 0.0, 0.9661051464753108, 0.9999779095030014], dtype=np.float) ref_result = np.array([-0.9999779095030014, -0.9661051464753108, 0.0, 0.9661051464753108, 0.9999779095030014], dtype=np.float)
cmp_result = [] cmp_result = []

1
tests/1d/numpy/universal_functions.py.exp
View file

@ -28,5 +28,6 @@ True
[True, True, True, True, True] [True, True, True, True, True]
[True, True, True, True, True] [True, True, True, True, True]
[True, True, True, True, True] [True, True, True, True, True]
[True, True, True, True, True]
[True, True, True, True] [True, True, True, True]
[True, True, True] [True, True, True]