From b68c394819cddfd1f1e5823cce642720c8408bcd Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= Date: Tue, 23 Jan 2018 09:36:57 +0200 Subject: [PATCH] [dipdip dirty] fk5q1n3 Former-commit-id: 52c789b790fb1b20e8ae199287e9b2e87816962e --- dipdip-dirty/lrhankel_recspace_dirty.c | 17 +++++++++++++++-- 1 file changed, 15 insertions(+), 2 deletions(-) diff --git a/dipdip-dirty/lrhankel_recspace_dirty.c b/dipdip-dirty/lrhankel_recspace_dirty.c index 829151f..72e85ce 100644 --- a/dipdip-dirty/lrhankel_recspace_dirty.c +++ b/dipdip-dirty/lrhankel_recspace_dirty.c @@ -42,7 +42,7 @@ complex double fk5q1n3l(double c, double k0, double k, const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) { double kk3 = 3*k*k; return ( - - d[0]*(kk3+4*a[0]*a[0]) + - FF*d[0]*(kk3+4*a[0]*a[0]) + 5*d[1]*(kk3+4*a[1]*a[1]) - 10*d[2]*(kk3+4*a[2]*a[2]) + 10*d[3]*(kk3+4*a[3]*a[3]) @@ -115,7 +115,20 @@ complex double fk5q1n2s(double c, double k0, double k, )/k0; } -lrhankelspec fk5q2n0s = fk5q2n0, fk5q2n0l = fk5q2n0; +complex double fk5q1n3s(double c, double k0, double k, + const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) { + double kk3 = 3*k*k; + return ( + - d[0]*(kk3+4*a[0]*a[0]) + + 5*d[1]*(kk3+4*a[1]*a[1]) + - 10*d[2]*(kk3+4*a[2]*a[2]) + + 10*d[3]*(kk3+4*a[3]*a[3]) + - 5*d[4]*(kk3+4*a[4]*a[4]) + + d[5]*(kk3+4*a[5]*a[5]) + )/(k0*k*k*k); +} + +const lrhankelspec fk5q2n0s = fk5q2n0, fk5q2n0l = fk5q2n0; complex double fk5q2n1s(double c, double k0, double k, const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {