RKMatrixEigenTransformations.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RKMatrixEigenTransformations.h

#pragma once

#include "EigenMatrixTypedefs.h"
#include "RKTools.h"

namespace genfit {

    template <unsigned int rows, unsigned int cols>
    Eigen::Matrix<double, rows, cols> RKMatrixToEigenMatrix(const RKMatrix<rows, cols>& rkMatrix) {
        Eigen::Matrix<double, rows, cols> eigenMatrix;

        for (unsigned int row=0; row < rows; ++row) {
            for (unsigned int col=0; col < cols; ++col) {
                eigenMatrix(row, col) = rkMatrix[cols*row + col];
            }
        }

        return eigenMatrix;
    }


    template <unsigned int rows, unsigned int cols>
    RKMatrix<rows, cols> eigenMatrixToRKMatrix(const Eigen::Matrix<double, rows, cols>& eigenMatrix) {
        RKMatrix<rows, cols> rkMatrix;

        for (unsigned int row=0; row < rows; ++row) {
            for (unsigned int col=0; col < cols; ++col) {
                rkMatrix(row, col) = eigenMatrix(row, col);
            }
        }

        return rkMatrix;
    }

}