TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2025-04-19 07:00:10 +02:00
Code Issues Projects Releases Wiki Activity

Checking context in SM

Browse Source
This commit is contained in:
Anthony Scemama 2021-10-06 23:44:06 +02:00
parent 71ad7abb7f
commit af54e7a7dc
1 changed files with 62 additions and 42 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
org/qmckl_sherman_morrison_woodbury.org
View File

@ -111,6 +111,10 @@ qmckl_exit_code qmckl_sherman_morrison_c_(const qmckl_context context,
const double* breakdown_p, const double* breakdown_p,
double* Slater_inv) { double* Slater_inv) {
if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) {
return QMCKL_NULL_CONTEXT;
}
const uint64_t Dim = *Dim_p; const uint64_t Dim = *Dim_p;
const uint64_t N_updates = *N_updates_p; const uint64_t N_updates = *N_updates_p;
const double breakdown = *breakdown_p; const double breakdown = *breakdown_p;
@ -364,6 +368,10 @@ qmckl_exit_code qmckl_woodbury_2_c_(const qmckl_context context,
D := V * S^{-1}, 2 x dim D := V * S^{-1}, 2 x dim
,*/ ,*/
if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) {
return QMCKL_NULL_CONTEXT;
}
const uint64_t Dim = *Dim_p; const uint64_t Dim = *Dim_p;
const double breakdown = *breakdown_p; const double breakdown = *breakdown_p;
@ -581,6 +589,10 @@ qmckl_exit_code qmckl_woodbury_3_c_(const qmckl_context context,
D := V * S^{-1}, 3 x dim D := V * S^{-1}, 3 x dim
,*/ ,*/
if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) {
return QMCKL_NULL_CONTEXT;
}
const uint64_t Dim = *Dim_p; const uint64_t Dim = *Dim_p;
const double breakdown = *breakdown_p; const double breakdown = *breakdown_p;
@ -814,6 +826,10 @@ qmckl_exit_code qmckl_sherman_morrison_splitting_c_(const qmckl_context context,
const double* breakdown, const double* breakdown,
double* Slater_inv) { double* Slater_inv) {
if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) {
return QMCKL_NULL_CONTEXT;
}
double later_updates[*Dim * *N_updates]; double later_updates[*Dim * *N_updates];
uint64_t later_index[*N_updates]; uint64_t later_index[*N_updates];
uint64_t later = 0; uint64_t later = 0;
@ -989,6 +1005,10 @@ qmckl_exit_code qmckl_sherman_morrison_smw32s_c_(const qmckl_context context,
const double* breakdown_p, const double* breakdown_p,
double* Slater_inv) { double* Slater_inv) {
if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) {
return QMCKL_NULL_CONTEXT;
}
const uint64_t Dim = *Dim_p; const uint64_t Dim = *Dim_p;
const uint64_t N_updates = *N_updates_p; const uint64_t N_updates = *N_updates_p;
const double breakdown = *breakdown_p; const double breakdown = *breakdown_p;