Angular Momentum#

group AngularMomentum

Functions

cuestStatus_t cuestAngularMomentumCompute(
cuestHandle_t handle,
const cuestOEIntPlan_t plan,
const cuestAngularMomentumComputeParameters_t parameters,
cuestWorkspace_t *temporaryWorkspace,
cuestAngularMomentumComputeParametersComponent_t component,
const double *origin,
double *outMatrix
)#

Compute the atomic orbital angular momentum matrix using the provided cuestOEIntPlan_t plan.

The output buffer must be preallocated and sized for nao × nao elements, where nao is the number of atomic orbitals in the basis. Temporary workspace must be allocated using size returned by cuestAngularMomentumComputeWorkspaceQuery.

Since

v0.2.0

Parameters:

  • handle – [in] cuEST handle. Must not be NULL.

  • plan – [in] One-electron integral computation plan (opaque handle). Must not be NULL.

  • parameters – [in] Compute parameters (cuestAngularMomentumComputeParameters_t). Must not be NULL

  • temporaryWorkspace – [in] Temporary workspace buffers for computation (preallocated). Must not be NULL.

  • component – [in] Enum (cuestAngularMomentumComputeParametersComponent_t) defining the angular momentum operator (i.e. \(L_x\), \(L_y\), \(L_z\)).

  • origin – [in] Cartesian coordinates defining the operator origin [size: 3]. CPU array. Must not be NULL.

  • outMatrix – [out] Output GPU buffer for computed angular momentum matrix elements [size: nao × nao]. Must not be NULL. Elements of this array are overwritten by the calculation results.

Returns:

  • CUEST_STATUS_SUCCESS on success;

  • CUEST_STATUS_INVALID_HANDLE if the cuEST handle is NULL;

  • CUEST_STATUS_NULL_POINTER if any required pointer is NULL;

  • CUEST_STATUS_INVALID_TYPE if opaque handles are not the correct type;

  • CUEST_STATUS_INVALID_ARGUMENT if component is not a valid cuestAngularMomentumComputeParametersComponent_t value;

  • CUEST_STATUS_EXCEPTION or CUEST_STATUS_UNKNOWN_ERROR otherwise.

cuestStatus_t cuestAngularMomentumComputeWorkspaceQuery(
cuestHandle_t handle,
const cuestOEIntPlan_t plan,
const cuestAngularMomentumComputeParameters_t parameters,
cuestWorkspaceDescriptor_t *temporaryWorkspaceDescriptor,
cuestAngularMomentumComputeParametersComponent_t component,
const double *origin,
double *outMatrix
)#

Query the temporary workspace required (in bytes) for angular momentum matrix calculation.

This function computes the workspace requirements for a requested angular momentum calculation. Fill the provided workspace descriptor with host and device buffer sizes required. The outMatrix pointer may be NULL.

Since

v0.2.0

Parameters:

  • handle – [in] cuEST handle. Must not be NULL.

  • plan – [in] One-electron integral computation plan (opaque handle). Must not be NULL.

  • parameters – [in] Compute parameters (cuestAngularMomentumComputeParameters_t). Must not be NULL

  • temporaryWorkspaceDescriptor – [out] Output descriptor for temporary workspace sizes. Must not be NULL.

  • component – [in] Enum (cuestAngularMomentumComputeParametersComponent_t) defining the angular momentum operator (i.e. \(L_x\), \(L_y\), \(L_z\)).

  • origin – [in] Cartesian coordinates defining the operator origin [size: 3]. CPU array. Must not be NULL.

  • outMatrix – User matrix buffer (may be NULL during the workspace query).

Returns:

  • CUEST_STATUS_SUCCESS on success;

  • CUEST_STATUS_INVALID_HANDLE if the cuEST handle is NULL;

  • CUEST_STATUS_NULL_POINTER if any required pointer is NULL;

  • CUEST_STATUS_INVALID_TYPE if opaque handles are not the correct type;

  • CUEST_STATUS_INVALID_ARGUMENT if component is not a valid cuestAngularMomentumComputeParametersComponent_t value;

  • CUEST_STATUS_EXCEPTION or CUEST_STATUS_UNKNOWN_ERROR otherwise.

group AngularMomentumGrad

Functions

