From cc417889e17b88b48f41368e41d69c067f04b525 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= Date: Thu, 18 Jan 2018 12:56:04 +0000 Subject: [PATCH] [dirty dipoles] removing unused functions Former-commit-id: f1203d5250a8934663fdd2600b89f9170127dbe3 --- dipdip-dirty/bessels.c | 31 ------------------------------- 1 file changed, 31 deletions(-) diff --git a/dipdip-dirty/bessels.c b/dipdip-dirty/bessels.c index 80df3ff..d1d08de 100644 --- a/dipdip-dirty/bessels.c +++ b/dipdip-dirty/bessels.c @@ -38,37 +38,6 @@ complex double * hankelcoefftable_init(size_t maxn) { return hct; } -void hankelLR2_fill(complex double *target, size_t maxn, complex double *hct, - unsigned kappa, double c, double x) { - memset(target, 0, (maxn+1)*sizeof(complex double)); - double regularisator = pow(1. - exp(-c * x), (double) kappa); - complex double expix = cexp(I * x); - double xfrac = 1.; // x ** (-1-k) - for (size_t k = 0; k < 2; ++k) { - xfrac /= x; - for(size_t n = k; n <= maxn; ++n) - target[n] += regularisator * xfrac * hankelcoeffs_get(hct,n)[k]; - } - for(size_t n = 0; n <= maxn; ++n) - target[n] *= expix; -} - -void hankelSR2_fill(complex double *target, size_t maxn, complex double *hct, - unsigned kappa, double c, double x) { - memset(target, 0, (maxn+1)*sizeof(complex double)); - double antiregularisator = 1 - pow(1. - exp(-c * x), (double) kappa); - complex double expix = cexp(I * x); - double xfrac = 1.; // x ** (-1-k) - for (size_t k = 0; k < maxn; ++k) { - xfrac /= x; - for(size_t n = k; n <= maxn; ++n) - target[n] += ((k<2) ? antiregularisator : 1) - * xfrac * hankelcoeffs_get(hct,n)[k]; - } - for(size_t n = 0; n <= maxn; ++n) - target[n] *= expix; -} - void hankelparts_fill(complex double *lrt, complex double *srt, size_t maxn, size_t lrk_cutoff, complex double *hct, unsigned kappa, double c, double x) {