QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

translations.c return accidentally deleted function, 

remove more legacy stuff


Former-commit-id: 6bb2033234a893e85a0242c10156a700bd2d3ceb
This commit is contained in:
Marek Nečada 2019-08-19 09:29:14 +03:00
parent 8b9c5a794c
commit 9ea487df84
2 changed files with 12 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
qpms/translations.c
View File

@ -290,6 +290,18 @@ complex double qpms_trans_single_B(qpms_normalisation_t norm,
return presum * sum; return presum * sum;
} }
void qpms_trans_calculator_free(qpms_trans_calculator *c) {
free(c->A_multipliers[0]);
free(c->A_multipliers);
free(c->B_multipliers[0]);
free(c->B_multipliers);
#ifdef LATTICESUMS
qpms_ewald3_constants_free(e3c);
#endif
free(c->legendre0);
free(c);
}
static inline size_t qpms_trans_calculator_index_mnmunu(const qpms_trans_calculator *c, static inline size_t qpms_trans_calculator_index_mnmunu(const qpms_trans_calculator *c,
int m, int n, int mu, int nu){ int m, int n, int mu, int nu){
return c->nelem * qpms_mn2y(m,n) + qpms_mn2y(mu,nu); return c->nelem * qpms_mn2y(m,n) + qpms_mn2y(mu,nu);
@ -380,8 +392,6 @@ void qpms_trans_calculator_multipliers_B(
// This is according to the Cruzan-type formula [Xu](59) // This is according to the Cruzan-type formula [Xu](59)
assert(Qmax == gauntB_Q_max(-m,n,mu,nu)); assert(Qmax == gauntB_Q_max(-m,n,mu,nu));
double normlogfac= qpms_trans_normlogfac(norm,m,n,mu,nu); double normlogfac= qpms_trans_normlogfac(norm,m,n,mu,nu);
double normfac = qpms_trans_normfac(norm,m,n,mu,nu); double normfac = qpms_trans_normfac(norm,m,n,mu,nu);
/// N<-M type coefficients w.r.t. Kristensson's convention. Csphase has been already taken into acct ^^^. /// N<-M type coefficients w.r.t. Kristensson's convention. Csphase has been already taken into acct ^^^.
@ -495,9 +505,6 @@ qpms_trans_calculator
} }
free(qmaxes); free(qmaxes);
#ifdef LATTICESUMS_OLD
c->hct = hankelcoefftable_init(2*lMax+1);
#endif
#ifdef LATTICESUMS32 #ifdef LATTICESUMS32
c->e3c = qpms_ewald3_constants_init(2 * lMax + 1, /*csphase*/ qpms_normalisation_t_csphase(normalisation)); c->e3c = qpms_ewald3_constants_init(2 * lMax + 1, /*csphase*/ qpms_normalisation_t_csphase(normalisation));
#endif #endif

7
qpms/translations.h
View File

@ -32,10 +32,6 @@
#include <stdbool.h> #include <stdbool.h>
#include <stddef.h> #include <stddef.h>
#ifdef LATTICESUMS_OLD
#include "bessels.h"
#endif
#if defined LATTICESUMS32 || defined LATTICESUMS31 #if defined LATTICESUMS32 || defined LATTICESUMS31
#include "ewald.h" #include "ewald.h"
#endif #endif
@ -83,9 +79,6 @@ typedef struct qpms_trans_calculator {
#if defined LATTICESUMS32 || defined LATTICESUMS31 #if defined LATTICESUMS32 || defined LATTICESUMS31
qpms_ewald3_constants_t *e3c; qpms_ewald3_constants_t *e3c;
#endif
#ifdef LATTICESUMS_OLD
complex double *hct; // Hankel function coefficient table
#endif #endif
double *legendre0; // Zero-argument Legendre functions – this might go outside #ifdef in the end... double *legendre0; // Zero-argument Legendre functions – this might go outside #ifdef in the end...
} qpms_trans_calculator; } qpms_trans_calculator;