src/math/simpson/lib.hpp
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | #pragma once | ||
| 2 | |||
| 3 | #include "../../util/consts/lib.hpp" | ||
| 4 | |||
| 5 | namespace tifa_libs { | ||
| 6 | |||
| 7 | template <std::floating_point FP, FP (*f)(FP)> | ||
| 8 | class simpson_impl { | ||
| 9 | 849050 | static CEXP FP sps(FP l, FP r, FP fl, FP fmid, FP fr) NE { return (fl + 4 * fmid + fr) * (r - l) / 6; } | |
| 10 | 419874 | CEXP FP asr(FP l, FP r, FP area, FP eps, i64 dep) CNE { | |
| 11 | 419874 | const FP mid = r - (r - l) * .5, lmid = mid - (mid - l) * .5, rmid = r - (r - mid) * .5; | |
| 12 | 419874 | const FP fl = f(l), flmid = f(lmid), fmid = f(mid), frmid = f(rmid), fr = f(r); | |
| 13 | 419874 | const FP ls = sps(l, mid, fl, flmid, fmid), rs = sps(mid, r, fmid, frmid, fr); | |
| 14 |
5/6✓ Branch 1 taken 214588 times.
✓ Branch 2 taken 205286 times.
✓ Branch 3 taken 214588 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 214588 times.
✓ Branch 6 taken 205286 times.
|
419874 | if (abs(ls + rs - area) <= 15 * eps && dep <= 0) return ls + rs + (ls + rs - area) / 15; |
| 15 | 205286 | return asr(l, mid, ls, eps * .5, dep - 1) + asr(mid, r, rs, eps * .5, dep - 1); | |
| 16 | } | ||
| 17 | |||
| 18 | public: | ||
| 19 | 9302 | CEXP FP operator()(FP l, FP r, FP eps = eps_v<FP>, i64 min_dep = -1) CNE { return asr(l, r, sps(l, r, f(l), f(r - (r - l) * .5), f(r)), eps, min_dep); } | |
| 20 | }; | ||
| 21 | |||
| 22 | } // namespace tifa_libs | ||
| 23 |