Esempio n. 1
0
 inline bool check_greater_or_equal(const char* function,
                                    const T_y& y,
                                    const T_low& low,
                                    const char* name,  
                                    T_result* result) {
   return check_greater_or_equal(function,y,low,name,result,
                                 default_policy());
 }                               
Esempio n. 2
0
 Eigen::Matrix<fvar<T>, Eigen::Dynamic, Eigen::Dynamic>
 qr_Q(const Eigen::Matrix<fvar<T>, Eigen::Dynamic, Eigen::Dynamic>& m) {
   typedef Eigen::Matrix<fvar<T>, Eigen::Dynamic, Eigen::Dynamic>
     matrix_fwd_t;
   check_nonzero_size("qr_Q", "m", m);
   check_greater_or_equal("qr_Q", "m.rows()", m.rows(), m.cols());
   Eigen::HouseholderQR< matrix_fwd_t > qr(m.rows(), m.cols());
   qr.compute(m);
   matrix_fwd_t Q = qr.householderQ();
   for (int i = 0; i < m.cols(); i++)
     if (qr.matrixQR()(i, i) < 0.0)
       Q.col(i) *= -1.0;
   return Q;
 }