cuestStatus_t cuestAngularMomentumDerivativeCompute(
cuestHandle_t handle,
const cuestOEIntPlan_t plan,
const cuestAngularMomentumDerivativeComputeParameters_t parameters,
cuestWorkspace_t *temporaryWorkspace,
cuestAngularMomentumComputeParametersComponent_t component,
const double *origin,
const double *densityMatrix,
double *outGradient
)#

Compute the electronic gradient (derivative) of the atomic orbital angular momentum matrix using the supplied integral plan and density matrix.

Both the input density matrix and the output gradient must be preallocated. The density matrix must be an (nao x nao) array and the output gradient must have space for (natom x 3) elements, where nao is the number of atomic orbitals and natom the number of atoms in the system. All required temporary workspace must be allocated by the user prior to calling this function.

Since

v0.2.0

Parameters:

  • handle – [in] cuEST handle. Must not be NULL.

  • plan – [in] One-electron integral computation plan (opaque handle). Must not be NULL.

  • parameters – [in] Compute parameters (cuestAngularMomentumDerivativeComputeParameters_t). Must not be NULL

  • temporaryWorkspace – [in] Temporary workspace buffers for computation (preallocated). Must not be NULL.

  • component – [in] Enum (cuestAngularMomentumComputeParametersComponent_t) defining the angular momentum operator (i.e. \(L_x\), \(L_y\), \(L_z\)).

  • origin – [in] Cartesian coordinates defining the operator origin [size: 3]. CPU array. Must not be NULL.

  • densityMatrix – [in] Pointer to input density matrix [size: nao × nao]. GPU array. Must not be NULL.

  • outGradient – [out] Pointer to GPU buffer for computed gradients [size: natom × 3]. Must not be NULL. Elements of this array are overwritten by the calculation results.

Returns:

  • CUEST_STATUS_SUCCESS on success;

  • CUEST_STATUS_INVALID_HANDLE if the cuEST handle is NULL;

  • CUEST_STATUS_NULL_POINTER if any required pointer is NULL;

  • CUEST_STATUS_INVALID_TYPE if opaque handles are not the correct type;

  • CUEST_STATUS_INVALID_ARGUMENT if component is not a valid cuestAngularMomentumComputeParametersComponent_t value;

  • CUEST_STATUS_EXCEPTION or CUEST_STATUS_UNKNOWN_ERROR otherwise.

cuestStatus_t cuestAngularMomentumDerivativeComputeWorkspaceQuery(
cuestHandle_t handle,
const cuestOEIntPlan_t plan,
const cuestAngularMomentumDerivativeComputeParameters_t parameters,
cuestWorkspaceDescriptor_t *temporaryWorkspaceDescriptor,
cuestAngularMomentumComputeParametersComponent_t component,
const double *origin,
const double *densityMatrix,
double *outGradient
)#

Query the temporary workspace required (in bytes) for angular momentum matrix gradient calculation.

This function computes the workspace requirements of a single gradient evaluation, based on the supplied plan and dimensions. The workspace descriptor is filled with host and device memory requirements. Output buffers may be NULL for the purposes of the query.

Since

v0.2.0

Parameters:

  • handle – [in] cuEST handle. Must not be NULL.

  • plan – [in] One-electron integral computation plan (opaque handle). Must not be NULL.

  • parameters – [in] Compute parameters (cuestAngularMomentumDerivativeComputeParameters_t). Must not be NULL

  • temporaryWorkspaceDescriptor – [out] Output descriptor for temporary workspace sizes. Must not be NULL.

  • component – [in] Enum (cuestAngularMomentumComputeParametersComponent_t) defining the angular momentum operator (i.e. \(L_x\), \(L_y\), \(L_z\)).

  • origin – [in] Cartesian coordinates defining the operator origin [size: 3]. CPU array. Must not be NULL.

  • densityMatrix – Optional density matrix pointer (may be NULL).

  • outGradient – Optional buffer for computed gradient (may be NULL).

Returns:

  • CUEST_STATUS_SUCCESS on success;

  • CUEST_STATUS_INVALID_HANDLE if the cuEST handle is NULL;

  • CUEST_STATUS_NULL_POINTER if any required pointer is NULL;

  • CUEST_STATUS_INVALID_TYPE if opaque handles are not the correct type;

  • CUEST_STATUS_INVALID_ARGUMENT if component is not a valid cuestAngularMomentumComputeParametersComponent_t value;

  • CUEST_STATUS_EXCEPTION or CUEST_STATUS_UNKNOWN_ERROR otherwise.