|
|
Boost-Commit : |
Subject: [Boost-commit] svn:boost r51659 - in sandbox/numeric_bindings: boost/numeric/bindings/blas/level1 boost/numeric/bindings/lapack/computational boost/numeric/bindings/lapack/driver libs/numeric/bindings/tools/templates/computational
From: rutger_at_[hidden]
Date: 2009-03-09 09:41:18
Author: rutger
Date: 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
New Revision: 51659
URL: http://svn.boost.org/trac/boost/changeset/51659
Log:
Added workspace info of geqp3
Added:
sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/computational/geqp3.hpp (contents, props changed)
Text files modified:
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotmg.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbcon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbtrf.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebak.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebal.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gecon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/geqp3.hpp | 16 ++++++++--------
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getrf.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getri.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hecon.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetrf.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetri.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/lacon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbcon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbequ.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbstf.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbtrf.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pocon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/poequ.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potrf.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potri.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sycon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytrf.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytri.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/tbcon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trcon.hpp | 2 --
sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trtri.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heev.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heevd.hpp | 1 -
sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/lacgv.hpp | 1 -
34 files changed, 8 insertions(+), 52 deletions(-)
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -48,7 +48,6 @@
// templated specialization
template< typename VectorX >
static return_type compute( VectorX& x ) {
-
return detail::asum( traits::vector_size(x),
traits::vector_storage(x), traits::vector_stride(x) );
}
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -48,7 +48,6 @@
// templated specialization
template< typename VectorX >
static return_type compute( VectorX& x ) {
-
return detail::nrm2( traits::vector_size(x),
traits::vector_storage(x), traits::vector_stride(x) );
}
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -59,7 +59,6 @@
template< >
static return_type compute( value_type& a, value_type& b, real_type& c,
value_type& s ) {
-
detail::rotg( a, b, c, s );
}
};
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotmg.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotmg.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotmg.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -51,7 +51,6 @@
template< typename VectorDPARAM >
static return_type compute( real_type& d1, real_type& d2, real_type& x1,
real_type const y1, VectorDPARAM& dparam ) {
-
detail::rotmg( d1, d2, x1, y1, traits::vector_storage(dparam) );
}
};
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -60,7 +60,6 @@
// templated specialization
template< typename VectorX >
static return_type compute( value_type const a, VectorX& x ) {
-
detail::scal( traits::vector_size(x), a,
traits::vector_storage(x), traits::vector_stride(x) );
}
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbcon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbcon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbcon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -85,7 +85,6 @@
integer_t const kl, integer_t const ku, MatrixAB& ab,
VectorIPIV& ipiv, real_type const anorm, real_type& rcond,
integer_t& info, detail::workspace2< WORK, IWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( n >= 0 );
@@ -150,7 +149,6 @@
integer_t const kl, integer_t const ku, MatrixAB& ab,
VectorIPIV& ipiv, real_type const anorm, real_type& rcond,
integer_t& info, detail::workspace2< WORK, RWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( n >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbtrf.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbtrf.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gbtrf.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -66,7 +66,6 @@
static void compute( integer_t const m, integer_t const n,
integer_t const kl, integer_t const ku, MatrixAB& ab,
VectorIPIV& ipiv, integer_t& info ) {
-
#ifndef NDEBUG
assert( m >= 0 );
assert( n >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebak.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebak.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebak.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -110,7 +110,6 @@
static void compute( char const job, char const side, integer_t const ilo,
integer_t const ihi, VectorSCALE& scale, MatrixV& v,
integer_t& info ) {
-
#ifndef NDEBUG
assert( job == 'N' || job == 'P' || job == 'S' || job == 'B' );
assert( side == 'R' || side == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebal.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebal.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gebal.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -98,7 +98,6 @@
template< typename MatrixA, typename VectorSCALE >
static void compute( char const job, MatrixA& a, integer_t& ilo,
integer_t& ihi, VectorSCALE& scale, integer_t& info ) {
-
#ifndef NDEBUG
assert( job == 'N' || job == 'P' || job == 'S' || job == 'B' );
assert( traits::matrix_num_columns(a) >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gecon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gecon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/gecon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -79,7 +79,6 @@
static void compute( char const norm, MatrixA& a, real_type const anorm,
real_type& rcond, integer_t& info, detail::workspace2< WORK,
IWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( traits::matrix_num_columns(a) >= 0 );
@@ -137,7 +136,6 @@
static void compute( char const norm, MatrixA& a, real_type const anorm,
real_type& rcond, integer_t& info, detail::workspace2< WORK,
RWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( traits::matrix_num_columns(a) >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/geqp3.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/geqp3.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/geqp3.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -92,7 +92,7 @@
std::min(traits::matrix_num_rows(a),
traits::matrix_num_columns(a)) );
assert( traits::vector_size(work.select(real_type())) >=
- min_size_work( $CALL_MIN_SIZE ));
+ min_size_work( traits::matrix_num_columns(a) ));
#endif
detail::geqp3( traits::matrix_num_rows(a),
traits::matrix_num_columns(a), traits::matrix_storage(a),
@@ -107,7 +107,7 @@
static void compute( MatrixA& a, VectorJPVT& jpvt, VectorTAU& tau,
integer_t& info, minimal_workspace work ) {
traits::detail::array< real_type > tmp_work( min_size_work(
- $CALL_MIN_SIZE ) );
+ traits::matrix_num_columns(a) ) );
compute( a, jpvt, tau, info, workspace( tmp_work ) );
}
@@ -125,8 +125,8 @@
compute( a, jpvt, tau, info, workspace( tmp_work ) );
}
- static integer_t min_size_work( $ARGUMENTS ) {
- $MIN_SIZE
+ static integer_t min_size_work( integer_t const n ) {
+ return 3*n+1;
}
};
@@ -154,7 +154,7 @@
std::min(traits::matrix_num_rows(a),
traits::matrix_num_columns(a)) );
assert( traits::vector_size(work.select(value_type())) >=
- min_size_work( $CALL_MIN_SIZE ));
+ min_size_work( traits::matrix_num_columns(a) ));
assert( traits::vector_size(work.select(real_type())) >=
min_size_rwork( traits::matrix_num_columns(a) ));
#endif
@@ -172,7 +172,7 @@
static void compute( MatrixA& a, VectorJPVT& jpvt, VectorTAU& tau,
integer_t& info, minimal_workspace work ) {
traits::detail::array< value_type > tmp_work( min_size_work(
- $CALL_MIN_SIZE ) );
+ traits::matrix_num_columns(a) ) );
traits::detail::array< real_type > tmp_rwork( min_size_rwork(
traits::matrix_num_columns(a) ) );
compute( a, jpvt, tau, info, workspace( tmp_work, tmp_rwork ) );
@@ -195,8 +195,8 @@
compute( a, jpvt, tau, info, workspace( tmp_work, tmp_rwork ) );
}
- static integer_t min_size_work( $ARGUMENTS ) {
- $MIN_SIZE
+ static integer_t min_size_work( integer_t const n ) {
+ return n+1;
}
static integer_t min_size_rwork( integer_t const n ) {
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getrf.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getrf.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getrf.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -60,7 +60,6 @@
// templated specialization
template< typename MatrixA, typename VectorIPIV >
static void compute( MatrixA& a, VectorIPIV& ipiv, integer_t& info ) {
-
#ifndef NDEBUG
assert( traits::matrix_num_rows(a) >= 0 );
assert( traits::matrix_num_columns(a) >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getri.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getri.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/getri.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -75,7 +75,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( MatrixA& a, VectorIPIV& ipiv, integer_t& info,
detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( traits::matrix_num_columns(a) >= 0 );
assert( traits::leading_dimension(a) >= std::max(1,
@@ -129,7 +128,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( MatrixA& a, VectorIPIV& ipiv, integer_t& info,
detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( traits::matrix_num_columns(a) >= 0 );
assert( traits::leading_dimension(a) >= std::max(1,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hecon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hecon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hecon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -60,7 +60,6 @@
static void compute( char const uplo, MatrixA& a, VectorIPIV& ipiv,
real_type const anorm, real_type& rcond, integer_t& info,
detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( traits::matrix_num_columns(a) >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetrf.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetrf.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetrf.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -57,7 +57,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( MatrixA& a, VectorIPIV& ipiv, integer_t& info,
detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetri.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetri.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/hetri.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -57,7 +57,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( char const uplo, MatrixA& a, VectorIPIV& ipiv,
integer_t& info, detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( traits::matrix_num_columns(a) >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/lacon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/lacon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/lacon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -70,7 +70,6 @@
template< typename VectorX, typename V, typename ISGN >
static void compute( integer_t const n, VectorX& x, real_type& est,
integer_t& kase, detail::workspace2< V, ISGN > work ) {
-
#ifndef NDEBUG
assert( n >= 1 );
assert( traits::vector_size(work.select(real_type())) >=
@@ -119,7 +118,6 @@
template< typename VectorX, typename V >
static void compute( integer_t const n, VectorX& x, real_type& est,
integer_t& kase, detail::workspace1< V > work ) {
-
#ifndef NDEBUG
assert( n >= 1 );
assert( traits::vector_size(work.select(value_type())) >=
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbcon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbcon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbcon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -80,7 +80,6 @@
integer_t const kd, MatrixAB& ab, real_type const anorm,
real_type& rcond, integer_t& info, detail::workspace2< WORK,
IWORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( n >= 0 );
@@ -138,7 +137,6 @@
integer_t const kd, MatrixAB& ab, real_type const anorm,
real_type& rcond, integer_t& info, detail::workspace2< WORK,
RWORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( n >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbequ.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbequ.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbequ.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -99,7 +99,6 @@
template< typename MatrixAB, typename VectorS >
static void compute( integer_t const n, integer_t const kd, MatrixAB& ab,
VectorS& s, real_type& scond, real_type& amax, integer_t& info ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbstf.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbstf.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbstf.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -59,7 +59,6 @@
template< typename MatrixAB >
static void compute( integer_t const n, integer_t const kd, MatrixAB& ab,
integer_t& info ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbtrf.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbtrf.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pbtrf.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -59,7 +59,6 @@
template< typename MatrixAB >
static void compute( integer_t const n, integer_t const kd, MatrixAB& ab,
integer_t& info ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pocon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pocon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/pocon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -78,7 +78,6 @@
template< typename MatrixA, typename WORK, typename IWORK >
static void compute( MatrixA& a, real_type const anorm, real_type& rcond,
integer_t& info, detail::workspace2< WORK, IWORK > work ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
@@ -135,7 +134,6 @@
template< typename MatrixA, typename WORK, typename RWORK >
static void compute( MatrixA& a, real_type const anorm, real_type& rcond,
integer_t& info, detail::workspace2< WORK, RWORK > work ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/poequ.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/poequ.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/poequ.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -95,7 +95,6 @@
template< typename MatrixA, typename VectorS >
static void compute( MatrixA& a, VectorS& s, real_type& scond,
real_type& amax, integer_t& info ) {
-
#ifndef NDEBUG
assert( traits::matrix_num_columns(a) >= 0 );
assert( traits::leading_dimension(a) >= std::max(1,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potrf.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potrf.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potrf.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -58,7 +58,6 @@
// templated specialization
template< typename MatrixA >
static void compute( MatrixA& a, integer_t& info ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potri.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potri.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/potri.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -58,7 +58,6 @@
// templated specialization
template< typename MatrixA >
static void compute( MatrixA& a, integer_t& info ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sycon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sycon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sycon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -80,7 +80,6 @@
static void compute( char const uplo, MatrixA& a, VectorIPIV& ipiv,
real_type const anorm, real_type& rcond, integer_t& info,
detail::workspace2< WORK, IWORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( traits::matrix_num_columns(a) >= 0 );
@@ -141,7 +140,6 @@
static void compute( char const uplo, MatrixA& a, VectorIPIV& ipiv,
real_type const anorm, real_type& rcond, integer_t& info,
detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( traits::matrix_num_columns(a) >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytrf.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytrf.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytrf.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -75,7 +75,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( MatrixA& a, VectorIPIV& ipiv, integer_t& info,
detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
@@ -131,7 +130,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( MatrixA& a, VectorIPIV& ipiv, integer_t& info,
detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( traits::matrix_uplo_tag(a) == 'U' ||
traits::matrix_uplo_tag(a) == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytri.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytri.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/sytri.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -74,7 +74,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( char const uplo, MatrixA& a, VectorIPIV& ipiv,
integer_t& info, detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( traits::matrix_num_columns(a) >= 0 );
@@ -122,7 +121,6 @@
template< typename MatrixA, typename VectorIPIV, typename WORK >
static void compute( char const uplo, MatrixA& a, VectorIPIV& ipiv,
integer_t& info, detail::workspace1< WORK > work ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( traits::matrix_num_columns(a) >= 0 );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/tbcon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/tbcon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/tbcon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -82,7 +82,6 @@
integer_t const n, integer_t const kd, MatrixAB& ab,
real_type& rcond, integer_t& info, detail::workspace2< WORK,
IWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( uplo == 'U' || uplo == 'L' );
@@ -143,7 +142,6 @@
integer_t const n, integer_t const kd, MatrixAB& ab,
real_type& rcond, integer_t& info, detail::workspace2< WORK,
RWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( uplo == 'U' || uplo == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trcon.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trcon.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trcon.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -79,7 +79,6 @@
static void compute( char const norm, char const uplo, char const diag,
MatrixA& a, real_type& rcond, integer_t& info, detail::workspace2<
WORK, IWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( uplo == 'U' || uplo == 'L' );
@@ -140,7 +139,6 @@
static void compute( char const norm, char const uplo, char const diag,
MatrixA& a, real_type& rcond, integer_t& info, detail::workspace2<
WORK, RWORK > work ) {
-
#ifndef NDEBUG
assert( norm == '1' || norm == 'O' || norm == 'I' );
assert( uplo == 'U' || uplo == 'L' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trtri.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trtri.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/computational/trtri.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -59,7 +59,6 @@
template< typename MatrixA >
static void compute( char const uplo, char const diag, MatrixA& a,
integer_t& info ) {
-
#ifndef NDEBUG
assert( uplo == 'U' || uplo == 'L' );
assert( diag == 'N' || diag == 'U' );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heev.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heev.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heev.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -61,7 +61,6 @@
typename RWORK >
static void compute( char const jobz, MatrixA& a, VectorW& w,
integer_t& info, detail::workspace2< WORK, RWORK > work ) {
-
#ifndef NDEBUG
assert( jobz == 'N' || jobz == 'V' );
assert( traits::matrix_uplo_tag(a) == 'U' ||
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heevd.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heevd.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/heevd.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -65,7 +65,6 @@
typename RWORK, typename IWORK >
static void compute( char const jobz, MatrixA& a, VectorW& w,
integer_t& info, detail::workspace3< WORK, RWORK, IWORK > work ) {
-
#ifndef NDEBUG
assert( jobz == 'N' || jobz == 'V' );
assert( traits::matrix_uplo_tag(a) == 'U' ||
Modified: sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/lacgv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/lacgv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/lacgv.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -51,7 +51,6 @@
template< typename VectorX >
static void compute( integer_t const n, VectorX& x,
integer_t const incx ) {
-
#ifndef NDEBUG
assert( n >= 0 );
#endif
Added: sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/computational/geqp3.hpp
==============================================================================
--- (empty file)
+++ sandbox/numeric_bindings/libs/numeric/bindings/tools/templates/computational/geqp3.hpp 2009-03-09 09:41:13 EDT (Mon, 09 Mar 2009)
@@ -0,0 +1,7 @@
+$TEMPLATE[geqp3.all.min_size_work.args]
+N
+$TEMPLATE[geqp3.real.min_size_work]
+return 3*n+1;
+$TEMPLATE[geqp3.complex.min_size_work]
+return n+1;
+$TEMPLATE[end]
Boost-Commit list run by bdawes at acm.org, david.abrahams at rcn.com, gregod at cs.rpi.edu, cpdaniel at pacbell.net, john at johnmaddock.co.uk