Main Conformance: Update dot product error bound
Update the error bound for dot product operations to
allow floating point accumulators where denormal
values are flushed to zero.
Change-Id: I1d5805f0ad613ca9d78c6a6da12e31d7df5f6cc4
Signed-off-by: Dominic Symes <dominic.symes@arm.com>
diff --git a/chapters/introduction.adoc b/chapters/introduction.adoc
index 26fef0e..d415e24 100644
--- a/chapters/introduction.adoc
+++ b/chapters/introduction.adoc
@@ -358,8 +358,8 @@
REQUIRE(out_ref == 0.0 && out_imp == 0.0);
out_err = 0.0;
} else { // 0.0 < out_bnd < infinity
- out_bnd = max(out_bnd, acc_min_normal);
- out_err = (static_cast<fp64_t>(out_imp) - out_ref) * acc_prec / out_bnd;
+ fp64_t out_err_bnd = max(out_bnd / acc_prec, acc_min_normal);
+ out_err = (static_cast<fp64_t>(out_imp) - out_ref) / out_err_bnd;
REQUIRE(abs(out_err) <= ksb);
}
out_err_sum += out_err;