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

Removed unnecessary __restrict keywords from public exposed functoin headers.

Browse Source
This commit is contained in:
François 2023-01-26 11:50:58 +01:00
parent 70241915aa
commit 4a9cecff61
1 changed files with 51 additions and 51 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

102
org/qmckl_sherman_morrison_woodbury.org
View File

@ -518,11 +518,11 @@ assert(rc == QMCKL_SUCCESS);
const qmckl_context context, const qmckl_context context,
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant); double* determinant);
#+end_src #+end_src
*** C source *** C source
@ -531,11 +531,11 @@ assert(rc == QMCKL_SUCCESS);
qmckl_exit_code qmckl_woodbury_2(const qmckl_context context, qmckl_exit_code qmckl_woodbury_2(const qmckl_context context,
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant) { double* determinant) {
/* /*
C := S^{-1} * U, dim x 2 C := S^{-1} * U, dim x 2
B := 1 + V * C, 2 x 2 B := 1 + V * C, 2 x 2
@ -590,8 +590,8 @@ qmckl_exit_code qmckl_woodbury_2(const qmckl_context context,
// tmp = B^{-1}D : 2 x LDS // tmp = B^{-1}D : 2 x LDS
double __attribute__((aligned(8))) tmp[2 * LDS]; double __attribute__((aligned(8))) tmp[2 * LDS];
double* __restrict r1dim = &(Slater_inv[row1 * LDS]); double* r1dim = &(Slater_inv[row1 * LDS]);
double* __restrict r2dim = &(Slater_inv[row2 * LDS]); double* r2dim = &(Slater_inv[row2 * LDS]);
IVDEP IVDEP
ALIGNED ALIGNED
for (uint64_t j = 0; j < LDS; j++) { for (uint64_t j = 0; j < LDS; j++) {
@ -733,11 +733,11 @@ assert(rc == QMCKL_SUCCESS);
const qmckl_context context, const qmckl_context context,
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant); double* determinant);
#+end_src #+end_src
*** C source *** C source
@ -746,11 +746,11 @@ assert(rc == QMCKL_SUCCESS);
qmckl_exit_code qmckl_woodbury_3(const qmckl_context context, qmckl_exit_code qmckl_woodbury_3(const qmckl_context context,
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant) { double* determinant) {
/* /*
C := S^{-1} * U, dim x 3 C := S^{-1} * U, dim x 3
B := 1 + V * C, 3 x 3 B := 1 + V * C, 3 x 3
@ -820,9 +820,9 @@ qmckl_exit_code qmckl_woodbury_3(const qmckl_context context,
// tmp = B^{-1}D : 3 x LDS // tmp = B^{-1}D : 3 x LDS
double __attribute__((aligned(8))) tmp[3 * LDS]; double __attribute__((aligned(8))) tmp[3 * LDS];
double* __restrict r1dim = &(Slater_inv[row1 * LDS]); double* r1dim = &(Slater_inv[row1 * LDS]);
double* __restrict r2dim = &(Slater_inv[row2 * LDS]); double* r2dim = &(Slater_inv[row2 * LDS]);
double* __restrict r3dim = &(Slater_inv[row3 * LDS]); double* r3dim = &(Slater_inv[row3 * LDS]);
IVDEP IVDEP
ALIGNED ALIGNED
for (uint64_t j = 0; j < LDS; j++) { for (uint64_t j = 0; j < LDS; j++) {
@ -975,11 +975,11 @@ assert(rc == QMCKL_SUCCESS);
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const uint64_t N_updates, const uint64_t N_updates,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant); double* determinant);
#+end_src #+end_src
*** C source *** C source
@ -989,11 +989,11 @@ qmckl_exit_code qmckl_sherman_morrison_splitting(const qmckl_context context,
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const uint64_t N_updates, const uint64_t N_updates,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant) { double* determinant) {
if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) { if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) {
return qmckl_failwith(context, return qmckl_failwith(context,
@ -1141,11 +1141,11 @@ assert(rc == QMCKL_SUCCESS);
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const uint64_t N_updates, const uint64_t N_updates,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant); double* determinant);
#+end_src #+end_src
*** C source *** C source
@ -1155,11 +1155,11 @@ qmckl_exit_code qmckl_sherman_morrison_smw32s(const qmckl_context context,
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const uint64_t N_updates, const uint64_t N_updates,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict determinant) { double* determinant) {
if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) { if (qmckl_context_check(context) == QMCKL_NULL_CONTEXT) {
return qmckl_failwith(context, return qmckl_failwith(context,
@ -1399,14 +1399,14 @@ These functions can only be used internally by the kernels in this module.
const uint64_t LDS, const uint64_t LDS,
const uint64_t Dim, const uint64_t Dim,
const uint64_t N_updates, const uint64_t N_updates,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict later_updates, double* later_updates,
uint64_t* __restrict later_index, uint64_t* later_index,
uint64_t* __restrict later, uint64_t* later,
double* __restrict determinant); double* determinant);
#+end_src #+end_src
*** C source *** C source
@ -1415,14 +1415,14 @@ These functions can only be used internally by the kernels in this module.
qmckl_exit_code qmckl_slagel_splitting(uint64_t LDS, qmckl_exit_code qmckl_slagel_splitting(uint64_t LDS,
uint64_t Dim, uint64_t Dim,
uint64_t N_updates, uint64_t N_updates,
const double* __restrict Updates, const double* Updates,
const uint64_t* __restrict Updates_index, const uint64_t* Updates_index,
const double breakdown, const double breakdown,
double* __restrict Slater_inv, double* Slater_inv,
double* __restrict later_updates, double* later_updates,
uint64_t* __restrict later_index, uint64_t* later_index,
uint64_t* __restrict later, uint64_t* later,
double* __restrict determinant) { double* determinant) {
double __attribute__((aligned(8))) C[LDS]; double __attribute__((aligned(8))) C[LDS];
double __attribute__((aligned(8))) D[LDS]; double __attribute__((aligned(8))) D[LDS];