[Boost-commit] svn:boost r58832 - sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/driver

Subject: [Boost-commit] svn:boost r58832 - sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/driver
From: rutger_at_[hidden]
Date: 2010-01-09 07:51:23

Author: rutger
Date: 2010-01-09 07:51:22 EST (Sat, 09 Jan 2010)
New Revision: 58832
URL: http://svn.boost.org/trac/boost/changeset/58832

Log:
fixed nlvl computation in gelsd

Text files modified: 
   sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/driver/gelsd.hpp |     2 +-                                      
   1 files changed, 1 insertions(+), 1 deletions(-)

Modified: sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/driver/gelsd.hpp
==============================================================================
--- sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/driver/gelsd.hpp	(original)
+++ sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/driver/gelsd.hpp	2010-01-09 07:51:22 EST (Sat, 09 Jan 2010)
@@ -9,7 +9,7 @@
 $TEMPLATE[gelsd.all.SMLSIZ.init]
 $INTEGER_TYPE smlsiz = ilaenv(9, "GELSD", "");
 $TEMPLATE[gelsd.all.NLVL.init]
-$INTEGER_TYPE nlvl = static_cast<$INTEGER_TYPE>(((std::log(static_cast<real_type>(minmn)) / std::log(static_cast<real_type>(2.))) / (smlsiz+1)) + 1);
+$INTEGER_TYPE nlvl = static_cast<$INTEGER_TYPE>(std::log(static_cast<real_type>(minmn)/static_cast<real_type>(smlsiz+1))/std::log(2.0)) + 1;
 $TEMPLATE[gelsd.complex.min_size_rwork]
 $INTEGER_TYPE smlsiz_plus_one = smlsiz + 1;
 return std::max< $INTEGER_TYPE >( 1, 10*minmn + 2*minmn*smlsiz + 8*minmn*nlvl + 3*smlsiz*nrhs + smlsiz_plus_one * smlsiz_plus_one );