From 1cd4733e851fcd405a6175163d6a267304c75f5b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= <marek@necada.org>
Date: Thu, 6 Sep 2018 11:52:54 +0300
Subject: [PATCH] Fix ewald32_sigma0()

Former-commit-id: 969764c516c0e0e8cc928e18447eb5327791d5bc
---
 qpms/ewald.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/qpms/ewald.c b/qpms/ewald.c
index 874701d..52f9a43 100644
--- a/qpms/ewald.c
+++ b/qpms/ewald.c
@@ -139,7 +139,7 @@ int ewald32_sigma0(complex double *result, double *err,
     const double eta, const double k)
 {
   qpms_csf_result gam;
-  int retval = complex_gamma_inc_e(-0.5, -k/sq(2*eta), &gam);
+  int retval = complex_gamma_inc_e(-0.5, -sq(k/(2*eta)), &gam);
   if (0 != retval)
     abort();
   *result = gam.val * c->legendre0[gsl_sf_legendre_array_index(0,0)] / 2 / M_SQRTPI;