TREX/Sherman-Morrison
1
0
Fork 0
mirror of https://github.com/TREX-CoE/Sherman-Morrison.git synced 2024-12-25 13:53:56 +01:00
Code Issues Projects Releases Wiki Activity

Cleaned-up main loop.

Browse Source
This commit is contained in:
François Coppens 2021-03-04 18:36:51 +01:00
parent 066e50914c
commit a9604b3838
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/SM_MaponiA3.cpp
View File

@ -76,18 +76,18 @@ void MaponiA3(double *Slater_inv, unsigned int Dim,
} tmp = p[l+1]; p[l+1] = p[lbar]; p[lbar] = tmp;
component = Updates_index[p[l+1] - 1];
beta = 1 + ylk[l+1 - 1][p[l+1]][component];
beta = 1 + ylk[l][p[l+1]][component];
if (beta == 0) {
cout << "Break-down occured. Exiting..." << endl;
exit(1);
}
for (k = l+1 + 1; k < N_updates + 1; k++) {
alpha = ylk[l+1 - 1][p[k]][component] / beta;
for (k = l+2; k < N_updates + 1; k++) {
alpha = ylk[l][p[k]][component] / beta;
cout << "( l, k, p[k], component ) = (" << l << ", " << k << ", " << p[k] << ", " << component << ")" << endl;
for (i = 1; i < Dim + 1; i++) {
cout << "ylk[" << l << "][p[" << k << "]][" << i << "] = ylk[" << l - 1 << "][p[" << k << "]][" << i << "] - alpha * ylk[" << l - 1 << "][p[" << l << "]][" << i << "]" << endl;
ylk[l+1][p[k]][i] = ylk[l+1 - 1][p[k]][i]
- alpha * ylk[l+1 - 1][p[l+1]][i];
ylk[l+1][p[k]][i] = ylk[l][p[k]][i]
- alpha * ylk[l][p[l+1]][i];
}
}
}