C++ (Cpp) RealMatrix::copy示例

示例#1

文件： AccpmLASolve.cpp 项目： manasdas17/scilab-mip

int
AccpmLACholeskyFactor(const RealMatrix &A, RealMatrix &L)
{
  L.copy(A);

  char uplo = 'L';
  long int n = A.size(0);
  long int lda = A.inc(0) * A.gdim(0);
  long int info;
  F77NAME(dpotrf) (&uplo, &n, &L(0,0), &lda, &info);
  if (info > 0) {
    std::cerr << "AccpmLACholeskyFactor: Matrix is not Symmetric-positive-definite." 
  	      << std::endl;
  } else if (info < 0) {
    std::cerr << "AccpmLACholeskyFactor: argument " << -info << " is invalid" 
  	      << std::endl;
  }
  return info;
}