Skip to content

Commit 12758ea

Browse files
lucascolleylucascolley
committed
STY: apply clang-format
1 parent d18e5a5 commit 12758ea

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

64 files changed

+16076
-13420
lines changed

‎include/xsf/airy.h

+1-1
Original file line numberDiff line numberDiff line change
@@ -212,7 +212,7 @@ inline void airyb(double x, double *ai, double *bi, double *ad, double *bd) {
212212
*ad = c1 * df - c2 * dg;
213213
*bd = sr3 * (c1 * df + c2 * dg);
214214
} else {
215-
km = (int) (24.5 - xa);
215+
km = (int)(24.5 - xa);
216216
if (xa < 6.0)
217217
km = 14;
218218
if (xa > 15.0)

‎include/xsf/alg.h

+6-2
Original file line numberDiff line numberDiff line change
@@ -4,8 +4,12 @@
44

55
namespace xsf {
66

7-
XSF_HOST_DEVICE inline double cbrt(double x) { return cephes::detail::cbrt(x); }
7+
XSF_HOST_DEVICE inline double cbrt(double x) {
8+
return cephes::detail::cbrt(x);
9+
}
810

9-
XSF_HOST_DEVICE inline float cbrt(float x) { return cbrt(static_cast<double>(x)); }
11+
XSF_HOST_DEVICE inline float cbrt(float x) {
12+
return cbrt(static_cast<double>(x));
13+
}
1014

1115
} // namespace xsf

‎include/xsf/amos/amos.h

+6,152-5,808
Large diffs are not rendered by default.

‎include/xsf/bessel.h

+81-27
Original file line numberDiff line numberDiff line change
@@ -638,13 +638,21 @@ T cyl_bessel_je(T v, T x) {
638638
return std::real(cyl_bessel_je(v, std::complex(x)));
639639
}
640640

641-
inline double cyl_bessel_y0(double x) { return cephes::y0(x); }
641+
inline double cyl_bessel_y0(double x) {
642+
return cephes::y0(x);
643+
}
642644

643-
inline float cyl_bessel_y0(float x) { return cyl_bessel_y0(static_cast<double>(x)); }
645+
inline float cyl_bessel_y0(float x) {
646+
return cyl_bessel_y0(static_cast<double>(x));
647+
}
644648

645-
inline double cyl_bessel_y1(double x) { return cephes::y1(x); }
649+
inline double cyl_bessel_y1(double x) {
650+
return cephes::y1(x);
651+
}
646652

647-
inline float cyl_bessel_y1(float x) { return cyl_bessel_y1(static_cast<double>(x)); }
653+
inline float cyl_bessel_y1(float x) {
654+
return cyl_bessel_y1(static_cast<double>(x));
655+
}
648656

649657
inline std::complex<double> cyl_bessel_ye(double v, std::complex<double> z) {
650658
int n = 1;
@@ -857,13 +865,21 @@ inline std::complex<float> cyl_hankel_2e(float v, std::complex<float> z) {
857865
);
858866
}
859867

860-
inline double cyl_bessel_j0(double x) { return cephes::j0(x); }
868+
inline double cyl_bessel_j0(double x) {
869+
return cephes::j0(x);
870+
}
861871

862-
inline float cyl_bessel_j0(float x) { return cyl_bessel_j0(static_cast<double>(x)); }
872+
inline float cyl_bessel_j0(float x) {
873+
return cyl_bessel_j0(static_cast<double>(x));
874+
}
863875

864-
inline double cyl_bessel_j1(double x) { return cephes::j1(x); }
876+
inline double cyl_bessel_j1(double x) {
877+
return cephes::j1(x);
878+
}
865879

866-
inline float cyl_bessel_j1(float x) { return cyl_bessel_j1(static_cast<double>(x)); }
880+
inline float cyl_bessel_j1(float x) {
881+
return cyl_bessel_j1(static_cast<double>(x));
882+
}
867883

868884
inline std::complex<double> cyl_bessel_j(double v, std::complex<double> z) {
869885
int n = 1;
@@ -990,25 +1006,45 @@ T cyl_bessel_y(T v, T x) {
9901006
return std::real(res);
9911007
}
9921008

993-
inline double cyl_bessel_i(double v, double x) { return cephes::iv(v, x); }
1009+
inline double cyl_bessel_i(double v, double x) {
1010+
return cephes::iv(v, x);
1011+
}
9941012

995-
inline float cyl_bessel_i(float v, float x) { return cyl_bessel_i(static_cast<double>(v), static_cast<double>(x)); }
1013+
inline float cyl_bessel_i(float v, float x) {
1014+
return cyl_bessel_i(static_cast<double>(v), static_cast<double>(x));
1015+
}
9961016

997-
inline double cyl_bessel_i0(double x) { return cephes::i0(x); }
1017+
inline double cyl_bessel_i0(double x) {
1018+
return cephes::i0(x);
1019+
}
9981020

999-
inline float cyl_bessel_i0(float x) { return cyl_bessel_i0(static_cast<double>(x)); }
1021+
inline float cyl_bessel_i0(float x) {
1022+
return cyl_bessel_i0(static_cast<double>(x));
1023+
}
10001024

1001-
inline double cyl_bessel_i0e(double x) { return cephes::i0e(x); }
1025+
inline double cyl_bessel_i0e(double x) {
1026+
return cephes::i0e(x);
1027+
}
10021028

1003-
inline float cyl_bessel_i0e(float x) { return cyl_bessel_i0e(static_cast<double>(x)); }
1029+
inline float cyl_bessel_i0e(float x) {
1030+
return cyl_bessel_i0e(static_cast<double>(x));
1031+
}
10041032

1005-
inline double cyl_bessel_i1(double x) { return cephes::i1(x); }
1033+
inline double cyl_bessel_i1(double x) {
1034+
return cephes::i1(x);
1035+
}
10061036

1007-
inline float cyl_bessel_i1(float x) { return cyl_bessel_i1(static_cast<double>(x)); }
1037+
inline float cyl_bessel_i1(float x) {
1038+
return cyl_bessel_i1(static_cast<double>(x));
1039+
}
10081040

1009-
inline double cyl_bessel_i1e(double x) { return cephes::i1e(x); }
1041+
inline double cyl_bessel_i1e(double x) {
1042+
return cephes::i1e(x);
1043+
}
10101044

1011-
inline float cyl_bessel_i1e(float x) { return cyl_bessel_i1e(static_cast<double>(x)); }
1045+
inline float cyl_bessel_i1e(float x) {
1046+
return cyl_bessel_i1e(static_cast<double>(x));
1047+
}
10121048

10131049
inline std::complex<double> cyl_bessel_i(double v, std::complex<double> z) {
10141050
int n = 1;
@@ -1112,21 +1148,37 @@ T cyl_bessel_k(T v, T z) {
11121148
return std::real(cyl_bessel_k(v, std::complex(z)));
11131149
}
11141150

1115-
inline double cyl_bessel_k0(double x) { return cephes::k0(x); }
1151+
inline double cyl_bessel_k0(double x) {
1152+
return cephes::k0(x);
1153+
}
11161154

1117-
inline float cyl_bessel_k0(float x) { return cyl_bessel_k0(static_cast<double>(x)); }
1155+
inline float cyl_bessel_k0(float x) {
1156+
return cyl_bessel_k0(static_cast<double>(x));
1157+
}
11181158

1119-
inline double cyl_bessel_k0e(double x) { return cephes::k0e(x); }
1159+
inline double cyl_bessel_k0e(double x) {
1160+
return cephes::k0e(x);
1161+
}
11201162

1121-
inline float cyl_bessel_k0e(float x) { return cyl_bessel_k0e(static_cast<double>(x)); }
1163+
inline float cyl_bessel_k0e(float x) {
1164+
return cyl_bessel_k0e(static_cast<double>(x));
1165+
}
11221166

1123-
inline double cyl_bessel_k1(double x) { return cephes::k1(x); }
1167+
inline double cyl_bessel_k1(double x) {
1168+
return cephes::k1(x);
1169+
}
11241170

1125-
inline float cyl_bessel_k1(float x) { return cyl_bessel_k1(static_cast<double>(x)); }
1171+
inline float cyl_bessel_k1(float x) {
1172+
return cyl_bessel_k1(static_cast<double>(x));
1173+
}
11261174

1127-
inline double cyl_bessel_k1e(double x) { return cephes::k1e(x); }
1175+
inline double cyl_bessel_k1e(double x) {
1176+
return cephes::k1e(x);
1177+
}
11281178

1129-
inline float cyl_bessel_k1e(float x) { return cyl_bessel_k1e(static_cast<double>(x)); }
1179+
inline float cyl_bessel_k1e(float x) {
1180+
return cyl_bessel_k1e(static_cast<double>(x));
1181+
}
11301182

11311183
inline std::complex<double> cyl_hankel_1(double v, std::complex<double> z) {
11321184
int n = 1;
@@ -1193,7 +1245,9 @@ inline std::complex<float> cyl_hankel_2(float v, std::complex<float> z) {
11931245
return static_cast<std::complex<float>>(cyl_hankel_2(static_cast<double>(v), static_cast<std::complex<double>>(z)));
11941246
}
11951247

1196-
inline double besselpoly(double a, double lambda, double nu) { return cephes::besselpoly(a, lambda, nu); }
1248+
inline double besselpoly(double a, double lambda, double nu) {
1249+
return cephes::besselpoly(a, lambda, nu);
1250+
}
11971251

11981252
inline float besselpoly(float a, float lambda, float nu) {
11991253
return besselpoly(static_cast<double>(a), static_cast<double>(lambda), static_cast<double>(nu));

‎include/xsf/beta.h

+12-4
Original file line numberDiff line numberDiff line change
@@ -4,12 +4,20 @@
44

55
namespace xsf {
66

7-
XSF_HOST_DEVICE inline double beta(double a, double b) { return cephes::beta(a, b); }
7+
XSF_HOST_DEVICE inline double beta(double a, double b) {
8+
return cephes::beta(a, b);
9+
}
810

9-
XSF_HOST_DEVICE inline float beta(float a, float b) { return beta(static_cast<double>(a), static_cast<double>(b)); }
11+
XSF_HOST_DEVICE inline float beta(float a, float b) {
12+
return beta(static_cast<double>(a), static_cast<double>(b));
13+
}
1014

11-
XSF_HOST_DEVICE inline double betaln(double a, double b) { return cephes::lbeta(a, b); }
15+
XSF_HOST_DEVICE inline double betaln(double a, double b) {
16+
return cephes::lbeta(a, b);
17+
}
1218

13-
XSF_HOST_DEVICE inline float betaln(float a, float b) { return betaln(static_cast<double>(a), static_cast<double>(b)); }
19+
XSF_HOST_DEVICE inline float betaln(float a, float b) {
20+
return betaln(static_cast<double>(a), static_cast<double>(b));
21+
}
1422

1523
} // namespace xsf

‎include/xsf/cdflib.h

+8-10
Original file line numberDiff line numberDiff line change
@@ -43,10 +43,10 @@ XSF_HOST_DEVICE inline double gdtrib(double a, double p, double x) {
4343
}
4444
double q = 1.0 - p;
4545
auto func = [a, p, q, x](double b) {
46-
if (p <= q) {
47-
return cephes::igam(b, a * x) - p;
48-
}
49-
return q - cephes::igamc(b, a * x);
46+
if (p <= q) {
47+
return cephes::igam(b, a * x) - p;
48+
}
49+
return q - cephes::igamc(b, a * x);
5050
};
5151
double lower_bound = std::numeric_limits<double>::min();
5252
double upper_bound = std::numeric_limits<double>::max();
@@ -71,9 +71,8 @@ XSF_HOST_DEVICE inline double gdtrib(double a, double p, double x) {
7171
* number of iterations needed in this bracket search to check all normalized
7272
* floating point values.
7373
*/
74-
auto [xl, xr, f_xl, f_xr, bracket_status] = detail::bracket_root_for_cdf_inversion(
75-
func, 1.0, lower_bound, upper_bound, -0.875, 7.0, 0.125, 8, false, 342
76-
);
74+
auto [xl, xr, f_xl, f_xr, bracket_status] =
75+
detail::bracket_root_for_cdf_inversion(func, 1.0, lower_bound, upper_bound, -0.875, 7.0, 0.125, 8, false, 342);
7776
if (bracket_status == 1) {
7877
set_error("gdtrib", SF_ERROR_UNDERFLOW, NULL);
7978
return 0.0;
@@ -86,9 +85,8 @@ XSF_HOST_DEVICE inline double gdtrib(double a, double p, double x) {
8685
set_error("gdtrib", SF_ERROR_OTHER, "Computational Error");
8786
return std::numeric_limits<double>::quiet_NaN();
8887
}
89-
auto [result, root_status] = detail::find_root_chandrupatla(
90-
func, xl, xr, f_xl, f_xr, std::numeric_limits<double>::epsilon(), 1e-100, 100
91-
);
88+
auto [result, root_status] =
89+
detail::find_root_chandrupatla(func, xl, xr, f_xl, f_xr, std::numeric_limits<double>::epsilon(), 1e-100, 100);
9290
if (root_status) {
9391
/* The root finding return should only fail if there's a bug in our code. */
9492
set_error("gdtrib", SF_ERROR_OTHER, "Computational Error, (%.17g, %.17g, %.17g)", a, p, x);

‎include/xsf/cephes/beta.h

+3-3
Original file line numberDiff line numberDiff line change
@@ -161,9 +161,9 @@ namespace cephes {
161161
b = Gamma(b);
162162
if (std::isinf(y)) {
163163
goto overflow;
164-
}
164+
}
165165

166-
if (std::abs(std::abs(a*y) - 1.0) > std::abs(std::abs(b*y) - 1.0)) {
166+
if (std::abs(std::abs(a * y) - 1.0) > std::abs(std::abs(b * y) - 1.0)) {
167167
y = b * y;
168168
y *= a;
169169
} else {
@@ -239,7 +239,7 @@ namespace cephes {
239239
return (sign * std::numeric_limits<double>::infinity());
240240
}
241241

242-
if (std::abs(std::abs(a*y) - 1.0) > std::abs(std::abs(b*y) - 1.0)) {
242+
if (std::abs(std::abs(a * y) - 1.0) > std::abs(std::abs(b * y) - 1.0)) {
243243
y = b * y;
244244
y *= a;
245245
} else {

‎include/xsf/cephes/chdtr.h

+1-1
Original file line numberDiff line numberDiff line change
@@ -165,7 +165,7 @@ namespace cephes {
165165
if (x < 0.0) {
166166
set_error("chdtr", SF_ERROR_DOMAIN, NULL);
167167
return (std::numeric_limits<double>::quiet_NaN());
168-
}
168+
}
169169
return (igamc(df / 2.0, x / 2.0));
170170
}
171171

0 commit comments

Comments
 (0)