TREX/Sherman-Morrison
1
0
Fork 0
mirror of https://github.com/TREX-CoE/Sherman-Morrison.git synced 2024-12-26 14:23:47 +01:00
Code Issues Projects Releases Wiki Activity

Replaced old matMul() function calls.

Browse Source
This commit is contained in:
François Coppens 2021-03-11 15:04:35 +01:00
parent d6169cc6ab
commit a1c77a22c8
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