diff options
| author | Kacper <kacper@mail.openlinux.dev> | 2025-12-25 19:24:38 +0100 |
|---|---|---|
| committer | Kacper <kacper@mail.openlinux.dev> | 2025-12-25 20:35:03 +0100 |
| commit | a984eb367c032dbe2577f01238c3d1268526be70 (patch) | |
| tree | 437fef40379b2758b129ccea39df3570fa2d145e /lib/libm/j0.c | |
| parent | 8834571b202cf4dc9c649cfb096c213b6ecf1566 (diff) | |
Clang-tidy fixes
Diffstat (limited to 'lib/libm/j0.c')
| -rw-r--r-- | lib/libm/j0.c | 45 |
1 files changed, 20 insertions, 25 deletions
diff --git a/lib/libm/j0.c b/lib/libm/j0.c index 52d3f483..b466c1ad 100644 --- a/lib/libm/j0.c +++ b/lib/libm/j0.c @@ -63,7 +63,7 @@ static double pzero(double), qzero(double); static const double invsqrtpi = 5.64189583547756279280e-01, /* 0x3FE20DD7, 0x50429B6D */ - tpi = 6.36619772367581382433e-01; /* 0x3FE45F30, 0x6DC9C883 */ + tpi = 6.36619772367581382433e-01; /* 0x3FE45F30, 0x6DC9C883 */ /* common method when |x|>=2 */ static double common(uint32_t ix, double x, int y0) @@ -103,13 +103,13 @@ static double common(uint32_t ix, double x, int y0) /* R0/S0 on [0, 2.00] */ static const double R02 = 1.56249999999999947958e-02, /* 0x3F8FFFFF, 0xFFFFFFFD */ - R03 = -1.89979294238854721751e-04, /* 0xBF28E6A5, 0xB61AC6E9 */ - R04 = 1.82954049532700665670e-06, /* 0x3EBEB1D1, 0x0C503919 */ - R05 = -4.61832688532103189199e-09, /* 0xBE33D5E7, 0x73D63FCE */ - S01 = 1.56191029464890010492e-02, /* 0x3F8FFCE8, 0x82C8C2A4 */ - S02 = 1.16926784663337450260e-04, /* 0x3F1EA6D2, 0xDD57DBF4 */ - S03 = 5.13546550207318111446e-07, /* 0x3EA13B54, 0xCE84D5A9 */ - S04 = 1.16614003333790000205e-09; /* 0x3E1408BC, 0xF4745D8F */ + R03 = -1.89979294238854721751e-04, /* 0xBF28E6A5, 0xB61AC6E9 */ + R04 = 1.82954049532700665670e-06, /* 0x3EBEB1D1, 0x0C503919 */ + R05 = -4.61832688532103189199e-09, /* 0xBE33D5E7, 0x73D63FCE */ + S01 = 1.56191029464890010492e-02, /* 0x3F8FFCE8, 0x82C8C2A4 */ + S02 = 1.16926784663337450260e-04, /* 0x3F1EA6D2, 0xDD57DBF4 */ + S03 = 5.13546550207318111446e-07, /* 0x3EA13B54, 0xCE84D5A9 */ + S04 = 1.16614003333790000205e-09; /* 0x3E1408BC, 0xF4745D8F */ double j0(double x) { @@ -148,16 +148,16 @@ double j0(double x) static const double u00 = -7.38042951086872317523e-02, /* 0xBFB2E4D6, 0x99CBD01F */ - u01 = 1.76666452509181115538e-01, /* 0x3FC69D01, 0x9DE9E3FC */ - u02 = -1.38185671945596898896e-02, /* 0xBF8C4CE8, 0xB16CFA97 */ - u03 = 3.47453432093683650238e-04, /* 0x3F36C54D, 0x20B29B6B */ - u04 = -3.81407053724364161125e-06, /* 0xBECFFEA7, 0x73D25CAD */ - u05 = 1.95590137035022920206e-08, /* 0x3E550057, 0x3B4EABD4 */ - u06 = -3.98205194132103398453e-11, /* 0xBDC5E43D, 0x693FB3C8 */ - v01 = 1.27304834834123699328e-02, /* 0x3F8A1270, 0x91C9C71A */ - v02 = 7.60068627350353253702e-05, /* 0x3F13ECBB, 0xF578C6C1 */ - v03 = 2.59150851840457805467e-07, /* 0x3E91642D, 0x7FF202FD */ - v04 = 4.41110311332675467403e-10; /* 0x3DFE5018, 0x3BD6D9EF */ + u01 = 1.76666452509181115538e-01, /* 0x3FC69D01, 0x9DE9E3FC */ + u02 = -1.38185671945596898896e-02, /* 0xBF8C4CE8, 0xB16CFA97 */ + u03 = 3.47453432093683650238e-04, /* 0x3F36C54D, 0x20B29B6B */ + u04 = -3.81407053724364161125e-06, /* 0xBECFFEA7, 0x73D25CAD */ + u05 = 1.95590137035022920206e-08, /* 0x3E550057, 0x3B4EABD4 */ + u06 = -3.98205194132103398453e-11, /* 0xBDC5E43D, 0x693FB3C8 */ + v01 = 1.27304834834123699328e-02, /* 0x3F8A1270, 0x91C9C71A */ + v02 = 7.60068627350353253702e-05, /* 0x3F13ECBB, 0xF578C6C1 */ + v03 = 2.59150851840457805467e-07, /* 0x3E91642D, 0x7FF202FD */ + v04 = 4.41110311332675467403e-10; /* 0x3DFE5018, 0x3BD6D9EF */ double y0(double x) { @@ -183,10 +183,7 @@ double y0(double x) if (ix >= 0x3e400000) { /* x >= 2**-27 */ /* large ulp error near the first zero, x ~= 0.89 */ z = x * x; - u = u00 + - z * (u01 + - z * (u02 + - z * (u03 + z * (u04 + z * (u05 + z * u06))))); + u = u00 + z * (u01 + z * (u02 + z * (u03 + z * (u04 + z * (u05 + z * u06))))); v = 1.0 + z * (v01 + z * (v02 + z * (v03 + z * v04))); return u / v + tpi * (j0(x) * log(x)); } @@ -401,8 +398,6 @@ static double qzero(double x) } z = 1.0 / (x * x); r = p[0] + z * (p[1] + z * (p[2] + z * (p[3] + z * (p[4] + z * p[5])))); - s = 1.0 + - z * (q[0] + - z * (q[1] + z * (q[2] + z * (q[3] + z * (q[4] + z * q[5]))))); + s = 1.0 + z * (q[0] + z * (q[1] + z * (q[2] + z * (q[3] + z * (q[4] + z * q[5]))))); return (-.125 + r / s) / x; } |