sum[( 2 +
n ) % 3] = sum1; \
{ \
s21 = s2 - s1; \
if ( s21 == s10 || ( fabs( s2.real() * heps ) >= fabs( s21.real() ) && \
fabs( s2.imag() * heps ) >= fabs( s21.imag() ) ) ) \
break; \
dv = sump; \
sump = s1 + 1. / ( 1. / s21 - 1. / s10 ); \
} \
if ( fabs( sump.real() * teps ) >= fabs( sump.real() - dv.real() ) && \
fabs( sump.imag() * teps ) >= fabs( sump.imag() - dv.imag() ) ) \
break;
std::complex< double > ncomplex
17#define stepWynn( n ) \
18 sum[( 2 + n ) % 3] = sum1; \
19 { \
20 const ncomplex s2 = sum[( 2 + n ) % 3]; \
21 const ncomplex s1 = sum[( 1 + n ) % 3]; \
22 const ncomplex s10 = s21; \
23 s21 = s2 - s1; \
24 if ( s21 == s10 || ( fabs( s2.real() * heps ) >= fabs( s21.real() ) && \
25 fabs( s2.imag() * heps ) >= fabs( s21.imag() ) ) ) \
26 break; \
27 dv = sump; \
28 sump = s1 + 1. / ( 1. / s21 - 1. / s10 ); \
29 } \
30 if ( fabs( sump.real() * teps ) >= fabs( sump.real() - dv.real() ) && \
31 fabs( sump.imag() * teps ) >= fabs( sump.imag() - dv.imag() ) ) \
32 break;