#include <math.h> 
#include <time.h> 
#include <sstream> 
#include <ostream> 
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/io.hpp> 
#include <boost/numeric/ublas/operation.hpp>
 
using namespace std;
using namespace boost::numeric::ublas;

int main() {

  int M= 100;
  compressed_matrix<double> A(M,M); 
  A(0,0)=-2; 
  A(0,1)=1; 
  A(M-1,M-1)=-2; 
  A(0,M-1)=1; 
  A(M-1,0)=1; 
  A(M-1,M-2)=1; 
     
  for (int p=1;p<M-1;p++){ 
        A(p,p-1)=1; 
        A(p,p)=-2; 
        A(p,p+1)=1; 
        } 
  cout <<"after A matrix generation..\n";    
  /*------------generation of the H matrix------------------*/ 

  clock_t start,end;
  double tempo;
  start=clock();
  compressed_matrix<double> H(M,M);
  axpy_prod(A,A,H,false);
  end=clock();
  tempo=((long double)(end-start))/CLOCKS_PER_SEC;      
  cout << "Sparse_matrix multiplication done. time: "<< tempo << " secs.\n";
  A.~compressed_matrix();
  H.~compressed_matrix();
  
}