Changeset a5b9608 in tol
- Timestamp:
- Jul 29, 2015, 8:03:06 PM (11 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- b44a3b0
- Parents:
- 4e57976
- Location:
- tol
- Files:
-
- 3 edited
-
bmath/blinalg/tol_lapack.cpp (modified) (2 diffs)
-
bmath/blinalg/tol_lapack.h (modified) (2 diffs)
-
btol/vmatrix_type/vmat_chol.cpp (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
tol/bmath/blinalg/tol_lapack.cpp
r4e57976 ra5b9608 68 68 int n = S.Columns(); 69 69 double* X_ = X.GetData().GetBuffer(); 70 int res = clapack_dpotrf( int(CblasRowMajor), LAPACK_UPLO(Uplo), n, X_, n);70 int res = clapack_dpotrf( CblasRowMajor, LAPACK_UPLO(Uplo), n, X_, n ); 71 71 replicateTriang(X,Uplo,false); 72 72 return(res); … … 103 103 int n = A.Columns(); 104 104 double* X_ = X.GetData().GetBuffer(); 105 int res = clapack_dpotri( int(CblasRowMajor), LAPACK_UPLO(Uplo), n, X_, n);105 int res = clapack_dpotri( CblasRowMajor, LAPACK_UPLO(Uplo), n, X_, n ); 106 106 replicateTriang(X,Uplo,true); 107 107 return(res); -
tol/bmath/blinalg/tol_lapack.h
r4e57976 ra5b9608 32 32 #undef I 33 33 #define LAPACK_UPLO(UpLo) (UpLo==CblasLower?'L':'U') 34 #define LAPACK_ORDER(Order) (int)(order) 34 35 #define clapack_dpotrf LAPACKE_dpotrf 35 36 #define clapack_dpotri LAPACKE_dpotri … … 38 39 #include <clapack.h> 39 40 #define LAPACK_UPLO(UpLo) UpLo 41 #define LAPACK_ORDER(Order) Order 40 42 END_DECLS 41 43 #endif -
tol/btol/vmatrix_type/vmat_chol.cpp
r4e57976 ra5b9608 39 39 double* x = (double*)L.s_.blasRdense_->x; 40 40 memcpy(x, y, n*n*sizeof(double)); 41 res = clapack_dpotrf( CblasColMajor, LAPACK_UPLO(CblasLower), n, x, n);41 res = clapack_dpotrf( CblasColMajor, LAPACK_UPLO(CblasLower), n, x, n ); 42 42 for(i=0; i<n; i++) for(j=i+1; j<n; j++) x[j*n+i] = 0; 43 43 isOk = (res==0);
Note: See TracChangeset
for help on using the changeset viewer.
