Initial commit

include/expr_compiler.h

@@ -0,0 +1,15 @@
+// Reads a human math expression from the command line or a file
+//   and converts it from infix to postfix. e.g; a + b -> a b +
+//   If it is already in postfix then it will skip the infix conversion
+//   Finally it will then convert it into a bytecode stack
+
+#include "expr_compiler.h"
+// For type conversion of human-readable formats
+//  Supports all float formats
+#include "typing.h"
+
+// Converts human readable linear expression into bytecode for the executor
+uint8_t* infix_compile_expression(const char* src);
+// Converts postfix linear expression into bytecode for the executor
+uint8_t* postfix_compile_expression(const char* src);
+

include/expr_eval.h

@@ -0,0 +1,183 @@
+// A virtual-machine to execute a bytecode program based on a mathematical expression
+#include <stdint.h>
+
+typedef double calc_real_t;
+typedef uint64_t calc_natural_t;
+typedef int64_t calc_integer_t;
+typedef calc_real_t[2] calc_complex_t;
+typedef struct {
+  uint32_t m_poffset;
+  uint32_t m_size;
+} calc_vector_t;
+typedef ptrdiff_t calc_idx_t;
+
+typedef uint8_t calc_type_t;
+enum CALC_TYPE
+{
+  CALC_TYPE_NIL = 0,
+  CALC_TYPE_F64 = 1,
+  CALC_TYPE_I64 = 2,
+  CALC_TYPE_U64 = 3,
+  CALC_MOD_VECTOR  = 4,
+
+  CALC_UNIT_RAD = 0
+  CALC_UNIT_DEG = 8,
+  CALC_UNIT_GRAD = 16,
+
+  CALC_FORMAT_RECT = 0
+  CALC_FORMAT_POLAR = 32,
+};
+
+
+struct calc_stack_value_t
+{  
+  union
+  {
+    calc_integer_t m_integer;
+    calc_natural_t m_natural;
+    calc_real_t m_real;
+    calc_complex_t m_complex;
+    calc_vector_t m_vector;
+  }
+};
+
+struct calc_state_t
+{
+  calc_stack_value_t* m_pstack;
+  calc_type_t* m_pstacktypes;
+  uint8_t* m_pmem;
+  uint8_t* m_pprogram;
+  size_t m_stacksize;
+  size_t m_memsize;
+  size_t m_programsize;
+};
+
+void calc_push_real(calc_state_t* state, double val);
+void calc_push_natural(calc_state_t* state, uint64_t val);
+void calc_push_integer(calc_state_t* state, int64_t val);
+void calc_push_complex(calc_state_t* state, calc_real_t re, calc_real_t im);
+
+void calc_push_real_vector(calc_state_t* state, size_t num_elements, calc_real_t* vals);
+void calc_push_natural_vector(calc_state_t* state, size_t num_elements, calc_real_t* vals);
+void calc_push_integer_vector(calc_state_t* state, size_t num_elements, calc_real_t* vals);
+void calc_push_complex_vector(calc_state_t* state, size_t num_elements, calc_stack_value_t* vals);
+
+calc_type_t calc_type(calc_state_t* state, calc_idx_t idx);
+const char* calc_typename(calc_state_t* state, calc_idx_t idx);
+calc_real_t calc_toreal(calc_state_t* state, calc_idx_t idx);
+calc_natural_t calc_tonatural(calc_state_t* state, calc_idx_t idx);
+calc_integer_t calc_tointeger(calc_state_t* state, calc_idx_t idx);
+calc_complex_t calc_tocomplex(calc_state_t* state, calc_idx_t idx);
+void calc_pop(calc_state_t* state, size_t num_items);
+void calc_remove(calc_state_t* state, calc_idx_t index);
+void calc_insert(calc_state_t* state, calc_idx_t dstidx, calc_idx_t srcidx);
+void calc_swap(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_push(calc_state_t* state, calc_idx_t srcidx);
+void calc_settype(calc_state_t* state, calc_idx_t idx, calc_type_t newtype);
+
+void calc_typemod_degrees(calc_state_t* state, calc_idx_t idx);
+void calc_typemod_radians(calc_state_t* state, calc_idx_t idx);
+void calc_typemod_gradient(calc_state_t* state, calc_idx_t idx);
+void calc_typemod_rectangular(calc_state_t* state, calc_idx_t idx);
+void calc_typemod_polar(calc_state_t* state, calc_idx_t idx);
+
+void calc_reinterpretreal(calc_state_t* state, calc_idx_t idx);
+void calc_reinterpretnatural(calc_state_t* state, calc_idx_t idx);
+void calc_reinterpretinteger(calc_state_t* state, calc_idx_t idx);
+void calc_reinterpretcomplex(calc_state_t* state, calc_idx_t idx);
+
+void calc_castreal(calc_state_t* state, calc_idx_t idx);
+void calc_castnatural(calc_state_t* state, calc_idx_t idx);
+void calc_castinteger(calc_state_t* state, calc_idx_t idx);
+void calc_castcomplex(calc_state_t* state, calc_idx_t idx);
+
+void calc_store_real(calc_state_t* state, size_t addr, calc_idx_t idx);
+void calc_store_natural(calc_state_t* state, size_t addr, calc_idx_t idx);
+void calc_store_integer(calc_state_t* state, size_t addr, calc_idx_t idx);
+void calc_store_complex(calc_state_t* state, size_t addr, calc_idx_t idx);
+void calc_store_vector(calc_state_t* state, size_t addr, calc_idx_t idx);
+
+void calc_load_real(calc_state_t* state, size_t addr);
+void calc_load_natural(calc_state_t* state, size_t addr);
+void calc_load_integer(calc_state_t* state, size_t addr);
+void calc_load_complex(calc_state_t* state, size_t addr);
+void calc_load_vector(calc_state_t* state, size_t addr);
+
+void calc_create_complex(calc_state_t* state, calc_idx_t idxre, calc_idx_t idxim);
+
+// Generic unary ops
+void calc_op_negate(calc_state_t* state, calc_idx_t idx);
+void calc_op_conjugate(calc_state_t* state, calc_idx_t idx);
+void calc_op_abs(calc_state_t* state, calc_idx_t idx);
+void calc_op_sqrt(calc_state_t* state, calc_idx_t idx);
+void calc_op_square(calc_state_t* state, calc_idx_t idx);
+void calc_op_recip(calc_state_t* state, calc_idx_t idx);
+void calc_op_ceil(calc_state_t* state, calc_idx_t idx);
+void calc_op_floor(calc_state_t* state, calc_idx_t idx);
+void calc_op_round(calc_state_t* state, calc_idx_t idx);
+void calc_op_real(calc_state_t* state, calc_idx_t idx);
+void calc_op_imag(calc_state_t* state, calc_idx_t idx);
+
+void calc_op_sin(calc_state_t* state, calc_idx_t idx);
+void calc_op_cos(calc_state_t* state, calc_idx_t idx);
+void calc_op_tan(calc_state_t* state, calc_idx_t idx);
+void calc_op_asin(calc_state_t* state, calc_idx_t idx);
+void calc_op_acos(calc_state_t* state, calc_idx_t idx);
+void calc_op_atan(calc_state_t* state, calc_idx_t idx);
+void calc_op_sinh(calc_state_t* state, calc_idx_t idx);
+void calc_op_cosh(calc_state_t* state, calc_idx_t idx);
+void calc_op_tanh(calc_state_t* state, calc_idx_t idx);
+void calc_op_asinh(calc_state_t* state, calc_idx_t idx);
+void calc_op_acosh(calc_state_t* state, calc_idx_t idx);
+void calc_op_atanh(calc_state_t* state, calc_idx_t idx);
+
+// Generic binary ops
+void calc_op_add(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_sub(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_mul(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_div(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_floor2(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_ceil2(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_round2(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_mod(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_atan2(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_pow(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_root(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+
+// Vector ops
+void calc_op_dot(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_cross(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_conv(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_fft(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_ifft(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_vset1(calc_state_t* state, calc_idx_t idxa);
+
+// Multi-arg ops
+void calc_op_reduce_add(calc_state_t* state, calc_idx_t idxbase, size_t nargs);
+void calc_op_reduce_rcpadd(calc_state_t* state, calc_idx_t idxbase, size_t nargs);
+void calc_op_reduce_addsub(calc_state_t* state, calc_idx_t idxbase, size_t nargs);
+void calc_op_reduce_rcpaddsub(calc_state_t* state, calc_idx_t idxbase, size_t nargs);
+void calc_op_reduce_product(calc_state_t* state, calc_idx_t idxbase, size_t nargs);
+
+// Takes a list of impedances and computes the equivalent series impedance
+void calc_op_ee_parallel(calc_state_t* state, calc_idx_t idxbase, size_t nargs);
+
+// Bitwise
+void calc_op_bnot(calc_state_t* state, calc_idx_t idxa);
+void calc_op_bpopcnt(calc_state_t* state, calc_idx_t idxa);
+void calc_op_band(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bor(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bxor(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bnand(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bnor(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bxnor(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bashl(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bashr(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_blshl(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_blshr(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_brotr(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_brotl(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bget(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_bset(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_breset(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);
+void calc_op_btoggle(calc_state_t* state, calc_idx_t idxa, calc_idx_t idxb);