コード例 #1
0
ファイル: tail.hpp プロジェクト: stan-dev/math
 inline
 Eigen::Matrix<T, 1, Eigen::Dynamic>
 tail(const Eigen::Matrix<T, 1, Eigen::Dynamic>& rv,
      size_t n) {
   if (n != 0)
     check_column_index("tail", "n", rv, n);
   return rv.tail(n);
 }
コード例 #2
0
ファイル: tail.hpp プロジェクト: actuariat/stan
 inline
 Eigen::Matrix<T,1,Eigen::Dynamic>
 tail(const Eigen::Matrix<T,1,Eigen::Dynamic>& rv,
      size_t n) {
   if (n != 0)
     stan::math::check_column_index("tail(%1%)",n,rv,"n",(double*)0);
   return rv.tail(n);
 }
コード例 #3
0
ファイル: tail.hpp プロジェクト: stan-dev/math
 inline
 Eigen::Matrix<T, Eigen::Dynamic, 1>
 tail(const Eigen::Matrix<T, Eigen::Dynamic, 1>& v,
      size_t n) {
   if (n != 0)
     check_row_index("tail", "n", v, n);
   return v.tail(n);
 }
コード例 #4
0
ファイル: tail.hpp プロジェクト: danstowell/stan
 inline
 Eigen::Matrix<T,1,Eigen::Dynamic>
 tail(const Eigen::Matrix<T,1,Eigen::Dynamic>& rv,
      size_t n) {
   if (n != 0)
     validate_column_index(rv, rv.size() - n + 1, "tail");
   return rv.tail(n);
 }
コード例 #5
0
ファイル: tail.hpp プロジェクト: actuariat/stan
 inline
 Eigen::Matrix<T,Eigen::Dynamic,1>
 tail(const Eigen::Matrix<T,Eigen::Dynamic,1>& v,
      size_t n) {
   if (n != 0)
     stan::math::check_row_index("tail(%1%)",n,v,"n",(double*)0);
   return v.tail(n);
 }
コード例 #6
0
ファイル: tail.hpp プロジェクト: danstowell/stan
 inline
 Eigen::Matrix<T,Eigen::Dynamic,1>
 tail(const Eigen::Matrix<T,Eigen::Dynamic,1>& v,
      size_t n) {
   if (n != 0)
     validate_row_index(v, v.size() - n + 1, "tail");
   return v.tail(n);
 }
コード例 #7
0
ファイル: Polynomial.hpp プロジェクト: caomw/polynomial
 void rootBounds( double &lb, double &ub )
 {
     Eigen::Matrix<double,deg,1> mycoef = coef.tail(deg).array().abs();
     mycoef /= fabs(coef(0));
     mycoef(0) += 1.;
     ub = mycoef.maxCoeff();
     lb = -ub;
 }
コード例 #8
0
ファイル: cholesky.hpp プロジェクト: caomw/slam-4
void cholesky_update (Eigen::Matrix<double, N, N>& L, Eigen::Matrix<double, N, 1> v) {

	Eigen::JacobiRotation<double> rot;

	for (int i = 0; i < N; ++i) {
		rot.makeGivens(L(i,i), -v(i), &L(i,i)), v(i) = 0;
		if (i < N-1) apply_jacobi_rotation(L.col(i).tail(N-i-1), v.tail(N-i-1), rot);
	}
}
コード例 #9
0
ファイル: Polynomial.hpp プロジェクト: caomw/polynomial
 void realRoots(std::vector<double> &roots) const
 {
     if ( coef[0] == 0 )
     {
         Internal::RootFinder<deg-1>::compute(coef.tail(deg),roots);
     } else {
         Internal::RootFinder<deg>::compute(coef,roots);
     }
 }
コード例 #10
0
ファイル: Polynomial.hpp プロジェクト: caomw/polynomial
 void realRootsSturm(const double lb, const double ub, std::vector<double> &roots) const
 {
     if ( coef[0] == 0 )
     {
         Internal::SturmRootFinder<deg-1> sturm( coef.tail(deg) );
         sturm.realRoots( lb, ub, roots );
     } else {
         Internal::SturmRootFinder<deg> sturm( coef );
         sturm.realRoots( lb, ub, roots );
     }
 }