MapAnalyser.h

Go to the documentation of this file.

#ifndef MAP_ANALYSER_H
#define MAP_ANALYSER_H

#include "Classic/FixedAlgebra/FMatrix.h"

#include <array>
#include <complex>

#include <Ippl.h>

class MapAnalyser
{
public:
    typedef FMatrix<double, 6, 6> fMatrix_t;
    typedef FMatrix<std::complex<double>, 6, 6> cfMatrix_t;

    MapAnalyser();

    void linTAnalyze(const fMatrix_t& tMatrix);


    void linSigAnalyze(fMatrix_t& sigMatrix);
private:

    void eigenDecomp_m(const fMatrix_t& M, cfMatrix_t& eigenVal, cfMatrix_t& eigenVec, cfMatrix_t& invEigenVec);

    cfMatrix_t getBlockDiagonal_m(const fMatrix_t& M, cfMatrix_t& eigenVecM, cfMatrix_t& invEigenVecM);


    void printPhaseShift_m(fMatrix_t& Sigma, fMatrix_t tM, cfMatrix_t& oldN);



    void setNMatrix_m(fMatrix_t& M, cfMatrix_t& N, cfMatrix_t& invN);




    fMatrix_t createRotMatrix_m(std::array<double, 3> phi);

    fMatrix_t createSkewMatrix_m();

    fMatrix_t realPartOfMatrix_m(cfMatrix_t cM);

    fMatrix_t imagPartOfMatrix_m(cfMatrix_t cM);

    cfMatrix_t complexTypeMatrix_m(fMatrix_t M);

    cfMatrix_t invertMatrix_m(const cfMatrix_t &M);

    void rearrangeEigen_m(cfMatrix_t& EigenVal, cfMatrix_t& EigenVec);

private:
    IpplTimings::TimerRef mapAnalysis_m;
    IpplTimings::TimerRef bunchAnalysis_m;
};

#endif