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

vswf computation (start writing) 

Former-commit-id: 8886648a52b94a21a17e2fede87941a65597d740
This commit is contained in:
Marek Nečada 2017-12-17 17:46:37 +02:00
parent bb289a46b9
commit ef8129d93c
1 changed files with 45 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

45
qpms/vswf_xu.c Normal file
View File

@ -0,0 +1,45 @@
#include "translations.h"
#include <math.h>
#include <stdlib.h> //abort();
/*
* References:
* [1] Yu-Lin Xu, Journal of Computational Physics 127, 285–298 (1996)
*/
/*
* The following value delimits the interval for which the computation
* of the pi and tau functions [1,(37)] actually takes place.
* For x in [-1, -1 + PITAU_0THRESHOLD] the x->-1 limit is taken,
* for x in [-1 + PITAU_0THRESHOLD, 1 - PITAU_0THRESHOLD] value for x is calculated,
* for x in [1 - PITAU_THRESHOLD, 1] the x->1 limit is taken.
*
* low-priority TODO: take more than 0th order expansion at x->±1
*/
#define PITAU_0THRESHOLD 1e-7
int taumncos_Xu_zerolim_fill(int lmax, double *where);
double * taumncos_Xu_zerolim_get(int lmax);
int taumncos_Xu_pilim_fill(int lmax, double *where);
double * taumncos_Xu_pilim_get(int lmax);
int pimncos_Xu_zerolim_fill(int lmax, double *where);
double * pimncos_Xu_zerolim_get(int lmax);
int pimncos_Xu_pilim_fill(int lmax, double *where);
double * pimncos_Xu_pilim_get(int lmax);
// [1] (37)
static inline double pimncos_Xu(double theta) {
abort();
}
static inline double taumncos_Xu(double theta) {
abort();
}
// [1] (36)
complex double qpms_vswf_single_mg_Xu(int m, int n, sph_t kdlj,
qpms_bessel_t btyp)
{
}