[Boost-commit] svn:boost r50297 - sandbox/libs/numeric/bindings/blas/test

Subject: [Boost-commit] svn:boost r50297 - sandbox/libs/numeric/bindings/blas/test
From: thomas.klimpel_at_[hidden]
Date: 2008-12-16 09:51:56

Author: klimpel
Date: 2008-12-16 09:51:55 EST (Tue, 16 Dec 2008)
New Revision: 50297
URL: http://svn.boost.org/trac/boost/changeset/50297

Log:
minor modifications, to let regression tests pass with gfortran

Text files modified: 
   sandbox/libs/numeric/bindings/blas/test/blas1.cpp |    14 ++++++++------                          
   1 files changed, 8 insertions(+), 6 deletions(-)

Modified: sandbox/libs/numeric/bindings/blas/test/blas1.cpp
==============================================================================
--- sandbox/libs/numeric/bindings/blas/test/blas1.cpp	(original)
+++ sandbox/libs/numeric/bindings/blas/test/blas1.cpp	2008-12-16 09:51:55 EST (Tue, 16 Dec 2008)
@@ -241,22 +241,23 @@
   template <typename W>
   int operator()(W& w) const {
      using namespace boost::numeric::bindings::blas ;
+     real_type safety_factor (1.5);
 
      copy_vector( w );
 
      // Test blas routines
      value_type prod = dotc( this->v_, w );
      if ( std::abs(prod - inner_prod( conj(this->v1_ref_), this->v2_ref_ ))
-          > std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
+          > safety_factor*std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
 
      prod = dotu( this->v_, w );
      if ( std::abs(prod - inner_prod( this->v1_ref_, this->v2_ref_ ))
-          > std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
+          > safety_factor*std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
 
      axpy( value_type(2.0), this->v_, w );
      for (size_t i=0; i<this->size(); ++i)
         if ( std::abs(w[i] - (this->v2_ref_(i) + value_type(2.0)*this->v1_ref_(i)))
-          > std::numeric_limits< real_type >::epsilon() * std::abs(w[i])) return 255 ;
+          > safety_factor*std::numeric_limits< real_type >::epsilon() * std::abs(w[i])) return 255 ;
 
      scal( value_type(0.0), w ) ;
      copy( this->v_, w ) ;
@@ -285,22 +286,23 @@
   template <typename W>
   int operator()(W& w) const {
      using namespace boost::numeric::bindings::blas ;
+     real_type safety_factor (1.5);
 
      copy_vector( w );
 
      // Test blas routines
      value_type prod = dotc( this->v_, w );
      if ( std::abs(prod - inner_prod( conj(this->v1_ref_), this->v2_ref_ ))
-          > std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
+          > safety_factor*std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
 
      prod = dotu( this->v_, w );
      if ( std::abs(prod - inner_prod( this->v1_ref_, this->v2_ref_ ))
-          > std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
+          > safety_factor*std::numeric_limits< real_type >::epsilon() * std::abs(prod)) return 255 ;
 
      axpy( value_type(2.0), this->v_, w );
      for (size_t i=0; i<this->size(); ++i)
         if ( std::abs(w[i] - (this->v2_ref_(i) + value_type(2.0)*this->v1_ref_(i)))
-          > std::numeric_limits< real_type >::epsilon() * std::abs(w[i])) return 255 ;
+          > safety_factor*std::numeric_limits< real_type >::epsilon() * std::abs(w[i])) return 255 ;
 
      scal( value_type(0.0), w ) ;
      copy( this->v_, w ) ;