{

// Test muParser.

// Evaluate "i+1" one million times.

{

mup::ParserX parser;//(mup::pckCOMMON | mup::pckNON_COMPLEX | mup::pckMATRIX | mup::pckUNIT);

mup::Value value1{0};

mup::Variable variable1{& value1};

parser.DefineVar("v1", variable1);

parser.SetExpr("v1 + 1");

const auto timeStart = std::chrono::high_resolution_clock::now();

for(int i = 0; i < 1000000; i++) {

value1 = i;

if(parser.Eval().GetInteger() != i + 1) {

std::cerr << __FILE__ << "(" << __LINE__<< ") "

<< "Wrong value for variable: " << parser.Eval().GetInteger()

<< ", should be: " << i + 1

<< std::endl;

abort();

}

}

const auto timeEnd = std::chrono::high_resolution_clock::now();

const auto total = std::chrono::duration_cast<std::chrono::duration<double> >(timeEnd - timeStart).count();

std::cout << "1e6 evaluations of `i+1` took muParserX (s): " << total << std::endl;

}

// Test EigenLab.

// Evaluate "i+1" one million times with expression caching off.

{

EigenLab::ParserXd parser;

double variable1;

parser.var("v1").setShared(& variable1);

EigenLab::ValueXd result;

parser.setCacheExpressions(false);

const auto timeStart = std::chrono::high_resolution_clock::now();

for(int i = 0; i < 1000000; i++) {

variable1 = i;

result = parser.eval("v1 + 1");

if(result.matrix()(0, 0) != i + 1) {

std::cerr << __FILE__ << "(" << __LINE__<< ") "

<< "Wrong value for variable: " << result.matrix()(0, 0)

<< ", should be: " << i + 1

<< std::endl;

abort();

}

}

const auto timeEnd = std::chrono::high_resolution_clock::now();

const auto total = std::chrono::duration_cast<std::chrono::duration<double> >(timeEnd - timeStart).count();

std::cout << "1e6 evaluations of `i+1` (expression caching off) took EigenLab (s): " << total << std::endl;

}

// Test EigenLab.

// Evaluate "i+1" one million times with expression caching on.

{

EigenLab::ParserXd parser;

double variable1;

parser.var("v1").setShared(& variable1);

EigenLab::ValueXd result;

parser.setCacheExpressions(true);

const auto timeStart = std::chrono::high_resolution_clock::now();

for(int i = 0; i < 1000000; i++) {

variable1 = i;

result = parser.eval("v1 + 1");

if(result.matrix()(0, 0) != i + 1) {

std::cerr << __FILE__ << "(" << __LINE__<< ") "

<< "Wrong value for variable: " << result.matrix()(0, 0)

<< ", should be: " << i + 1

<< std::endl;

abort();

}

}

const auto timeEnd = std::chrono::high_resolution_clock::now();

const auto total = std::chrono::duration_cast<std::chrono::duration<double> >(timeEnd - timeStart).count();

std::cout << "1e6 evaluations of `i+1` (expression caching on) took EigenLab (s): " << total << std::endl;

}

// Test EigenLab.

// Evaluate "i+1" on an array of one million values.

{

EigenLab::ParserXd parser;

const auto timeStart = std::chrono::high_resolution_clock::now();

Eigen::MatrixXd variable1(1, 1000000);

parser.var("v1").setShared(variable1);

for(int i = 0; i < 1000000; i++) {

variable1(0, i) = i;

}

parser.eval("v1 = v1 + 1");

for(int i = 0; i < 1000000; i++) {

if(variable1(i) != i + 1) {

std::cerr << __FILE__ << "(" << __LINE__<< ") "

<< "Wrong value for variable: " << variable1(i)

<< ", should be: " << i + 1

<< std::endl;

abort();

}

}

const auto timeEnd = std::chrono::high_resolution_clock::now();

const auto total = std::chrono::duration_cast<std::chrono::duration<double> >(timeEnd - timeStart).count();

std::cout << "Creation of array i=[0:1e6-1] and evaluation of `i+1` took EigenLab (s): " << total << std::endl;

}

return 0;