BoostMatrix.h

Go to the documentation of this file.

//
// This header file is made to facilitate boost's matrix and vector_t operations in OPAL.
//
// Copyright (c) 2023, Paul Scherrer Institute, Villigen PSI, Switzerland
// All rights reserved
//
// This file is part of OPAL.
//
// OPAL is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
 
#ifndef OPAL_MATRIX_HH
#define OPAL_MATRIX_HH
 
#include <boost/numeric/ublas/matrix.hpp>
 
typedef boost::numeric::ublas::matrix<double> matrix_t;
 
template <class T>
T prod_boost_vector(boost::numeric::ublas::matrix<double> rotation, const T& vect) {
    boost::numeric::ublas::vector<double> boostVector(3);
 
    boostVector(0) = vect[0];
    boostVector(1) = vect[1];
    boostVector(2) = vect[2];
 
    boostVector = boost::numeric::ublas::prod(rotation, boostVector);
 
    T prodVector(0.0);
 
    prodVector[0] = boostVector(0);
    prodVector[1] = boostVector(1);
    prodVector[2] = boostVector(2);
 
    return prodVector;
}
 
#endif