Transformation.hpp

Go to the documentation of this file.

#ifndef HOPS_TRANSFORMATION_HPP
#define HOPS_TRANSFORMATION_HPP
 
#include <utility>
 
namespace hops {
 
    template<typename MatrixType, typename VectorType>
    class Transformation {
    public:
        Transformation() = default;
 
        Transformation(MatrixType matrix, VectorType shift) : matrix(matrix), shift(shift) {}
 
        VectorType apply(VectorType vector) {
            return matrix * vector + shift;
        }
 
        VectorType applyInverse(VectorType vector) {
            if (!matrix.isLowerTriangular()) {
                return matrix.inverse() * (vector - shift);
            }
            return matrix.template triangularView<Eigen::Lower>().solve(vector - shift);
        }
 
        MatrixType matrix;
        VectorType shift;
    };
}
 
#endif //HOPS_TRANSFORMATION_HPP