(original) (raw)

/* ************************************************************** * C++ Mathematical Expression Toolkit Library * * * * Function Call Benchmark (Custom vs Composite vs Native) * * Author: Arash Partow (1999-2024) * * URL: https://www.partow.net/programming/exprtk/index.html * * * * Copyright notice: * * Free use of the Mathematical Expression Toolkit Library is * * permitted under the guidelines and in accordance with the * * most current version of the MIT License. * * https://www.opensource.org/licenses/MIT * * SPDX-License-Identifier: MIT * * * ************************************************************** */ #include #include #include "exprtk.hpp" template struct func0 : public exprtk::ifunction{ using exprtk::ifunction::operator(); func0() : exprtk::ifunction(2) {} inline T operator()(const T& v1, const T& v2) { return T(1) + std::cos(v1 * v2) / T(3); } }; template inline T func1(T v1, T v2) { return T(1) + std::cos(v1 * v2) / T(3); } template void function_call_benchmark01() { typedef exprtk::symbol_table symbol_table_t; typedef exprtk::expression expression_t; typedef exprtk::parser parser_t; typedef exprtk::function_compositor compositor_t; typedef typename compositor_t::function function_t; T v1 = T(1); T v2 = T(2); func0 f0; symbol_table_t symbol_table; symbol_table.add_constants(); symbol_table.add_variable("v1" , v1 ); symbol_table.add_variable("v2" , v2 ); symbol_table.add_function("func0", f0 ); symbol_table.add_function("func1", func1); compositor_t compositor(symbol_table); compositor.add( function_t("func2") .var("v1").var("v2") .expression(" 1 + cos(v1 * v2) / 3;")); const std::string program0 = "func0(v1,v2);"; const std::string program1 = "func1(v1,v2);"; const std::string program2 = "func2(v1,v2);"; expression_t expression0(symbol_table); expression_t expression1(symbol_table); expression_t expression2(symbol_table); parser_t parser; if (!parser.compile(program0,expression0)) { printf("Error: %s\tExpression: %s\n", parser.error().c_str(), program0.c_str()); return; } if (!parser.compile(program1,expression1)) { printf("Error: %s\tExpression: %s\n", parser.error().c_str(), program1.c_str()); return; } if (!parser.compile(program2,expression2)) { printf("Error: %s\tExpression: %s\n", parser.error().c_str(), program1.c_str()); return; } const std::size_t rounds = 100000000; { T total = T(0); exprtk::timer timer; timer.start(); for (std::size_t r = 0; r < rounds; ++r) { total += expression0.value(); std::swap(v1,v2); } timer.stop(); printf("[custom function] Total time: %8.3fsec\tRate: %15.3f\tTotal:%20.5f\n", timer.time(), rounds / timer.time(), total); } { T total = T(0); exprtk::timer timer; timer.start(); for (std::size_t r = 0; r < rounds; ++r) { total += expression1.value(); std::swap(v1,v2); } timer.stop(); printf("[free function ] Total time: %8.3fsec\tRate: %15.3f\tTotal:%20.5f\n", timer.time(), rounds / timer.time(), total); } { T total = T(0); exprtk::timer timer; timer.start(); for (std::size_t r = 0; r < rounds; ++r) { total += expression2.value(); std::swap(v1,v2); } timer.stop(); printf("[compositor ] Total time: %8.3fsec\tRate: %15.3f\tTotal:%20.5f\n", timer.time(), rounds / timer.time(), total); } { T total = T(0); exprtk::timer timer; timer.start(); for (std::size_t r = 0; r < rounds; ++r) { total += func1(v1,v2); std::swap(v1,v2); } timer.stop(); printf("[native ] Total time: %8.3fsec\tRate: %15.3f\tTotal:%20.5f\n", timer.time(), rounds / timer.time(), total); } } int main() { function_call_benchmark01(); return 0; }