|
6 | 6 | // SCOREC/core apf/apfShape.cc @ 7cd76473 |
7 | 7 |
|
8 | 8 | namespace { |
9 | | -template <typename Array> KOKKOS_INLINE_FUNCTION bool sumsToOne(Array &xi, double tol=10*MeshField::MachinePrecision) { |
10 | | - const bool sums_to_one = [](){ |
| 9 | +template <typename Array> |
| 10 | +KOKKOS_INLINE_FUNCTION bool |
| 11 | +sumsToOne(Array &xi, double tol = 10 * MeshField::MachinePrecision) { |
| 12 | + const bool sums_to_one = []() { |
11 | 13 | auto sum = 0.0; |
12 | 14 | for (size_t i = 0; i < xi.size(); i++) { |
13 | 15 | sum += xi[i]; |
14 | 16 | } |
15 | 17 | return (Kokkos::fabs(sum - 1) <= tol); |
16 | 18 | }(); |
17 | | - if(!sums_to_one) { |
18 | | - for (int i = 0; i < xi.size(); i++) { |
| 19 | + if (!sums_to_one) { |
| 20 | + for (int i = 0; i < xi.size(); i++) { |
19 | 21 | printf("%e ", xi[i]); |
20 | | - } |
21 | | - printf("\n"); |
22 | | - printf("sum: %e tol: %e \n", std::fabs(sum-1), tol); |
| 22 | + } |
| 23 | + printf("\n"); |
| 24 | + printf("sum: %e tol: %e \n", std::fabs(sum - 1), tol); |
23 | 25 | } |
24 | 26 | return sums_to_one; |
25 | 27 | } |
26 | 28 |
|
27 | 29 | template <typename Array> |
28 | | -KOKKOS_INLINE_FUNCTION bool greaterThanOrEqualZero(Array &xi, double tol=1E-12) { |
| 30 | +KOKKOS_INLINE_FUNCTION bool greaterThanOrEqualZero(Array &xi, |
| 31 | + double tol = 1E-12) { |
29 | 32 | for (size_t i = 0; i < xi.size(); i++) { |
30 | | - if(xi[i] < -tol){ |
31 | | - printf("failure %d, %e, %e\n", i, xi[i], tol); |
| 33 | + if (xi[i] < -tol) { |
| 34 | + printf("failure %d, %e, %e\n", i, xi[i], tol); |
32 | 35 | return false; |
33 | 36 | } |
34 | 37 | } |
|
0 commit comments