void doit(const Expr& e,
const Expr& tests,
const Expr& unks,
const Expr& u0,
const Expr& unkParams,
const Expr& paramVals,
const EvalContext& region)
{
TimeMonitor t0(doitTimer());
EvalManager mgr;
mgr.setRegion(region);
mgr.setVerb(5);
static RCP<AbstractEvalMediator> mediator
= rcp(new StringEvalMediator());
mgr.setMediator(mediator);
const EvaluatableExpr* ev
= dynamic_cast<const EvaluatableExpr*>(e[0].ptr().get());
Expr fixed;
Expr fixed0;
Out::os() << "params = " << unkParams << std::endl;
Out::os() << "param vals = " << paramVals << std::endl;
DerivSet d = SymbPreprocessor::setupSensitivities(e[0],
tests,
unks,
u0,
unkParams,
paramVals,
fixed,
fixed0,
fixed0,
fixed0,
region,
Sensitivities);
Tabs tab;
Out::os() << tab << "done setup" << std::endl;
Out::os() << tab << *ev->sparsitySuperset(region) << std::endl;
// ev->showSparsity(Out::os(), region);
// RCP<EvalVectorArray> results;
Array<double> constantResults;
Array<RCP<EvalVector> > vectorResults;
Out::os() << tab << "starting eval" << std::endl;
ev->evaluate(mgr, constantResults, vectorResults);
ev->sparsitySuperset(region)->print(Out::os(), vectorResults, constantResults);
// results->print(Out::os(), ev->sparsitySuperset(region).get());
}
