
/* Test program to test find functions of triagular matrices
 *
 * author: Gunter Winkler ( guwi17 at gmx dot de )
 */


// ublas headers

#include <boost/numeric/ublas/experimental/sparse_view.hpp>

#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/io.hpp>

#include <boost/numeric/ublas/traits/c_array.hpp>

// other boost headers

// headers for testcase

#define BOOST_TEST_MODULE SparseMatrixErasureTest
#include <boost/test/included/unit_test.hpp>

// standard and system headers

#include <iostream>
#include <string>

namespace ublas = boost::numeric::ublas;

    /*
      sparse input matrix:

      1 2 0 0
      0 3 9 0
      0 1 4 0
    */

    static const std::string inputMatrix = "[3,4]((1,2,0,0),(0,3,9,0),(0,1,4,0))\n";

    const unsigned int NNZ  = 6;
    const unsigned int IB   = 1;
    const double VA[]       = { 1.0, 2.0, 3.0, 9.0, 1.0, 4.0 };
    const unsigned int IA[] = { 1, 3, 5, 7 };
    const unsigned int JA[] = { 1, 2, 2, 3, 2, 3 };

BOOST_AUTO_TEST_CASE( test_construction_and_basic_operations )
{

    typedef ublas::matrix<double> DENSE_MATRIX;
    
    // prepare data

    DENSE_MATRIX A;

    std::istringstream iss(inputMatrix);
    iss >> A;

    std::cout << A << std::endl;

    std::cout << ( ublas::make_compressed_matrix_view<ublas::row_major,IB>(3,4,NNZ,IA,JA,VA) ) << std::endl;

    typedef ublas::compressed_matrix_view<ublas::row_major, IB, unsigned int [4], unsigned int [NNZ], double[NNZ]> COMPMATVIEW;

    COMPMATVIEW viewA(3,4,NNZ,IA,JA,VA);

    std::cout << viewA << std::endl;

}


BOOST_AUTO_TEST_CASE( test_construction_from_pointers )
{

    std::cout << ( ublas::make_compressed_matrix_view<ublas::column_major,IB>(4,3,NNZ
                                                                              , ublas::c_array_view<const unsigned int>(4,&(IA[0]))
                                                                              , ublas::c_array_view<const unsigned int>(6,&(JA[0]))
                                                                              , ublas::c_array_view<const double>(6,&(VA[0]))) ) << std::endl;

    unsigned int * ia = new unsigned int[4]();
    unsigned int * ja = new unsigned int[6]();
    double * va = new double[6]();

    std::copy(&(IA[0]),&(IA[4]),ia);
    std::copy(&(JA[0]),&(JA[6]),ja);
    std::copy(&(VA[0]),&(VA[6]),va);

    typedef ublas::compressed_matrix_view<ublas::column_major
      , IB
      , ublas::c_array_view<unsigned int>
      , ublas::c_array_view<unsigned int>
      , ublas::c_array_view<double> > COMPMATVIEW;

    COMPMATVIEW viewA(4,3,NNZ
                      , ublas::c_array_view<unsigned int>(4,ia)
                      , ublas::c_array_view<unsigned int>(6,ja)
                      , ublas::c_array_view<double>(6,va));

    std::cout << viewA << std::endl;

    delete[] va;
    delete[] ja;
    delete[] ia;

}