Revision | a54239bb8579c29e73aaa703ab63b2c60b16bb5d (tree) |
---|---|
Zeit | 2012-12-05 16:55:12 |
Autor | Katsuhiko Nishimra <ktns.87@gmai...> |
Commiter | Katsuhiko Nishimra |
Copy matrix using dcopy. #28764
git-svn-id: https://svn.sourceforge.jp/svnroot/molds/trunk@1178 1136aad2-a195-0410-b898-f5ea1d11b9d8
@@ -496,12 +496,7 @@ void BFGS::ShiftHessianRedundantMode(double** matrixHessian, | ||
496 | 496 | // Diagonalize hessian |
497 | 497 | MallocerFreer::GetInstance()->Malloc(&vectorHessianEigenValues, dimension); |
498 | 498 | MallocerFreer::GetInstance()->Malloc(&vectorsHessianModes, dimension, dimension); |
499 | -#pragma omp parallel for schedule(auto) | |
500 | - for(int i=0; i<dimension; i++){ | |
501 | - for(int j=0; j<dimension; j++){ | |
502 | - vectorsHessianModes[i][j] = matrixHessian[i][j]; | |
503 | - } | |
504 | - } | |
499 | + MolDS_wrappers::Blas::GetInstance()->Dcopy(dimension*dimension, &matrixHessian[0][0], &vectorsHessianModes[0][0]); | |
505 | 500 | bool calcEigenVectors = true; |
506 | 501 | MolDS_wrappers::Lapack::GetInstance()->Dsyevd(&vectorsHessianModes[0], |
507 | 502 | &vectorHessianEigenValues[0], |