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

Small bugfix in qmckl_slagel_splitting()

Browse Source
This commit is contained in:
Francois Coppens 2022-07-21 08:16:25 +02:00
parent 0a083e2875
commit f35ad6a777
2 changed files with 3 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
independent_test_harness/sm.c
View File

@ -322,9 +322,7 @@ uint32_t qmckl_woodbury_k(const uint64_t vLDS,
j += min(abs(ipiv[i] - i), 1);
det *= B[(N_updates + 1) * i];
}
if (j & 1 == 0) {
det = -det;
}
if (j & 1 == 0) det = -det; // multiply det with -1 if j is even
// Check if determinant of B is not too close to zero
if (fabs(det) < breakdown) {
@ -368,8 +366,7 @@ uint32_t qmckl_slagel_splitting(
const uint32_t LDS = 24;
const uint32_t Dim = 21;
// double __attribute__((aligned(8))) C[N_updates * Dim];
double __attribute__((aligned(8))) C[Dim];
double __attribute__((aligned(8))) C[LDS];
double __attribute__((aligned(8))) D[LDS];
uint32_t l = 0;

2
independent_test_harness/test.c
View File

@ -4,7 +4,7 @@
#define DATASET "dataset_329d_zeropadded_cm.hdf5"
// #define DATASET "dataset_15784d_zeropadded_cm.hdf5"
#define REPETITIONS 100000
#define REPETITIONS 1000000
uint64_t n_splits;
uint64_t block_fail;