diff options
Diffstat (limited to 'pseudocode')
-rw-r--r-- | pseudocode/library/generic_helpers.tosac | 6 | ||||
-rw-r--r-- | pseudocode/library/numeric_accuracy_helpers.tosac | 14 |
2 files changed, 20 insertions, 0 deletions
diff --git a/pseudocode/library/generic_helpers.tosac b/pseudocode/library/generic_helpers.tosac index 6dc2755..cffd55c 100644 --- a/pseudocode/library/generic_helpers.tosac +++ b/pseudocode/library/generic_helpers.tosac @@ -69,3 +69,9 @@ in_out_t maximum_u<in_out_t>(); // return the minimum value when interpreting type in_out_t as an unsigned value as returned by the make_unsigned helper. in_out_t minimum_u<in_out_t>(); + +// return true if the given value is a NaN. Only valid for floating-point types +bool is_a_NaN(fp64_t value); + +// return true if value is a normal fp64 value (Not zero, subnormal, infinite or NaN) +bool is_normal_fp64(fp64_t value);
\ No newline at end of file diff --git a/pseudocode/library/numeric_accuracy_helpers.tosac b/pseudocode/library/numeric_accuracy_helpers.tosac index b89d898..dbff2dd 100644 --- a/pseudocode/library/numeric_accuracy_helpers.tosac +++ b/pseudocode/library/numeric_accuracy_helpers.tosac @@ -56,3 +56,17 @@ int normal_frac<in_t> () { case fp8e5m2_t: return 2; } } + +double calcAbsErrorBound<in_t>(double bound_magnitude, double bounds_value, + double lower_bound, double normal_divisor) { + double error_bound = 0.0; + // Avoid cases where we generate an error_bound of NaN by multiplying inf * 0 + if (is_finite(bounds_value) || abs(bound_magnitude) != 0.0) { + double value_bound = abs(bound_magnitude) * bounds_value; + if (lower_bound > 0) { + value_bound = max(lower_bound, value_bound); + } + error_bound = exp2(-normal_frac<in_t> / normal_divisor) * value_bound; + } + return error_bound; +}
\ No newline at end of file |