From a0da08b637212d34924f4ab1ac2003ee4ffa6dd7 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= Date: Thu, 13 Sep 2018 05:32:07 +0300 Subject: [PATCH] Fix normalisation of test sigma0 factor. Ewald sum working now! Former-commit-id: e65a370e83c68c6c3553f13312b9280a3413b200 --- tests/ewalds.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/tests/ewalds.c b/tests/ewalds.c index bfd3e2e..76672d9 100644 --- a/tests/ewalds.c +++ b/tests/ewalds.c @@ -213,7 +213,7 @@ ewaldtest_triang_results *ewaldtest_triang(const ewaldtest_triang_params p) { // Now calculate the reference values [LT(4.16)] results->regsigmas_416 = calloc(nelem_sc, sizeof(complex double)); - results->regsigmas_416[0] = -1/M_SQRTPI; + results->regsigmas_416[0] = -2 * c->legendre0[gsl_sf_legendre_array_index(0,0)]; { double legendres[gsl_sf_legendre_array_n(p.lMax)];