promote_scalar.hpp

Go to the documentation of this file.

#ifndef STAN__MATH__MATRIX__PROMOTE_SCALAR_HPP
#define STAN__MATH__MATRIX__PROMOTE_SCALAR_HPP

#include <stan/math/functions/promote_scalar.hpp>
#include <stan/math/matrix/promote_scalar_type.hpp>
#include <stan/math/matrix/Eigen.hpp>

namespace stan {

  namespace math {

    template <typename T, typename S>
    struct promote_scalar_struct<T, Eigen::Matrix<S,-1,-1> > {

      static Eigen::Matrix<typename promote_scalar_type<T,S>::type, -1,-1>
      apply(const Eigen::Matrix<S, -1,-1>& x) {
        Eigen::Matrix<typename promote_scalar_type<T,S>::type, -1,-1>
          y(x.rows(), x.cols());
        for (size_t i = 0; i < x.size(); ++i)
          y(i) = promote_scalar_struct<T,S>::apply(x(i));
        return y;
      }

    };


    template <typename T, typename S>
    struct promote_scalar_struct<T, Eigen::Matrix<S,1,-1> > {

      static Eigen::Matrix<typename promote_scalar_type<T,S>::type, 1,-1>
      apply(const Eigen::Matrix<S, 1,-1>& x) {
        Eigen::Matrix<typename promote_scalar_type<T,S>::type, 1,-1>
          y(x.rows(), x.cols());
        for (size_t i = 0; i < x.size(); ++i)
          y(i) = promote_scalar_struct<T,S>::apply(x(i));
        return y;
      }

    };


    template <typename T, typename S>
    struct promote_scalar_struct<T, Eigen::Matrix<S,-1,1> > {

      static Eigen::Matrix<typename promote_scalar_type<T,S>::type, -1,1>
      apply(const Eigen::Matrix<S, -1,1>& x) {
        Eigen::Matrix<typename promote_scalar_type<T,S>::type, -1,1>
          y(x.rows(), x.cols());
        for (size_t i = 0; i < x.size(); ++i)
          y(i) = promote_scalar_struct<T,S>::apply(x(i));
        return y;
      }

    };


  }
}


#endif