sample68.cpp

Go to the documentation of this file.

#include <iostream>

// #include <boost/numeric/ublas/config.hpp>
// #undef BOOST_UBLAS_STRICT_VECTOR_SPARSE
// #undef BOOST_UBLAS_STRICT_MATRIX_SPARSE
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/io.hpp>

using std::cout;
using std::endl;

using namespace boost::numeric::ublas;

int main(int argc, char *argv[])
{

  // typedef sparse_matrix<double>  MY_STIMA;
  typedef compressed_matrix<double>  MY_STIMA;

  MY_STIMA M(4,4,16);

  for (int i=0; i<M.size1(); ++i) M(i,i)=1.0;
  M(0,0)=3;
  M(0,1)=-1;
  M(0,2)=-2;
  M(1,0)=-1;
  M(1,1)=4;
  M(1,2)=-2;
  M(2,0)=-2;
  M(2,1)=-2;
  M(2,2)=9;

  cout << "input M = " << M << "\n";

  MY_STIMA::iterator1 r_it1 = M.begin1();
  MY_STIMA::iterator2 c_it1 = M.begin2();
  while ( (r_it1 != M.end1()) && (c_it1 != M.end2()) ) {
#ifndef BOOST_UBLAS_NO_NESTED_CLASS_RELATION
        MY_STIMA::iterator1  r_it2 = c_it1.begin();
        MY_STIMA::iterator2  c_it2 = r_it1.begin();
#else
        MY_STIMA::iterator1  r_it2 = begin(c_it1, iterator2_tag());
        MY_STIMA::iterator2  c_it2 = begin(r_it1, iterator1_tag());
#endif
        MY_STIMA::iterator2  c;

        c = c_it2;
#ifndef BOOST_UBLAS_NO_NESTED_CLASS_RELATION
        while ( (r_it2 != c_it1.end()) && (c != r_it1.end()) ) {
#else
        while ( (r_it2 != end(c_it1, iterator2_tag())) && (c != end(r_it1, iterator1_tag())) ) {
#endif
          while (c.index2() < r_it2.index1()) {
#ifndef BOOST_UBLAS_NO_NESTED_CLASS_RELATION
              assert (c != r_it1.end());
#else
              assert (c != end(r_it1, iterator1_tag()));
#endif
              ++c;
          }
#ifndef BOOST_UBLAS_NO_NESTED_CLASS_RELATION
          assert (c != r_it1.end());
#else
          assert (c != end(r_it1, iterator1_tag()));
#endif
          while (c.index2() > r_it2.index1()) {
#ifndef BOOST_UBLAS_NO_NESTED_CLASS_RELATION
              assert (r_it2 != c_it1.end());
#else
              assert (r_it2 != end(c_it1, iterator2_tag()));
#endif
              ++r_it2;
          }
#ifndef BOOST_UBLAS_NO_NESTED_CLASS_RELATION
          assert (r_it2 != c_it1.end());
#else
          assert (r_it2 != end(c_it1, iterator2_tag()));
#endif
          if (c.index2() == r_it2.index1()) {
                cout << r_it2.index1() << ", " << c.index2() << ": " << (*c) << " += " << (*r_it2);
                double temp = (*r_it2);
                (*c) += temp;
                cout << " --> " << (*c) << " ? \n";
                ++ r_it2;
                ++ c;
          }
        }
        ++ r_it1;
        ++ c_it1;
  }

  cout << "output M = " << M << endl;

  return EXIT_SUCCESS;
}

---