blob: 22c9ab9a54013dea66738644b31c148c7dd31c6b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
|
/* ========================================
Copyright (C) 2024-2025 Amlal El Mahrouss, licensed under the Apache 2.0 license.
======================================== */
#pragma once
#include <defines>
#ifndef NAN
#define NAN (__builtin_nanf(""))
#endif // !NAN
/// @file base_math.h
/// @brief Base Mathematic functions.
#ifdef __LIBCXX_USE_DOUBLE__
typedef double real_type;
#else
typedef float real_type;
#endif
namespace std::base_math {
inline constexpr static auto not_a_number = NAN;
/// =========================================================== ///
/// @brief Power of Exponent function.
/// =========================================================== ///
template <size_t Exponent>
inline real_type pow(real_type in) {
if (Exponent == 0) return 1; // Any number to the power of 0 is 1.
if (Exponent == 1) return in; // Any number to the power of 1 is itself.
real_type result = 1;
for (auto i = 0UL; i < Exponent; ++i) result *= in;
return result;
}
/// =========================================================== ///
/// @brief Square root function.
/// =========================================================== ///
inline real_type sqrt(real_type in) {
if (in == 0) return 0;
if (in == not_a_number) return not_a_number;
auto constexpr const static Base = 2;
auto x = in / Base;
for (int i = 0; i < 10; ++i) {
x = (x + in / x) / Base;
}
return x;
}
/// =========================================================== ///
/// @brief Square of function, with Base template argument.
/// @param of Base argument to find the square of.
/// =========================================================== ///
template <size_t Base>
inline real_type surd(real_type in) {
if (in == 0) return 0;
if (in == 1) return 1;
if (Base == 1) return in;
if (Base == 2) return sqrt(in);
return not_a_number;
}
/// =========================================================== ///
/// @brief Linear interpolation equation solver.
/// @param from where?
/// @param to to?
/// @param Updated diff value according to difference.
/// =========================================================== ///
inline real_type lerp(real_type to, real_type from, real_type stat) {
real_type diff = (to - from);
return from + (diff * stat);
}
using real_domain = double;
struct complex_domain final {
double Re;
double Im;
};
typename<class Result> using callable_type = Result (*)(size_t n, ...);
} // namespace std::base_math
#ifdef __cpp_lib_base_math
#define __cpp_lib_base_math 1
#endif
|
