/*
 **************************************************************
 *         C++ Mathematical Expression Toolkit Library        *
 *                                                            *
 * Examples and Unit-Tests                                    *
 * Author: Arash Partow (1999-2013)                           *
 * 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 <cmath>
#include <cstddef>
#include <cstdio>
#include <deque>
#include <fstream>
#include <iostream>
#include <string>
#include <vector>

#include "exprtk.hpp"


typedef std::pair<std::string,double> test_t;

static const test_t test_list[] =
                  {
                     // Note: Each of following tests should all
                     // compile down to a single literal node.
                     test_t("0",0.0),
                     test_t("1",1.0),
                     test_t("2",2.0),
                     test_t("3",3.0),
                     test_t("4",4.0),
                     test_t("5",5.0),
                     test_t("6",6.0),
                     test_t("7",7.0),
                     test_t("8",8.0),
                     test_t("9",9.0),
                     test_t("12.12",12.12),
                     test_t("123.123",123.123),
                     test_t("1234.1234",1234.1234),
                     test_t("12345.12345",12345.12345),
                     test_t("123456.123456",123456.123456),
                     test_t("0.0",0.0),
                     test_t("1.0",1.0),
                     test_t("2.0",2.0),
                     test_t("3.0",3.0),
                     test_t("4.0",4.0),
                     test_t("5.0",5.0),
                     test_t("6.0",6.0),
                     test_t("7.0",7.0),
                     test_t("8.0",8.0),
                     test_t("9.0",9.0),
                     test_t("0.0",0.0),
                     test_t("1.1",1.1),
                     test_t("2.2",2.2),
                     test_t("3.3",3.3),
                     test_t("4.4",4.4),
                     test_t("5.5",5.5),
                     test_t("6.6",6.6),
                     test_t("7.7",7.7),
                     test_t("8.8",8.8),
                     test_t("9.9",9.9),
                     test_t("+0",0.0),
                     test_t("+1",1.0),
                     test_t("+2",2.0),
                     test_t("+3",3.0),
                     test_t("+4",4.0),
                     test_t("+5",5.0),
                     test_t("+6",6.0),
                     test_t("+7",7.0),
                     test_t("+8",8.0),
                     test_t("+9",9.0),
                     test_t("+0.0",0.0),
                     test_t("+1.0",1.0),
                     test_t("+2.0",2.0),
                     test_t("+3.0",3.0),
                     test_t("+4.0",4.0),
                     test_t("+5.0",5.0),
                     test_t("+6.0",6.0),
                     test_t("+7.0",7.0),
                     test_t("+8.0",8.0),
                     test_t("+9.0",9.0),
                     test_t("+0.0",0.0),
                     test_t("+1.1",1.1),
                     test_t("+2.2",2.2),
                     test_t("+3.3",3.3),
                     test_t("+4.4",4.4),
                     test_t("+5.5",5.5),
                     test_t("+6.6",6.6),
                     test_t("+7.7",7.7),
                     test_t("+8.8",8.8),
                     test_t("+9.9",9.9),
                     test_t("-0",-0.0),
                     test_t("-1",-1.0),
                     test_t("-2",-2.0),
                     test_t("-3",-3.0),
                     test_t("-4",-4.0),
                     test_t("-5",-5.0),
                     test_t("-6",-6.0),
                     test_t("-7",-7.0),
                     test_t("-8",-8.0),
                     test_t("-9",-9.0),
                     test_t("-0.0",-0.0),
                     test_t("-1.0",-1.0),
                     test_t("-2.0",-2.0),
                     test_t("-3.0",-3.0),
                     test_t("-4.0",-4.0),
                     test_t("-5.0",-5.0),
                     test_t("-6.0",-6.0),
                     test_t("-7.0",-7.0),
                     test_t("-8.0",-8.0),
                     test_t("-9.0",-9.0),
                     test_t("-0.0",-0.0),
                     test_t("-1.1",-1.1),
                     test_t("-2.2",-2.2),
                     test_t("-3.3",-3.3),
                     test_t("-4.4",-4.4),
                     test_t("-5.5",-5.5),
                     test_t("-6.6",-6.6),
                     test_t("-7.7",-7.7),
                     test_t("-8.8",-8.8),
                     test_t("-9.9",-9.9),
                     test_t("0.0e+0",+0.0e+0),
                     test_t("1.1e+1",+1.1e+1),
                     test_t("2.2e+2",+2.2e+2),
                     test_t("3.3e+3",+3.3e+3),
                     test_t("4.4e+4",+4.4e+4),
                     test_t("5.5e+5",+5.5e+5),
                     test_t("6.6e+6",+6.6e+6),
                     test_t("7.7e+7",+7.7e+7),
                     test_t("8.8e+8",+8.8e+8),
                     test_t("9.9e+9",+9.9e+9),
                     test_t("-0.0e+0",-0.0e+0),
                     test_t("-1.1e+1",-1.1e+1),
                     test_t("-2.2e+2",-2.2e+2),
                     test_t("-3.3e+3",-3.3e+3),
                     test_t("-4.4e+4",-4.4e+4),
                     test_t("-5.5e+5",-5.5e+5),
                     test_t("-6.6e+6",-6.6e+6),
                     test_t("-7.7e+7",-7.7e+7),
                     test_t("-8.8e+8",-8.8e+8),
                     test_t("-9.9e+9",-9.9e+9),
                     test_t("0.0E+0",+0.0E+0),
                     test_t("1.1E+1",+1.1E+1),
                     test_t("2.2E+2",+2.2E+2),
                     test_t("3.3E+3",+3.3E+3),
                     test_t("4.4E+4",+4.4E+4),
                     test_t("5.5E+5",+5.5E+5),
                     test_t("6.6E+6",+6.6E+6),
                     test_t("7.7E+7",+7.7E+7),
                     test_t("8.8E+8",+8.8E+8),
                     test_t("9.9E+9",+9.9E+9),
                     test_t("-0.0E+0",-0.0E+0),
                     test_t("-1.1E+1",-1.1E+1),
                     test_t("-2.2E+2",-2.2E+2),
                     test_t("-3.3E+3",-3.3E+3),
                     test_t("-4.4E+4",-4.4E+4),
                     test_t("-5.5E+5",-5.5E+5),
                     test_t("-6.6E+6",-6.6E+6),
                     test_t("-7.7E+7",-7.7E+7),
                     test_t("-8.8E+8",-8.8E+8),
                     test_t("-9.9E+9",-9.9E+9),
                     test_t("(0)",0.0),
                     test_t("(1)",1.0),
                     test_t("(2)",2.0),
                     test_t("(3)",3.0),
                     test_t("(4)",4.0),
                     test_t("(5)",5.0),
                     test_t("(6)",6.0),
                     test_t("(7)",7.0),
                     test_t("(8)",8.0),
                     test_t("(9)",9.0),
                     test_t("(0.0)",0.0),
                     test_t("(1.0)",1.0),
                     test_t("(2.0)",2.0),
                     test_t("(3.0)",3.0),
                     test_t("(4.0)",4.0),
                     test_t("(5.0)",5.0),
                     test_t("(6.0)",6.0),
                     test_t("(7.0)",7.0),
                     test_t("(8.0)",8.0),
                     test_t("(9.0)",9.0),
                     test_t("(0.0)",0.0),
                     test_t("(1.1)",1.1),
                     test_t("(2.2)",2.2),
                     test_t("(3.3)",3.3),
                     test_t("(4.4)",4.4),
                     test_t("(5.5)",5.5),
                     test_t("(6.6)",6.6),
                     test_t("(7.7)",7.7),
                     test_t("(8.8)",8.8),
                     test_t("(9.9)",9.9),
                     test_t("(+0)",0.0),
                     test_t("(+1)",1.0),
                     test_t("(+2)",2.0),
                     test_t("(+3)",3.0),
                     test_t("(+4)",4.0),
                     test_t("(+5)",5.0),
                     test_t("(+6)",6.0),
                     test_t("(+7)",7.0),
                     test_t("(+8)",8.0),
                     test_t("(+9)",9.0),
                     test_t("(+0.0)",0.0),
                     test_t("(+1.0)",1.0),
                     test_t("(+2.0)",2.0),
                     test_t("(+3.0)",3.0),
                     test_t("(+4.0)",4.0),
                     test_t("(+5.0)",5.0),
                     test_t("(+6.0)",6.0),
                     test_t("(+7.0)",7.0),
                     test_t("(+8.0)",8.0),
                     test_t("(+9.0)",9.0),
                     test_t("(+0.0)",0.0),
                     test_t("(+1.1)",1.1),
                     test_t("(+2.2)",2.2),
                     test_t("(+3.3)",3.3),
                     test_t("(+4.4)",4.4),
                     test_t("(+5.5)",5.5),
                     test_t("(+6.6)",6.6),
                     test_t("(+7.7)",7.7),
                     test_t("(+8.8)",8.8),
                     test_t("(+9.9)",9.9),
                     test_t("(-0)",-0.0),
                     test_t("(-1)",-1.0),
                     test_t("(-2)",-2.0),
                     test_t("(-3)",-3.0),
                     test_t("(-4)",-4.0),
                     test_t("(-5)",-5.0),
                     test_t("(-6)",-6.0),
                     test_t("(-7)",-7.0),
                     test_t("(-8)",-8.0),
                     test_t("(-9)",-9.0),
                     test_t("(-0.0)",-0.0),
                     test_t("(-1.0)",-1.0),
                     test_t("(-2.0)",-2.0),
                     test_t("(-3.0)",-3.0),
                     test_t("(-4.0)",-4.0),
                     test_t("(-5.0)",-5.0),
                     test_t("(-6.0)",-6.0),
                     test_t("(-7.0)",-7.0),
                     test_t("(-8.0)",-8.0),
                     test_t("(-9.0)",-9.0),
                     test_t("(-0.0)",-0.0),
                     test_t("(-1.1)",-1.1),
                     test_t("(-2.2)",-2.2),
                     test_t("(-3.3)",-3.3),
                     test_t("(-4.4)",-4.4),
                     test_t("(-5.5)",-5.5),
                     test_t("(-6.6)",-6.6),
                     test_t("(-7.7)",-7.7),
                     test_t("(-8.8)",-8.8),
                     test_t("(-9.9)",-9.9),
                     test_t("1234567890",1234567890),
                     test_t("123456789.0",123456789.0),
                     test_t("+1234567890",1234567890),
                     test_t("+123456789.0",123456789.0),
                     test_t("-1234567890",-1234567890),
                     test_t("-123456789.0",-123456789.0),
                     test_t("1234.567890",1234.567890),
                     test_t("-1234.567890",-1234.567890),
                     test_t("0+9",9.0),
                     test_t("1+8",9.0),
                     test_t("2+7",9.0),
                     test_t("3+6",9.0),
                     test_t("4+5",9.0),
                     test_t("5+4",9.0),
                     test_t("6+3",9.0),
                     test_t("7+2",9.0),
                     test_t("8+1",9.0),
                     test_t("9+0",9.0),
                     test_t(" 0 + 9 ",9.0),
                     test_t(" 1 + 8 ",9.0),
                     test_t(" 2 + 7 ",9.0),
                     test_t(" 3 + 6 ",9.0),
                     test_t(" 4 + 5 ",9.0),
                     test_t(" 5 + 4 ",9.0),
                     test_t(" 6 + 3 ",9.0),
                     test_t(" 7 + 2 ",9.0),
                     test_t(" 8 + 1 ",9.0),
                     test_t(" 9 + 0 ",9.0),
                     test_t("( 0 + 9 )",9.0),
                     test_t("( 1 + 8 )",9.0),
                     test_t("( 2 + 7 )",9.0),
                     test_t("( 3 + 6 )",9.0),
                     test_t("( 4 + 5 )",9.0),
                     test_t("( 5 + 4 )",9.0),
                     test_t("( 6 + 3 )",9.0),
                     test_t("( 7 + 2 )",9.0),
                     test_t("( 8 + 1 )",9.0),
                     test_t("( 9 + 0 )",9.0),
                     test_t("1+2",+3.0),
                     test_t("1-2",-1.0),
                     test_t("1*2",+2.0),
                     test_t("1/2",+0.5),
                     test_t("1.1+2.2", +3.3),
                     test_t("1.1-2.2", -1.1),
                     test_t("1.1*2.2",+2.42),
                     test_t("1.1/2.2", +0.5),
                     test_t("0-9",-9.0),
                     test_t("1-8",-7.0),
                     test_t("2-7",-5.0),
                     test_t("3-6",-3.0),
                     test_t("4-5",-1.0),
                     test_t("5-4",+1.0),
                     test_t("6-3",+3.0),
                     test_t("7-2",+5.0),
                     test_t("8-1",+7.0),
                     test_t("9-0",+9.0),
                     test_t(" 0 - 9 ",-9.0),
                     test_t(" 1 - 8 ",-7.0),
                     test_t(" 2 - 7 ",-5.0),
                     test_t(" 3 - 6 ",-3.0),
                     test_t(" 4 - 5 ",-1.0),
                     test_t(" 5 - 4 ",+1.0),
                     test_t(" 6 - 3 ",+3.0),
                     test_t(" 7 - 2 ",+5.0),
                     test_t(" 8 - 1 ",+7.0),
                     test_t(" 9 - 0 ",+9.0),
                     test_t("( 0 - 9 )",-9.0),
                     test_t("( 1 - 8 )",-7.0),
                     test_t("( 2 - 7 )",-5.0),
                     test_t("( 3 - 6 )",-3.0),
                     test_t("( 4 - 5 )",-1.0),
                     test_t("( 5 - 4 )",+1.0),
                     test_t("( 6 - 3 )",+3.0),
                     test_t("( 7 - 2 )",+5.0),
                     test_t("( 8 - 1 )",+7.0),
                     test_t("( 9 - 0 )",+9.0),
                     test_t("-(1+2)",-3.0),
                     test_t("+(1+2)",+3.0),
                     test_t("+(1-2)",-1.0),
                     test_t("-(1-2)",+1.0),
                     test_t("1.1+2.2+3.3",+6.6),
                     test_t("+1.1+2.2+3.3",+6.6),
                     test_t("-1.1-2.2-3.3",-6.6),
                     test_t("1.1*2.2*3.3",+7.986),
                     test_t("+1.1*2.2*3.3",+7.986),
                     test_t("-1.1*-2.2*-3.3",-7.986),
                     test_t("1 + 1/2",+1.5),
                     test_t("1 + (1/2)",+1.5),
                     test_t("1.1 + 1.1/2.2",+1.6),
                     test_t("1.1 + (1.1/2.2)",+1.6),
                     test_t("2 * 1/2",+1.0),
                     test_t("2 * (1/2)",+1.0),
                     test_t("2.2 * 1.1/2.2",+1.1),
                     test_t("2.2 * (1.1/2.2)",+1.1),
                     test_t("1^2",1.0),
                     test_t("2^1",2.0),
                     test_t("2^3",8.0),
                     test_t("-2^3",-8.0),
                     test_t("-2^4",-16.0),
                     test_t("(-2)^3",-8.0),
                     test_t("(-2)^4",+16.0),
                     test_t("3^2^1",9.0),
                     test_t("1.1^2.2",1.23328630055466251099),
                     test_t("2.2^1.1",2.3804822576003541627),
                     test_t("2.2^3.3",13.48946876053338489127),
                     test_t("3.3^2.2^1.1",17.15193942371376191362),
                     test_t("1.1^(1.1 * 2.2)", 1.25941916576299080582),
                     test_t("2.2^(1.1 * 3.3)",17.49823848953534759743),
                     test_t("3.3^(1.1 * 2.2)",17.98058156638874965269),
                     test_t("equal(1.23^3,(1.23 * 1.23 * 1.23))",1.0),
                     test_t("equal(1.23^-3,1/(1.23 * 1.23 * 1.23))",1.0),
                     test_t("equal((2.1 + 1.23^3),(2.1 + [1.23 * 1.23 * 1.23]))",1.0),
                     test_t("equal((2.1 - 1.23^3),(2.1 - [1.23 * 1.23 * 1.23]))",1.0),
                     test_t("equal((2.1 * 1.23^3),(2.1 * [1.23 * 1.23 * 1.23]))",1.0),
                     test_t("equal((2.1 / 1.23^3),(2.1 / [1.23 * 1.23 * 1.23]))",1.0),
                     test_t("equal((1.23^3 + 2.1),({1.23 * 1.23 * 1.23} + 2.1))",1.0),
                     test_t("equal((1.23^3 - 2.1),({1.23 * 1.23 * 1.23} - 2.1))",1.0),
                     test_t("equal((1.23^3 * 2.1),({1.23 * 1.23 * 1.23} * 2.1))",1.0),
                     test_t("equal((1.23^3 / 2.1),({1.23 * 1.23 * 1.23} / 2.1))",1.0),
                     test_t("equal(1.0^(1.0/2.0),sqrt(1.0))",1.0),
                     test_t("equal(1.0^(1.0/2.0),root(1.0,2.0))",1.0),
                     test_t("equal(1.0^(1.0/3.0),root(1.0,3.0))",1.0),
                     test_t("equal(1.0^(1.0/4.0),root(1.0,4.0))",1.0),
                     test_t("equal(1.0^(1.0/5.0),root(1.0,5.0))",1.0),
                     test_t("equal(1.0^(1.0/6.0),root(1.0,6.0))",1.0),
                     test_t("equal(1.0^(1.0/7.0),root(1.0,7.0))",1.0),
                     test_t("equal(1.0^(1.0/8.0),root(1.0,8.0))",1.0),
                     test_t("equal(1.0^(1.0/9.0),root(1.0,9.0))",1.0),
                     test_t("equal(2.0^(1.0/2.0),sqrt(2.0))",1.0),
                     test_t("equal(2.0^(1.0/2.0),root(2.0,2.0))",1.0),
                     test_t("equal(3.0^(1.0/3.0),root(3.0,3.0))",1.0),
                     test_t("equal(4.0^(1.0/4.0),root(4.0,4.0))",1.0),
                     test_t("equal(5.0^(1.0/5.0),root(5.0,5.0))",1.0),
                     test_t("equal(6.0^(1.0/6.0),root(6.0,6.0))",1.0),
                     test_t("equal(7.0^(1.0/7.0),root(7.0,7.0))",1.0),
                     test_t("equal(8.0^(1.0/8.0),root(8.0,8.0))",1.0),
                     test_t("equal(9.0^(1.0/9.0),root(9.0,9.0))",1.0),
                     test_t("1 <  2", 1.0),
                     test_t("1 <= 2", 1.0),
                     test_t("1.1 <= 2.2", 1.0),
                     test_t("(1.0 + 0.1) <= (2.0 + 0.2)", 1.0),
                     test_t("1 >  2", 0.0),
                     test_t("1 >= 2", 0.0),
                     test_t("1.1 >= 2.2", 0.0),
                     test_t("(1.0 + 0.1) >= (2.0 + 0.2)", 0.0),
                     test_t("1 <> 2", 1.0),
                     test_t("1 != 2", 1.0),
                     test_t("1.1 <> 2.2", 1.0),
                     test_t("1.1 != 2.2", 1.0),
                     test_t("(1.0 + 0.1) <> (2.0 + 0.2)", 1.0),
                     test_t("(1.0 + 0.1) != (2.0 + 0.2)", 1.0),
                     test_t("1 == 1", 1.0),
                     test_t("1.1 == 1.1", 1.0),
                     test_t("1 = 1", 1.0),
                     test_t("1.1 = 1.1", 1.0),
                     test_t("1 <> 1", 0.0),
                     test_t("1 != 1", 0.0),
                     test_t("1.1 <> 1.1", 0.0),
                     test_t("1.1 != 1.1", 0.0),
                     test_t("(1.0 + 0.1) <> (1.0 + 0.1)", 0.0),
                     test_t("(1.0 + 0.1) != (1.0 + 0.1)", 0.0),
                     test_t("equal(1.1,1.1)",1.0),
                     test_t("equal(1.1,2.2)",0.0),
                     test_t("not_equal(1.1,1.1)",0.0),
                     test_t("not_equal(1.1,2.2)",1.0),
                     test_t("1 and 1",1.0),
                     test_t("1 and 0",0.0),
                     test_t("0 and 1",0.0),
                     test_t("0 and 0",0.0),
                     test_t("1.0 and 1.0",1.0),
                     test_t("1.0 and 0.0",0.0),
                     test_t("0.0 and 1.0",0.0),
                     test_t("0.0 and 0.0",0.0),
                     test_t("(1 and 1)",1.0),
                     test_t("(1 and 0)",0.0),
                     test_t("(0 and 1)",0.0),
                     test_t("(0 and 0)",0.0),
                     test_t("(1.0 and 1.0)",1.0),
                     test_t("(1.0 and 0.0)",0.0),
                     test_t("(0.0 and 1.0)",0.0),
                     test_t("(0.0 and 0.0)",0.0),
                     test_t("1 or 1",1.0),
                     test_t("1 or 0",1.0),
                     test_t("0 or 1",1.0),
                     test_t("0 or 0",0.0),
                     test_t("1.0 or 1.0",1.0),
                     test_t("1.0 or 0.0",1.0),
                     test_t("0.0 or 1.0",1.0),
                     test_t("0.0 or 0.0",0.0),
                     test_t("(1 or 1)",1.0),
                     test_t("(1 or 0)",1.0),
                     test_t("(0 or 1)",1.0),
                     test_t("(0 or 0)",0.0),
                     test_t("(1.0 or 1.0)",1.0),
                     test_t("(1.0 or 0.0)",1.0),
                     test_t("(0.0 or 1.0)",1.0),
                     test_t("(0.0 or 0.0)",0.0),
                     test_t("1 nand 1",0.0),
                     test_t("1 nand 0",1.0),
                     test_t("0 nand 1",1.0),
                     test_t("0 nand 0",1.0),
                     test_t("1.0 nand 1.0",0.0),
                     test_t("1.0 nand 0.0",1.0),
                     test_t("0.0 nand 1.0",1.0),
                     test_t("0.0 nand 0.0",1.0),
                     test_t("(1 nand 1)",0.0),
                     test_t("(1 nand 0)",1.0),
                     test_t("(0 nand 1)",1.0),
                     test_t("(0 nand 0)",1.0),
                     test_t("(1.0 nand 1.0)",0.0),
                     test_t("(1.0 nand 0.0)",1.0),
                     test_t("(0.0 nand 1.0)",1.0),
                     test_t("(0.0 nand 0.0)",1.0),
                     test_t("1 nor 1",0.0),
                     test_t("1 nor 0",0.0),
                     test_t("0 nor 1",0.0),
                     test_t("0 nor 0",1.0),
                     test_t("1.0 nor 1.0",0.0),
                     test_t("1.0 nor 0.0",0.0),
                     test_t("0.0 nor 1.0",0.0),
                     test_t("0.0 nor 0.0",1.0),
                     test_t("(1 nor 1)",0.0),
                     test_t("(1 nor 0)",0.0),
                     test_t("(0 nor 1)",0.0),
                     test_t("(0 nor 0)",1.0),
                     test_t("(1.0 nor 1.0)",0.0),
                     test_t("(1.0 nor 0.0)",0.0),
                     test_t("(0.0 nor 1.0)",0.0),
                     test_t("(0.0 nor 0.0)",1.0),
                     test_t("0 xor 0",0.0),
                     test_t("0 xor 1",1.0),
                     test_t("1 xor 0",1.0),
                     test_t("1 xor 1",0.0),
                     test_t("0.0 xor 0.0",0.0),
                     test_t("0.0 xor 1.0",1.0),
                     test_t("1.0 xor 0.0",1.0),
                     test_t("1.0 xor 1.0",0.0),
                     test_t("(0 xor 0)",0.0),
                     test_t("(0 xor 1)",1.0),
                     test_t("(1 xor 0)",1.0),
                     test_t("(1 xor 1)",0.0),
                     test_t("(0.0 xor 0.0)",0.0),
                     test_t("(0.0 xor 1.0)",1.0),
                     test_t("(1.0 xor 0.0)",1.0),
                     test_t("(1.0 xor 1.0)",0.0),
                     test_t("1 & 1",1.0),
                     test_t("1 & 0",0.0),
                     test_t("0 & 1",0.0),
                     test_t("0 & 0",0.0),
                     test_t("1.0 & 1.0",1.0),
                     test_t("1.0 & 0.0",0.0),
                     test_t("0.0 & 1.0",0.0),
                     test_t("0.0 & 0.0",0.0),
                     test_t("(1 & 1)",1.0),
                     test_t("(1 & 0)",0.0),
                     test_t("(0 & 1)",0.0),
                     test_t("(0 & 0)",0.0),
                     test_t("(1.0 & 1.0)",1.0),
                     test_t("(1.0 & 0.0)",0.0),
                     test_t("(0.0 & 1.0)",0.0),
                     test_t("(0.0 & 0.0)",0.0),
                     test_t("1 | 1",1.0),
                     test_t("1 | 0",1.0),
                     test_t("0 | 1",1.0),
                     test_t("0 | 0",0.0),
                     test_t("1.0 | 1.0",1.0),
                     test_t("1.0 | 0.0",1.0),
                     test_t("0.0 | 1.0",1.0),
                     test_t("0.0 | 0.0",0.0),
                     test_t("(1 | 1)",1.0),
                     test_t("(1 | 0)",1.0),
                     test_t("(0 | 1)",1.0),
                     test_t("(0 | 0)",0.0),
                     test_t("(1.0 | 1.0)",1.0),
                     test_t("(1.0 | 0.0)",1.0),
                     test_t("(0.0 | 1.0)",1.0),
                     test_t("(0.0 | 0.0)",0.0),
                     test_t("(1 nand 1) == not(1 and 1)",1.0),
                     test_t("(1 nand 0) == not(1 and 0)",1.0),
                     test_t("(0 nand 1) == not(0 and 1)",1.0),
                     test_t("(0 nand 0) == not(0 and 0)",1.0),
                     test_t("(1 nor 1) == not(1 or 1)",1.0),
                     test_t("(1 nor 0) == not(1 or 0)",1.0),
                     test_t("(0 nor 1) == not(0 or 1)",1.0),
                     test_t("(0 nor 0) == not(0 or 0)",1.0),
                     test_t("(1.0 nand 1.0) == not(1.0 and 1.0)",1.0),
                     test_t("(1.0 nand 0.0) == not(1.0 and 0.0)",1.0),
                     test_t("(0.0 nand 1.0) == not(0.0 and 1.0)",1.0),
                     test_t("(0.0 nand 0.0) == not(0.0 and 0.0)",1.0),
                     test_t("(1.0 nor 1.0) == not(1.0 or 1.0)",1.0),
                     test_t("(1.0 nor 0.0) == not(1.0 or 0.0)",1.0),
                     test_t("(0.0 nor 1.0) == not(0.0 or 1.0)",1.0),
                     test_t("(0.0 nor 0.0) == not(0.0 or 0.0)",1.0),
                     test_t("(1 nand 1) == not(1 & 1)",1.0),
                     test_t("(1 nand 0) == not(1 & 0)",1.0),
                     test_t("(0 nand 1) == not(0 & 1)",1.0),
                     test_t("(0 nand 0) == not(0 & 0)",1.0),
                     test_t("(1 nor 1) == not(1 | 1)",1.0),
                     test_t("(1 nor 0) == not(1 | 0)",1.0),
                     test_t("(0 nor 1) == not(0 | 1)",1.0),
                     test_t("(0 nor 0) == not(0 | 0)",1.0),
                     test_t("(1.0 nand 1.0) == not(1.0 & 1.0)",1.0),
                     test_t("(1.0 nand 0.0) == not(1.0 & 0.0)",1.0),
                     test_t("(0.0 nand 1.0) == not(0.0 & 1.0)",1.0),
                     test_t("(0.0 nand 0.0) == not(0.0 & 0.0)",1.0),
                     test_t("(1.0 nor 1.0) == not(1.0 | 1.0)",1.0),
                     test_t("(1.0 nor 0.0) == not(1.0 | 0.0)",1.0),
                     test_t("(0.0 nor 1.0) == not(0.0 | 1.0)",1.0),
                     test_t("(0.0 nor 0.0) == not(0.0 | 0.0)",1.0),
                     test_t("mand(1,1)",1.0),
                     test_t("mand(1,0)",0.0),
                     test_t("mand(0,1)",0.0),
                     test_t("mand(0,0)",0.0),
                     test_t("mand(1.0,1.0)",1.0),
                     test_t("mand(1.0,0.0)",0.0),
                     test_t("mand(0.0,1.0)",0.0),
                     test_t("mand(0.0,0.0)",0.0),
                     test_t("mor(1,1)",1.0),
                     test_t("mor(1,0)",1.0),
                     test_t("mor(0,1)",1.0),
                     test_t("mor(0,0)",0.0),
                     test_t("mor(1.0,1.0)",1.0),
                     test_t("mor(1.0,0.0)",1.0),
                     test_t("mor(0.0,1.0)",1.0),
                     test_t("mor(0.0,0.0)",0.0),
                     test_t("(1 nand 1) == not(mand(1,1))",1.0),
                     test_t("(1 nand 0) == not(mand(1,0))",1.0),
                     test_t("(0 nand 1) == not(mand(0,1))",1.0),
                     test_t("(0 nand 0) == not(mand(0,0))",1.0),
                     test_t("(1 nor 1) == not(mor(1,1))",1.0),
                     test_t("(1 nor 0) == not(mor(1,0))",1.0),
                     test_t("(0 nor 1) == not(mor(0,1))",1.0),
                     test_t("(0 nor 0) == not(mor(0,0))",1.0),
                     test_t("(1.0 nand 1.0) == not(mand(1.0,1.0))",1.0),
                     test_t("(1.0 nand 0.0) == not(mand(1.0,0.0))",1.0),
                     test_t("(0.0 nand 1.0) == not(mand(0.0,1.0))",1.0),
                     test_t("(0.0 nand 0.0) == not(mand(0.0,0.0))",1.0),
                     test_t("(1.0 nor 1.0) == not(mor(1.0,1.0))",1.0),
                     test_t("(1.0 nor 0.0) == not(mor(1.0,0.0))",1.0),
                     test_t("(0.0 nor 1.0) == not(mor(0.0,1.0))",1.0),
                     test_t("(0.0 nor 0.0) == not(mor(0.0,0.0))",1.0),
                     test_t("abs(1)",1.0),
                     test_t("abs(-1)",1.0),
                     test_t("abs(1.0)",1.0),
                     test_t("abs(-1.0)",1.0),
                     test_t("min(1,2)",1.0),
                     test_t("min(1,2,3)",1.0),
                     test_t("min(1,2,3,4)",1.0),
                     test_t("min(1,2,3,4,5)",1.0),
                     test_t("min(1,2,3,4,5,6)",1.0),
                     test_t("min(1.1,2.2)",1.1),
                     test_t("min(1.1,2.2,3.3)",1.1),
                     test_t("min(1.1,2.2,3.3,4.4)",1.1),
                     test_t("min(1.1,2.2,3.3,4.4,5.5)",1.1),
                     test_t("min(1.1,2.2,3.3,4.4,5.5,6.6)",1.1),
                     test_t("min(min(1,2),min(3,4))",1.0),
                     test_t("max(1,2)",2.0),
                     test_t("max(1,2,3)",3.0),
                     test_t("max(1,2,3,4)",4.0),
                     test_t("max(1,2,3,4,5)",5.0),
                     test_t("max(1,2,3,4,5,6)",6.0),
                     test_t("max(1.1,2.2)",2.2),
                     test_t("max(1.1,2.2,3.3)",3.3),
                     test_t("max(1.1,2.2,3.3,4.4)",4.4),
                     test_t("max(1.1,2.2,3.3,4.4,5.5)",5.5),
                     test_t("max(1.1,2.2,3.3,4.4,5.5,6.6)",6.6),
                     test_t("max(max(1,2),max(3,4))",4.0),
                     test_t("avg(1,2)",1.5),
                     test_t("avg(1,2,3)",2.0),
                     test_t("avg(1,2,3,4)",2.5),
                     test_t("avg(1,2,3,4,5)",3.0),
                     test_t("avg(1.1,2.2)",1.65),
                     test_t("avg(1.1,2.2,3.3)",2.2),
                     test_t("avg(1.1,2.2,3.3,4.4)",2.75),
                     test_t("avg(1.1,2.2,3.3,4.4,5.5)",3.3),
                     test_t("avg(1.1,2.2,3.3,4.4,5.5,6.6)",3.85),
                     test_t("sum(1,2)",3.0),
                     test_t("sum(1,2,3)",6.0),
                     test_t("sum(1,2,3,4)",10),
                     test_t("sum(1,2,3,4,5)",15.0),
                     test_t("sum(1,2,3,4,5,6)",21),
                     test_t("sum(1.1,2.2)",3.3),
                     test_t("sum(1.1,2.2,3.3)",6.6),
                     test_t("sum(1.1,2.2,3.3,4.4)",11.0),
                     test_t("sum(1.1,2.2,3.3,4.4,5.5)",16.5),
                     test_t("sum(1.1,2.2,3.3,4.4,5.5,6.6)",23.1),
                     test_t("mul(1,2)",2.0),
                     test_t("mul(1,2,3)",6.0),
                     test_t("mul(1,2,3,4)",24.0),
                     test_t("mul(1,2,3,4,5)",120.0),
                     test_t("mul(1,2,3,4,5,6)",720.0),
                     test_t("mul(1.1,2.2)",2.42),
                     test_t("mul(1.1,2.2,3.3)",7.986),
                     test_t("mul(1.1,2.2,3.3,4.4)",35.1384),
                     test_t("mul(1.1,2.2,3.3,4.4,5.5)",193.2612),
                     test_t("mul(1.1,2.2,3.3,4.4,5.5,6.6)",1275.52392),
                     test_t("equal(sum(1.1,2.2,3.3,4.4,5.5,6.6,7.7,8.8,9.9),(1.1+2.2+3.3+4.4+5.5+6.6+7.7+8.8+9.9))",1.0),
                     test_t("equal(mul(1.1,2.2,3.3,4.4,5.5,6.6,7.7,8.8,9.9),(1.1*2.2*3.3*4.4*5.5*6.6*7.7*8.8*9.9))",1.0),
                     test_t("equal(min(1.1,2.2,3.3,4.4,5.5,6.6,7.7,8.8,9.9),1.1)",1.0),
                     test_t("equal(max(1.1,2.2,3.3,4.4,5.5,6.6,7.7,8.8,9.9),9.9)",1.0),
                     test_t("equal(avg(1.1,2.2,3.3,4.4,5.5,6.6,7.7,8.8,9.9),5.5)",1.0),
                     test_t("floor(1.0)",1.0),
                     test_t("floor(1.1)",1.0),
                     test_t("floor(-1.0)",-1.0),
                     test_t("floor(-1.1)",-2.0),
                     test_t("ceil(1.0)",1.0),
                     test_t("ceil(1.1)",2.0),
                     test_t("ceil(-1.0)",-1.0),
                     test_t("ceil(-1.1)",-1.0),
                     test_t("round(1.1)",1.0),
                     test_t("round(1.49)",1.0),
                     test_t("round(1.5)",2.0),
                     test_t("round(1.9)",2.0),
                     test_t("roundn(1/3,2)",0.33),
                     test_t("roundn(1/3,5)",0.33333),
                     test_t("roundn(2/3,2)",0.67),
                     test_t("roundn(2/3,5)",0.66667),
                     test_t("roundn(1.0/3.0,2.0)",0.33),
                     test_t("roundn(1.0/3.0,5.0)",0.33333),
                     test_t("roundn(2.0/3.0,2.0)",0.67),
                     test_t("roundn(2.0/3.0,5.0)",0.66667),
                     test_t("cos(0.0)",1.0),
                     test_t("sin(0.0)",0.0),
                     test_t("equal(sin(pi/4.0),cos(pi/4.0))",1.0),
                     test_t("equal(sin(pi/6.0),cos(pi/3.0))",1.0),
                     test_t("(sin(pi/4.0) - cos(pi/4.0)) <= epsilon",1.0),
                     test_t("(cos(pi/3.0) - sin(pi/6.0)) <= epsilon",1.0),
                     test_t("sin(deg2rad(30))",0.5),
                     test_t("cos(deg2rad(60))",0.5),
                     test_t("sin(deg2rad(30)) + cos(deg2rad(60))",1.0),
                     test_t("equal(sin(deg2rad(30))/cos(deg2rad(30)),tan(deg2rad(30)))",1.0),
                     test_t("exp(1.0)",2.71828182845904523536028747135266249775724),
                     test_t("exp(0.0)",1.0),
                     test_t("log(2.7182818284590451)",1.0),
                     test_t("log10(10.0)",1.0),
                     test_t("frac(12.34) + trunc(12.34)",12.34),
                     test_t("hypot(3.0,4.0)",5.0),
                     test_t("hypot(1.0,sqrt(3.0))",2.0),
                     test_t("if(1 < 2, 3, 4)",3.0),
                     test_t("if(1.1 < 2.2, 3.3, 4.4)",3.3),
                     test_t("if((1.0+1.1) < (2.0+1.2), 3.3, 4.4)",3.3),
                     test_t("if(1 = 2, 3, 4)",4.0),
                     test_t("if(1.1 = 2.2, 3.3, 4.4)",4.4),
                     test_t("if((1.0+1.1) = (2.0+1.2), 3.3, 4.4)",4.4),
                     test_t("if(1 == 2, 3, 4)",4.0),
                     test_t("if(1.1 == 2.2, 3.3, 4.4)",4.4),
                     test_t("if((1.0+1.1) == (2.0+1.2), 3.3, 4.4)",4.4),
                     test_t("if(1 >= 2, 3, 4)",4.0),
                     test_t("if(1.1 >= 2.2, 3.3, 4.4)",4.4),
                     test_t("if((1.0+1.1) >= (2.0+1.2), 3.3, 4.4)",4.4),
                     test_t("if(((1.0 + 2.0) == 3.0) and ((4.0 + 5.0) < 9.0),1,2)",2.0),
                     test_t("(3.0 - 1.0 - 2.0) == ((3.0 - 1.0) - 2.0)",1.0),
                     test_t("true == true",1.0),
                     test_t("false == false",1.0),
                     test_t("true != false",1.0),
                     test_t("false != true",1.0),
                     test_t("(1 < 2) == true",1.0),
                     test_t("(1 > 2) == false",1.0),
                     test_t("true == (1 < 2)",1.0),
                     test_t("false == (1 > 2)",1.0),
                     test_t("(1 > 2) != true",1.0),
                     test_t("(1 < 2) != false",1.0),
                     test_t("true != (1 > 2)",1.0),
                     test_t("false != (1 < 2)",1.0),
                     test_t("(true and true) == true",1.0),
                     test_t("(false and false) == false",1.0),
                     test_t("(true or true) == true",1.0),
                     test_t("(false or false) == false",1.0),
                     test_t("(true and false) == false",1.0),
                     test_t("(false and true) == false",1.0),
                     test_t("(true or false) == true",1.0),
                     test_t("(false or true) == true",1.0),
                     test_t("(true & true) == true",1.0),
                     test_t("(false & false) == false",1.0),
                     test_t("(true | true) == true",1.0),
                     test_t("(false | false) == false",1.0),
                     test_t("(true & false) == false",1.0),
                     test_t("(false & true) == false",1.0),
                     test_t("(true | false) == true",1.0),
                     test_t("(false | true) == true",1.0),
                     test_t("clamp(-1,1,+1)",1.0),
                     test_t("clamp(-1,-1.5,+1.0)",-1.0),
                     test_t("clamp(-1,+1.5,+1.0)",+1.0),
                     test_t("clamp(-1,-1.5,+1.0) + clamp(-1,+1.5,+1.0)",0.0),
                     test_t("inrange(-2,1,+2) == ((-2 <= 1) and (1 <= +2))",1.0),
                     test_t("inrange(-2,1,+2) == if(({-2 <= 1} and [1 <= +2]),1.0,0.0)",1.0),
                     test_t("sgn( 0)", 0.0),
                     test_t("sgn(+3)",+1.0),
                     test_t("sgn(-3)",-1.0),
                     test_t("equal($f00(1.1,2.2,3.3),(1.1+2.2)/3.3)",1.0),
                     test_t("equal($f01(1.1,2.2,3.3),(1.1+2.2)*3.3)",1.0),
                     test_t("equal($f02(1.1,2.2,3.3),(1.1+2.2)-3.3)",1.0),
                     test_t("equal($f03(1.1,2.2,3.3),(1.1+2.2)+3.3)",1.0),
                     test_t("equal($f04(1.1,2.2,3.3),(1.1-2.2)/3.3)",1.0),
                     test_t("equal($f05(1.1,2.2,3.3),(1.1-2.2)*3.3)",1.0),
                     test_t("equal($f06(1.1,2.2,3.3),(1.1*2.2)+3.3)",1.0),
                     test_t("equal($f07(1.1,2.2,3.3),(1.1*2.2)-3.3)",1.0),
                     test_t("equal($f08(1.1,2.2,3.3),(1.1*2.2)/3.3)",1.0),
                     test_t("equal($f09(1.1,2.2,3.3),(1.1*2.2)*3.3)",1.0),
                     test_t("equal($f10(1.1,2.2,3.3),(1.1/2.2)+3.3)",1.0),
                     test_t("equal($f11(1.1,2.2,3.3),(1.1/2.2)-3.3)",1.0),
                     test_t("equal($f12(1.1,2.2,3.3),(1.1/2.2)/3.3)",1.0),
                     test_t("equal($f13(1.1,2.2,3.3),(1.1/2.2)*3.3)",1.0),
                     test_t("equal($f14(1.1,2.2,3.3),1.1/(2.2+3.3))",1.0),
                     test_t("equal($f15(1.1,2.2,3.3),1.1/(2.2-3.3))",1.0),
                     test_t("equal($f16(1.1,2.2,3.3),1.1/(2.2*3.3))",1.0),
                     test_t("equal($f17(1.1,2.2,3.3),1.1/(2.2/3.3))",1.0),
                     test_t("equal($f18(1.1,2.2,3.3),1.1*(2.2+3.3))",1.0),
                     test_t("equal($f19(1.1,2.2,3.3),1.1*(2.2-3.3))",1.0),
                     test_t("equal($f20(1.1,2.2,3.3),1.1*(2.2*3.3))",1.0),
                     test_t("equal($f21(1.1,2.2,3.3),1.1*(2.2/3.3))",1.0),
                     test_t("equal($f22(1.1,2.2,3.3),1.1-(2.2+3.3))",1.0),
                     test_t("equal($f23(1.1,2.2,3.3),1.1-(2.2-3.3))",1.0),
                     test_t("equal($f24(1.1,2.2,3.3),1.1-(2.2/3.3))",1.0),
                     test_t("equal($f25(1.1,2.2,3.3),1.1-(2.2*3.3))",1.0),
                     test_t("equal($f26(1.1,2.2,3.3),1.1+(2.2*3.3))",1.0),
                     test_t("equal($f27(1.1,2.2,3.3),1.1+(2.2/3.3))",1.0),
                     test_t("equal($f28(1.1,2.2,3.3),1.1+(2.2+3.3))",1.0),
                     test_t("equal($f29(1.1,2.2,3.3),1.1+(2.2-3.3))",1.0),
                     test_t("equal($f30(1.1,2.2,3.3),1.1*2.2^2+3.3)",1.0),
                     test_t("equal($f31(1.1,2.2,3.3),1.1*2.2^3+3.3)",1.0),
                     test_t("equal($f32(1.1,2.2,3.3),1.1*2.2^4+3.3)",1.0),
                     test_t("equal($f33(1.1,2.2,3.3),1.1*2.2^5+3.3)",1.0),
                     test_t("equal($f34(1.1,2.2,3.3),1.1*2.2^6+3.3)",1.0),
                     test_t("equal($f35(1.1,2.2,3.3),1.1*2.2^7+3.3)",1.0),
                     test_t("equal($f36(1.1,2.2,3.3),1.1*2.2^8+3.3)",1.0),
                     test_t("equal($f37(1.1,2.2,3.3),1.1*2.2^9+3.3)",1.0),
                     test_t("equal($f38(1.1,2.2,3.3),1.1*log(2.2)+3.3)",1.0),
                     test_t("equal($f39(1.1,2.2,3.3),1.1*log(2.2)-3.3)",1.0),
                     test_t("equal($f40(1.1,2.2,3.3),1.1*log10(2.2)+3.3)",1.0),
                     test_t("equal($f41(1.1,2.2,3.3),1.1*log10(2.2)-3.3)",1.0),
                     test_t("equal($f42(1.1,2.2,3.3),1.1*sin(2.2)+3.3)",1.0),
                     test_t("equal($f43(1.1,2.2,3.3),1.1*sin(2.2)-3.3)",1.0),
                     test_t("equal($f44(1.1,2.2,3.3),1.1*cos(2.2)+3.3)",1.0),
                     test_t("equal($f45(1.1,2.2,3.3),1.1*cos(2.2)-3.3)",1.0),
                     test_t("equal($f46(1.1,2.2,3.3),if(0!=1.1,2.2,3.3))",1.0),
                     test_t("equal($f47(1.1,2.2,3.3,4.4),1.1+((2.2+3.3)/4.4))",1.0),
                     test_t("equal($f48(1.1,2.2,3.3,4.4),1.1+((2.2+3.3)*4.4))",1.0),
                     test_t("equal($f49(1.1,2.2,3.3,4.4),1.1+((2.2-3.3)/4.4))",1.0),
                     test_t("equal($f50(1.1,2.2,3.3,4.4),1.1+((2.2-3.3)*4.4))",1.0),
                     test_t("equal($f51(1.1,2.2,3.3,4.4),1.1+((2.2*3.3)/4.4))",1.0),
                     test_t("equal($f52(1.1,2.2,3.3,4.4),1.1+((2.2*3.3)*4.4))",1.0),
                     test_t("equal($f53(1.1,2.2,3.3,4.4),1.1+((2.2/3.3)+4.4))",1.0),
                     test_t("equal($f54(1.1,2.2,3.3,4.4),1.1+((2.2/3.3)/4.4))",1.0),
                     test_t("equal($f55(1.1,2.2,3.3,4.4),1.1+((2.2/3.3)*4.4))",1.0),
                     test_t("equal($f56(1.1,2.2,3.3,4.4),1.1-((2.2+3.3)/4.4))",1.0),
                     test_t("equal($f57(1.1,2.2,3.3,4.4),1.1-((2.2+3.3)*4.4))",1.0),
                     test_t("equal($f58(1.1,2.2,3.3,4.4),1.1-((2.2-3.3)/4.4))",1.0),
                     test_t("equal($f59(1.1,2.2,3.3,4.4),1.1-((2.2-3.3)*4.4))",1.0),
                     test_t("equal($f60(1.1,2.2,3.3,4.4),1.1-((2.2*3.3)/4.4))",1.0),
                     test_t("equal($f61(1.1,2.2,3.3,4.4),1.1-((2.2*3.3)*4.4))",1.0),
                     test_t("equal($f62(1.1,2.2,3.3,4.4),1.1-((2.2/3.3)/4.4))",1.0),
                     test_t("equal($f63(1.1,2.2,3.3,4.4),1.1-((2.2/3.3)*4.4))",1.0),
                     test_t("equal($f64(1.1,2.2,3.3,4.4),((1.1+2.2)*3.3)-4.4)",1.0),
                     test_t("equal($f65(1.1,2.2,3.3,4.4),((1.1-2.2)*3.3)-4.4)",1.0),
                     test_t("equal($f66(1.1,2.2,3.3,4.4),((1.1*2.2)*3.3)-4.4)",1.0),
                     test_t("equal($f67(1.1,2.2,3.3,4.4),((1.1/2.2)*3.3)-4.4)",1.0),
                     test_t("equal($f68(1.1,2.2,3.3,4.4),((1.1+2.2)/3.3)-4.4)",1.0),
                     test_t("equal($f69(1.1,2.2,3.3,4.4),((1.1-2.2)/3.3)-4.4)",1.0),
                     test_t("equal($f70(1.1,2.2,3.3,4.4),((1.1*2.2)/3.3)-4.4)",1.0),
                     test_t("equal($f71(1.1,2.2,3.3,4.4),((1.1/2.2)/3.3)-4.4)",1.0),
                     test_t("equal($f72(1.1,2.2,3.3,4.4),(1.1*2.2)+(3.3*4.4))",1.0),
                     test_t("equal($f73(1.1,2.2,3.3,4.4),(1.1*2.2)-(3.3*4.4))",1.0),
                     test_t("equal($f74(1.1,2.2,3.3,4.4),(1.1*2.2)+(3.3/4.4))",1.0),
                     test_t("equal($f75(1.1,2.2,3.3,4.4),(1.1*2.2)-(3.3/4.4))",1.0),
                     test_t("equal($f76(1.1,2.2,3.3,4.4),(1.1/2.2)+(3.3/4.4))",1.0),
                     test_t("equal($f77(1.1,2.2,3.3,4.4),(1.1/2.2)-(3.3/4.4))",1.0),
                     test_t("equal($f78(1.1,2.2,3.3,4.4),(1.1/2.2)-(3.3*4.4))",1.0),
                     test_t("equal($f79(1.1,2.2,3.3,4.4),1.1/(2.2+(3.3*4.4)))",1.0),
                     test_t("equal($f80(1.1,2.2,3.3,4.4),1.1/(2.2-(3.3*4.4)))",1.0),
                     test_t("equal($f81(1.1,2.2,3.3,4.4),1.1*(2.2+(3.3*4.4)))",1.0),
                     test_t("equal($f82(1.1,2.2,3.3,4.4),1.1*(2.2-(3.3*4.4)))",1.0),
                     test_t("equal($f83(1.1,2.2,3.3,4.4),1.1*2.2^2+3.3*4.4^2)",1.0),
                     test_t("equal($f84(1.1,2.2,3.3,4.4),1.1*2.2^3+3.3*4.4^3)",1.0),
                     test_t("equal($f85(1.1,2.2,3.3,4.4),1.1*2.2^4+3.3*4.4^4)",1.0),
                     test_t("equal($f86(1.1,2.2,3.3,4.4),1.1*2.2^5+3.3*4.4^5)",1.0),
                     test_t("equal($f87(1.1,2.2,3.3,4.4),1.1*2.2^6+3.3*4.4^6)",1.0),
                     test_t("equal($f88(1.1,2.2,3.3,4.4),1.1*2.2^7+3.3*4.4^7)",1.0),
                     test_t("equal($f89(1.1,2.2,3.3,4.4),1.1*2.2^8+3.3*4.4^8)",1.0),
                     test_t("equal($f90(1.1,2.2,3.3,4.4),1.1*2.2^9+3.3*4.4^9)",1.0),
                     test_t("equal($f91(1.1,2.2,3.3,4.4),if(1.1 and 2.2,3.3,4.4))",1.0),
                     test_t("equal($f92(1.1,2.2,3.3,4.4),if(1.1 or 2.2,3.3,4.4))",1.0),
                     test_t("equal($f93(1.1,2.2,3.3,4.4),if(1.1 < 2.2,3.3,4.4))",1.0),
                     test_t("equal($f94(1.1,2.2,3.3,4.4),if(1.1 <= 2.2,3.3,4.4))",1.0),
                     test_t("equal($f95(1.1,2.2,3.3,4.4),if(1.1 > 2.2,3.3,4.4))",1.0),
                     test_t("equal($f96(1.1,2.2,3.3,4.4),if(1.1 >= 2.2,3.3,4.4))",1.0),
                     test_t("equal($f97(1.1,2.2,3.3,4.4),if(equal(1.1,2.2),3.3,4.4))",1.0),
                     test_t("equal($f98(1.1,2.2,3.3,4.4),1.1*sin(2.2)+3.3*cos(4.4))",1.0),
                     test_t("equal((48.0/2.0*(9.0+3.0)),288.0)",1.0),
                     test_t("equal((48.0/2.0(9.0+3.0)),288.0)",1.0),
                     test_t("equal((6.0/2.0(1.0+2.0)),9.0)",1.0),
                     test_t("1+2+3+4+5+6+7+8+9+0",45.0),
                     test_t("1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 9 + 0",45.0),
                     test_t("1.0 + 2.0 + 3.0 + 4.0 + 5.0 + 6.0 + 7.0 + 8.0 + 9.0 + 0.0",45.0),
                     test_t("(1+2)+(3+4)+(5+6)+(7+8)+(9+0)",45.0),
                     test_t("(1-2)+(3-4)+(5-6)+(7-8)+(9-0)",+5.0),
                     test_t("(1+2)-(3+4)-(5+6)-(7+8)-(9+0)",-39.0),
                     test_t("(1.0+2.0)+(3.0+4.0)+(5.0+6.0)+(7.0+8.0)+(9.0+0.0)",45.0),
                     test_t("(1.0-2.0)+(3.0-4.0)+(5.0-6.0)+(7.0-8.0)+(9.0-0.0)",+5.0),
                     test_t("(1.0+2.0)-(3.0+4.0)-(5.0+6.0)-(7.0+8.0)-(9.0+0.0)",-39.0),
                     test_t("[(1.0+2.0)+[3.0+4.0]+(5.0+6.0)]+([7.0+8.0]+(9.0+0.0))",45.0),
                     test_t("([1.0-2.0]+(3.0-4.0)+[5.0-6.0])+[(7.0-8.0)+[9.0-0.0]]",+5.0),
                     test_t("((1.0+2.0))-[(3.0+4.0)]-((5.0+6.0))-[(7.0+8.0)]-((9.0+0.0))",-39.0),
                     test_t("{[(1.0+2.0)+[3.0+4.0]+({5.0+6.0})]}+({[7.0+8.0]+(9.0+0.0)})",45.0),
                     test_t("{([1.0-2.0]+(3.0-4.0)+[5.0-6.0])}+[({+7.0}-{+8.0})+[{+9.0-0.0}]]",+5.0),
                     test_t("((+1.0+2.0))-[({+3.0+4.0})]-(({+5.0+6.0}))-[({+7.0}+8.0)]-(({+9.0}+{0.0}))",-39.0),
                     test_t("1+2-3*4/5+6-7*8/9+0",0.37777777777777777778),
                     test_t("1.1+2.2-3.3*4.4/5.5+6.6-7.7*8.8/9.9+0.0",0.41555555555555555556),
                     test_t("(1+2)-(3*4)/(5+6)-(7*8)/(9+0)",-4.31313131313131313131),
                     test_t("1/1+1/2+1/3+1/4+1/5+1/6+1/7+1/8+1/9",2.82896825396825396825),
                     test_t("(1/1)+(1/2)+(1/3)+(1/4)+(1/5)+(1/6)+(1/7)+(1/8)+(1/9)",2.82896825396825396825),
                     test_t("1.0/1.0+1.0/2.0+1.0/3.0+1.0/4.0+1.0/5.0+1.0/6.0+1.0/7.0+1.0/8.0+1.0/9",2.82896825396825396825),
                     test_t("(1.0/1.0)+(1.0/2.0)+(1.0/3.0)+(1.0/4.0)+(1.0/5.0)+(1.0/6.0)+(1.0/7.0)+(1.0/8.0)+(1.0/9)",2.82896825396825396825),
                     test_t("1/1*1/2*1/3*1/4*1/5*1/6*1/7*1/8*1/9",0.00000275573192239859),
                     test_t("(1/1)*(1/2)*(1/3)*(1/4)*(1/5)*(1/6)*(1/7)*(1/8)*(1/9)",0.00000275573192239859),
                     test_t("1.0/1.0*1.0/2.0*1.0/3.0*1.0/4.0*1.0/5.0*1.0/6.0*1.0/7.0*1.0/8.0*1.0/9",0.00000275573192239859),
                     test_t("(1.0/1.0)*(1.0/2.0)*(1.0/3.0)*(1.0/4.0)*(1.0/5.0)*(1.0/6.0)*(1.0/7.0)*(1.0/8.0)*(1.0/9)",0.00000275573192239859),
                     test_t("equal(poly01(1.2345,2.2,1.1),(2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly02(1.2345,3.3,2.2,1.1),(3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly03(1.2345,4.4,3.3,2.2,1.1),(4.4*1.2345^3+3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly04(1.2345,5.5,4.4,3.3,2.2,1.1),(5.5*1.2345^4+4.4*1.2345^3+3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly05(1.2345,6.6,5.5,4.4,3.3,2.2,1.1),(6.6*1.2345^5+5.5*1.2345^4+4.4*1.2345^3+3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly06(1.2345,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(7.7*1.2345^6+6.6*1.2345^5+5.5*1.2345^4+4.4*1.2345^3+3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly07(1.2345,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(8.8*1.2345^7+7.7*1.2345^6+6.6*1.2345^5+5.5*1.2345^4+4.4*1.2345^3+3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly08(1.2345,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(9.9*1.2345^8+8.8*1.2345^7+7.7*1.2345^6+6.6*1.2345^5+5.5*1.2345^4+4.4*1.2345^3+3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly09(1.2345,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(1.1*1.2345^9+9.9*1.2345^8+8.8*1.2345^7+7.7*1.2345^6+6.6*1.2345^5+5.5*1.2345^4+4.4*1.2345^3+3.3*1.2345^2+2.2*1.2345^1+1.1))",1.0),
                     test_t("equal(poly10(1.37,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(2.2*1.37^10+1.1*1.37^9+9.9*1.37^8+8.8*1.37^7+7.7*1.37^6+6.6*1.37^5+5.5*1.37^4+4.4*1.37^3+3.3*1.37^2+2.2*1.37^1+1.1))",1.0),
                     test_t("equal(poly11(1.37,3.3,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(3.3*1.37^11+2.2*1.37^10+1.1*1.37^9+9.9*1.37^8+8.8*1.37^7+7.7*1.37^6+6.6*1.37^5+5.5*1.37^4+4.4*1.37^3+3.3*1.37^2+2.2*1.37^1+1.1))",1.0),
                     test_t("equal(poly12(1.37,4.4,3.3,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(4.4*1.37^12+3.3*1.37^11+2.2*1.37^10+1.1*1.37^9+9.9*1.37^8+8.8*1.37^7+7.7*1.37^6+6.6*1.37^5+5.5*1.37^4+4.4*1.37^3+3.3*1.37^2+2.2*1.37^1+1.1))",1.0),
                     test_t("equal(\t \n(\n \r1.1\t\t - \n\n 2.2\n\n/\r3.3\t),(1.1-2.2/3.3))",1.0),
                     test_t("equal((pi^2^3),(pi^8))",1.0),
                     test_t("equal((pi^(2^3)),(pi^8))",1.0),
                     test_t("equal(pi^2^3-pi^8,0)",1.0),
                     test_t("equal((2*pi^2^3),2*(pi^8))",1.0),
                     test_t("equal((pi^2^3*2),2*(pi^8))",1.0),
                     test_t("equal((pi^2^3/2),(pi^8)/2)",1.0),
                     test_t("equal((pi^2.2^3.3),(pi^13.4894687605338489))",1.0),
                     test_t("equal((pi^(2.2^3.3)),(pi^13.4894687605338489))",1.0),
                     test_t("equal((2.2*pi^2.2^3.3),2.2*(pi^13.4894687605338489))",1.0),
                     test_t("equal((pi^2.2^3.3*2),2*(pi^13.4894687605338489))",1.0),
                     test_t("equal((pi^2.2^3.3/2.2),(pi^13.4894687605338489)/2.2)",1.0),
                     test_t("equal((pi^-2^3),1/(pi^8))",1.0),
                     test_t("equal((pi^(-2^3)),1/(pi^8))",1.0),
                     test_t("equal((pi^2^-3),(pi^(1/8)))",1.0),
                     test_t("equal((pi^(2^-3)),(pi^(1/8)))",1.0),
                     test_t("equal((pi^-2^-3),1/(pi^(1/8)))",1.0),
                     test_t("equal((pi^(-2^-3)),1/(pi^(1/8)))",1.0),
                     test_t("equal((-pi^2^3),(-pi^8))",1.0),
                     test_t("equal((-pi^(2^3)),(-pi^8))",1.0),
                     test_t("equal(-pi^2^3--pi^8,0)",1.0),
                     test_t("equal((2*-pi^2^3),2*(-pi^8))",1.0),
                     test_t("equal((-pi^2^3*2),2*(-pi^8))",1.0),
                     test_t("equal((-pi^2^3/2),(-pi^8)/2)",1.0),
                     test_t("equal((-pi^2.2^3.3),(-pi^13.4894687605338489))",1.0),
                     test_t("equal((-pi^(2.2^3.3)),(-pi^13.4894687605338489))",1.0),
                     test_t("equal((2.2*-pi^2.2^3.3),2.2*(-pi^13.4894687605338489))",1.0),
                     test_t("equal((-pi^2.2^3.3*2),2*(-pi^13.4894687605338489))",1.0),
                     test_t("equal((-pi^2.2^3.3/2.2),(-pi^13.4894687605338489)/2.2)",1.0),
                     test_t("equal((-pi^-2^3),1/(-pi^8))",1.0),
                     test_t("equal((-pi^(-2^3)),1/(-pi^8))",1.0),
                     test_t("equal((-pi^2^-3),(-pi^(1/8)))",1.0),
                     test_t("equal((-pi^(2^-3)),(-pi^(1/8)))",1.0),
                     test_t("equal((-pi^-2^-3),1/(-pi^(1/8)))",1.0),
                     test_t("equal((-pi^(-2^-3)),1/(-pi^(1/8)))",1.0),
                     test_t("equal((+pi^+2^+3),(+pi^+8))",1.0),
                     test_t("equal((+pi^(2^3)),(+pi^+8))",1.0),
                     test_t("equal(+pi^+2^+3-+pi^+8,0)",1.0),
                     test_t("equal((+2*+pi^+2^+3),+2*(+pi^+8))",1.0),
                     test_t("equal((+pi^+2^+3*+2),+2*(+pi^+8))",1.0),
                     test_t("equal((+pi^+2^+3/+2),(+pi^+8)/+2)",1.0),
                     test_t("equal((+pi^+2.2^+3.3),(+pi^+13.4894687605338489))",1.0),
                     test_t("equal((+pi^(+2.2^+3.3)),(+pi^+13.4894687605338489))",1.0),
                     test_t("equal((+2.2*+pi^+2.2^+3.3),+2.2*(+pi^+13.4894687605338489))",1.0),
                     test_t("equal((+pi^+2.2^+3.3*+2),+2*(+pi^+13.4894687605338489))",1.0),
                     test_t("equal((+pi^+2.2^+3.3/+2.2),(+pi^+13.4894687605338489)/+2.2)",1.0),
                     test_t("equal((+pi^-2^3),1/(+pi^+8))",1.0),
                     test_t("equal((+pi^(-2^3)),1/(+pi^+8))",1.0),
                     test_t("equal((+pi^2^-3),(+pi^(+1/+8)))",1.0),
                     test_t("equal((+pi^(2^-3)),(+pi^(+1/+8)))",1.0),
                     test_t("equal((+pi^-2^-3),1/(+pi^(+1/+8)))",1.0),
                     test_t("equal((+pi^(-2^-3)),1/(+pi^(+1/+8)))",1.0),
                     test_t("equal((-pi^+2^+3),(-pi^+8))",1.0),
                     test_t("equal((-pi^(2^3)),(-pi^+8))",1.0),
                     test_t("equal(-pi^+2^+3--pi^+8,0)",1.0),
                     test_t("equal((+2*-pi^+2^+3),2*(-pi^+8))",1.0),
                     test_t("equal((-pi^+2^+3*2),2*(-pi^+8))",1.0),
                     test_t("equal((-pi^+2^+3/+2),(-pi^+8)/+2)",1.0),
                     test_t("equal((-pi^+2.2^+3.3),(-pi^+13.4894687605338489))",1.0),
                     test_t("equal((-pi^(2.2^3.3)),(-pi^+13.4894687605338489))",1.0),
                     test_t("equal((+2.2*-pi^+2.2^+3.3),2.2*(-pi^+13.4894687605338489))",1.0),
                     test_t("equal((-pi^+2.2^+3.3*2),2*(-pi^+13.4894687605338489))",1.0),
                     test_t("equal((-pi^+2.2^+3.3/+2.2),(-pi^+13.4894687605338489)/+2.2)",1.0),
                     test_t("equal((-pi^-2^3),1/(-pi^+8))",1.0),
                     test_t("equal((-pi^(-2^3)),1/(-pi^+8))",1.0),
                     test_t("equal((-pi^2^-3),(-pi^(+1/+8)))",1.0),
                     test_t("equal((-pi^(2^-3)),(-pi^(+1/+8)))",1.0),
                     test_t("equal((-pi^-2^-3),1/(-pi^(+1/+8)))",1.0),
                     test_t("equal((-pi^(-2^-3)),1/(-pi^(+1/+8)))",1.0)
                  };

static const std::size_t test_list_size = sizeof(test_list) / sizeof(test_t);

template <typename T>
inline bool not_equal_impl(const T& t1,
                           const T& t2,
                           const T& epsilon = 0.0000000001/*std::numeric_limits<T>::epsilon()*/)
{
   if (t1 != t1) return true;
   if (t2 != t2) return true;
   T diff = std::abs(t1 - t2);
   T eps_norm = (std::max(T(1.0),std::max(std::abs(t1),std::abs(t2))) * epsilon);
   return diff > eps_norm;
}

