Skip to content
This repository has been archived by the owner on Oct 19, 2020. It is now read-only.

Commit

Permalink
Browse files Browse the repository at this point in the history
Fix asserts in gbsv.hpp as suggested by Vardan.
  • Loading branch information
inducer committed Nov 16, 2008
1 parent bdda751 commit fbbc7c9
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions boost/numeric/bindings/lapack/gbsv.hpp
Expand Up @@ -54,13 +54,13 @@ namespace boost { namespace numeric { namespace bindings {
int const m = traits::matrix_size2 (a);
assert (traits::vector_size (ipiv) == (m < n ? m : n));

// if the matrix has kl lower and ku upper diagonals, then we should have
// if the matrix has kl lower and ku upper diagonals, then we should have
// allocated kl lower and kl+ku upper diagonals
int const kl = traits::matrix_lower_bandwidth (a);
int const ku = traits::matrix_upper_bandwidth (a) - kl;
int const ld = traits::leading_dimension (a);

assert(ku > 0);
assert(ku >= 0);

int info;
detail::gbtrf (n, m, kl, ku,
Expand Down Expand Up @@ -101,13 +101,13 @@ namespace boost { namespace numeric { namespace bindings {
assert (n == traits::matrix_size1 (b));
assert (n == traits::vector_size (ipiv));

// if the matrix has kl lower and ku upper diagonals, then we should have
// if the matrix has kl lower and ku upper diagonals, then we should have
// allocated kl lower and kl+ku upper diagonals
int const kl = traits::matrix_lower_bandwidth (a);
int const ku = traits::matrix_upper_bandwidth (a) - kl;
int const ld = traits::leading_dimension (a);

assert(ku > 0);
assert(ku >= 0);

int info;
detail::gbtrs (trans, n, kl, ku, traits::matrix_size2 (b),
Expand Down

0 comments on commit fbbc7c9

Please sign in to comment.