TREX/Sherman-Morrison
1
0
Fork 0
mirror of https://github.com/TREX-CoE/Sherman-Morrison.git synced 2025-01-12 22:18:36 +01:00
Code Issues Projects Releases Wiki Activity

Cleanup

Browse Source
This commit is contained in:
Francois Coppens 2021-03-02 10:24:14 +01:00 committed by François Coppens
parent 983f87d504
commit 3f60779700
2 changed files with 4 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
include/Helpers.hpp
View File

@ -79,7 +79,7 @@ T *matMul(T *A, T *B, unsigned int M) {
template<typename T> template<typename T>
void matMul2(T *A, T *B, T *C, unsigned int M) { void matMul2(T *A, T *B, T *C, unsigned int M) {
memset(C,0,M*M*sizeof(T)); memset(C, 0, M*M*sizeof(T));
for (unsigned int i = 0; i < M; i++) { for (unsigned int i = 0; i < M; i++) {
for (unsigned int j = 0; j < M; j++) { for (unsigned int j = 0; j < M; j++) {
for (unsigned int k = 0; k < M; k++) { for (unsigned int k = 0; k < M; k++) {

6
src/SM_MaponiA3.cpp
View File

@ -68,10 +68,9 @@ void MaponiA3(double *Slater_inv, unsigned int Dim,
} }
} }
double *next = new double[Dim*Dim] {0};
double *last;
last = Slater_inv;
// Construct A-inverse from A0-inverse and the ylk // Construct A-inverse from A0-inverse and the ylk
double *last = Slater_inv;
double *next = new double[Dim*Dim] {0};
for (l = 0; l < N_updates; l++) { for (l = 0; l < N_updates; l++) {
k = l + 1; k = l + 1;
component = Updates_index[p[k] - 1]; component = Updates_index[p[k] - 1];
@ -82,7 +81,6 @@ void MaponiA3(double *Slater_inv, unsigned int Dim,
* ylk[l][p[k]][i + 1] / beta; * ylk[l][p[k]][i + 1] / beta;
} }
} }
matMul2(Al, last, next, Dim); matMul2(Al, last, next, Dim);
double *tmp = next; double *tmp = next;
next = last; next = last;