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

Dudom 

Former-commit-id: ce1f24b4b792b797985113d2c76f2191847c91ac
This commit is contained in:
Marek Nečada 2018-09-02 06:19:50 +03:00
parent 6c09121a5d
commit 83253bfc0c
3 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
qpms/qpms_types.h
View File

@ -19,7 +19,8 @@ typedef size_t qpms_y_t;
typedef enum { typedef enum {
QPMS_SUCCESS = 0, QPMS_SUCCESS = 0,
QPMS_ERROR = 1 QPMS_ERROR = 1,
QPMS_ENOMEM = 8
} qpms_errno_t; } qpms_errno_t;
// Normalisations // Normalisations

4
qpms/translations.c
View File

@ -495,8 +495,8 @@ void qpms_trans_calculator_free(qpms_trans_calculator *c) {
free(c->B_multipliers); free(c->B_multipliers);
#ifdef LATTICESUMS_OLD #ifdef LATTICESUMS_OLD
free(c->hct); free(c->hct);
free(c->legendre0);
#endif #endif
free(c->legendre0);
free(c); free(c);
} }
@ -907,10 +907,10 @@ qpms_trans_calculator
free(qmaxes); free(qmaxes);
#ifdef LATTICESUMS_OLD #ifdef LATTICESUMS_OLD
c->hct = hankelcoefftable_init(2*lMax+1); c->hct = hankelcoefftable_init(2*lMax+1);
#endif
c->legendre0 = malloc(gsl_sf_legendre_array_n(2*lMax+1) * sizeof(double)); c->legendre0 = malloc(gsl_sf_legendre_array_n(2*lMax+1) * sizeof(double));
if (gsl_sf_legendre_array_e(GSL_SF_LEGENDRE_NONE,2*lMax+1, if (gsl_sf_legendre_array_e(GSL_SF_LEGENDRE_NONE,2*lMax+1,
0,-1,c->legendre0)) abort(); // TODO maybe use some "precise" analytical formula instead? 0,-1,c->legendre0)) abort(); // TODO maybe use some "precise" analytical formula instead?
#endif
return c; return c;
} }

2
qpms/translations.h
View File

@ -47,8 +47,8 @@ typedef struct qpms_trans_calculator {
#endif #endif
#ifdef LATTICESUMS_OLD #ifdef LATTICESUMS_OLD
complex double *hct; // Hankel function coefficient table complex double *hct; // Hankel function coefficient table
double *legendre0; // Zero-argument Legendre functions – this might go outside #ifdef in the end...
#endif #endif
double *legendre0; // Zero-argument Legendre functions – this might go outside #ifdef in the end...
} qpms_trans_calculator; } qpms_trans_calculator;