diff --git a/Core/FormFactors/src/FormFactorPyramid.cpp b/Core/FormFactors/src/FormFactorPyramid.cpp index fbcb76c71132f2a5ae78f88989cd9fcf7167ab75..ed781b61966edde50ae6de98b3ab9d441c55a3f0 100644 --- a/Core/FormFactors/src/FormFactorPyramid.cpp +++ b/Core/FormFactors/src/FormFactorPyramid.cpp @@ -118,10 +118,10 @@ complex_t FormFactorPyramid::fullPyramidPrimitive(complex_t a, complex_t b, comp { const complex_t im(0, 1); complex_t phase = std::exp(im * c * z); - complex_t nominator = std::sin(a * z) * (b * (a * a - b * b + c * c) * std::cos(b * z) + complex_t numerator = std::sin(a * z) * (b * (a * a - b * b + c * c) * std::cos(b * z) + im * c * (a * a + b * b - c * c) * std::sin(b * z)) + a * std::cos(a * z) * ((-a * a + b * b + c * c) * std::sin(b * z) + 2.0 * im * b * c * std::cos(b * z)); complex_t denominator = a * b * (a - b - c) * (a + b - c) * (a - b + c) * (a + b + c); - return -4.0 * phase * nominator / denominator; + return -4.0 * phase * numerator / denominator; }