当前位置: 首页>>代码示例>>C++>>正文


C++ AbstractMatrix::rows方法代码示例

本文整理汇总了C++中AbstractMatrix::rows方法的典型用法代码示例。如果您正苦于以下问题:C++ AbstractMatrix::rows方法的具体用法?C++ AbstractMatrix::rows怎么用?C++ AbstractMatrix::rows使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在AbstractMatrix的用法示例。


在下文中一共展示了AbstractMatrix::rows方法的3个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。

示例1: COException

typename Cholesky<Matrix>::DMatrix Cholesky<Matrix>::doSolve( const AbstractMatrix& b )
{
	if ( this->rowNum() != b.rows() )
	{
		std::cerr<<"The order of matrix is "<<this->rowNum()<<" and the dimension of vector is "<<b.rows()<<std::endl;
		throw COException("Cholesky solving error: the size if not consistent!");
	}
	DMatrix result(b);
	copt_lapack_potrs('U',this->rowNum(),b.cols(),__a,this->lda(),result.dataPtr(),result.lda(),&__info);
	if ( __info != 0 )
		std::cerr<<"Warning in Cholesky solver: solving is wrong!"<<std::endl;
	return result;
}
开发者ID:Evanzhouwen,项目名称:copt,代码行数:13,代码来源:CholeskySolver.hpp

示例2: domain_error

void QRDecomposition<T>::operator()(const AbstractMatrix<T>& aMatrix)
{
  q = GenericMatrix<T>(aMatrix.rows());
  r = TriangularMatrix<T>(aMatrix.rows(), TRIANGLE_TYPE::UPPER);
  T normResult;
  TwoNorm<T> norm;

  /*for(int i = 0; i < aMatrix.rows(); i++)
  {
    for(int j = 0; j < aMatrix.rows(); j++)
    {
      q[i][j] = 0;
      r[i][j] = 0;
      //q(i, j) = 0;
      //r(i, j) = 0;
    }
  }*/

  int rows = static_cast<int>(aMatrix.rows());
  for(int i = 0; i < rows; i++)
  {
    if(i == 0)
    {
      normResult = norm(aMatrix.getColumn(0));
      if(normResult == 0)
      {
        throw std::domain_error("Error in QRDecomposition: Attempted division by zero");
      }

      for(int j = 0; j < rows; j++)
      {
        q[j][0] = aMatrix[j][0] * (1 / normResult);
      }
    }
    else
    {
      AlgebraVector<T> temp(aMatrix[0]);

      temp = r[0][i] * q.getColumn(0);
      //temp = r(0, i) * q.getColumn(0);

      for(int k = 1; k < i; k++)
      {
        temp += r[k][i] * q.getColumn(k);
        //temp += r(k, i) * q.getColumn(k);
      }

      temp = aMatrix.getColumn(i) - temp;

      r[i][i] = norm(temp);
      //r(i, i) = norm(temp);

      //if(r(i, i) == 0)
      if(r[i][i] == 0)
      {
        throw std::domain_error("QRDecomposition: Attempted division by zero");
      }

      for(int j = 0; j < rows; j++)
      {
        q[j][i] = (1 / r[i][i]) * temp[j];
        //q(j, i) = (1 / r(i, i)) * temp[j];
      }
    }

    for(int j = 0; j < rows; j++)
    {
      r[i][j] = aMatrix.getColumn(j) * q.getColumn(i);
      //r(i, j) = aMatrix.getColumn(j) * q.getColumn(i);
    }
  }
}
开发者ID:Jmgiacone,项目名称:CS5201,代码行数:72,代码来源:qrdecomposition.hpp

示例3: operator

 void operator()(AbstractMatrix<M>& m) const {
     for (int i=0; i<m.rows(); ++i)
         for (int j=0; j<m.cols(); ++j)
             m(i,j) = rand() * a - b;
 }
开发者ID:qqchen,项目名称:LATL,代码行数:5,代码来源:expr.hpp


注:本文中的AbstractMatrix::rows方法示例由纯净天空整理自Github/MSDocs等开源代码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。