bool doTest(const char *equation, double x, double result, const double tol = 0.00000000001) {
yy_scan_string(equation);
int rc = yyparse();
if (rc == 0) {
vectorsUsed.clear();
Equation::Node *eq = static_cast<Equation::Node*>(ParsedEquation);
assert(eq);
ParsedEquation = 0L;
Equation::Context ctx;
ctx.sampleCount = 2;
ctx.noPoint = NOPOINT;
ctx.x = x;
ctx.xVector = xVector;
if (xVector) {
ctx.sampleCount = xVector->length();
}
Equation::FoldVisitor vis(&ctx, &eq);
if (eq->isConst() && !dynamic_cast<Equation::Number*>(eq)) {
if (!optimizerFailed) {
optimizerFailed = true;
::rc--;
printf("Optimizer bug: found an unoptimized const equation. Optimizing for coverage purposes.\n");
}
double v = eq->value(&ctx);
delete eq;
eq = new Equation::Number(v);
}
KstScalarMap scm;
KstStringMap stm;
eq->collectObjects(vectorsUsed, scm, stm);
//
// do not use -1 in the call to eq->update( ... )
// as this will simply return without updating the plugin,
// unless an underlying data object was updated...
//
eq->update(0, &ctx);
double v = eq->value(&ctx);
delete eq;
if (fabs(v - result) < tol || (result != result && v != v) || (result == INF && v == INF) || (result == -INF && v == -INF)) {
return true;
} else {
printf("Result: %.16f\n", v);
return false;
}
} else {
// Parse error
printf("Failures on [%s] -------------------------\n", equation);
for (QStringList::ConstIterator i = Equation::errorStack.constBegin(); i != Equation::errorStack.constEnd(); ++i) {
printf("%s\n", (*i).latin1());
}
printf("------------------------------------------\n");
delete (Equation::Node*)ParsedEquation;
ParsedEquation = 0L;
return false;
}
}
double Equation::interpret(const char *txt, bool *ok) {
if (!txt || !*txt) {
if (ok) {
*ok = false;
}
return 0.0;
}
mutex().lock();
yy_scan_string(txt);
int rc = yyparse();
if (rc == 0) {
Equation::Node *eq = static_cast<Equation::Node*>(ParsedEquation);
ParsedEquation = 0L;
mutex().unlock();
Equation::Context ctx;
ctx.sampleCount = 2;
ctx.noPoint = KST::NOPOINT;
ctx.x = 0.0;
ctx.xVector = 0L;
Equation::FoldVisitor vis(&ctx, &eq);
double v = eq->value(&ctx);
delete eq;
if (ok) {
*ok = true;
}
return v;
} else {
ParsedEquation = 0L;
mutex().unlock();
if (ok) {
*ok = false;
}
return 0.0;
}
}