SourceForge has been redesigned. Learn more.

matrix size - blas2.cpp

  • Vish108

    Vish108 - 2013-10-04

    I changed the sizes in blas2.cpp provided in the tutorial. Please see below:

    ublas::vector<ScalarType> rhs(100000);
      for (unsigned int i = 0; i < rhs.size(); ++i)
        rhs(i) = random<ScalarType>();
      ublas::vector<ScalarType> rhs2 = rhs;
      ublas::vector<ScalarType> result = ublas::zero_vector<ScalarType>(60000);

    And I got the following error:

    Check failed in file /usr/include/boost/numeric/ublas/functional.hpp at line 1362:
    size_j == 0 || size_i <= (std::numeric_limits<size_type>::max) () / size_j
    terminate called after throwing an instance of 'boost::numeric::ublas::bad_size'
      what():  bad size
    Aborted (core dumped)

    Why did this not pass the numeric limits test? Thank you.

  • Karl Rupp

    Karl Rupp - 2013-10-10


    sorry for the delay, I still haven't figured out how to get notifications for new threads on this forum :-(

    Later in the blas2-example there is a matrix created from the sizes provided to rhs and result. If you increase the size, then the matrix size blows up and hence the error.

    Hope that helps :-)

    Best regards,


Log in to post a comment.