@@ -536,17 +536,20 @@ static @ctype@ _npy_scaled_cexp@c@(@type@ x, @type@ y, npy_int expt)
536536 return npy_cpack @c @(npy_cosh @c @(x ), x * y );
537537 absx = npy_fabs @c @(x );
538538 if (absx < CCOSH_BIG ) /* small x: normal case */
539- return npy_cpack @c @(npy_cosh @c @(x ) * npy_cos @c @(y ), npy_sinh @c @(x ) * npy_sin @c @(y ));
539+ return npy_cpack @c @(npy_cosh @c @(x ) * npy_cos @c @(y ),
540+ npy_sinh @c @(x ) * npy_sin @c @(y ));
540541
541542 /* |x| >= 22, so cosh(x) ~= exp(|x|) */
542543 if (absx < SCALED_CEXP_LOWER @C @) {
543544 /* x < 710: exp(|x|) won't overflow */
544545 h = npy_exp @c @(absx ) * 0.5 ;
545- return npy_cpack @c @(h * npy_cos @c @(y ), npy_copysign @c @(h , x ) * npy_sin @c @(y ));
546+ return npy_cpack @c @(h * npy_cos @c @(y ),
547+ npy_copysign @c @(h , x ) * npy_sin @c @(y ));
546548 } else if (absx < SCALED_CEXP_UPPER @C @) {
547549 /* x < 1455: scale to avoid overflow */
548550 z = _npy_scaled_cexp @c @(absx , y , -1 );
549- return npy_cpack @c @(npy_creal @c @(z ), npy_cimag @c @(z ) * npy_copysign @c @(1 , x ));
551+ return npy_cpack @c @(npy_creal @c @(z ),
552+ npy_cimag @c @(z ) * npy_copysign @c @(1 , x ));
550553 } else {
551554 /* x >= 1455: the result always overflows */
552555 h = CCOSH_HUGE * x ;
@@ -661,17 +664,20 @@ static @ctype@ _npy_scaled_cexp@c@(@type@ x, @type@ y, npy_int expt)
661664 return npy_cpack @c @(npy_sinh @c @(x ), y );
662665 absx = npy_fabs @c @(x );
663666 if (absx < CSINH_BIG ) /* small x: normal case */
664- return npy_cpack @c @(npy_sinh @c @(x ) * npy_cos @c @(y ), npy_cosh @c @(x ) * npy_sin @c @(y ));
667+ return npy_cpack @c @(npy_sinh @c @(x ) * npy_cos @c @(y ),
668+ npy_cosh @c @(x ) * npy_sin @c @(y ));
665669
666670 /* |x| >= 22, so cosh(x) ~= exp(|x|) */
667671 if (absx < SCALED_CEXP_LOWER @C @) {
668672 /* x < 710: exp(|x|) won't overflow */
669673 h = npy_exp @c @(npy_fabs @c @(x )) * 0.5 ;
670- return npy_cpack @c @(npy_copysign @c @(h , x ) * npy_cos @c @(y ), h * npy_sin @c @(y ));
674+ return npy_cpack @c @(npy_copysign @c @(h , x ) * npy_cos @c @(y ),
675+ h * npy_sin @c @(y ));
671676 } else if (x < SCALED_CEXP_UPPER @C @) {
672677 /* x < 1455: scale to avoid overflow */
673678 z = _npy_scaled_cexp @c @(absx , y , -1 );
674- return npy_cpack @c @(npy_creal @c @(z ) * npy_copysign @c @(1 , x ), npy_cimag @c @(z ));
679+ return npy_cpack @c @(npy_creal @c @(z ) * npy_copysign @c @(1 , x ),
680+ npy_cimag @c @(z ));
675681 } else {
676682 /* x >= 1455: the result always overflows */
677683 h = CSINH_HUGE * x ;
@@ -1058,7 +1064,8 @@ static inline void _do_hard_work@c@(@type@ x, @type@ y, @type@ *rx,
10581064/*
10591065 * Optimized version of clog() for |z| finite and larger than ~RECIP_EPSILON.
10601066 */
1061- static inline void _clog_for_large_values @c @(@type @ x , @type @ y , @type @ * rr , @type @ * ri )
1067+ static inline void _clog_for_large_values @c @(@type @ x , @type @ y ,
1068+ @type @ * rr , @type @ * ri )
10621069{
10631070#if @precision @ == 1
10641071 const npy_float QUARTER_SQRT_MAX = 0x1p61f ;
@@ -1387,7 +1394,7 @@ static inline npy_float _real_part_reciprocalf(npy_float x, npy_float y)
13871394#endif
13881395#if @precision @ == 2
13891396#define BIAS (DBL_MAX_EXP - 1)
1390- /* XXX more guard digits are useful iff there is extra precision. */
1397+ /* more guard digits are useful iff there is extra precision. */
13911398#define CUTOFF (DBL_MANT_DIG / 2 + 1) /* just half or 1 guard digit */
13921399static inline npy_double _real_part_reciprocal (npy_double x , npy_double y )
13931400{
@@ -1421,7 +1428,8 @@ static inline npy_double _real_part_reciprocal(npy_double x, npy_double y)
14211428#if @precision @ == 3
14221429#define BIAS (LDBL_MAX_EXP - 1)
14231430#define CUTOFF (LDBL_MANT_DIG / 2 + 1)
1424- static inline npy_longdouble _real_part_reciprocall (npy_longdouble x , npy_longdouble y )
1431+ static inline npy_longdouble _real_part_reciprocall (npy_longdouble x ,
1432+ npy_longdouble y )
14251433{
14261434 npy_longdouble scale ;
14271435 union IEEEl2bitsrep ux , uy , us ;
0 commit comments