matrix_banded.hpp

#ifndef MATRIX_BANDED_H
#define MATRIX_BANDED_H
#include <memory>
 
#include "sll/matrix.hpp"
 
class Matrix_Banded : public Matrix
{
public:
    Matrix_Banded(int n, int kl, int ku);
    virtual double get_element(int i, int j) const override;
    virtual void set_element(int i, int j, double a_ij) override;
 
protected:
    virtual int factorize_method() override;
    virtual int solve_inplace_method(double* b, char transpose, int n_equations) const override;
    int const kl; // no. of subdiagonals
    int const ku; // no. of superdiagonals
    int const c; // no. of columns in q
    std::unique_ptr<int[]> ipiv; // pivot indices
    std::unique_ptr<double[]> q; // banded matrix representation
};
 
#endif // MATRIX_BANDED_H