The qmckl.h header file installed in the ${prefix}/include directory
@@ -385,12 +385,12 @@ Both files are located in the include/ directory.
In a traditional source code, most of the lines of source files of a program @@ -435,8 +435,8 @@ interactively, in the same spirit as Jupyter notebooks.
For a tutorial on literate programming with org-mode, follow this link. @@ -467,8 +467,8 @@ org-mode.
Most of the codes of the TREX CoE are written in Fortran with some scripts in @@ -516,8 +516,8 @@ For more guidelines on using Fortran to generate a C interface, see
The authors should follow the recommendations of the C99
@@ -535,8 +535,8 @@ Compliance can be checked with cppcheck as:
The proposed API should allow the library to: deal with memory transfers @@ -547,8 +547,8 @@ functions (see below).
To avoid namespace collisions, we use qmckl_ as a prefix for all exported
@@ -573,8 +573,8 @@ form is allowed.
In the C language, the number of bits used by the integer types can change @@ -606,15 +606,15 @@ bindings in other languages in other repositories.
Global variables should be avoided in the library, because it is
possible that one single program needs to use multiple instances
of the library. To solve this problem we propose to use a pointer
to a context variable, built by the library with the
-qmckl_context_create function. The =context= contains the global
+qmckl_context_create function. The =context= contains the global
state of the library, and is used as the first argument of many
QMCkl functions.
A single qmckl.h header to be distributed by the library
@@ -717,8 +717,8 @@ and the types definitions should be written in the *_f_type.f90 fil
Low-level functions are very simple functions which are leaves of @@ -727,14 +727,14 @@ the function call tree (they don't call any other QMCkl function).
These functions are pure, and unaware of the QMCkl
-context. They are not allowed to allocate/deallocate memory, and
+context. They are not allowed to allocate/deallocate memory, and
if they need temporary memory it should be provided in input.
High-level functions are at the top of the function call tree. @@ -747,27 +747,27 @@ temporary storage, to simplify the use of accelerators.
The high-level functions should be pure, unless the introduction
of non-purity is justified. All the side effects should be made in
-the context variable.
+the context variable.
The number of bits of precision required for a function should be
given as an input of low-level computational functions. This input
will be used to define the values of the different thresholds that
might be used to avoid computing unnecessary noise. High-level
-functions will use the precision specified in the context
+functions will use the precision specified in the context
variable.
Reducing the scaling of an algorithm usually implies also reducing @@ -783,7 +783,7 @@ implemented adapted to different problem sizes.
The following arrays are stored in the context: @@ -503,8 +503,8 @@ coefficient = [ 0.006068, 0.045308, 0.202822, 0.503903, 0.383421,
typedef struct qmckl_ao_basis_struct {
@@ -533,8 +533,8 @@ struct is then initialized and provided == true.
When all the data for the AOs have been provided, the following
@@ -548,8 +548,8 @@ function returns true.
To set the basis set, all the following functions need to be @@ -572,17 +572,17 @@ called. When
The qmckl_ao_power function computes all the powers of the n
@@ -594,7 +594,7 @@ the \(n\) points:
\[ P_{ik} = X_i^k \]
context is not QMCKL_NULL_CONTEXTqmckl_exit_code qmckl_ao_power ( @@ -682,8 +682,8 @@ the \(n\) points:
integer function qmckl_ao_power_f(context, n, X, LMAX, P, ldp) result(info) @@ -734,15 +734,15 @@ the \(n\) points:
integer(c_int32_t) function test_qmckl_ao_power(context) bind(C) @@ -793,8 +793,8 @@ the \(n\) points:
A polynomial is centered on a nucleus \(\mathbf{R}_i\)
@@ -839,7 +839,7 @@ Laplacians at a given point in space, of all polynomials with an
angular momentum up to lmax.
context is not QMCKL_NULL_CONTEXTqmckl_exit_code qmckl_ao_polynomial_vgl ( @@ -963,8 +963,8 @@ For example, with a=0, b=2 and c=1 the string is "yyz"
integer function qmckl_ao_polynomial_vgl_f(context, X, R, lmax, n, L, ldl, VGL, ldv) result(info) @@ -1099,16 +1099,16 @@ For example, with a=0, b=2 and c=1 the string is "yyz"
integer(c_int32_t) function test_qmckl_ao_polynomial_vgl(context) bind(C) @@ -1213,12 +1213,12 @@ assert(0 == test_qmckl_ao_polynomial_vgl(context));
qmckl_ao_gaussian_vgl computes the values, gradients and
@@ -1450,21 +1450,21 @@ Requirements
The context variable is a handle for the state of the library,
@@ -338,7 +338,7 @@ A value of QMCKL_NULL_CONTEXT for the context is equivalent to a
typedef int64_t qmckl_context ; +typedef int64_t qmckl_context ; #define QMCKL_NULL_CONTEXT (qmckl_context) 0
ctx is a qmckl_context_struct* pointer.
The context keeps a ``date'' that allows to check which data needs @@ -367,7 +367,7 @@ coordinates are updated.
When a new element is added to the context, the functions -qmcklcontextcreate, qmcklcontextdestroy and qmcklcontextcopy +qmcklcontextcreate, qmcklcontextdestroy and qmcklcontextcopy should be updated inorder to make deep copies.
@@ -416,8 +416,8 @@ if the context is valid,QMCKL_NULL_CONTEXT otherwise.
To create a new context, qmckl_context_create() should be used.
@@ -491,8 +491,8 @@ To create a new context, qmckl_context_create() should be used.
For thread safety, the context may be locked/unlocked. The lock is @@ -537,8 +537,8 @@ number of times the thread has locked it is saved in the
qmckl_context_copy makes a deep copy of a context. It returns
@@ -586,8 +586,8 @@ number of times the thread has locked it is saved in the
The context is destroyed with qmckl_context_destroy, leaving the ancestors untouched.
@@ -641,7 +641,7 @@ It frees the context, and returns the previous context.
qmckl_distance_sqqmckl_distance_sq
qmckl_distance_sq computes the matrix of the squared distances
@@ -379,7 +379,7 @@ between all pairs of points in two sets, one point within each set:
\]
context is not QMCKL_NULL_CONTEXTqmckl_exit_code qmckl_distance ( @@ -512,8 +512,8 @@ between all pairs of points in two sets, one point within each set:
integer function qmckl_distance_sq_f(context, transa, transb, m, n, & @@ -648,8 +648,8 @@ between all pairs of points in two sets, one point within each set:
This function might be more efficient when A and B are
@@ -659,12 +659,12 @@ transposed.
qmckl_distanceqmckl_distance
qmckl_distance computes the matrix of the distances between all
@@ -677,7 +677,7 @@ pairs of points in two sets, one point within each set:
\]
context is not QMCKL_NULL_CONTEXTqmckl_exit_code qmckl_distance ( @@ -810,8 +810,8 @@ pairs of points in two sets, one point within each set:
integer function qmckl_distance_f(context, transa, transb, m, n, & @@ -950,8 +950,8 @@ pairs of points in two sets, one point within each set:
This function might be more efficient when A and B are
@@ -964,7 +964,7 @@ transposed.
The following data stored in the context: @@ -421,8 +421,8 @@ The following data stored in the context:
typedef struct qmckl_electron_struct {
@@ -450,8 +450,8 @@ struct is then initialized and provided == true.
When all the data relative to electrons have been set, the
@@ -465,8 +465,8 @@ following function returns true.
To set the data relative to the electrons in the context, the @@ -496,8 +496,8 @@ electrons have been set.
/* Reference input data */
@@ -542,8 +542,8 @@ rc = qmckl_set_electron_coord (context, coord);
The computed data is stored in the context so that it can be reused @@ -556,12 +556,12 @@ current date is stored.
qmckl_exit_code qmckl_get_electron_ee_distance(qmckl_context context, double* distance); @@ -570,10 +570,10 @@ current date is stored.