From d203d0e438cc0ea99f2eb6aabb45c9c82edfde5b Mon Sep 17 00:00:00 2001
From: Erik Schnetter <schnetter@gmail.com>
Date: Wed, 3 Jul 2013 14:40:41 -0400
Subject: Take maximum vector element into account when testing vector
 reductions

---
 test.cc | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

(limited to 'test.cc')

diff --git a/test.cc b/test.cc
index 308b0ab..58346c3 100644
--- a/test.cc
+++ b/test.cc
@@ -467,7 +467,11 @@ struct vecmathlib_test {
                          const real_t accuracy)
   {
     const real_t dr = rstd - rvml;
-    const real_t scale = fabs(rstd) + fabs(rvml) + R(1.0);
+    real_t maxabs = 0.0;
+    for (int i=0; i<realvec_t::size; ++i) {
+      maxabs = fmax(maxabs, fabs(x[i]));
+    }
+    const real_t scale = fabs(rstd) + fabs(rvml) + fabs(maxabs) + R(1.0);
     const bool isbad = fabs(dr) > accuracy * scale;
     if (isbad) {
       ++ num_errors;
-- 
cgit v1.1