TREX/Sherman-Morrison
1
0
Fork 0
mirror of https://github.com/TREX-CoE/Sherman-Morrison.git synced 2024-12-27 06:43:55 +01:00
Code Issues Projects Releases Wiki Activity

Merge pull request #14 from fmgjcoppens/fix/repl_matMul_calls

Browse Source 
Replaced old matMul() function calls.
This commit is contained in:
François Coppens 2021-03-11 15:09:07 +01:00 committed by GitHub
commit 4c7a567750
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
tests/test_external_h5.cpp
View File

@ -74,7 +74,8 @@ int test_cycle(H5File file, int cycle) {
showMatrix(slater_inverse, dim, "NEW Inverse"); showMatrix(slater_inverse, dim, "NEW Inverse");
#endif #endif
double * res = matMul(slater_matrix, slater_inverse, dim); double * res = new double[dim*dim] {0};
matMul(slater_matrix, slater_inverse, res, dim);
bool ok = is_identity(res, dim, 0.5e-4); bool ok = is_identity(res, dim, 0.5e-4);
#ifdef DEBUG #ifdef DEBUG

3
tests/test_internal_h5.cpp
View File

@ -74,7 +74,8 @@ int test_cycle(H5File file, int cycle) {
showMatrix(slater_inverse, dim, "NEW Inverse"); showMatrix(slater_inverse, dim, "NEW Inverse");
#endif #endif
double * res = matMul(slater_matrix, slater_inverse, dim); double * res = new double[dim*dim] {0};
matMul(slater_matrix, slater_inverse, res, dim);
bool ok = is_identity(res, dim, 0.5e-4); bool ok = is_identity(res, dim, 0.5e-4);
#ifdef DEBUG #ifdef DEBUG