8
#include <libciomr/libciomr.h>
11
#include <libpsio/psio.h>
14
** chkpt_rd_contr_full(): Reads in the normalized contraction coefficients.
16
** takes no arguments.
18
** returns: double **contr Normalized contraction coefficients are
19
** returned as a matrix of doubles.
23
double **chkpt_rd_contr_full(void)
25
double **contr, *temp_contr;
26
int nprim, i, j, ij = 0;
28
keyword = chkpt_build_keyword("Contraction coefficients");
30
nprim = chkpt_rd_nprim();
32
temp_contr = init_array(MAXANGMOM*nprim);
33
contr = block_matrix(nprim,MAXANGMOM);
35
psio_read_entry(PSIF_CHKPT, keyword, (char *) temp_contr,
36
MAXANGMOM*nprim*sizeof(double));
38
/* Picking non-zero coefficients to the "master" array contr */
39
for(i=0,ij=0; i < MAXANGMOM; i++)
40
for(j=0; j < nprim; j++, ij++) {
41
contr[j][i] = temp_contr[ij];