Fix relrod calculation

("k" was renamed to "wavenumber")

1 files changed, 2 insertions, 2 deletions

diff --git a/src/relrod.c b/src/relrod.c
index c55432ec..ec8ad579 100644
--- a/src/relrod.c
+++ b/src/relrod.c
@@ -109,8 +109,8 @@ void get_reflections(struct image *image, UnitCell *cell)
 		/* Next, solve the relrod equation to calculate
 		 * the excitation error */
 		a = 1.0;
-		b = 2.0*(gn - wavenumber);
-		c = -2.0*gn*k + g_sq;
+		b = 2.0*(wavenumber + gn);
+		c = -2.0*gn*wavenumber + g_sq;
 		t = -0.5*(b + sign(b)*sqrt(b*b - 4.0*a*c));
 		s1 = t/a;
 		s2 = c/t;