// misc. boost includes
#include <boost/timer.hpp>

#include "trisolve.h"

// boost ublas includes
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/matrix_proxy.hpp>
#include <boost/numeric/ublas/triangular.hpp>
#include <boost/numeric/ublas/io.hpp>

// system includes
#include <iostream>
#include <limits>

using namespace boost::numeric::ublas;

int main ( int argc, char* argv[])
{
    int iterations = atoi (argv[1]);
    int gatesperpath = atoi (argv[2]);

        boost::timer t;
    {
      t.restart();
      unsigned int i = 0;
      double j;
      identity_matrix<double> ident (gatesperpath);
      compressed_matrix<double>   lembda1 (gatesperpath, gatesperpath);
      lembda1 (3, 1) = -0.1;
      vector <double> omegas(gatesperpath);
      vector<double> slopes;
    
      int iter = iterations;
      while (iter--) {
        // (I+lembda) slopes = omegas
        slopes = solve (lembda1, omegas, unit_lower_tag ()); 
      }
      std::cout << t.elapsed () << std::endl;
    }

    {
      t.restart();
      unsigned int i = 0;
      double j;
      identity_matrix<double> ident (gatesperpath);
      compressed_matrix<double>   lembda1 (gatesperpath, gatesperpath);
      lembda1 (3, 1) = -0.1;
      vector <double> omegas(gatesperpath);
      vector<double> slopes;
      
      int iter = iterations;
      while (iter--) {
        // (I+lembda) slopes = omegas
        slopes = omegas;
        inplace_solve_lower<true>(lembda1, slopes, row_major_tag() );
      }
      std::cout << t.elapsed () << std::endl;
    }

    return 0;
}