diff --git a/scrambling.cpp b/scrambling.cpp index b0dc18f..4930cd7 100644 --- a/scrambling.cpp +++ b/scrambling.cpp @@ -87,12 +87,12 @@ std::complex crossRatio (std::complex x1, std::complex x double Fitzpatrick(std::complex z, std::complex zbar, double phase) { // phase is necessary to take into account nontrivial monodromy - std::complex c1(1, 0); // one as a complex number + std::complex one = 1; // one as a complex number std::complex i(0,1); // imaginary unit double alphaExpr = 0.5-alpha/2; std::complex exponent1 = exp(phase*i*alphaExpr); std::complex exponent2 = exp(phase*i*alpha); return real(exponent1 * pow(z, alphaExpr) * pow(zbar, alphaExpr) - * (c1 - exponent2 * pow(z, alpha)) * (c1 - pow(zbar, alpha)) - / ( alpha*alpha * (c1-z) * (c1-zbar) ) ); + * (one - exponent2 * pow(z, alpha)) * (one - pow(zbar, alpha)) + / ( alpha*alpha * (one-z) * (one-zbar) ) ); }