lil_sparse_matrix.h

// -*- mode: c++ -*-
#ifndef _LIL_SPARSE_MATRIX_H_
#define _LIL_SPARSE_MATRIX_H_

#include <vector>
#include <string>
#include <fstream>
#include <limits>


using std::vector;

template<class el_type>
class lil_sparse_matrix 
{

public:

    typedef vector<int> idx_vector_type;
    typedef vector<el_type>  elt_vector_type;
    
    friend std::ostream & operator<<(std::ostream& os, const lil_sparse_matrix& A) 
    {
        os << A.to_string();
        return os;
    };

    int m_n_rows;
    int m_n_cols;
    int nnz_count;
    el_type eps;

    vector<idx_vector_type> m_idx;
    vector<elt_vector_type> m_x;

    lil_sparse_matrix (int n_rows, int n_cols) : m_n_rows(n_rows), m_n_cols (n_cols)
    {
        nnz_count = 0;
        eps = 1e-8;
    }
    
    int n_rows() const
    {
        return m_n_rows;
    }

    int n_cols() const
    {
        return m_n_cols;
    }

    int nnz() const 
    {
        return nnz_count;
    };
    
    virtual el_type coeff(const int& i, const int& j, int offset = 0) const = 0;
    
    virtual std::string to_string() const = 0;

    virtual ~lil_sparse_matrix()
    {
    }
};

#endif // _LIL_SPARSE_MATRIX_H_