[Boost-commit] svn:boost r59083 - sandbox/numeric_bindings/libs/numeric/bindings/lapack/test

Subject: [Boost-commit] svn:boost r59083 - sandbox/numeric_bindings/libs/numeric/bindings/lapack/test
From: thomas.klimpel_at_[hidden]
Date: 2010-01-16 18:07:23

Author: klimpel
Date: 2010-01-16 18:07:23 EST (Sat, 16 Jan 2010)
New Revision: 59083
URL: http://svn.boost.org/trac/boost/changeset/59083

Log:
Merge from numeric_bindings-v1
Properties modified: 
   sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_steqr.cpp   (contents, props changed)
Text files modified: 
   sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_steqr.cpp |     7 +++++--                                 
   1 files changed, 5 insertions(+), 2 deletions(-)

Modified: sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_steqr.cpp
==============================================================================
--- sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_steqr.cpp	(original)
+++ sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_steqr.cpp	2010-01-16 18:07:23 EST (Sat, 16 Jan 2010)
@@ -31,6 +31,7 @@
 
    typedef ublas::matrix<T, ublas::column_major> matrix_type ;
    typedef ublas::vector<T>                      vector_type ;
+   real_type safety_factor (1.5);
 
    // Set matrix
    matrix_type z( n, n );
@@ -47,14 +48,14 @@
      for (int j=0; j<d.size(); ++j) {
        sum += z(i,j)*z(i,j) * d(j) ;
      }
-     if (std::abs( sum - 2.0 ) > 10 * std::numeric_limits<T>::epsilon() ) return 1 ;
+     if (std::abs( sum - 2.0 ) > safety_factor*10 * std::numeric_limits<T>::epsilon() ) return 1 ;
 
      if (i>0) {
        sum = 0.0 ;
        for (int j=0; j<d.size(); ++j) {
          sum += z(i-1,j)*z(i,j) * d(j) ;
        }
-       if (std::abs( sum + 1.0 ) > 10 * std::numeric_limits<T>::epsilon() ) return 1 ;
+       if (std::abs( sum + 1.0 ) > safety_factor*10 * std::numeric_limits<T>::epsilon() ) return 1 ;
      }
    }
 
@@ -65,7 +66,9 @@
 
 int main() {
    // Run tests for different value_types
+   std::cout << "float\n" ;
    if (do_value_type<float>()) return 255;
+   std::cout << "double\n" ;
    if (do_value_type<double>()) return 255;
 
    std::cout << "Regression test succeeded\n" ;