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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
// license:GPL-2.0+
// copyright-holders:Couriersud
#ifndef PFUNCTION_H_
#define PFUNCTION_H_
///
/// \file pfunction.h
///
#include "pmath.h"
#include "pstring.h"
#include <vector>
namespace plib {
//============================================================
// function evaluation
//============================================================
/// \brief Class providing support for evaluating expressions
///
/// \tparam NT Number type, should be float or double
///
template <typename NT>
class pfunction
{
enum rpn_cmd
{
ADD,
MULT,
SUB,
DIV,
EQ,
NE,
LT,
GT,
LE,
GE,
IF,
POW,
LOG,
SIN,
COS,
MIN,
MAX,
RAND, /// random number between 0 and 1
TRUNC,
PUSH_CONST,
PUSH_INPUT
};
struct rpn_inst
{
rpn_inst() : m_cmd(ADD), m_param(plib::constants<NT>::zero()) { }
rpn_cmd m_cmd;
NT m_param;
};
public:
using value_type = NT;
using inputs_container = std::vector<pstring>;
using values_container = std::vector<value_type>;
/// \brief Constructor
///
pfunction()
: m_lfsr(0xace1U) // NOLINT
{
}
/// \brief Compile an expression
///
/// \param expr infix or postfix expression. default is infix, postrix
/// to be prefixed with rpn, e.g. "rpn:A B + 1.3 /"
/// \param inputs Vector of input variables, e.g. {"A","B"}
///
void compile(const pstring &expr, const inputs_container &inputs) noexcept(false);
/// \brief Compile a rpn expression
///
/// \param expr Reverse polish notation expression, e.g. "A B + 1.3 /"
/// \param inputs Vector of input variables, e.g. {"A","B"}
///
void compile_postfix(const pstring &expr, const inputs_container &inputs) noexcept(false);
/// \brief Compile an infix expression
///
/// \param expr Infix expression, e.g. "(A+B)/1.3"
/// \param inputs Vector of input variables, e.g. {"A","B"}
///
void compile_infix(const pstring &expr, const inputs_container &inputs) noexcept(false);
/// \brief Evaluate the expression
///
/// \param values for input variables, e.g. {1.1, 2.2}
/// \return value of expression
///
value_type evaluate(const values_container &values = values_container()) noexcept;
template <typename ST>
void save_state(ST &st)
{
st.save_item(m_lfsr, "m_lfsr");
}
private:
void compile_postfix(const inputs_container &inputs,
const std::vector<pstring> &cmds, const pstring &expr);
std::vector<rpn_inst> m_precompiled; //!< precompiled expression
std::uint16_t m_lfsr; //!< lfsr used for generating random numbers
};
extern template class pfunction<float>;
extern template class pfunction<double>;
extern template class pfunction<long double>;
#if (PUSE_FLOAT128)
extern template class pfunction<FLOAT128>;
#endif
} // namespace plib
#endif // PEXCEPTION_H_
|