Bruce/tomo
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Activity

Fix some minor floating point promotion issues

Browse Source
This commit is contained in:
Bruce Hill 2025-03-28 13:46:35 -04:00
parent a79b3c2216
commit 9bc18cb617
3 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/stdlib/nums.c
View File

@ -125,11 +125,11 @@ public CONSTFUNC bool Num32$near(float a, float b, float ratio, float absolute)
if (a == b) return true;
float diff = fabs(a - b);
float diff = fabsf(a - b);
if (diff < absolute) return true;
else if (isnan(diff)) return false;
float epsilon = fabs(a * ratio) + fabs(b * ratio);
float epsilon = fabsf(a * ratio) + fabsf(b * ratio);
if (isinf(epsilon)) epsilon = FLT_MAX;
return (diff < epsilon);
}

2
src/stdlib/nums.h
View File

@ -100,7 +100,7 @@ MACROLIKE CONSTFUNC float Num32$from_int(Int_t i, bool truncate) {
float ret = (float)mpz_get_d(*i.big);
if (!truncate) {
mpz_t roundtrip;
mpz_init_set_d(roundtrip, ret);
mpz_init_set_d(roundtrip, (double)ret);
if unlikely (mpz_cmp(*i.big, roundtrip) != 0)
fail("Could not convert integer to 32-bit floating point without losing precision: ", i);
}

4
src/stdlib/optionals.c
View File

@ -81,9 +81,9 @@ public void Optional$deserialize(FILE *in, void *outval, Array_t *pointers, cons
else if (nonnull->tag == TableInfo)
*(Table_t*)outval = (Table_t){.entries={.length=-1}};
else if (nonnull == &Num$info)
*(double*)outval = NAN;
*(double*)outval = (double)NAN;
else if (nonnull == &Num32$info)
*(float*)outval = NAN;
*(float*)outval = (float)NAN;
else if (nonnull->tag == StructInfo || (nonnull->tag == OpaqueInfo && type->size > nonnull->size))
memset(outval + type->size, -1, (size_t)(type->size - nonnull->size));
else