8#ifndef GRID_DGEMM_COEFFICIENTS_H
9#define GRID_DGEMM_COEFFICIENTS_H
13#include "../common/grid_common.h"
18 const double rp[3],
const int *lmax,
22 const int *
const lmax,
const int lp,
23 const double prefactor,
29 const tensor *coef,
const int *coef_offset_,
const int *lmin,
30 const int *lmax,
const double *ra,
const double *rb,
const double *rab,
31 const double prefactor,
tensor *co);
41 double *
const coef_xyz);
43 double *
const coef_xyz);
46 const double *__restrict coef_ijk,
47 double *__restrict coef_xyz);
51 const double prefactor,
void grid_transform_coef_ijk_to_xyz_cp2k(const int lp, const double dh[3][3], const double *__restrict coef_ijk, double *__restrict coef_xyz)
void grid_transform_coef_xyz_to_ijk(const double dh[3][3], const tensor *coef_xyz)
void grid_compute_coefficients_dgemm(const int *const lmin, const int *const lmax, const int lp, const double prefactor, const tensor *const alpha, const tensor *const pab, tensor *coef_xyz)
void compute_compact_polynomial_coefficients(const tensor *coef, const int *coef_offset_, const int *lmin, const int *lmax, const double *ra, const double *rb, const double *rab, const double prefactor, tensor *co)
void transform_xyz_to_triangular(const tensor *const coef, double *const coef_xyz)
void transform_triangular_to_xyz(const double *const coef_xyz, tensor *const coef)
void transform_yxz_to_triangular(const tensor *const coef, double *const coef_xyz)
void grid_transform_coef_jik_to_yxz(const double dh[3][3], const tensor *coef_xyz)
void grid_compute_vab(const int *const lmin, const int *const lmax, const int lp, const double prefactor, const tensor *const alpha, const tensor *const coef_xyz, tensor *pab)
void grid_prepare_alpha_dgemm(const double ra[3], const double rb[3], const double rp[3], const int *lmax, tensor *alpha)