Line data    Source code

       1             : /*----------------------------------------------------------------------------*/
       2             : /*  CP2K: A general program to perform molecular dynamics simulations         */
       3             : /*  Copyright 2000-2024 CP2K developers group <https://cp2k.org>              */
       4             : /*                                                                            */
       5             : /*  SPDX-License-Identifier: BSD-3-Clause                                     */
       6             : /*----------------------------------------------------------------------------*/
       7             : 
       8             : #include <assert.h>
       9             : #include <stdlib.h>
      10             : #include <string.h>
      11             : 
      12             : #include "grid_basis_set.h"
      13             : 
      14             : /*******************************************************************************
      15             :  * \brief Allocates a basis set which can be passed to grid_create_task_list.
      16             :  *        See grid_task_list.h for details.
      17             :  * \author Ole Schuett
      18             :  ******************************************************************************/
      19       14401 : void grid_create_basis_set(const int nset, const int nsgf, const int maxco,
      20             :                            const int maxpgf, const int lmin[nset],
      21             :                            const int lmax[nset], const int npgf[nset],
      22             :                            const int nsgf_set[nset], const int first_sgf[nset],
      23       14401 :                            const double sphi[nsgf][maxco],
      24       14401 :                            const double zet[nset][maxpgf],
      25             :                            grid_basis_set **basis_set_out) {
      26             : 
      27       14401 :   grid_basis_set *basis_set = malloc(sizeof(grid_basis_set));
      28       14401 :   assert(basis_set != NULL);
      29             : 
      30       14401 :   basis_set->nset = nset;
      31       14401 :   basis_set->nsgf = nsgf;
      32       14401 :   basis_set->maxco = maxco;
      33       14401 :   basis_set->maxpgf = maxpgf;
      34             : 
      35       14401 :   size_t size = nset * sizeof(int);
      36       14401 :   basis_set->lmin = malloc(size);
      37       14401 :   assert(basis_set->lmin != NULL);
      38       14401 :   memcpy(basis_set->lmin, lmin, size);
      39       14401 :   basis_set->lmax = malloc(size);
      40       14401 :   assert(basis_set->lmax != NULL);
      41       14401 :   memcpy(basis_set->lmax, lmax, size);
      42       14401 :   basis_set->npgf = malloc(size);
      43       14401 :   assert(basis_set->npgf != NULL);
      44       14401 :   memcpy(basis_set->npgf, npgf, size);
      45       14401 :   basis_set->nsgf_set = malloc(size);
      46       14401 :   assert(basis_set->nsgf_set != NULL);
      47       14401 :   memcpy(basis_set->nsgf_set, nsgf_set, size);
      48       14401 :   basis_set->first_sgf = malloc(size);
      49       14401 :   assert(basis_set->first_sgf != NULL);
      50       14401 :   memcpy(basis_set->first_sgf, first_sgf, size);
      51       14401 :   size = nsgf * maxco * sizeof(double);
      52       14401 :   basis_set->sphi = malloc(size);
      53       14401 :   assert(basis_set->sphi != NULL);
      54       14401 :   memcpy(basis_set->sphi, sphi, size);
      55       14401 :   size = nset * maxpgf * sizeof(double);
      56       14401 :   basis_set->zet = malloc(size);
      57       14401 :   assert(basis_set->zet != NULL);
      58       14401 :   memcpy(basis_set->zet, zet, size);
      59             : 
      60       14401 :   *basis_set_out = basis_set;
      61       14401 : }
      62             : 
      63             : /*******************************************************************************
      64             :  * \brief Deallocates given basis set.
      65             :  * \author Ole Schuett
      66             :  ******************************************************************************/
      67       14401 : void grid_free_basis_set(grid_basis_set *basis_set) {
      68       14401 :   free(basis_set->lmin);
      69       14401 :   free(basis_set->lmax);
      70       14401 :   free(basis_set->npgf);
      71       14401 :   free(basis_set->nsgf_set);
      72       14401 :   free(basis_set->first_sgf);
      73       14401 :   free(basis_set->sphi);
      74       14401 :   free(basis_set->zet);
      75       14401 :   free(basis_set);
      76       14401 : }
      77             : 
      78             : // EOF