exprtk_pgo.cpp

Go to the documentation of this file.

/*
 **************************************************************
 *         C++ Mathematical Expression Toolkit Library        *
 *                                                            *
 * ExprTk PGO Suite                                           *
 * Author: Arash Partow (1999-2016)                           *
 * URL: http://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 Common Public License.         *
 * http://www.opensource.org/licenses/cpl1.0.php              *
 *                                                            *
 **************************************************************
*/
 
 
#include <cstdio>
#include <fstream>
#include <iostream>
#include <string>
#include <deque>
 
#include "exprtk.hpp"
 
 
static const double lower_bound_x = -100.0;
static const double lower_bound_y = -100.0;
static const double upper_bound_x = +100.0;
static const double upper_bound_y = +100.0;
static const double delta         = 0.0111;
 
 
template <typename Allocator,
          template <typename,typename> class Sequence>
inline std::size_t load_expressions(const std::string& file_name,
                                    Sequence<std::string,Allocator>& sequence)
{
   std::ifstream stream(file_name.c_str());
   if (!stream) return 0;
   std::string buffer;
   buffer.reserve(1024);
   std::size_t line_count = 0;
   while (std::getline(stream,buffer))
   {
      ++line_count;
      sequence.push_back(buffer);
   }
   return line_count;
}
 
int main()
{
   std::deque<std::string> expr_str_list;
 
   double x = 0.0;
   double y = 0.0;
   double z = 0.0;
   double w = 0.0;
 
   exprtk::symbol_table<double> symbol_table;
   symbol_table.add_constants();
   symbol_table.add_variable("x",x);
   symbol_table.add_variable("y",y);
   symbol_table.add_variable("z",z);
   symbol_table.add_variable("w",w);
 
   std::deque<exprtk::expression<double> > expression_list;
 
   exprtk::parser<double> parser;
   exprtk::expression<double> temp_expr;
 
   temp_expr.register_symbol_table(symbol_table);
 
   {
      static const std::size_t rounds = 10;
      for (std::size_t r = 0; r < rounds; r++)
      {
         load_expressions("pgo_expression_list.txt",expr_str_list);
 
         for (std::size_t i = 0; i < expr_str_list.size(); ++i)
         {
            if (!parser.compile(expr_str_list[i],temp_expr))
            {
               printf("[load_expression] - [%04d] Parser Error: %s\tExpression: %s\n",
                      i,
                      parser.error().c_str(),
                      expr_str_list[i].c_str());
            }
            else
               expression_list.push_back(temp_expr);
         }
 
         std::cout << "Expression load: " << r << "  Loaded: " << expression_list.size() << std::endl;
         expr_str_list.clear();
         expression_list.clear();
      }
   }
 
   {
      load_expressions("pgo_expression_list.txt",expr_str_list);
 
      for (std::size_t i = 0; i < expr_str_list.size(); ++i)
      {
         if (!parser.compile(expr_str_list[i],temp_expr))
         {
               printf("[load_expression] - [%04d] Parser Error: %s\tExpression: %s\n",
                      i,
                      parser.error().c_str(),
                      expr_str_list[i].c_str());
         }
         else
            expression_list.push_back(temp_expr);
      }
 
      static const std::size_t rounds = 100000;
      double total = 0.0;
 
      exprtk::timer timer;
      timer.start();
 
      x = 0.001;
      y = 0.002;
      z = 0.003;
      w = 0.004;
 
      for (std::size_t r = 0; r < rounds; r++)
      {
         for (std::size_t i = 0; i < expression_list.size(); ++i)
         {
            total += expression_list[i].value();
         }
 
         x += 0.01;
         y += 0.02;
         z += 0.03;
         w += 0.04;
      }
 
      timer.stop();
 
      if (0.0 != total)
      {
         printf("Total Time:%12.8fsec Rate:%15.5fevals/sec\n",
                timer.time(),
                (rounds * expression_list.size()) / timer.time());
      }
   }
 
   return 0;
}