template <typename T>
inline bool not_equal(const T& t0, const T& t1,
                      const T& epsilon = T(0.0000000001))
{
   return not_equal_impl(t0,t1,epsilon);
}

inline bool not_equal(const float& t0, const float& t1, const float& epsilon = 0.000001f)
{
   return not_equal_impl(t0,t1,epsilon);
}

template <typename T>
inline bool test_expression(const std::string& expression_string, const T& expected_result)
{
   exprtk::symbol_table<T> symbol_table;
   symbol_table.add_constants();

   exprtk::polynomial<T, 1> poly01;
   exprtk::polynomial<T, 2> poly02;
   exprtk::polynomial<T, 3> poly03;
   exprtk::polynomial<T, 4> poly04;
   exprtk::polynomial<T, 5> poly05;
   exprtk::polynomial<T, 6> poly06;
   exprtk::polynomial<T, 7> poly07;
   exprtk::polynomial<T, 8> poly08;
   exprtk::polynomial<T, 9> poly09;
   exprtk::polynomial<T,10> poly10;
   exprtk::polynomial<T,11> poly11;
   exprtk::polynomial<T,12> poly12;

   symbol_table.add_function("poly01", poly01);
   symbol_table.add_function("poly02", poly02);
   symbol_table.add_function("poly03", poly03);
   symbol_table.add_function("poly04", poly04);
   symbol_table.add_function("poly05", poly05);
   symbol_table.add_function("poly06", poly06);
   symbol_table.add_function("poly07", poly07);
   symbol_table.add_function("poly08", poly08);
   symbol_table.add_function("poly09", poly09);
   symbol_table.add_function("poly10", poly10);
   symbol_table.add_function("poly11", poly11);
   symbol_table.add_function("poly12", poly12);

   exprtk::expression<T> expression;
   expression.register_symbol_table(symbol_table);

   {
      exprtk::parser<T> parser;

      if (!parser.compile(expression_string,expression))
      {
         printf("test_expression() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expression_string.c_str());
         return false;
      }
   }

   if (!exprtk::expression_helper<T>::is_head_constant(expression))
   {
      printf("test_expression() - Error: Expression did not compile to a constant!   Expression: %s\n",
             expression_string.c_str());
      return false;
   }

   T result = expression.value();

   if (not_equal(result,expected_result))
   {
      printf("test_expression() - Computation Error:  Expression: [%s]\tExpected: %19.15f\tResult: %19.15f\n",
             expression_string.c_str(),
             expected_result,
             result);
      return false;
   }

   return true;
}

template <typename T>
inline bool run_test00()
{
   const std::size_t rounds = 10;
   for (std::size_t r = 0; r < rounds; ++r)
   {
      for (std::size_t i = 0; i < test_list_size; ++i)
      {
         if (!test_expression<T>(test_list[i].first,T(test_list[i].second)))
         {
            return false;
         }
      }
   }
   return true;
}

template <typename T>
struct test_xy
{
   test_xy(std::string e, const T& v0, const T& v1, const T& r)
   : expr(e),
     x(v0),
     y(v1),
     result(r)
   {}

   std::string expr;
   T x;
   T y;
   T result;
};

template <typename T>
inline bool run_test01()
{
   static const test_xy<T> test_list[] =
                           {
                              test_xy<T>("x + y"             ,T(2.2),T(3.3),T(5.5  )),
                              test_xy<T>("x - y"             ,T(3.3),T(2.2),T(1.1  )),
                              test_xy<T>("x * y"             ,T(3.3),T(2.2),T(7.26 )),
                              test_xy<T>("x / y"             ,T(3.3),T(2.2),T(1.5  )),
                              test_xy<T>("(x + y) * (x + y)" ,T(2.2),T(3.3),T(30.25)),
                              test_xy<T>("(x + y) / (x + y)" ,T(2.2),T(3.3),T(1.0  )),
                              test_xy<T>("x + y > x and  x + y > y" ,T(2.2),T(3.3),T(1.0)),
                              test_xy<T>("1 + (x + y)"       ,T(2.2),T(3.3),T(6.5  )),
                              test_xy<T>("(x + y) - 1"       ,T(2.2),T(3.3),T(4.5  )),
                              test_xy<T>("1 + (x + y) * 2"   ,T(2.2),T(3.3),T(12.0  )),
                              test_xy<T>("2 * (x + y) - 1"   ,T(2.2),T(3.3),T(10.0  )),
                              test_xy<T>("y + (x + 1)"       ,T(2.2),T(3.3),T(6.5  )),
                              test_xy<T>("(x + 1) + y"       ,T(2.2),T(3.3),T(6.5  )),
                              test_xy<T>("2 * x"  ,T(2.2),T(0.0),T(4.4)),
                              test_xy<T>("x * 2"  ,T(2.2),T(0.0),T(4.4)),
                              test_xy<T>("1.1 + x",T(2.2),T(0.0),T(3.3)),
                              test_xy<T>("x + 1.1",T(2.2),T(0.0),T(3.3)),
                              test_xy<T>("x * 1 == x"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("1 * x == x"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("y * 1 == y"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("1 * y == y"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x * 0 == 0"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("0 * x == 0"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("y * 0 == 0"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("0 * y == 0"    ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x + 1 == 1 + x",T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("y + 1 == 1 + y",T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x + y == y + x",T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x * y == y * x",T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x < y"         ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("y > x"         ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x <= y"        ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("y >= x"        ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x + y > y"     ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x + y > x"     ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x * y > y"     ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("x * y > x"     ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("(x + y) > y"   ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("(x + y) > x"   ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("(x * y) > y"   ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("(x * y) > x"   ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("(2x + 3y)  == (2*x + 3*y)" ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("2(x +  y)  == (2*x + 2*y)" ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>(" (x +  y)3 == (3*x + 3*y)" ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("2x + 3y    ==   2*x + 3*y" ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("2(x +  y)  ==   2*x + 2*y" ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>(" (x +  y)3 ==   3*x + 3*y" ,T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("equal(x+y^3/7,x+(y*y*y)/7)",T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("equal(1-x^3+y^2*7,1-(x*x*x)+(y*y)*7)",T(2.0),T(3.0),T(1.0)),
                              test_xy<T>("equal( x^0,1)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^1,x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^2,x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^3,x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^4,x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^5,x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^6,x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^7,x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^8,x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^9,x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^10,x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^11,x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^12,x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^13,x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^14,x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^15,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^16,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^17,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^18,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^19,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^20,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^21,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^22,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^23,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^24,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^25,x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( y^0,1)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^1,y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^2,y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^3,y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^4,y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^5,y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^6,y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^7,y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^8,y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^9,y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^10,y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^11,y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^12,y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^13,y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^14,y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^15,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^16,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^17,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^18,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^19,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^20,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^21,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^22,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^23,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^24,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^25,y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( x^-0,1/1)",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-1,1/(x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-2,1/(x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-3,1/(x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-4,1/(x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-5,1/(x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-6,1/(x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-7,1/(x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-8,1/(x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( x^-9,1/(x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-10,1/(x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-11,1/(x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-12,1/(x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-13,1/(x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-14,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-15,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-16,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-17,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-18,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-19,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-20,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-21,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-22,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-23,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-24,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal(x^-25,1/(x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x*x))",T(12.34),T(0.0),T(1.0)),
                              test_xy<T>("equal( y^-0,1/1)",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-1,1/(y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-2,1/(y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-3,1/(y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-4,1/(y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-5,1/(y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-6,1/(y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-7,1/(y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-8,1/(y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal( y^-9,1/(y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-10,1/(y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-11,1/(y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-12,1/(y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-13,1/(y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-14,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-15,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-16,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-17,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-18,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-19,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-20,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-21,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-22,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-23,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-24,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("equal(y^-25,1/(y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y*y))",T(0.0),T(12.34),T(1.0)),
                              test_xy<T>("0 * (abs  (x) + acos (y) + asin (x) + atan (y))",T(1.0),T(1.0),T(0.0)),
                              test_xy<T>("0 * (ceil (x) + cos  (y) + cosh (x) + exp  (y))",T(1.0),T(1.0),T(0.0)),
                              test_xy<T>("0 * (floor(x) + log  (y) + log10(x) + round(y))",T(1.0),T(1.0),T(0.0)),
                              test_xy<T>("0 * (sin  (x) + sinh (y) + sqrt (x) + tan  (y))",T(1.0),T(1.0),T(0.0)),
                              test_xy<T>("0 * (sec  (x) + csc  (y) + tanh (x) + cot  (y))",T(1.0),T(1.0),T(0.0)),
                              test_xy<T>("0 * (erf  (x) + erfc (y) + sgn  (y) + frac (y))",T(1.0),T(1.0),T(0.0)),
                              test_xy<T>("0 * (deg2grad(x) + grad2deg(y) + rad2deg(x) + deg2rad(y))",T(1.0),T(1.0),T(0.0)),
                           };

   static const std::size_t test_list_size = sizeof(test_list) / sizeof(test_xy<T>);

   const std::size_t rounds = 60;

   for (std::size_t r = 0; r < rounds; ++r)
   {
      for (std::size_t i = 0; i < test_list_size; ++i)
      {
         test_xy<T>& test = const_cast<test_xy<T>&>(test_list[i]);

         exprtk::symbol_table<T> symbol_table;
         symbol_table.add_variable("x",test.x);
         symbol_table.add_variable("y",test.y);

         exprtk::expression<T> expression;
         expression.register_symbol_table(symbol_table);

         {
            exprtk::parser<T> parser;
            if (!parser.compile(test.expr,expression))
            {
               printf("run_test01() - Error: %s   Expression: %s\n",
                      parser.error().c_str(),
                      test.expr.c_str());
               return false;
            }
         }

         T result = expression.value();

         if (not_equal(result,test.result))
         {
            printf("run_test01() - Computation Error:  Expression: [%s]\tExpected: %19.15f\tResult: %19.15f\n",
                   test.expr.c_str(),
                   test.result,
                   result);
            return false;
         }
      }
   }

   return true;
}

template <typename T>
struct test_ab
{
   test_ab(std::string e, const std::string& v0, const std::string& v1, const T& r)
   : expr(e),
     a(v0),
     b(v1),
     c("ccc"),
     result(r)
   {}

   std::string expr;
   std::string a;
   std::string b;
   std::string c;
   T result;
};

template <typename T>
inline bool run_test02()
{
   static const test_ab<T> test_list[] =
                           {
                              test_ab<T>("'aaa' == 'aaa'"              ,"","",T(1.0)),
                              test_ab<T>("'aaa'  < 'bbb'"              ,"","",T(1.0)),
                              test_ab<T>("'aaa' <= 'bbb'"              ,"","",T(1.0)),
                              test_ab<T>("'bbb' >  'aaa'"              ,"","",T(1.0)),
                              test_ab<T>("'bbb' >= 'aaa'"              ,"","",T(1.0)),
                              test_ab<T>("'aaa' != 'aaa'"              ,"","",T(0.0)),
                              test_ab<T>("'aaa' != 'bbb'"              ,"","",T(1.0)),
                              test_ab<T>("'aaa' + '123' == 'aaa123'"   ,"","",T(1.0)),
                              test_ab<T>("'aaa123' == 'aaa' + '123'"   ,"","",T(1.0)),
                              test_ab<T>("('aaa' + '123') == 'aaa123'" ,"","",T(1.0)),
                              test_ab<T>("'aaa123' == ('aaa' + '123')" ,"","",T(1.0)),
                              test_ab<T>("'aaa' in 'aaa123'"           ,"","",T(1.0)),
                              test_ab<T>("'123' in 'aaa123'"           ,"","",T(1.0)),
                              test_ab<T>("'a123b' like '*123*'"        ,"","",T(1.0)),
                              test_ab<T>("'a123b' like '*123?'"        ,"","",T(1.0)),
                              test_ab<T>("'1XYZ2' ilike '*xyz*'"       ,"","",T(1.0)),
                              test_ab<T>("'1XYZ2' ilike '*xyz?'"       ,"","",T(1.0)),
                              test_ab<T>("inrange('aaa','bbb','ccc')"  ,"","",T(1.0)),
                              test_ab<T>("a == b"                      ,"aaa","aaa",T(1.0)),
                              test_ab<T>("a != b"                      ,"aaa","bbb",T(1.0)),
                              test_ab<T>("a <  b"                      ,"aaa","bbb",T(1.0)),
                              test_ab<T>("a <= b"                      ,"aaa","bbb",T(1.0)),
                              test_ab<T>("b >  a"                      ,"aaa","bbb",T(1.0)),
                              test_ab<T>("b >= a"                      ,"aaa","bbb",T(1.0)),
                              test_ab<T>("a in b"                      ,"aaa","aaa123",T(1.0)),
                              test_ab<T>("a in b"                      ,"123","aaa123",T(1.0)),
                              test_ab<T>("a == 'aaa'"                  ,"aaa","aaa",T(1.0)),
                              test_ab<T>("'aaa' == a"                  ,"aaa","aaa",T(1.0)),
                              test_ab<T>("a != 'bbb'"                  ,"aaa","bbb",T(1.0)),
                              test_ab<T>("'bbb' != a"                  ,"aaa","bbb",T(1.0)),
                              test_ab<T>("a <  'bbb'"                  ,"aaa","bbb",T(1.0)),
                              test_ab<T>("a <= 'bbb'"                  ,"aaa","bbb",T(1.0)),
                              test_ab<T>("'bbb' >  a"                  ,"aaa","bbb",T(1.0)),
                              test_ab<T>("'bbb' >= a"                  ,"aaa","bbb",T(1.0)),
                              test_ab<T>("a in 'aaa123'"               ,"aaa","aaa123",T(1.0)),
                              test_ab<T>("a in 'aaa123'"               ,"123","aaa123",T(1.0)),
                              test_ab<T>("'aaa' in b"                  ,"aaa","aaa123",T(1.0)),
                              test_ab<T>("'123' in b"                  ,"aaa","aaa123",T(1.0)),
                              test_ab<T>("(a < b) or (a == b)"         ,"aaa","bbb",T(1.0)),
                              test_ab<T>("(a == b) or (a < b)"         ,"aaa","bbb",T(1.0)),
                              test_ab<T>("(b > a) or (b == a)"         ,"aaa","bbb",T(1.0)),
                              test_ab<T>("(b == a) or (b > a)"         ,"aaa","bbb",T(1.0)),
                              test_ab<T>("(a < b) and (b > a)"         ,"aaa","bbb",T(1.0)),
                              test_ab<T>("a like '*123*'"              ,"a123b","",T(1.0)),
                              test_ab<T>("a like '*123?'"              ,"a123b","",T(1.0)),
                              test_ab<T>("'a123b' like b"              ,"a123b","*123*",T(1.0)),
                              test_ab<T>("'a123b' like b"              ,"a123b","*123?",T(1.0)),
                              test_ab<T>("a ilike '*xyz*'"             ,"1XYZ2","",T(1.0)),
                              test_ab<T>("a ilike '*xyz?'"             ,"1XYZ2","",T(1.0)),
                              test_ab<T>("'1XYZ2' ilike b"             ,"","*xyz*",T(1.0)),
                              test_ab<T>("'1XYZ2' ilike b"             ,"","*xyz?",T(1.0)),
                              test_ab<T>("inrange(a,'bbb',c)"          ,"aaa","bbb",T(1.0)),
                              test_ab<T>("inrange('aaa',b,'ccc')"      ,"aaa","bbb",T(1.0)),
                              test_ab<T>("inrange(a,b,c)"              ,"aaa","bbb",T(1.0)),
                              test_ab<T>("inrange(a,b,'ccc')"          ,"aaa","bbb",T(1.0)),
                              test_ab<T>("inrange('aaa',b,c)"          ,"aaa","bbb",T(1.0)),
                              test_ab<T>("inrange('aaa',b,c)"          ,"aaa","bbb",T(1.0)),
                              test_ab<T>("'!@#$%^&*([{}])-=' != ')]}{[(*&^%$#@!'","","",T(1.0)),
                              test_ab<T>("('!@#$%^&*([{}])-=') != (')]}{[(*&^%$#@!')","","",T(1.0)),
                              test_ab<T>("{[('a')]} == [{('a')}]","","",T(1.0)),
                              test_ab<T>("{[('!@#$%^&*([{}])-=')]} != [{(')]}{[(*&^%$#@!')}]","","",T(1.0)),
                              test_ab<T>("'!@#$%^&*([{}])-=' == '!@#$%^&*([{}])-='","","",T(1.0)),
                              test_ab<T>("('!@#$%^&*([{}])-=') == ('!@#$%^&*([{}])-=')","","",T(1.0)),
                              test_ab<T>("{[('!@#$%^&*([{}])-=')]} == [{('!@#$%^&*([{}])-=')}]","","",T(1.0))
                           };

   static const std::size_t test_list_size = sizeof(test_list) / sizeof(test_ab<T>);

   const std::size_t rounds = 50;
   for (std::size_t r = 0; r < rounds; ++r)
   {
      for (std::size_t i = 0; i < test_list_size; ++i)
      {
         test_ab<T>& test = const_cast<test_ab<T>&>(test_list[i]);

         exprtk::symbol_table<T> symbol_table;
         symbol_table.add_stringvar("a",test.a);
         symbol_table.add_stringvar("b",test.b);
         symbol_table.add_stringvar("c",test.c);

         exprtk::expression<T> expression;
         expression.register_symbol_table(symbol_table);

         {
            exprtk::parser<T> parser;

            if (!parser.compile(test.expr,expression))
            {
               printf("run_test02() - Error: %s   Expression: %s\n",
                      parser.error().c_str(),
                      test.expr.c_str());
               return false;
            }
         }

         T result = expression.value();

         if (not_equal(result,test.result))
         {
            printf("run_test02() - Computation Error:  Expression: [%s]\tExpected: %19.15f\tResult: %19.15f\n",
                   test.expr.c_str(),
                   test.result,
                   result);
            return false;
         }
      }

   }

   return true;
}

template <typename T>
inline bool run_test03()
{
   std::string expression_string = "a+b+c+d+e+f+g+h+i+j+k+l+m+n+o+p+q+r+s+t+u+v+w+x+y+z+"
                                   "A+B+C+D+E+F+G+H+I+J+K+L+M+N+O+P+Q+R+S+T+U+V+W+X+Y+Z+"
                                   "aa+bb+cc+dd+ee+ff+gg+hh+ii+jj+kk+ll+mm+nn+oo+pp+qq+rr+ss+tt+uu+vv+ww+xx+yy+zz+"
                                   "AA+BB+CC+DD+EE+FF+GG+HH+II+JJ+KK+LL+MM+NN+OO+PP+QQ+RR+SS+TT+UU+VV+WW+XX+YY+ZZ+"
                                   "Aa+Bb+Cc+Dd+Ee+Ff+Gg+Hh+Ii+Jj+Kk+Ll+Mm+Nn+Oo+Pp+Qq+Rr+Ss+Tt+Uu+Vv+Ww+Xx+Yy+Zz+"
                                   "a0+b1+c2+d3+e4+f5+g6+h7+i8+j9+k0+l1+m2+n3+o4+p5+q6+r7+s8+t9+u0+v1+w2+x3+y4+z5+"
                                   "A0+B1+C2+D3+E4+F5+G6+H7+I8+J9+K0+L1+M2+N3+O4+P5+Q6+R7+S8+T9+U0+V1+W2+X3+Y4+Z5+"
                                   "aa0+bb1+cc2+dd3+ee4+ff5+gg6+hh7+ii8+jj9+kk0+ll1+mm2+nn3+oo4+pp5+qq6+rr7+ss8+tt9+uu0+vv1+ww2+xx3+yy4+zz5+"
                                   "AA0+BB1+CC2+DD3+EE4+FF5+GG6+HH7+II8+JJ9+KK0+LL1+MM2+NN3+OO4+PP5+QQ6+RR7+SS8+TT9+UU0+VV1+WW2+XX3+YY4+ZZ5+"
                                   "Aa0+Bb1+Cc2+Dd3+Ee4+Ff5+Gg6+Hh7+Ii8+Jj9+Kk0+Ll1+Mm2+Nn3+Oo4+Pp5+Qq6+Rr7+Ss8+Tt9+Uu0+Vv1+Ww2+Xx3+Yy4+Zz5";

   static const std::string variable_list[] =
                               {
                                  "a", "b", "c", "d", "e", "f", "g", "h", "i", "j", "k", "l", "m", "n", "o", "p", "q",
                                  "r", "s", "t", "u", "v", "w", "x", "y", "z", "A", "B", "C", "D", "E", "F", "G", "H",
                                  "I", "J", "K", "L", "M", "N", "O", "P", "Q", "R", "S", "T", "U", "V", "W", "X", "Y",
                                  "Z", "aa", "bb", "cc", "dd", "ee", "ff", "gg", "hh", "ii", "jj", "kk", "ll", "mm",
                                  "nn", "oo", "pp", "qq", "rr", "ss", "tt", "uu", "vv", "ww", "xx", "yy", "zz", "AA",
                                  "BB", "CC", "DD", "EE", "FF", "GG", "HH", "II", "JJ", "KK", "LL", "MM", "NN", "OO",
                                  "PP", "QQ", "RR", "SS", "TT", "UU", "VV", "WW", "XX", "YY", "ZZ", "Aa", "Bb", "Cc",
                                  "Dd", "Ee", "Ff", "Gg", "Hh", "Ii", "Jj", "Kk", "Ll", "Mm", "Nn", "Oo", "Pp", "Qq",
                                  "Rr", "Ss", "Tt", "Uu", "Vv", "Ww", "Xx", "Yy", "Zz", "a0", "b1", "c2", "d3", "e4",
                                  "f5", "g6", "h7", "i8", "j9", "k0", "l1", "m2", "n3", "o4", "p5", "q6", "r7", "s8",
                                  "t9", "u0", "v1", "w2", "x3", "y4", "z5", "A0", "B1", "C2", "D3", "E4", "F5", "G6",
                                  "H7", "I8", "J9", "K0", "L1", "M2", "N3", "O4", "P5", "Q6", "R7", "S8", "T9", "U0",
                                  "V1", "W2", "X3", "Y4", "Z5", "aa0", "bb1", "cc2", "dd3", "ee4", "ff5", "gg6", "hh7",
                                  "ii8", "jj9", "kk0", "ll1", "mm2", "nn3", "oo4", "pp5", "qq6", "rr7", "ss8", "tt9",
                                  "uu0", "vv1", "ww2", "xx3", "yy4", "zz5", "AA0", "BB1", "CC2", "DD3", "EE4", "FF5",
                                  "GG6", "HH7", "II8", "JJ9", "KK0", "LL1", "MM2", "NN3", "OO4", "PP5", "QQ6", "RR7",
                                  "SS8", "TT9", "UU0", "VV1", "WW2", "XX3", "YY4", "ZZ5", "Aa0", "Bb1", "Cc2", "Dd3",
                                  "Ee4", "Ff5", "Gg6", "Hh7", "Ii8", "Jj9", "Kk0", "Ll1", "Mm2", "Nn3", "Oo4", "Pp5",
                                  "Qq6", "Rr7", "Ss8", "Tt9", "Uu0", "Vv1", "Ww2", "Xx3", "Yy4", "Zz5"
                               };

   static const std::size_t variable_list_size = sizeof(variable_list) / sizeof(std::string);

   static const std::size_t rounds = 300;

   for (std::size_t r = 0; r < rounds; ++r)
   {
      exprtk::symbol_table<T> symbol_table;
      exprtk::expression<T> expression;

      std::vector<T> v;
      v.resize(variable_list_size);

      for (std::size_t i = 0; i < variable_list_size; ++i)
      {
         v[i] = T(i);
         symbol_table.add_variable(variable_list[i],v[i]);
      }

      symbol_table.add_constants();

      expression.register_symbol_table(symbol_table);

      exprtk::parser<T> parser;

      if (!parser.compile(expression_string,expression))
      {
         printf("run_test03() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expression_string.c_str());
         return false;
      }

      expression.value();
   }

   return true;
}

template <typename T>
inline T clamp(const T& l, const T& v, const T& u)
{
   return (v < l) ? l : ((v > u) ? u : v);
}

template <typename T>
inline bool run_test04()
{
   std::string expression_string = "clamp(-1.0,sin(2 * pi * x) + cos(y / 2 * pi),+1.0)";

   exprtk::symbol_table<T> symbol_table;
   exprtk::expression<T> expression;

   T x = T(-1000.0);
   T y = T(-1000.0);

   symbol_table.add_variable("x",x);
   symbol_table.add_variable("y",y);
   symbol_table.add_constants();

   expression.register_symbol_table(symbol_table);

   {
      exprtk::parser<T> parser;

      if (!parser.compile(expression_string,expression))
      {
         printf("run_test04() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expression_string.c_str());
         return false;
      }
   }

   const T pi = T(3.14159265358979323846);
   const T increment = T(0.0001);

   while ((x <= T(+1000.0)) && (y <= T(+1000.0)))
   {
      T result1 = expression.value();
      T result2 = clamp<T>(-1.0,std::sin(2 * pi * x) + std::cos(y / 2 * pi),+1.0);

      if (not_equal(result1,result2))
      {
         printf("run_test04() - Computation Error:  Expression: [%s]\tExpected: %19.15f\tResult: %19.15f x:%19.15f\ty:%19.15f\n",
                expression_string.c_str(),
                result1,
                result2,
                x,
                y);
         return false;
      }

      x += increment;
      y += increment;
   }

   return true;
}

template <typename T>
inline bool run_test05()
{
   typedef exprtk::expression<T> expression_t;
   std::string expression_string = "clamp(-1.0,sin(2 * pi * x_var123) + cos(y_var123 / 2 * pi),+1.0)";

   exprtk::symbol_table<T> symbol_table;
   std::deque<expression_t> expression_list;

   T x = T(-1000.0);
   T y = T(-1000.0);

   symbol_table.add_variable("x_var123",x);
   symbol_table.add_variable("y_var123",y);
   symbol_table.add_constants();

   const std::size_t expression_count = 10;
   for (std::size_t i = 0; i < expression_count; ++i)
   {
      expression_t e;
      exprtk::parser<T> parser;
      e.register_symbol_table(symbol_table);

      if (!parser.compile(expression_string,e))
      {
         printf("run_test05() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expression_string.c_str());
         return false;
      }

      expression_list.push_back(e);
   }

   const T pi = T(3.14159265358979323846);
   const T increment = T(0.001);

   while ((x <= T(+1000.0)) && (y <= T(+1000.0)))
   {
      T real_result = clamp<T>(-1.0,std::sin(2 * pi * x) + std::cos(y / 2 * pi),+1.0);

      for (std::size_t i = 0; i < expression_list.size(); ++i)
      {
         expression_t& expr = expression_list[i];

         T result = expr.value();

         if (not_equal(result,real_result))
         {
            printf("run_test05() - Computation Error:  Expression: [%s]\tExpected: %19.15f\tResult: %19.15f x:%19.15f\ty:%19.15f\tIndex:%d\n",
                   expression_string.c_str(),
                   real_result,
                   result,
                   x,
                   y,
                   static_cast<unsigned int>(i));
            return false;
         }
      }

      x += increment;
      y += increment;
   }

   return true;
}

template <typename T>
inline bool run_test06()
{
   typedef exprtk::expression<T> expression_t;
   std::string expression_string = "sqrt(1 - (x^2))";

   T x = T(0.0);

   exprtk::symbol_table<T> symbol_table;
   symbol_table.add_variable("x",x);

   expression_t expression;
   expression.register_symbol_table(symbol_table);

   exprtk::parser<T> parser;

   if (!parser.compile(expression_string,expression))
   {
      printf("run_test06() - Error: %s   Expression: %s\n",
             parser.error().c_str(),
             expression_string.c_str());
      return false;
   }

   T total_area1 = exprtk::integrate(expression,x,T(-1.0),T(1.0));
   T total_area2 = exprtk::integrate(expression,"x",T(-1.0),T(1.0));
   const T pi = T(3.14159265358979323846);

   if (not_equal(total_area1,total_area2,T(0.000001)))
   {
      printf("run_test06() - Integration Error:  area1 != area2\n");
      return false;
   }

   if (not_equal(total_area1,T(pi)/T(2.0),T(0.000001)))
   {
      printf("run_test06() - Integration Error:  Expected: %19.15f\tResult: %19.15f\n",
             pi/T(2.0),
             total_area1);
      return false;
   }

   return true;
}

template <typename T>
inline bool run_test07()
{
   typedef exprtk::expression<T> expression_t;
   std::string expression_string = "sin(2x+1/3)";

   T x = T(0.0);

   exprtk::symbol_table<T> symbol_table;
   symbol_table.add_variable("x",x);

   expression_t expression;
   expression.register_symbol_table(symbol_table);

   exprtk::parser<T> parser;

   if (!parser.compile(expression_string,expression))
   {
      printf("run_test07() - Error: %s   Expression: %s\n",
             parser.error().c_str(),
             expression_string.c_str());

      return false;
   }

   for (x = T(-200.0); x < T(200); x += T(0.0001))
   {
      {
         T deriv1_real_result = T(2.0) * std::cos(T(2.0) * x + T(1.0/3.0));
         T deriv1_result1 = exprtk::derivative(expression,x);
         T deriv1_result2 = exprtk::derivative(expression,"x");

         if (not_equal(deriv1_result1,deriv1_result2,T(0.00001)))
         {
            printf("run_test07() - 1st Derivative Error:  result1 != result2\n");
            return false;
         }

         if (not_equal(deriv1_result1,deriv1_real_result,T(0.00001)))
         {
            printf("run_test07() - 1st Derivative Error:  x: %19.15f\tExpected: %19.15f\tResult: %19.15f\n",
                   x,
                   deriv1_real_result,
                   deriv1_result1);
            return false;
         }
      }

      {
         T deriv2_real_result = T(-4.0) * std::sin(T(2.0) * x + T(1.0/3.0));
         T deriv2_result1 = exprtk::second_derivative(expression,x);
         T deriv2_result2 = exprtk::second_derivative(expression,"x");

         if (not_equal(deriv2_result1,deriv2_result2,T(0.0000001)))
         {
            printf("run_test07() - 2nd Derivative Error:  result1 != result2\n");
            return false;
         }

         if (not_equal(deriv2_result1,deriv2_real_result,T(0.01)))
         {
            printf("run_test07() - 2nd Derivative Error:  x: %19.15f\tExpected: %19.15f\tResult: %19.15f\n",
                   x,
                   deriv2_real_result,
                   deriv2_result1);
            return false;
         }
      }

      {
         T deriv3_real_result = T(-8.0) * std::cos(T(2.0) * x + T(1.0/3.0));
         T deriv3_result1 = exprtk::third_derivative(expression,x);
         T deriv3_result2 = exprtk::third_derivative(expression,"x");

         if (not_equal(deriv3_result1,deriv3_result2,T(0.0000001)))
         {
            printf("run_test07() - 3rd Derivative Error:  result1 != result2\n");
            return false;
         }

         if (not_equal(deriv3_result1,deriv3_real_result,T(0.01)))
         {
            printf("run_test07() - 3rd Derivative Error:  x: %19.15f\tExpected: %19.15f\tResult: %19.15f\n",
                   x,
                   deriv3_real_result,
                   deriv3_result1);
            return false;
         }
      }
   }

   return true;
}

template <typename T>
inline bool run_test08()
{

   static const std::string expr_str[] =
                              {
                                 "x", "y", "z", "w", "u",
                                 "x + y + z + w + u",
                                 "x + y / z * w ^ u",
                                 "x:=1.1", "y:=2.2", "z:=3.3", "w:=4.4", "u:=5.5",
                                 "x:=x+1.1", "y:=y+2.2", "z:=z+3.3", "w:=w+4.4", "u:=u+5.5",
                                 "x:=1.1+x", "y:=2.2+y", "z:=3.3+z", "w:=4.4+w", "u:=5.5+u",
                                 "x:=(x <= 1.1)",
                                 "y:=(2.2 >= y)",
                                 "z:=(3.3 and z)",
                                 "w:=(4.4 or w)",
                                 "u:=(u xor 5.5)",
                                 "min(x,y) + min(x,y,z) + min(x,y,z,w) + min(x,y,z,w,y)",
                                 "max(x,y) + max(x,y,z) + max(x,y,z,w) + max(x,y,z,w,y)",
                                 "avg(x,y)",
                                 "avg(x,y,z)",
                                 "avg(x,y,z,w)",
                                 "avg(x,y,z,w,u)",
                                 "(u := u := min(x:=1,y:=2,z:=3)) == 1",
                                 "(2x+3y+4z+5w)==(2*x+3*y+4*z+5*w)",
                                 "(3(x+y)/2+1)==(3*(x+y)/2+1)",
                                 "((x+y)3+1/4)==((x+y)*3+1/4)",
                                 "((x+y)z+1/2)==((x+y)*z+1/2)",
                                 "(x+y^3/z) == (x+(y*y*y)/z)",
                                 "(z-x^3+y^2*7) == (z-(x*x*x)+(y*y)*7)",
                                 "(3min(x,y))==(3*min(x,y))",
                                 "(sin(x)y)==(sin(x)*y)",
                                 "(sin(x)cos(y)+1)==(sin(x)*cos(y)+1)",
                                 "(sgn(sin(x))cos(sgn(y))+1)==(sgn(sin(x))*cos(sgn(y))+1)",
                                 "equal($f00(x,y,z),(x+y)/z)",
                                 "equal($f01(x,y,z),(x+y)*z)",
                                 "equal($f02(x,y,z),(x+y)-z)",
                                 "equal($f03(x,y,z),(x+y)+z)",
                                 "equal($f04(x,y,z),(x-y)/z)",
                                 "equal($f05(x,y,z),(x-y)*z)",
                                 "equal($f06(x,y,z),(x*y)+z)",
                                 "equal($f07(x,y,z),(x*y)-z)",
                                 "equal($f08(x,y,z),(x*y)/z)",
                                 "equal($f09(x,y,z),(x*y)*z)",
                                 "equal($f10(x,y,z),(x/y)+z)",
                                 "equal($f11(x,y,z),(x/y)-z)",
                                 "equal($f12(x,y,z),(x/y)/z)",
                                 "equal($f13(x,y,z),(x/y)*z)",
                                 "equal($f14(x,y,z),x/(y+z))",
                                 "equal($f15(x,y,z),x/(y-z))",
                                 "equal($f16(x,y,z),x/(y*z))",
                                 "equal($f17(x,y,z),x/(y/z))",
                                 "equal($f18(x,y,z),x*(y+z))",
                                 "equal($f19(x,y,z),x*(y-z))",
                                 "equal($f20(x,y,z),x*(y*z))",
                                 "equal($f21(x,y,z),x*(y/z))",
                                 "equal($f22(x,y,z),x-(y+z))",
                                 "equal($f23(x,y,z),x-(y-z))",
                                 "equal($f24(x,y,z),x-(y/z))",
                                 "equal($f25(x,y,z),x-(y*z))",
                                 "equal($f26(x,y,z),x+(y*z))",
                                 "equal($f27(x,y,z),x+(y/z))",
                                 "equal($f28(x,y,z),x+(y+z))",
                                 "equal($f29(x,y,z),x+(y-z))",
                                 "equal($f30(x,y,z),x*y^2+z)",
                                 "equal($f31(x,y,z),x*y^3+z)",
                                 "equal($f32(x,y,z),x*y^4+z)",
                                 "equal($f33(x,y,z),x*y^5+z)",
                                 "equal($f34(x,y,z),x*y^6+z)",
                                 "equal($f35(x,y,z),x*y^7+z)",
                                 "equal($f36(x,y,z),x*y^8+z)",
                                 "equal($f37(x,y,z),x*y^9+z)",
                                 "equal($f38(x,y,z),x*log(y)+z)",
                                 "equal($f39(x,y,z),x*log(y)-z)",
                                 "equal($f40(x,y,z),x*log10(y)+z)",
                                 "equal($f41(x,y,z),x*log10(y)-z)",
                                 "equal($f42(x,y,z),x*sin(y)+z)",
                                 "equal($f43(x,y,z),x*sin(y)-z)",
                                 "equal($f44(x,y,z),x*cos(y)+z)",
                                 "equal($f45(x,y,z),x*cos(y)-z)",
                                 "equal($f46(x,y,z),if(0!=x,y,z))",
                                 "equal($f47(x,y,z,w),x+((y+z)/w))",
                                 "equal($f48(x,y,z,w),x+((y+z)*w))",
                                 "equal($f49(x,y,z,w),x+((y-z)/w))",
                                 "equal($f50(x,y,z,w),x+((y-z)*w))",
                                 "equal($f51(x,y,z,w),x+((y*z)/w))",
                                 "equal($f52(x,y,z,w),x+((y*z)*w))",
                                 "equal($f53(x,y,z,w),x+((y/z)+w))",
                                 "equal($f54(x,y,z,w),x+((y/z)/w))",
                                 "equal($f55(x,y,z,w),x+((y/z)*w))",
                                 "equal($f56(x,y,z,w),x-((y+z)/w))",
                                 "equal($f57(x,y,z,w),x-((y+z)*w))",
                                 "equal($f58(x,y,z,w),x-((y-z)/w))",
                                 "equal($f59(x,y,z,w),x-((y-z)*w))",
                                 "equal($f60(x,y,z,w),x-((y*z)/w))",
                                 "equal($f61(x,y,z,w),x-((y*z)*w))",
                                 "equal($f62(x,y,z,w),x-((y/z)/w))",
                                 "equal($f63(x,y,z,w),x-((y/z)*w))",
                                 "equal($f64(x,y,z,w),((x+y)*z)-w)",
                                 "equal($f65(x,y,z,w),((x-y)*z)-w)",
                                 "equal($f66(x,y,z,w),((x*y)*z)-w)",
                                 "equal($f67(x,y,z,w),((x/y)*z)-w)",
                                 "equal($f68(x,y,z,w),((x+y)/z)-w)",
                                 "equal($f69(x,y,z,w),((x-y)/z)-w)",
                                 "equal($f70(x,y,z,w),((x*y)/z)-w)",
                                 "equal($f71(x,y,z,w),((x/y)/z)-w)",
                                 "equal($f72(x,y,z,w),(x*y)+(z*w))",
                                 "equal($f73(x,y,z,w),(x*y)-(z*w))",
                                 "equal($f74(x,y,z,w),(x*y)+(z/w))",
                                 "equal($f75(x,y,z,w),(x*y)-(z/w))",
                                 "equal($f76(x,y,z,w),(x/y)+(z/w))",
                                 "equal($f77(x,y,z,w),(x/y)-(z/w))",
                                 "equal($f78(x,y,z,w),(x/y)-(z*w))",
                                 "equal($f79(x,y,z,w),x/(y+(z*w)))",
                                 "equal($f80(x,y,z,w),x/(y-(z*w)))",
                                 "equal($f81(x,y,z,w),x*(y+(z*w)))",
                                 "equal($f82(x,y,z,w),x*(y-(z*w)))",
                                 "equal($f83(x,y,z,w),x*y^2+z*w^2)",
                                 "equal($f84(x,y,z,w),x*y^3+z*w^3)",
                                 "equal($f85(x,y,z,w),x*y^4+z*w^4)",
                                 "equal($f86(x,y,z,w),x*y^5+z*w^5)",
                                 "equal($f87(x,y,z,w),x*y^6+z*w^6)",
                                 "equal($f88(x,y,z,w),x*y^7+z*w^7)",
                                 "equal($f89(x,y,z,w),x*y^8+z*w^8)",
                                 "equal($f90(x,y,z,w),x*y^9+z*w^9)",
                                 "equal($f91(x,y,z,w),if(x and y,z,w))",
                                 "equal($f92(x,y,z,w),if(x or y,z,w))",
                                 "equal($f93(x,y,z,w),if(x < y,z,w))",
                                 "equal($f94(x,y,z,w),if(x <= y,z,w))",
                                 "equal($f95(x,y,z,w),if(x > y,z,w))",
                                 "equal($f96(x,y,z,w),if(x >= y,z,w))",
                                 "equal($f97(x,y,z,w),if(equal(x,y),z,w))",
                                 "equal($f98(x,y,z,w),x*sin(y)+z*cos(w))"
                              };
   static const std::size_t expr_str_size = sizeof(expr_str) / sizeof(std::string);

   static const std::size_t rounds = 25;

   for (std::size_t i = 0; i < rounds; ++i)
   {
      for (std::size_t j = 0; j < expr_str_size; ++j)
      {
         typedef exprtk::expression<T> expression_t;

         T x = T(1.12345);
         T y = T(2.12345);
         T z = T(3.12345);
         T w = T(4.12345);
         T u = T(5.12345);

         exprtk::symbol_table<T> symbol_table;
         symbol_table.add_variable("x",x);
         symbol_table.add_variable("y",y);
         symbol_table.add_variable("z",z);
         symbol_table.add_variable("w",w);
         symbol_table.add_variable("u",u);

         expression_t expression;
         expression.register_symbol_table(symbol_table);

         exprtk::parser<T> parser;

         if (!parser.compile(expr_str[j],expression))
         {
            printf("run_test08() - Error: %s   Expression: %s\n",
                   parser.error().c_str(),
                   expr_str[j].c_str());
            return false;
         }

         expression.value();
      }
   }

   return true;
}

template <typename T>
struct myfunc : public exprtk::ifunction<T>
{
   myfunc() : exprtk::ifunction<T>(2) {}

   inline T operator()(const T& v1, const T& v2)
   {
      return T(1) + (v1 * v2) / T(3);
   }
};

template <typename T>
inline bool run_test09()
{
   static const std::size_t rounds = 1000;
   for (std::size_t i = 0; i < rounds; ++i)
   {
      typedef exprtk::expression<T> expression_t;
      std::string expression_string = "myfunc0(sin(x*pi),y/2)+myfunc1(sin(x*pi),y/2)+"
                                      "myfunc2(sin(x*pi),y/2)+myfunc3(sin(x*pi),y/2)+"
                                      "myfunc4(sin(x*pi),y/2)+myfunc5(sin(x*pi),y/2)+"
                                      "myfunc6(sin(x*pi),y/2)+myfunc7(sin(x*pi),y/2)+"
                                      "myfunc8(sin(x*pi),y/2)+myfunc9(sin(x*pi),y/2)+"
                                      "myfunc0(sin(x*pi),y/2)+myfunc1(sin(x*pi),y/2)+"
                                      "myfunc2(sin(x*pi),y/2)+myfunc3(sin(x*pi),y/2)+"
                                      "myfunc4(sin(x*pi),y/2)+myfunc5(sin(x*pi),y/2)+"
                                      "myfunc6(sin(x*pi),y/2)+myfunc7(sin(x*pi),y/2)+"
                                      "myfunc8(sin(x*pi),y/2)+myfunc9(sin(x*pi),y/2)+"
                                      "myfunc0(sin(x*pi),y/2)+myfunc1(sin(x*pi),y/2)+"
                                      "myfunc2(sin(x*pi),y/2)+myfunc3(sin(x*pi),y/2)+"
                                      "myfunc4(sin(x*pi),y/2)+myfunc5(sin(x*pi),y/2)+"
                                      "myfunc6(sin(x*pi),y/2)+myfunc7(sin(x*pi),y/2)+"
                                      "myfunc8(sin(x*pi),y/2)+myfunc9(sin(x*pi),y/2)+"
                                      "myfunc0(sin(x*pi),y/2)+myfunc1(sin(x*pi),y/2)+"
                                      "myfunc2(sin(x*pi),y/2)+myfunc3(sin(x*pi),y/2)+"
                                      "myfunc4(sin(x*pi),y/2)+myfunc5(sin(x*pi),y/2)+"
                                      "myfunc6(sin(x*pi),y/2)+myfunc7(sin(x*pi),y/2)+"
                                      "myfunc8(sin(x*pi),y/2)+myfunc9(sin(x*pi),y/2)";

      T x = T(1.0) + (i/T(10000.0));
      T y = T(2.0) + (i/T(10000.0));
      myfunc<T> mf;

      exprtk::symbol_table<T> symbol_table;
      symbol_table.add_variable("x",x);
      symbol_table.add_variable("y",y);
      symbol_table.add_function("myfunc0",mf);
      symbol_table.add_function("myfunc1",mf);
      symbol_table.add_function("myfunc2",mf);
      symbol_table.add_function("myfunc3",mf);
      symbol_table.add_function("myfunc4",mf);
      symbol_table.add_function("myfunc5",mf);
      symbol_table.add_function("myfunc6",mf);
      symbol_table.add_function("myfunc7",mf);
      symbol_table.add_function("myfunc8",mf);
      symbol_table.add_function("myfunc9",mf);
      symbol_table.add_constants();

      expression_t expression;
      expression.register_symbol_table(symbol_table);

      exprtk::parser<T> parser;

      if (!parser.compile(expression_string,expression))
      {
         printf("run_test09() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expression_string.c_str());
         return false;
      }

      const T pi = T(3.141592653589793238462);

      T result = expression.value();

      T expected = T(4.0) *
                   (
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0)) +
                      mf(sin(x*pi),y/T(2.0))
                   );

      if (not_equal(result,expected,T(0.0000001)))
      {
         printf("run_test09() - Error Expected: %19.15f\tResult: %19.15f\n",
                expected,
                result);
         return false;
      }
   }

   return true;
}

template <typename T>
inline bool run_test10()
{
   T  x = T(1.1);
   T  y = T(2.2);
   T xx = T(3.3);
   T yy = T(4.4);

   std::string i  = "A String";
   std::string j  = "Another String";
   std::string ii = "A String";
   std::string jj = "Another String";

   exprtk::symbol_table<T> symbol_table;

   typedef exprtk::expression<T> expression_t;

   struct test
   {
      static inline bool variable(exprtk::symbol_table<T>& symbol_table, const std::string& variable_name, const T& value)
      {
         exprtk::details::variable_node<T>* var = symbol_table.get_variable(variable_name);
         if (var)
            return (!not_equal(var->ref(),value));
         else
            return false;
      }

      static inline bool string(exprtk::symbol_table<T>& symbol_table, const std::string& string_name, const std::string& str)
      {
         exprtk::details::stringvar_node<T>* str_node = symbol_table.get_stringvar(string_name);
         if (str_node)
            return (str_node->ref() == str);
         else
            return false;
      }
   };

   static const std::size_t rounds = 10;

   for (std::size_t r = 0; r < rounds; ++r)
   {
      symbol_table.add_variable("x",  x);
      symbol_table.add_variable("y",  y);
      symbol_table.add_variable("xx",xx);
      symbol_table.add_variable("yy",yy);

      if (!symbol_table.symbol_exists("x"))
      {
         printf("run_test10() - Symbol 'x' does not exist!\n");
         return false;
      }
      else if (!symbol_table.symbol_exists("y"))
      {
         printf("run_test10() - Symbol 'y' does not exist!\n");
         return false;
      }
      else if (!symbol_table.symbol_exists("xx"))
      {
         printf("run_test10() - Symbol 'xx' does not exist!\n");
         return false;
      }
      else if (!symbol_table.symbol_exists("yy"))
      {
         printf("run_test10() - Symbol 'yy' does not exist!\n");
         return false;
      }
      else if (!test::variable(symbol_table,"x",x))
      {
         printf("run_test10() - Symbol 'x' value failure!\n");
         return false;
      }
      else if (!test::variable(symbol_table,"y",y))
      {
         printf("run_test10() - Symbol 'y' value failure!\n");
         return false;
      }
      else if (!test::variable(symbol_table,"xx",xx))
      {
         printf("run_test10() - Symbol 'xx' value failure!\n");
         return false;
      }
      else if (!test::variable(symbol_table,"yy",yy))
      {
         printf("run_test10() - Symbol 'yy' value failure!\n");
         return false;
      }

      if (!symbol_table.remove_variable("x"))
      {
         printf("run_test10() - Failed to remove symbol 'x'!\n");
         return false;
      }
      else if (!symbol_table.remove_variable("y"))
      {
         printf("run_test10() - Failed to remove symbol 'y'!\n");
         return false;
      }
      else if (!symbol_table.remove_variable("xx"))
      {
         printf("run_test10() - Failed to remove symbol 'xx'!\n");
         return false;
      }
      else if (!symbol_table.remove_variable("yy"))
      {
         printf("run_test10() - Failed to remove symbol 'yy'!\n");
         return false;
      }
   }

   for (std::size_t r = 0; r < rounds; ++r)
   {
      myfunc<T> mf;

      symbol_table.add_function("f",mf);
      symbol_table.add_function("f1",mf);

      if (!symbol_table.symbol_exists("f"))
      {
         printf("run_test10() - function 'f' does not exist!\n");
         return false;
      }
      else if (!symbol_table.symbol_exists("f1"))
      {
         printf("run_test10() - function 'f1' does not exist!\n");
         return false;
      }

      if (!symbol_table.remove_function("f"))
      {
         printf("run_test10() - Failed to remove function 'f'!\n");
         return false;
      }
      else if (!symbol_table.remove_function("f1"))
      {
         printf("run_test10() - Failed to remove function 'f1'!\n");
         return false;
      }
   }

   for (std::size_t r = 0; r < rounds; ++r)
   {
      symbol_table.add_stringvar("i",i);
      symbol_table.add_stringvar("j",j);

      symbol_table.add_stringvar("ii",ii);
      symbol_table.add_stringvar("jj",jj);

      if (!symbol_table.symbol_exists("i"))
      {
         printf("run_test10() - String 'i' does not exist!\n");
         return false;
      }
      else if (!symbol_table.symbol_exists("j"))
      {
         printf("run_test10() - String 'j' does not exist!\n");
         return false;
      }
      else if (!symbol_table.symbol_exists("ii"))
      {
         printf("run_test10() - String 'ii' does not exist!\n");
         return false;
      }
      else if (!symbol_table.symbol_exists("jj"))
      {
         printf("run_test10() - String 'jj' does not exist!\n");
         return false;
      }
      else if (!test::string(symbol_table,"i",i))
      {
         printf("run_test10() - String 'i' value failure!\n");
         return false;
      }
      else if (!test::string(symbol_table,"j",j))
      {
         printf("run_test10() - String 'j' value failure!\n");
         return false;
      }
      else if (!test::string(symbol_table,"ii",ii))
      {
         printf("run_test10() - String 'ii' value failure!\n");
         return false;
      }
      else if (!test::string(symbol_table,"jj",jj))
      {
         printf("run_test10() - String 'jj' value failure!\n");
         return false;
      }
      else if (!symbol_table.remove_stringvar("i"))
      {
         printf("run_test10() - Failed to remove String 'i'!\n");
         return false;
      }
      else if (!symbol_table.remove_stringvar("j"))
      {
         printf("run_test10() - Failed to remove String 'j'!\n");
         return false;
      }
      else if (!symbol_table.remove_stringvar("ii"))
      {
         printf("run_test10() - Failed to remove String 'ii'!\n");
         return false;
      }
      else if (!symbol_table.remove_stringvar("jj"))
      {
         printf("run_test10() - Failed to remove String 'jj'!\n");
         return false;
      }
   }

   for (std::size_t r = 0; r < rounds; ++r)
   {
      symbol_table.add_variable("x",  x);
      symbol_table.add_variable("y",  y);
      symbol_table.add_variable("xx",xx);
      symbol_table.add_variable("yy",yy);

      std::vector<std::string> expected_var_list;

      expected_var_list.push_back( "x");
      expected_var_list.push_back( "y");
      expected_var_list.push_back("xx");
      expected_var_list.push_back("yy");

      std::deque<std::pair<std::string,T> > variable_list;

      symbol_table.get_variable_list(variable_list);

      if (variable_list.size() != expected_var_list.size())
      {
         printf("run_test10() - Failed to get variable list (1)\n");
         return false;
      }

      std::size_t found_count = 0;

      for (std::size_t i = 0; i < variable_list.size(); ++i)
      {
         for (std::size_t j = 0; j < expected_var_list.size(); ++j)
         {
            if (variable_list[i].first == expected_var_list[j])
            {
               ++found_count;
               break;
            }
         }
      }

      if (found_count != expected_var_list.size())
      {
         printf("run_test10() - Failed to get variable list (2)\n");
         return false;
      }
   }

   for (std::size_t r = 0; r < rounds; ++r)
   {
      symbol_table.add_variable("x",  x);
      symbol_table.add_variable("y",  y);
      symbol_table.add_variable("xx",xx);
      symbol_table.add_variable("yy",yy);

      std::vector<std::string> expected_var_list;

      expected_var_list.push_back( "x");
      expected_var_list.push_back( "y");
      expected_var_list.push_back("xx");
      expected_var_list.push_back("yy");

      std::deque<std::string> variable_list;

      symbol_table.get_variable_list(variable_list);

      if (variable_list.size() != expected_var_list.size())
      {
         printf("run_test10() - Failed to get variable list (3)\n");
         return false;
      }

      std::size_t found_count = 0;

      for (std::size_t i = 0; i < variable_list.size(); ++i)
      {
         for (std::size_t j = 0; j < expected_var_list.size(); ++j)
         {
            if (variable_list[i] == expected_var_list[j])
            {
               ++found_count;
               break;
            }
         }
      }

      if (found_count != expected_var_list.size())
      {
         printf("run_test10() - Failed to get variable list (4)\n");
         return false;
      }
   }

   for (std::size_t r = 0; r < rounds; ++r)
   {
      symbol_table.add_stringvar( "i", i);
      symbol_table.add_stringvar( "j", j);
      symbol_table.add_stringvar("ii",ii);
      symbol_table.add_stringvar("jj",jj);

      std::vector<std::string> expected_var_list;

      expected_var_list.push_back( "i");
      expected_var_list.push_back( "j");
      expected_var_list.push_back("ii");
      expected_var_list.push_back("jj");

      std::deque<std::pair<std::string,std::string> > stringvar_list;

      symbol_table.get_stringvar_list(stringvar_list);

      if (stringvar_list.size() != expected_var_list.size())
      {
         printf("run_test10() - Failed to get stringvar list (1)\n");
         return false;
      }

      std::size_t found_count = 0;

      for (std::size_t i = 0; i < stringvar_list.size(); ++i)
      {
         for (std::size_t j = 0; j < expected_var_list.size(); ++j)
         {
            if (stringvar_list[i].first == expected_var_list[j])
            {
               ++found_count;
               break;
            }
         }
      }

      if (found_count != expected_var_list.size())
      {
         printf("run_test10() - Failed to get stringvar list (2)\n");
         return false;
      }
   }

   for (std::size_t r = 0; r < rounds; ++r)
   {
      symbol_table.add_stringvar("i",  i);
      symbol_table.add_stringvar("j",  j);
      symbol_table.add_stringvar("ii",ii);
      symbol_table.add_stringvar("jj",jj);

      std::vector<std::string> expected_var_list;

      expected_var_list.push_back( "i");
      expected_var_list.push_back( "j");
      expected_var_list.push_back("ii");
      expected_var_list.push_back("jj");

      std::deque<std::string> stringvar_list;

      symbol_table.get_stringvar_list(stringvar_list);

      if (stringvar_list.size() != expected_var_list.size())
      {
         printf("run_test10() - Failed to get stringvar list (3)\n");
         return false;
      }

      std::size_t found_count = 0;

      for (std::size_t i = 0; i < stringvar_list.size(); ++i)
      {
         for (std::size_t j = 0; j < expected_var_list.size(); ++j)
         {
            if (stringvar_list[i] == expected_var_list[j])
            {
               ++found_count;
               break;
            }
         }
      }

      if (found_count != expected_var_list.size())
      {
         printf("run_test10() - Failed to get stringvar list (4)\n");
         return false;
      }
   }

   {
      T x0 = T(0);
      T y0 = T(0);
      T z0 = T(0);

      std::string expression_string = "(x0 + y0) / z0";

      static const std::size_t rounds = 100;

      for (std::size_t i = 0; i < rounds; ++i)
      {
         expression_t expression0;

         x0 = T(i + 1.11);
         y0 = T(i + 2.22);
         z0 = T(i + 3.33);

         exprtk::symbol_table<T> st0;

         st0.add_variable("x0",x0);
         st0.add_variable("y0",y0);
         st0.add_variable("z0",z0);

         expression0.register_symbol_table(st0);

         {
            exprtk::parser<T> parser;
            if (!parser.compile(expression_string,expression0))
            {
               printf("run_test10() - Error: %s   Expression: %s\n",
                      parser.error().c_str(),
                      expression_string.c_str());
               return false;
            }
         }

         {
            expression_t expression1;
            exprtk::symbol_table<T> st1 = st0;
            expression1.register_symbol_table(st1);

            {
               exprtk::parser<T> parser;
               if (!parser.compile(expression_string,expression1))
               {
                  printf("run_test10() - Error: %s   Expression: %s\n",
                         parser.error().c_str(),
                         expression_string.c_str());
                  return false;
               }
            }

            st1.remove_variable("x0");
            st1.remove_variable("y0");
            st1.remove_variable("z0");
         }
      }
   }

   return true;
}

template <typename T>
inline bool run_test11()
{
   typedef exprtk::expression<T> expression_t;
   std::string expression_string = "(x + y) / 3";

   T x = T(1.0);
   T y = T(2.0);

   exprtk::symbol_table<T> symbol_table;
   symbol_table.add_variable("x",x);
   symbol_table.add_variable("y",y);

   expression_t expression;
   expression.register_symbol_table(symbol_table);

   static const std::size_t rounds = 500;

   for (std::size_t i = 0; i < rounds; ++i)
   {
      {
         exprtk::parser<T> parser;
         if (!parser.compile(expression_string,expression))
         {
            printf("run_test11() - Error: %s   Expression: %s\n",
                   parser.error().c_str(),
                   expression_string.c_str());
            return false;
         }
      }

      if (not_equal(expression.value(),(x + y)/T(3.0),T(0.000001)))
      {
         printf("run_test11() - Error in evaluation!(1)\n");
         return false;
      }

      expression.release();

      if (false == (!expression))
      {
         printf("run_test11() - Error in evaluation!(2)\n");
         return false;
      }

      {
         exprtk::parser<T> parser;
         if (!parser.compile(expression_string,expression))
         {
            printf("run_test11() - Error: %s   Expression: %s\n",
                   parser.error().c_str(),
                   expression_string.c_str());

            return false;
         }
      }

      expression.value();

      if (not_equal(expression.value(),(x + y)/T(3.0),T(0.000001)))
      {
         printf("run_test11() - Error in evaluation!(3)\n");
         return false;
      }
   }

   if (!exprtk::pgo_primer<T>())
   {
      printf("run_test11() - Failed PGO primer\n");
      return false;
   }

   return true;
}

template <typename T>
inline bool run_test12()
{
   typedef exprtk::expression<T> expression_t;
   static const std::string expression_string[] =
                            {
                               "equal(poly01(x,2.2,1.1),(2.2x^1+1.1))",
                               "equal(poly02(x,3.3,2.2,1.1),(3.3x^2+2.2x^1+1.1))",
                               "equal(poly03(x,4.4,3.3,2.2,1.1),(4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly04(x,5.5,4.4,3.3,2.2,1.1),(5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly05(x,6.6,5.5,4.4,3.3,2.2,1.1),(6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly06(x,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly07(x,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly08(x,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly09(x,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly10(x,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(2.2x^10+1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly11(x,3.3,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(3.3x^11+2.2x^10+1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equal(poly12(x,4.4,3.3,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(4.4x^12+3.3x^11+2.2x^10+1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "EquaL(Poly01(x,2.2,1.1),(2.2x^1+1.1))",
                               "eQuAl(pOly02(x,3.3,2.2,1.1),(3.3x^2+2.2x^1+1.1))",
                               "eqUal(poLy03(x,4.4,3.3,2.2,1.1),(4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "eQuAl(polY04(x,5.5,4.4,3.3,2.2,1.1),(5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "EquAl(pOLy05(x,6.6,5.5,4.4,3.3,2.2,1.1),(6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "EqUaL(pOly06(x,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "Equal(Poly07(x,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "eQual(PoLy08(x,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "eqUal(pOlY09(x,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equAl(POLY10(x,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(2.2x^10+1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "equaL(PolY11(x,3.3,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(3.3x^11+2.2x^10+1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))",
                               "EQUAL(pOLy12(x,4.4,3.3,2.2,1.1,9.9,8.8,7.7,6.6,5.5,4.4,3.3,2.2,1.1),(4.4x^12+3.3x^11+2.2x^10+1.1x^9+9.9x^8+8.8x^7+7.7x^6+6.6x^5+5.5x^4+4.4x^3+3.3x^2+2.2x^1+1.1))"
                            };
   static const std::size_t expression_string_size = sizeof(expression_string) / sizeof(std::string);

   T x = T(1.23456);

   exprtk::polynomial<T, 1> poly01;
   exprtk::polynomial<T, 2> poly02;
   exprtk::polynomial<T, 3> poly03;
   exprtk::polynomial<T, 4> poly04;
   exprtk::polynomial<T, 5> poly05;
   exprtk::polynomial<T, 6> poly06;
   exprtk::polynomial<T, 7> poly07;
   exprtk::polynomial<T, 8> poly08;
   exprtk::polynomial<T, 9> poly09;
   exprtk::polynomial<T,10> poly10;
   exprtk::polynomial<T,11> poly11;
   exprtk::polynomial<T,12> poly12;

   exprtk::symbol_table<T> symbol_table;

   symbol_table.add_variable("x",x);
   symbol_table.add_function("poly01", poly01);
   symbol_table.add_function("poly02", poly02);
   symbol_table.add_function("poly03", poly03);
   symbol_table.add_function("poly04", poly04);
   symbol_table.add_function("poly05", poly05);
   symbol_table.add_function("poly06", poly06);
   symbol_table.add_function("poly07", poly07);
   symbol_table.add_function("poly08", poly08);
   symbol_table.add_function("poly09", poly09);
   symbol_table.add_function("poly10", poly10);
   symbol_table.add_function("poly11", poly11);
   symbol_table.add_function("poly12", poly12);

   expression_t expression;
   expression.register_symbol_table(symbol_table);

   static const std::size_t rounds = 500;

   for (std::size_t i = 0; i < rounds; ++i)
   {
      for (std::size_t j = 0; j < expression_string_size; ++j)
      {
         const std::string& expr_str = expression_string[j];

         {
            exprtk::parser<T> parser;
            if (!parser.compile(expr_str,expression))
            {
               printf("run_test12() - Error: %s   Expression: %s\n",
                      parser.error().c_str(),
                      expr_str.c_str());
               return false;
            }
         }

         if (T(1.0) != expression.value())
         {
            printf("run_test12() - Error in evaluation! Expression: %s\n",expr_str.c_str());
            return false;
         }
      }
   }

   return true;
}

template <typename T>
struct sine_deg : public exprtk::ifunction<T>
{
   sine_deg() : exprtk::ifunction<T>(1) {}

   inline T operator()(const T& v)
   {
      return std::sin((v * T(exprtk::details::numeric::constant::pi))/T(180.0));
   }
};

template <typename T>
struct cosine_deg : public exprtk::ifunction<T>
{
   cosine_deg() : exprtk::ifunction<T>(1) {}

   inline T operator()(const T& v)
   {
      return std::cos((v * T(exprtk::details::numeric::constant::pi))/T(180.0));
   }
};

template <typename T>
inline bool run_test13()
{
   typedef exprtk::expression<T> expression_t;
   static const std::string expression_string[] =
                            {
                               "equal(sin(30),0.5)",
                               "equal(cos(60),0.5)",
                               "equal(sin(60),sqrt(3)/2)",
                               "equal(cos(30),sqrt(3)/2)",
                               "equal(sin(x_deg),0.5)",
                               "equal(cos(y_deg),0.5)",
                               "equal(sin(y_deg),sqrt(3)/2)",
                               "equal(cos(x_deg),sqrt(3)/2)",
                            };
   static const std::size_t expression_string_size = sizeof(expression_string) / sizeof(std::string);

   T x_deg = T(30);
   T y_deg = T(60);

   sine_deg<T> sine;
   cosine_deg<T> cosine;

   exprtk::symbol_table<T> symbol_table;

   symbol_table.add_variable("x_deg",x_deg);
   symbol_table.add_variable("y_deg",y_deg);

   symbol_table.add_function("sine_deg",sine);
   symbol_table.add_function("cosine_deg",cosine);

   expression_t expression;
   expression.register_symbol_table(symbol_table);

   static const std::size_t rounds = 100;

   for (std::size_t i = 0; i < rounds; ++i)
   {
      for (std::size_t j = 0; j < expression_string_size; ++j)
      {
         const std::string& expr_str = expression_string[j];

         {
            exprtk::parser<T> parser;

            parser.replace_symbol("sin","sine_deg");
            parser.replace_symbol("cos","cosine_deg");

            if (!parser.compile(expr_str,expression))
            {
               printf("run_test13() - Error: %s   Expression: %s\n",
                      parser.error().c_str(),
                      expr_str.c_str());
               return false;
            }
         }

         if (T(1.0) != expression.value())
         {
            printf("run_test13() - Error in evaluation! Expression: %s\n",expr_str.c_str());
            return false;
         }
      }
   }

   return true;
}

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))
   {
      if (buffer.empty())
         continue;
      else if ('#' == buffer[0])
         continue;
      ++line_count;
      sequence.push_back(buffer);
   }
   return line_count;
}

template <typename T>
inline bool run_test14()
{
   typedef exprtk::expression<T> expression_t;

   T x = T(0);
   T y = T(0);
   T z = T(0);
   T w = T(0);

   exprtk::polynomial<T, 1> poly01;
   exprtk::polynomial<T, 2> poly02;
   exprtk::polynomial<T, 3> poly03;
   exprtk::polynomial<T, 4> poly04;
   exprtk::polynomial<T, 5> poly05;
   exprtk::polynomial<T, 6> poly06;
   exprtk::polynomial<T, 7> poly07;
   exprtk::polynomial<T, 8> poly08;
   exprtk::polynomial<T, 9> poly09;
   exprtk::polynomial<T,10> poly10;
   exprtk::polynomial<T,11> poly11;
   exprtk::polynomial<T,12> poly12;

   exprtk::symbol_table<T> 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);
   symbol_table.add_function("poly01", poly01);
   symbol_table.add_function("poly02", poly02);
   symbol_table.add_function("poly03", poly03);
   symbol_table.add_function("poly04", poly04);
   symbol_table.add_function("poly05", poly05);
   symbol_table.add_function("poly06", poly06);
   symbol_table.add_function("poly07", poly07);
   symbol_table.add_function("poly08", poly08);
   symbol_table.add_function("poly09", poly09);
   symbol_table.add_function("poly10", poly10);
   symbol_table.add_function("poly11", poly11);
   symbol_table.add_function("poly12", poly12);

   expression_t expression;
   expression.register_symbol_table(symbol_table);

   exprtk::parser<T> parser;

   std::deque<std::string> expr_str_list;

   if (0 == load_expressions("exprtk_functional_test.txt",expr_str_list))
   {
      return true;
   }

   std::deque<exprtk::expression<T> > expression_list;
   bool failure = false;
   static const std::size_t rounds = 5;

   for (std::size_t r = 0; r < rounds; ++r)
   {
      for (std::size_t i = 0; i < expr_str_list.size(); ++i)
      {
         exprtk::expression<T> current_expression;

         current_expression.register_symbol_table(symbol_table);

         if (!parser.compile(expr_str_list[i],current_expression))
         {
            printf("run_test14() - Error: %s   Expression: %s\n",
                   parser.error().c_str(),
                   expr_str_list[i].c_str());
            return false;
         }
         else
            expression_list.push_back(current_expression);
      }

      for (std::size_t i = 0; i < expression_list.size(); ++i)
      {
         T result = expression_list[i].value();
         if (result != T(1))
         {
            failure = true;
            printf("run_test14() - Error with expression: %s\n",
                   expr_str_list[i].c_str());
         }
      }

      expression_list.clear();

      if (failure)
         break;
   }

   return !failure;
}

template <typename T>
inline bool run_test15()
{
   typedef exprtk::expression<T> expression_t;

   T x = T(1.1);
   T y = T(2.2);
   T z = T(3.3);

   exprtk::symbol_table<T> symbol_table;
   symbol_table.add_constants();
   symbol_table.add_variable("x",x);
   symbol_table.add_variable("y",y);
   symbol_table.add_variable("z",z);

   static const std::string expr_str_list[] =
                               {
                                  "2 - (x + y) / z//Comment 01   ",
                                  "2 - (x + y) / z#Comment 02    ",
                                  "2 - (x + y) / z //Comment 03  ",
                                  "2 - (x + y) / z #Comment 04   ",
                                  "2 - (x + y) / z//Comment 05 \n",
                                  "2 - (x + y) / z#Comment 06  \n",
                                  "2 - (x + y) / z //Comment 07\n",
                                  "2 - (x + y) / z #Comment 08 \n",
                                  "/* Comment 09*/2 - (x + y) / z",
                                  "/* Comment 10*/2 - (x + y) / z\n",
                                  "/* Comment 11*/2 - (x + y) / z/* Comment 12*/",
                                  "/* Comment 13*/2 - (x + y) / z/* Comment 14*/\n",
                                  "2 - /* Comment 15 */(x + y) / z",
                                  "2 - /* Comment 15 */(x + y) /* Comment 16 *// z \n",
                                  "2 - /* Comment 17 */(x + y) /* Comment 18 */ / z //Comment 19\n",
                                  "2 - /* Comment 20 */(x + y) /* Comment 21 */ / z #Comment 22\n",
                                  "2 - /* Comment 23 */(x + y) /* Comment 24 */ / z //Comment 25",
                                  "2 - /* Comment 26 */(x + y) /* Comment 27 */ / z #Comment 28"
                               };
   static const std::size_t expr_str_list_size = sizeof(expr_str_list) / sizeof(std::string);

   std::deque<expression_t> expression_list;

   for (std::size_t i = 0; i < expr_str_list_size; ++i)
   {
      expression_t expression;
      expression.register_symbol_table(symbol_table);

      exprtk::parser<T> parser;

      if (!parser.compile(expr_str_list[i],expression))
      {
         printf("run_test15() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expr_str_list[i].c_str());
         return false;
      }
      else
         expression_list.push_back(expression);
   }

   expression_t base_expression;
   const std::string base_expr_str = "2 - (x + y) / z";

   {
      base_expression.register_symbol_table(symbol_table);

      exprtk::parser<T> parser;

      if (!parser.compile(base_expr_str,base_expression))
      {
         printf("run_test15() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                base_expr_str.c_str());
         return false;
      }
   }

   bool failure = false;

   for (std::size_t i = 0; i < expression_list.size(); ++i)
   {
      T base_result = base_expression.value();
      T result      = expression_list[i].value();

      if (not_equal(base_result,result))
      {
         printf("run_test15() - Error in evaluation! (1)  Base: %20.10f\tResult: %20.10f\tExpression: %s\n",
                base_result,
                result,
                expr_str_list[i].c_str());
         failure = true;
      }
   }

   return !failure;
}

template <typename T>
struct base_func : public exprtk::ifunction<T>
{
   typedef const T& type;
   base_func(const std::size_t& n) : exprtk::ifunction<T>(n) {}
   inline T operator()(type v0, type v1, type v2, type v3, type v4) { return (v0 + v1 + v2 + v3 + v4); }
   inline T operator()(type v0, type v1, type v2, type v3) { return (v0 + v1 + v2 + v3); }
   inline T operator()(type v0, type v1, type v2) { return (v0 + v1 + v2); }
   inline T operator()(type v0, type v1) { return (v0 + v1); }
   inline T operator()(type v0) { return v0; }
   inline T operator()() { return T(1.1234); }
};

template <typename T> struct test_func5 : public base_func<T> { test_func5() : base_func<T>(5){} };
template <typename T> struct test_func4 : public base_func<T> { test_func4() : base_func<T>(4){} };
template <typename T> struct test_func3 : public base_func<T> { test_func3() : base_func<T>(3){} };
template <typename T> struct test_func2 : public base_func<T> { test_func2() : base_func<T>(2){} };
template <typename T> struct test_func1 : public base_func<T> { test_func1() : base_func<T>(1){} };
template <typename T> struct test_func0 : public base_func<T> { test_func0() : base_func<T>(0){} };

template <typename T>
inline bool run_test16()
{
   typedef exprtk::expression<T> expression_t;

   T x = T(1.1);
   T y = T(2.2);
   T z = T(3.3);
   T w = T(4.4);
   T u = T(5.5);

   test_func0<T> test_func00;
   test_func1<T> test_func01;
   test_func2<T> test_func02;
   test_func3<T> test_func03;
   test_func4<T> test_func04;
   test_func5<T> test_func05;

   exprtk::symbol_table<T> 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);
   symbol_table.add_variable("u",u);

   symbol_table.add_function("test_func0",test_func00);
   symbol_table.add_function("test_func1",test_func01);
   symbol_table.add_function("test_func2",test_func02);
   symbol_table.add_function("test_func3",test_func03);
   symbol_table.add_function("test_func4",test_func04);
   symbol_table.add_function("test_func5",test_func05);

   static const std::string expr_str_list[] =
                               {
                                  "test_func0 = test_func0()",
                                  "test_func0 == test_func0()",
                                  "equal(1.1 + test_func0,test_func0() + 1.1)",
                                  "equal(test_func0 + 1.1,1.1 + test_func0())",
                                  "equal((1.1 + test_func0),(test_func0() + 1.1))",
                                  "equal((test_func0 + 1.1),(1.1 + test_func0()))",
                                  "equal(test_func0,test_func0())",
                                  "equal(test_func0,1.1234)",
                                  "equal(test_func0(),1.1234)",
                                  "equal(test_func0 + test_func0(),2 * 1.1234)",
                                  "equal((test_func0 + test_func0()),2 * 1.1234)",
                                  "equal((test_func0) + (test_func0()),2 * 1.1234)",
                                  "equal(test_func1(x),(x))",
                                  "equal(test_func2(x,y),(x + y))",
                                  "equal(test_func3(x,y,z),(x + y + z))",
                                  "equal(test_func4(x,y,z,w),(x + y + z + w))",
                                  "equal(test_func5(x,y,z,w,u),(x + y + z + w + u))",
                                  "equal(1.1 + test_func0,1.1234 + 1.1)",
                                  "equal(1.1 + test_func0(),1.1234 + 1.1)",
                                  "equal(1.1 + test_func1(x),(x + 1.1))",
                                  "equal(1.1 + test_func2(x,y),(x + y + 1.1))",
                                  "equal(1.1 + test_func3(x,y,z),(x + y + z + 1.1))",
                                  "equal(1.1 + test_func4(x,y,z,w),(x + y + z + w + 1.1))",
                                  "equal(1.1 + test_func5(x,y,z,w,u),(x + y + z + w + u + 1.1))",
                                  "equal(test_func0 + 2.2,1.1234 + 2.2)",
                                  "equal(test_func0() + 2.2,1.1234 + 2.2)",
                                  "equal(test_func1(x) + 2.2,(x + 2.2))",
                                  "equal(test_func2(x,y) + 2.2,(x + y + 2.2))",
                                  "equal(test_func3(x,y,z) + 2.2,(x + y + z + 2.2))",
                                  "equal(test_func4(x,y,z,w) + 2.2,(x + y + z + w + 2.2))",
                                  "equal(test_func5(x,y,z,w,u) + 2.2,(x + y + z + w + u + 2.2))",
                                  "equal({[test_func1(x)]},{[(x)]})",
                                  "equal({[test_func2(x,y)]},{[(x + y)]})",
                                  "equal({[test_func3(x,y,z)]},{[(x + y + z)]})",
                                  "equal({[test_func4(x,y,z,w)]},{[(x + y + z + w)]})",
                                  "equal({[test_func5(x,y,z,w,u)]},{[(x + y + z + w + u)]})",
                                  "equal(test_func1(2.12),(2.12))",
                                  "equal(test_func2(2.12,3.12),(2.12 + 3.12))",
                                  "equal(test_func3(2.12,3.12,4.12),(2.12 + 3.12 + 4.12))",
                                  "equal(test_func4(2.12,3.12,4.12,5.12),(2.12 + 3.12 + 4.12 + 5.12))",
                                  "equal(test_func5(2.12,3.12,4.12,5.12,6.12),(2.12 + 3.12 + 4.12 + 5.12 + 6.12))",
                                  "equal(1.1 + test_func1(2.12),(2.12 + 1.1))",
                                  "equal(1.1 + test_func2(2.12,3.12),(2.12 + 3.12 + 1.1))",
                                  "equal(1.1 + test_func3(2.12,3.12,4.12),(2.12 + 3.12 + 4.12 + 1.1))",
                                  "equal(1.1 + test_func4(2.12,3.12,4.12,5.12),(2.12 + 3.12 + 4.12 + 5.12 + 1.1))",
                                  "equal(1.1 + test_func5(2.12,3.12,4.12,5.12,6.12),(2.12 + 3.12 + 4.12 + 5.12 + 6.12 + 1.1))",
                                  "equal(test_func1(2.12) + 2.2,(2.12 + 2.2))",
                                  "equal(test_func2(2.12,3.12) + 2.2,(2.12 + 3.12 + 2.2))",
                                  "equal(test_func3(2.12,3.12,4.12) + 2.2,(2.12 + 3.12 + 4.12 + 2.2))",
                                  "equal(test_func4(2.12,3.12,4.12,5.12) + 2.2,(2.12 + 3.12 + 4.12 + 5.12 + 2.2))",
                                  "equal(test_func5(2.12,3.12,4.12,5.12,6.12) + 2.2,(2.12 + 3.12 + 4.12 + 5.12 + 6.12 + 2.2))",
                                  "equal({[test_func1(2.12)]},{[(2.12)]})",
                                  "equal({[test_func2(2.12,3.12)]},{[(2.12 + 3.12)]})",
                                  "equal({[test_func3(2.12,3.12,4.12)]},{[(2.12 + 3.12 + 4.12)]})",
                                  "equal({[test_func4(2.12,3.12,4.12,5.12)]},{[(2.12 + 3.12 + 4.12 + 5.12)]})",
                                  "equal({[test_func5(2.12,3.12,4.12,5.12,6.12)]},{[(2.12 + 3.12 + 4.12 + 5.12 + 6.12)]})",
                                  "TeSt_FuNc0 = tEsT_fUnC0()",
                                  "TEst_fuNC0 == tESt_fUNc0()",
                                  "EquaL(tEsT_FuNC1(x),(x))",
                                  "eQuAl(teSt_fUnc2(x,y),(x + y))",
                                  "EqUaL(tEsT_fUNc3(x,y,z),(x + y + z))",
                                  "EQUal(TEst_FunC4(x,y,z,w),(x + y + z + w))",
                                  "eqUAL(tEsT_FuNc5(x,y,z,w,u),(x + y + z + w + u))"
                               };
   static const std::size_t expr_str_list_size = sizeof(expr_str_list) / sizeof(std::string);

   std::deque<expression_t> expression_list;

   for (std::size_t i = 0; i < expr_str_list_size; ++i)
   {
      expression_t expression;
      expression.register_symbol_table(symbol_table);

      exprtk::parser<T> parser;

      if (!parser.compile(expr_str_list[i],expression))
      {
         printf("run_test16() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expr_str_list[i].c_str());
         return false;
      }
      else
         expression_list.push_back(expression);
   }

   bool failure = false;

   for (std::size_t i = 0; i < expression_list.size(); ++i)
   {
      if (T(1) != expression_list[i].value())
      {
         printf("run_test16() - Error in evaluation! (1) Expression: %s\n",
                expr_str_list[i].c_str());
         failure = true;
      }
   }

   return !failure;
}

template <typename T>
inline bool run_test17()
{
   typedef exprtk::expression<T> expression_t;

   T x = T(1.1);
   T y = T(2.2);
   T z = T(3.3);
   T w = T(4.4);
   T u = T(5.5);
   T v = T(6.6);
   T t = T(7.7);

   T one  = T(1.0);
   T zero = T(0.0);

   exprtk::symbol_table<T> 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);
   symbol_table.add_variable("u",u);
   symbol_table.add_variable("v",v);
   symbol_table.add_variable("t",t);

   symbol_table.add_variable("one",one);
   symbol_table.add_variable("zero",zero);

   static const std::string expr_str_list[] =
                               {
                                  "equal(mand(one,one),1.0)",
                                  "equal(mand(one,zero),0.0)",
                                  "equal(mand(zero,one),0.0)",
                                  "equal(mand(zero,zero),0.0)",
                                  "equal(mand(one,one),1.0)",
                                  "equal(mand(one,one,one),1.0)",
                                  "equal(mand(one,one,one,one),1.0)",
                                  "equal(mand(one,one,one,one,one),1.0)",
                                  "equal(mand(one,one,zero),0.0)",
                                  "equal(mand(one,one,one,zero),0.0)",
                                  "equal(mand(one,one,one,one,zero),0.0)",
                                  "equal(mand(one,one,one,one,one,zero),0.0)",
                                  "equal(mor(one,one),1.0)",
                                  "equal(mor(one,zero),1.0)",
                                  "equal(mor(zero,one),1.0)",
                                  "equal(mor(zero,zero),0.0)",
                                  "equal(mor(one,one),1.0)",
                                  "equal(mor(one,one,zero),1.0)",
                                  "equal(mor(one,zero,one),1.0)",
                                  "equal(mor(one,zero,one,zero),1.0)",
                                  "equal(mor(zero,one),1.0)",
                                  "equal(mor(zero,zero,one),1.0)",
                                  "equal(mor(zero,zero,zero,zero,one),1.0)",
                                  "equal(mor(zero,zero,zero,zero,zero,one),1.0)",
                                  "equal(mor(zero,zero,zero,zero,zero,zero,zero,zero),0.0)",
                                  "equal((one nand one),not(mand(one,one)))",
                                  "equal((one nand zero),not(mand(one,zero)))",
                                  "equal((zero nand one),not(mand(zero,one)))",
                                  "equal((zero nand zero),not(mand(zero,zero)))",
                                  "equal((one nor one),not(mor(one,one)))",
                                  "equal((one nor zero),not(mor(one,zero)))",
                                  "equal((zero nor one),not(mor(zero,one)))",
                                  "equal((zero nor zero),not(mor(zero,zero)))",
                                  "equal(sum(x,y,z,w,u,v,t),(x+y+z+w+u+v+t))",
                                  "equal(sum(x+t,y+v,z+u,w+w,u+z,v+y,t+x),2*(x+y+z+w+u+v+t))",
                                  "equal(mul(x,y,z,w,u,v,t),(x*y*z*w*u*v*t))",
                                  "equal(mul(x*t,y*v,z*u,w*w,u*z,v*y,t*x),(x*y*z*w*u*v*t)^2)",
                                  "equal(sum(x,1.0,y,2.0,z,3.0,w,4.0,u,5.0,v,6.0,t),(x+y+z+w+u+v+t+21.0))",
                                  "equal(sum(x+1.0,y+2.0,z+3.0,w+4.0,u+5.0,v+6.0,t),(x+y+z+w+u+v+t+21.0))",
                                  "equal(mul(x,1.0,y,2.0,z,3.0,w,4.0,u,5.0,v,6.0,t),(x*y*z*w*u*v*t*720.0))",
                                  "equal(mul(x*1.0,y*2.0,z*3.0,w*4.0,u*5.0,v*6.0,t),(x*y*z*w*u*v*t*720.0))",
                                  "equal(min(x,y,z,w,u,v,t,zero),zero)",
                                  "equal(min(x+y,z+w,u+v,t,zero),zero)",
                                  "equal(max(one,x,y,z,w,u,v,t),t)",
                                  "equal(max(x+y,z+w,u+v,t,one),u+v)",
                                  "equal(avg(x,y,z,w,u,v,t),(x+y+z+w+u+v+t)/7.0)",
                                  "equal(avg(x+t,y+v,z+u,w+w,u+z,v+y,t+x),2/7*(x+y+z+w+u+v+t))"
                               };
   static const std::size_t expr_str_list_size = sizeof(expr_str_list) / sizeof(std::string);

   std::deque<expression_t> expression_list;

   for (std::size_t i = 0; i < expr_str_list_size; ++i)
   {
      expression_t expression;
      expression.register_symbol_table(symbol_table);

      exprtk::parser<T> parser;

      if (!parser.compile(expr_str_list[i],expression))
      {
         printf("run_test17() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expr_str_list[i].c_str());
         return false;
      }
      else
         expression_list.push_back(expression);
   }

   bool failure = false;

   for (std::size_t i = 0; i < expression_list.size(); ++i)
   {
      if (T(1) != expression_list[i].value())
      {
         printf("run_test17() - Error in evaluation! (1) Expression: %s\n",
                expr_str_list[i].c_str());
         failure = true;
      }
   }

   return !failure;
}

template <typename T>
struct va_func : public exprtk::ivararg_function<T>
{
   inline T operator()(const std::vector<T>& arglist)
   {
      T result = T(0);
      for (std::size_t i = 0; i < arglist.size(); ++i)
      {
         result += arglist[i];
      }
      return result;
   }
};

template <typename T>
inline bool run_test18()
{
   typedef exprtk::expression<T> expression_t;

   T x = T(1.1);
   T y = T(2.2);
   T z = T(3.3);
   T w = T(4.4);
   T u = T(5.5);
   T v = T(6.6);
   T t = T(7.7);

   va_func<T> vaf;

   exprtk::symbol_table<T> 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);
   symbol_table.add_variable("u",u);
   symbol_table.add_variable("v",v);
   symbol_table.add_variable("t",t);
   symbol_table.add_vararg_function("va_func",vaf);

   static const std::string expr_str_list[] =
                               {
                                  "equal(va_func(1,2,3,4,5,6,7,8,9),(1+2+3+4+5+6+7+8+9))",
                                  "equal(va_func(1,x,3,y,5,z,7,w,9),(1+x+3+y+5+z+7+w+9))",
                                  "equal(va_func(x,2,y,4,z,6,w,8,u),(x+2+y+4+z+6+w+8+u))",
                                  "equal(va_func(x,y,z,w,u,v,t,1,2,3),(x+y+z+w+u+v+t+1+2+3))",
                                  "equal(va_func(x,y,z,w,u,v,t),(x+y+z+w+u+v+t))",
                                  "equal(va_func(x+t,y+v,z+u,w+w,u+z,v+y,t+x),2*(x+y+z+w+u+v+t))",
                                  "equal(1+va_func(1,x,3,y,5,z,7,w,9),(1+x+3+y+5+z+7+w+9)+1)",
                                  "equal(va_func(va_func(x,y,z,w,u,v,t),va_func(x,y,z,w,u,v,t)),2*(x+y+z+w+u+v+t))",
                                  "equal(va_func(va_func(x),va_func(y),va_func(z)),va_func(x,y,z))",
                                  "equal(va_func(va_func(va_func(va_func(va_func(va_func(va_func(va_func(x)))))))),x)",
                                  "equal(va_func(va_func(va_func(va_func(va_func(va_func(va_func(va_func(123.456)))))))),123.456)",
                                  "equal(va_func(va_func(va_func(va_func(va_func(va_func(va_func(va_func(x+1)))))))),x+1)",
                                  "equal(va_func(va_func(va_func(va_func(va_func(va_func(va_func(va_func(x+y)))))))),x+y)"
                               };
   static const std::size_t expr_str_list_size = sizeof(expr_str_list) / sizeof(std::string);

   std::deque<expression_t> expression_list;

   for (std::size_t i = 0; i < expr_str_list_size; ++i)
   {
      expression_t expression;
      expression.register_symbol_table(symbol_table);

      exprtk::parser<T> parser;

      if (!parser.compile(expr_str_list[i],expression))
      {
         printf("run_test18() - Error: %s   Expression: %s\n",
                parser.error().c_str(),
                expr_str_list[i].c_str());
         return false;
      }
      else
         expression_list.push_back(expression);
   }

   bool failure = false;

   for (std::size_t i = 0; i < expression_list.size(); ++i)
   {
      if (T(1) != expression_list[i].value())
      {
         printf("run_test18() - Error in evaluation! (1) Expression: %s\n",
                expr_str_list[i].c_str());
         failure = true;
      }
   }

   return !failure;
}

int main()
{
   #define perform_test(Type,Number) \
   { \
      exprtk::timer timer; \
      timer.start(); \
      if (!run_test##Number<Type>()) \
      { \
         printf("run_test"#Number" ("#Type") *** FAILED! ***\n"); \
      } \
      else \
      { \
         timer.stop(); \
         printf("run_test"#Number" ("#Type") - Result: SUCCESS     Time: %8.4fsec\n",timer.time()); \
      } \
   } \

   perform_test(double,00)
   perform_test(double,01)
   perform_test(double,02)
   perform_test(double,03)
   perform_test(double,04)
   perform_test(double,05)
   perform_test(double,06)
   perform_test(double,07)
   perform_test(double,08)
   perform_test(double,09)
   perform_test(double,10)
   perform_test(double,11)
   perform_test(double,12)
   perform_test(double,13)
   perform_test(double,14)
   perform_test(double,15)
   perform_test(double,16)
   perform_test(double,17)
   perform_test(double,18)

   #undef perform_test

   return 0;
}