2018-01-18 13:22:27 +02:00
|
|
|
#ifndef BESSELS_H
|
|
|
|
|
#define BESSELS_H
|
|
|
|
|
|
|
|
|
|
#include <stddef.h>
|
|
|
|
|
#include <complex.h>
|
|
|
|
|
|
|
|
|
|
complex double *hankelcoefftable_init(size_t maxn);
|
|
|
|
|
|
|
|
|
|
// general, gives the offset such that result[k] is
|
|
|
|
|
// the coefficient corresponding to the e**(I * x) * x**(-k-1)
|
|
|
|
|
// term of the Hankel function; no boundary checks!
|
|
|
|
|
static inline complex double *
|
|
|
|
|
hankelcoeffs_get(complex double *hankelcoefftable, size_t n){
|
2018-01-18 14:55:15 +02:00
|
|
|
return hankelcoefftable + n*(n+1)/2;
|
2018-01-18 13:22:27 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
2018-01-18 14:55:15 +02:00
|
|
|
// general; target_longrange and target_shortrange are of size (maxn+1)
|
|
|
|
|
void hankelparts_fill(complex double *target_longrange, complex double *target_shortrange,
|
|
|
|
|
size_t maxn, size_t longrange_k_cutoff, // x**(-(k+1)-1) terms go completely to short-range part
|
|
|
|
|
complex double *hankelcoefftable,
|
|
|
|
|
unsigned kappa, double c, double x);
|
|
|
|
|
|
2018-01-18 13:22:27 +02:00
|
|
|
#endif //BESSELS_H
|
