mdivide_left_tri_low.hpp

Go to the documentation of this file.

#ifndef __STAN__MATH__MATRIX__MDIVIDE_LEFT_TRI_LOW_HPP__
#define __STAN__MATH__MATRIX__MDIVIDE_LEFT_TRI_LOW_HPP__

#include <boost/math/tools/promotion.hpp>
#include <stan/math/matrix/Eigen.hpp>
#include <stan/math/matrix/mdivide_left_tri.hpp>

namespace stan {
  namespace math {

    template <typename T1, typename T2, int R1,int C1,int R2,int C2>
    inline 
    Eigen::Matrix<typename boost::math::tools::promote_args<T1,T2>::type,
    R1,C2>
    mdivide_left_tri_low(const Eigen::Matrix<T1,R1,C1> &A,
                         const Eigen::Matrix<T2,R2,C2> &b) {
//      stan::math::validate_square(A,"mdivide_left_tri_low/2");
//      stan::math::validate_multiplicable(A,b,"mdivide_left_tri_low");
//      return promote_common<Eigen::Matrix<T1,R1,C1>,
//      Eigen::Matrix<T2,R1,C1> >(A)
//      .template triangularView<Eigen::Lower>()
//      .solve( promote_common<Eigen::Matrix<T1,R2,C2>,
//             Eigen::Matrix<T2,R2,C2> >(b) );
      return mdivide_left_tri<Eigen::Lower,T1,T2,R1,C1,R2,C2>(A,b);
    }
    template <typename T,int R1, int C1>
    inline 
    Eigen::Matrix<T,R1,C1>
    mdivide_left_tri_low(const 
                         Eigen::Matrix<T,R1,C1> &A) {
//      stan::math::validate_square(A,"mdivide_left_tri_low/1");
//      int n = A.rows();
//      Eigen::Matrix<T,Eigen::Dynamic,Eigen::Dynamic> b;
//      b.setIdentity(n,n);
//      A.template triangularView<Eigen::Lower>().solveInPlace(b);
//      return b;
      return mdivide_left_tri<Eigen::Lower,T,R1,C1>(A);
    }

  }
}
#endif