summaryrefslogtreecommitdiffstatshomepage
path: root/src/lib/netlist/plib/pfunction.h
blob: 643a8b1ae2638e16b33bdb299f7610f32ba462a5 (plain) (blame)
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
128
129
// license:GPL-2.0+
// copyright-holders:Couriersud

#ifndef PFUNCTION_H_
#define PFUNCTION_H_

///
/// \file pfunction.h
///

#include "pmath.h"
#include "pstate.h"
#include "pstring.h"
#include "putil.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,
			POW,
			SIN,
			COS,
			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:

		/// \brief Constructor with state saving support
		///
		/// \param name Name of this object
		/// \param owner Owner of this object
		/// \param state_manager State manager to handle saving object state
		///
		///
		pfunction(const pstring &name, const void *owner, state_manager_t &state_manager)
		: m_lfsr(0xACE1u)
		{
			state_manager.save_item(owner, m_lfsr, name + ".lfsr");
		}

		/// \brief Constructor without state saving support
		///
		pfunction()
		: m_lfsr(0xACE1u)
		{
		}

		/// \brief Compile an expression
		///
		/// \param inputs Vector of input variables, e.g. {"A","B"}
		/// \param expr infix or postfix expression. default is infix, postrix
		///          to be prefixed with rpn, e.g. "rpn:A B + 1.3 /"
		///
		void compile(const std::vector<pstring> &inputs, const pstring &expr);

		/// \brief Compile a rpn expression
		///
		/// \param inputs Vector of input variables, e.g. {"A","B"}
		/// \param expr Reverse polish notation expression, e.g. "A B + 1.3 /"
		///
		void compile_postfix(const std::vector<pstring> &inputs, const pstring &expr);

		/// \brief Compile an infix expression
		///
		/// \param inputs Vector of input variables, e.g. {"A","B"}
		/// \param expr Infix expression, e.g. "(A+B)/1.3"
		///
		void compile_infix(const std::vector<pstring> &inputs, const pstring &expr);

		/// \brief Evaluate the expression
		///
		/// \param values for input variables, e.g. {1.1, 2.2}
		/// \return value of expression
		///
		NT evaluate(const std::vector<NT> &values) noexcept;

	private:

		void compile_postfix(const std::vector<pstring> &inputs,
				const std::vector<pstring> &cmds, const pstring &expr);

		NT lfsr_random() noexcept
		{
			std::uint16_t lsb = m_lfsr & 1;
			m_lfsr >>= 1;
			if (lsb)
				m_lfsr ^= 0xB400u; // taps 15, 13, 12, 10
			return static_cast<NT>(m_lfsr) / static_cast<NT>(0xffffu);
		}

		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_