Skip to content

Commit a320224

Browse files
lucascolleylucascolley
committed
run format
1 parent e16e499 commit a320224

37 files changed

+256
-751
lines changed

‎include/xsf/alg.h

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

55
namespace xsf {
66

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

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

1511
} // namespace xsf

‎include/xsf/bessel.h

+27-81
Original file line numberDiff line numberDiff line change
@@ -638,21 +638,13 @@ 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) {
642-
return cephes::y0(x);
643-
}
641+
inline double cyl_bessel_y0(double x) { return cephes::y0(x); }
644642

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

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

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

657649
inline std::complex<double> cyl_bessel_ye(double v, std::complex<double> z) {
658650
int n = 1;
@@ -865,21 +857,13 @@ inline std::complex<float> cyl_hankel_2e(float v, std::complex<float> z) {
865857
);
866858
}
867859

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

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

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

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

884868
inline std::complex<double> cyl_bessel_j(double v, std::complex<double> z) {
885869
int n = 1;
@@ -1006,45 +990,25 @@ T cyl_bessel_y(T v, T x) {
1006990
return std::real(res);
1007991
}
1008992

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

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

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

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

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

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

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

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

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

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

10491013
inline std::complex<double> cyl_bessel_i(double v, std::complex<double> z) {
10501014
int n = 1;
@@ -1148,37 +1112,21 @@ T cyl_bessel_k(T v, T z) {
11481112
return std::real(cyl_bessel_k(v, std::complex(z)));
11491113
}
11501114

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

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

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

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

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

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

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

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

11831131
inline std::complex<double> cyl_hankel_1(double v, std::complex<double> z) {
11841132
int n = 1;
@@ -1245,9 +1193,7 @@ inline std::complex<float> cyl_hankel_2(float v, std::complex<float> z) {
12451193
return static_cast<std::complex<float>>(cyl_hankel_2(static_cast<double>(v), static_cast<std::complex<double>>(z)));
12461194
}
12471195

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

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

‎include/xsf/beta.h

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

55
namespace xsf {
66

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

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

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

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

2315
} // namespace xsf

‎include/xsf/binom.h

+1-3
Original file line numberDiff line numberDiff line change
@@ -82,8 +82,6 @@ XSF_HOST_DEVICE inline double binom(double n, double k) {
8282
return 1 / (n + 1) / cephes::beta(1 + n - k, 1 + k);
8383
}
8484

85-
XSF_HOST_DEVICE inline float binom(float n, float k) {
86-
return binom(static_cast<double>(n), static_cast<double>(k));
87-
}
85+
XSF_HOST_DEVICE inline float binom(float n, float k) { return binom(static_cast<double>(n), static_cast<double>(k)); }
8886

8987
} // namespace xsf

‎include/xsf/cephes/dd_real.h

+12-34
Original file line numberDiff line numberDiff line change
@@ -112,24 +112,16 @@ namespace cephes {
112112
double hi, lo;
113113

114114
double_double() = default;
115-
double_double(double high, double low) : hi(high), lo(low) {
116-
}
117-
explicit double_double(double high) : hi(high), lo(0.0) {
118-
}
115+
double_double(double high, double low) : hi(high), lo(low) {}
116+
explicit double_double(double high) : hi(high), lo(0.0) {}
119117

120-
XSF_HOST_DEVICE explicit operator double() const {
121-
return hi;
122-
}
123-
XSF_HOST_DEVICE explicit operator int() const {
124-
return static_cast<int>(hi);
125-
}
118+
XSF_HOST_DEVICE explicit operator double() const { return hi; }
119+
XSF_HOST_DEVICE explicit operator int() const { return static_cast<int>(hi); }
126120
};
127121

128122
// Arithmetic operations
129123

130-
XSF_HOST_DEVICE inline double_double operator-(const double_double &x) {
131-
return double_double(-x.hi, -x.lo);
132-
}
124+
XSF_HOST_DEVICE inline double_double operator-(const double_double &x) { return double_double(-x.hi, -x.lo); }
133125

134126
XSF_HOST_DEVICE inline double_double operator+(const double_double &lhs, const double_double &rhs) {
135127
/* This one satisfies IEEE style error bound,
@@ -280,13 +272,9 @@ namespace cephes {
280272
return rhs < lhs;
281273
}
282274

283-
XSF_HOST_DEVICE inline bool operator<(const double lhs, const double_double &rhs) {
284-
return rhs > lhs;
285-
}
275+
XSF_HOST_DEVICE inline bool operator<(const double lhs, const double_double &rhs) { return rhs > lhs; }
286276

287-
XSF_HOST_DEVICE inline bool operator>(const double lhs, const double_double &rhs) {
288-
return rhs < lhs;
289-
}
277+
XSF_HOST_DEVICE inline bool operator>(const double lhs, const double_double &rhs) { return rhs < lhs; }
290278

291279
XSF_HOST_DEVICE inline bool operator<=(const double_double &lhs, const double_double &rhs) {
292280
if (lhs.hi < rhs.hi) {
@@ -313,13 +301,9 @@ namespace cephes {
313301
return rhs <= lhs;
314302
}
315303

316-
XSF_HOST_DEVICE inline bool operator>=(const double lhs, const double_double &rhs) {
317-
return rhs <= lhs;
318-
}
304+
XSF_HOST_DEVICE inline bool operator>=(const double lhs, const double_double &rhs) { return rhs <= lhs; }
319305

320-
XSF_HOST_DEVICE inline bool operator<=(const double lhs, const double_double &rhs) {
321-
return rhs >= lhs;
322-
}
306+
XSF_HOST_DEVICE inline bool operator<=(const double lhs, const double_double &rhs) { return rhs >= lhs; }
323307

324308
// Math functions
325309

@@ -328,13 +312,9 @@ namespace cephes {
328312
return double_double(lhs.hi * rhs, lhs.lo * rhs);
329313
}
330314

331-
XSF_HOST_DEVICE inline bool isfinite(const double_double &a) {
332-
return std::isfinite(a.hi);
333-
}
315+
XSF_HOST_DEVICE inline bool isfinite(const double_double &a) { return std::isfinite(a.hi); }
334316

335-
XSF_HOST_DEVICE inline bool isinf(const double_double &a) {
336-
return std::isinf(a.hi);
337-
}
317+
XSF_HOST_DEVICE inline bool isinf(const double_double &a) { return std::isinf(a.hi); }
338318

339319
XSF_HOST_DEVICE inline double_double round(const double_double &a) {
340320
double hi = two_nint(a.hi);
@@ -388,9 +368,7 @@ namespace cephes {
388368
return (a.hi >= 0.0) ? floor(a) : ceil(a);
389369
}
390370

391-
XSF_HOST_DEVICE inline double_double abs(const double_double &a) {
392-
return (a.hi < 0.0 ? -a : a);
393-
}
371+
XSF_HOST_DEVICE inline double_double abs(const double_double &a) { return (a.hi < 0.0 ? -a : a); }
394372

395373
XSF_HOST_DEVICE inline double_double fmod(const double_double &lhs, const double_double &rhs) {
396374
double_double n = trunc(lhs / rhs);

‎include/xsf/cephes/incbet.h

+1-3
Original file line numberDiff line numberDiff line change
@@ -82,9 +82,7 @@ namespace cephes {
8282
#pragma GCC push_options
8383
#pragma GCC optimize("00")
8484
#endif
85-
XSF_HOST_DEVICE inline double triple_product(double u, double v, double w) {
86-
return (u * v) * w;
87-
}
85+
XSF_HOST_DEVICE inline double triple_product(double u, double v, double w) { return (u * v) * w; }
8886
#if defined(__GNUC__) && defined(__i386__)
8987
#pragma GCC pop_options
9088
#endif

‎include/xsf/cephes/lanczos.h

+1-3
Original file line numberDiff line numberDiff line change
@@ -85,9 +85,7 @@ namespace cephes {
8585
-9.8907772644920670589288081640128194231, 6.565936202082889535528455955485877361223
8686
};
8787

88-
XSF_HOST_DEVICE double lanczos_sum(double x) {
89-
return ratevl(x, lanczos_num, 12, lanczos_denom, 12);
90-
}
88+
XSF_HOST_DEVICE double lanczos_sum(double x) { return ratevl(x, lanczos_num, 12, lanczos_denom, 12); }
9189

9290
XSF_HOST_DEVICE double lanczos_sum_near_1(double dx) {
9391
double result = 0;

‎include/xsf/cephes/owens_t.h

+2-6
Original file line numberDiff line numberDiff line change
@@ -116,13 +116,9 @@ namespace cephes {
116116
return owens_t_SELECT_METHOD[iaint * 15 + ihint];
117117
}
118118

119-
XSF_HOST_DEVICE inline double owens_t_norm1(double x) {
120-
return xsf::cephes::erf(x / std::sqrt(2)) / 2;
121-
}
119+
XSF_HOST_DEVICE inline double owens_t_norm1(double x) { return xsf::cephes::erf(x / std::sqrt(2)) / 2; }
122120

123-
XSF_HOST_DEVICE inline double owens_t_norm2(double x) {
124-
return xsf::cephes::erfc(x / std::sqrt(2)) / 2;
125-
}
121+
XSF_HOST_DEVICE inline double owens_t_norm2(double x) { return xsf::cephes::erfc(x / std::sqrt(2)) / 2; }
126122

127123
XSF_HOST_DEVICE inline double owensT1(double h, double a, double m) {
128124
int j = 1;

‎include/xsf/cephes/struve.h

+2-6
Original file line numberDiff line numberDiff line change
@@ -374,13 +374,9 @@ namespace cephes {
374374
}
375375
} // namespace detail
376376

377-
XSF_HOST_DEVICE inline double struve_h(double v, double z) {
378-
return detail::struve_hl(v, z, 1);
379-
}
377+
XSF_HOST_DEVICE inline double struve_h(double v, double z) { return detail::struve_hl(v, z, 1); }
380378

381-
XSF_HOST_DEVICE inline double struve_l(double v, double z) {
382-
return detail::struve_hl(v, z, 0);
383-
}
379+
XSF_HOST_DEVICE inline double struve_l(double v, double z) { return detail::struve_hl(v, z, 0); }
384380

385381
} // namespace cephes
386382
} // namespace xsf

‎include/xsf/cephes/tandg.h

+2-6
Original file line numberDiff line numberDiff line change
@@ -131,13 +131,9 @@ namespace cephes {
131131

132132
} // namespace detail
133133

134-
XSF_HOST_DEVICE inline double tandg(double x) {
135-
return (detail::tancot(x, 0));
136-
}
134+
XSF_HOST_DEVICE inline double tandg(double x) { return (detail::tancot(x, 0)); }
137135

138-
XSF_HOST_DEVICE inline double cotdg(double x) {
139-
return (detail::tancot(x, 1));
140-
}
136+
XSF_HOST_DEVICE inline double cotdg(double x) { return (detail::tancot(x, 1)); }
141137

142138
} // namespace cephes
143139
} // namespace xsf

‎include/xsf/cephes/zetac.h

+1-3
Original file line numberDiff line numberDiff line change
@@ -197,9 +197,7 @@ namespace cephes {
197197
* Compute zetac for small negative x. We can't use the reflection
198198
* formula because to double precision 1 - x = 1 and zetac(1) = inf.
199199
*/
200-
XSF_HOST_DEVICE inline double zetac_smallneg(double x) {
201-
return xsf::cephes::polevl(x, zetac_TAYLOR0, 9);
202-
}
200+
XSF_HOST_DEVICE inline double zetac_smallneg(double x) { return xsf::cephes::polevl(x, zetac_TAYLOR0, 9); }
203201

204202
/*
205203
* Compute zetac using the reflection formula (see DLMF 25.4.2) plus

0 commit comments

Comments
 (0)