TREX/qmckl
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmckl.git synced 2025-01-03 10:06:09 +01:00
Code Issues Projects Releases Wiki Activity

Fixed Fortran function call bug.

Browse Source
This commit is contained in:
Francois Coppens 2023-02-23 15:57:35 +01:00
parent 8216f682b3
commit 7a97aa4a77
1 changed files with 15 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
org/qmckl_sherman_morrison_woodbury.org
View File

@ -29,7 +29,7 @@ int main() {
This is the range that determines the how many high performance kernel instantces will be generated, using the C-function templates defined in the sections below. If the name of the C-function template is called ~qmckl_kernel_{Dim}~, then ~range(K, L+1)~ will results in kernel instances from ~qmckl_kernel_K~ to ~qmckl_kernel_L~. This is the range that determines the how many high performance kernel instantces will be generated, using the C-function templates defined in the sections below. If the name of the C-function template is called ~qmckl_kernel_{Dim}~, then ~range(K, L+1)~ will results in kernel instances from ~qmckl_kernel_K~ to ~qmckl_kernel_L~.
#+NAME:kernel_generator_range #+NAME:kernel_generator_range
#+begin_src python :noweb yes :exports none #+begin_src python :noweb yes :exports none
range(2, 3) range(2, 22)
#+end_src #+end_src
@ -244,7 +244,7 @@ integer function qmckl_sm_naive_doc_f(context, &
end do end do
! Copy updated inverse back to s_inv ! Copy updated inverse back to s_inv
call copy_back(Inverse, s_inv, lds, dim) call copy_back_inv(Inverse, s_inv, lds, dim)
info = QMCKL_SUCCESS info = QMCKL_SUCCESS
@ -1451,6 +1451,8 @@ integer recursive function qmckl_sm_splitting_doc_f( &
real*8 , intent(inout) :: s_inv(dim * lds) real*8 , intent(inout) :: s_inv(dim * lds)
real*8 , intent(inout) :: determinant real*8 , intent(inout) :: determinant
integer , external :: qmckl_sm_splitting_core_doc_f
integer*8 :: Later integer*8 :: Later
integer*8 , dimension(nupdates) :: Later_index integer*8 , dimension(nupdates) :: Later_index
real*8 , dimension(nupdates * lds) :: Later_updates real*8 , dimension(nupdates * lds) :: Later_updates
@ -1690,17 +1692,7 @@ qmckl_exit_code qmckl_sm_splitting(
Slater_inv, Slater_inv,
determinant); determinant);
#else #else
// return qmckl_sm_splitting_doc( return qmckl_sm_splitting_doc(
// context,
// LDS,
// Dim,
// N_updates,
// Updates,
// Updates_index,
// breakdown,
// Slater_inv,
// determinant);
return qmckl_sm_splitting_hpc(
context, context,
LDS, LDS,
Dim, Dim,
@ -1710,6 +1702,16 @@ qmckl_exit_code qmckl_sm_splitting(
breakdown, breakdown,
Slater_inv, Slater_inv,
determinant); determinant);
// return qmckl_sm_splitting_hpc(
// context,
// LDS,
// Dim,
// N_updates,
// Updates,
// Updates_index,
// breakdown,
// Slater_inv,
// determinant);
#endif #endif
return QMCKL_SUCCESS; return QMCKL_SUCCESS;