From 8bf6e713d5e9985886df581b2b00c57752d18022 Mon Sep 17 00:00:00 2001 From: Dominic Symes Date: Thu, 7 Sep 2023 11:10:59 +0100 Subject: Main Conformance: Add local_bound to RFFT2D Also express the fp32_t acc_min_normal calculation without a double negation. Signed-off-by: Dominic Symes Change-Id: I34ba274264bc511e60f0e01601d5d0a3f983e594 --- chapters/introduction.adoc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'chapters') diff --git a/chapters/introduction.adoc b/chapters/introduction.adoc index d415e24..c8cfc30 100644 --- a/chapters/introduction.adoc +++ b/chapters/introduction.adoc @@ -336,10 +336,10 @@ fp64_t out_err_sum = 0.0; fp64_t out_err_sumsq = 0.0; fp64_t acc_prec; // 1<<(M+1) where M is the number of mantissa bits fp64_t acc_min_normal; // accumulator minimum normal greater than zero -fp64_t two_m63 = -1.0/(fp64)((int64_t)-1<<63); // pow(2, -63) +fp64_t two_m42 = 1.0/static_cast((int64_t)1<<42); // pow(2, -42) switch (acc_t) { case fp32_t: acc_prec = static_cast(1<<24); // pow(2, 24) - acc_min_normal = two_m63 * two_m63; // pow(2, -126) + acc_min_normal = two_m42 * two_m42 * two_m42; // pow(2, -126) break; case fp16_t: acc_prec = static_cast(1<<11); // pow(2, 11) acc_min_normal = 1.0/static_cast(1<<14); // pow(2, -14) -- cgit v1.2.1