group.h

#include <iostream>
#include <typeinfo>
#include <complex>
#include <cmath>
#include <Eigen/Eigenvalues>
#include <unsupported/Eigen/CXX11/Tensor>

using namespace Eigen;
using namespace std;

class Group
{
    public:
        MatrixXd group;
        int order;
    // void cyclicGroup(int dim); //Deprecated -- Use CyclicGroup class
    bool checkIfGroup();
        MatrixXd regularRepresentationOfElement(int ele);
    void regularRepresentation();
    void printGroup();
    void correctedCyclicGroup();
    Tensor<complex<double>, 3> pauliMatrices();
        Matrix<complex<double>,2,2> pauliMatrix(int n);
        void printPauliMatrix(int n);
    template <typename T>
    void checkHermiticity(T mat);
};

class PauliGroup:public Group
{
  public:
    PauliGroup();
    // vector <Matrix<complex<double>,2,2> > PauliElements; <-- Initialization of Elements as arrays
    Tensor <complex<double>,3> PauliTensor;
    Matrix <complex<double>,2,2> PauliMatrix(int subScript);

    void printPauliMatrix(int subScript);
};

class CyclicGroup:public Group
{
  public:
    //Cyclic Group constructor with dimensions
     CyclicGroup(int dim);
         MatrixXd HDecomp(int ele);
    //Inherits other group public methods
};

#include "group.tpp"