main.h

#ifndef MAIN_H
#define MAIN_H

#include <iostream>
#include <vector>
#include <Eigen/Dense>
#include <Eigen/StdVector>
#include <unsupported/Eigen/KroneckerProduct>
#include "GlobalHamiltonianParameters.h"

#define kp kroneckerProduct

#ifdef realHamiltonian
    typedef double hamScalar;
#elif defined(complexHamiltonian)
    typedef std::complex<double> hamScalar;
#endif
#if defined(realHamiltonian) && defined(realObservables)
    #define obsRe
    typedef double obsScalar;
#elif defined(complexObservables) || defined(complexHamiltonian)
    #define obsRe std::real
    typedef std::complex<double> obsScalar;
#endif
#if defined(realHamiltonian) && defined(complexObservables)
    #define differentScalars
#endif

typedef Eigen::Matrix<hamScalar, d, d> MatrixD_t;
typedef Eigen::Matrix<hamScalar, Eigen::Dynamic, Eigen::Dynamic> MatrixX_t;
typedef Eigen::Matrix<hamScalar, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor>
    rmMatrixX_t;
typedef Eigen::Matrix<hamScalar, Eigen::Dynamic, 1> VectorX_t;
typedef Eigen::Matrix<obsScalar, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor>
    obsMatrixX_t;
typedef Eigen::Matrix<obsScalar, d, d> obsMatrixD_t;
typedef std::vector<std::pair<obsMatrixD_t, int>,
                    Eigen::aligned_allocator<std::pair<obsMatrixD_t, int>>>
                    opsVec;

#endif