linear_solver.h

#pragma once
#include "sparse_block_matrix.h"
#include "sparse_block_matrix_ccs.h"
 
namespace NDEVR {


    template <typename MatrixType>
    class LinearSolverCCS
    {
    public:
        LinearSolverCCS()
            : _ccsMatrix(0)
        {}
 
    protected:
        SparseBlockMatrixCCS<MatrixType> _ccsMatrix;    

        void initMatrixStructure(const SparseBlockMatrix<MatrixType>& A)
        {
            _ccsMatrix = SparseBlockMatrixCCS<MatrixType>(A.rowBlockIndices(), A.colBlockIndices());
            A.fillSparseBlockMatrixCCS(*_ccsMatrix);
        }
        EIGEN_MAKE_ALIGNED_OPERATOR_NEW
    };
 
}