#ifndef PETSCMATRIX_HH
#define PETSCMATRIX_HH

#include <stdexcept>

#include "petscksp.h"

class petscMatrix
{
    public:
        petscMatrix();
        petscMatrix(PetscInt, PetscInt, char) throw(std::runtime_error);
        Mat* assignSubMatrix();
        int createMatrix();
    private:
        Mat A;        // Operation matrix
        PetscInt m,n; // the row and column count 
                      // of the matrix
        char matType;
        // utility function to create the matrix

};
#endif