multiply.hpp

Go to the documentation of this file.

#ifndef __STAN__MATH__MATRIX__MULTIPLY_HPP__
#define __STAN__MATH__MATRIX__MULTIPLY_HPP__

#include <stdexcept>
#include <stan/math/matrix/Eigen.hpp>
#include <stan/math/matrix/validate_multiplicable.hpp>
#include <stan/math/matrix/validate_matching_sizes.hpp>

namespace stan {
  namespace math {

    template <int R, int C>
    inline
    Eigen::Matrix<double,R,C>
    multiply(const Eigen::Matrix<double,R,C>& m,
             double c) {
      return c * m;
    }
    template <int R, int C>
    inline
    Eigen::Matrix<double,R,C>
    multiply(double c,
             const Eigen::Matrix<double,R,C>& m) {
      return c * m;
    }


    template<int R1,int C1,int R2,int C2>
    inline Eigen::Matrix<double,R1,C2> multiply(const Eigen::Matrix<double,R1,C1>& m1,
                                                const Eigen::Matrix<double,R2,C2>& m2) {
      
      validate_multiplicable(m1,m2,"multiply");
      return m1*m2;
    }

    template<int C1,int R2>
    inline double multiply(const Eigen::Matrix<double,1,C1>& rv,
                           const Eigen::Matrix<double,R2,1>& v) {
      stan::math::validate_matching_sizes(rv,v,"multiply");
      if (rv.size() != v.size()) 
        throw std::domain_error("rv.size() != v.size()");
      return rv.dot(v);
    }

  }
}
#endif