Deriv DiffOpEvaluator::remainder(const MultipleDeriv& big,
const MultipleDeriv& little, int verb) const
{
Tabs tab;
SUNDANCE_MSG5(verb, tab << "computing remainder: big=" << big << ", little="
<< little);
TEUCHOS_TEST_FOR_EXCEPT(big.order()-little.order() != 1);
MultipleDeriv r;
if (little.order()==0) r = big;
else r = big.factorOutDeriv(little);
SUNDANCE_MSG5(verb, tab << "remainder = " << r);
TEUCHOS_TEST_FOR_EXCEPT(r.order() != 1);
return *(r.begin());
}
Set<MultipleDeriv> DiffOpEvaluator
::backedDerivs(const MultipleDeriv& mu,
const Set<MultipleDeriv>& W1, int verb) const
{
Tabs tabs;
SUNDANCE_MSG3(verb, tabs << "computing backed-out derivs for mu= " << mu
<< ", W1=" << W1);
Set<MultipleDeriv> rtn;
if (mu.order() != 0)
{
const MultiIndex& alpha = expr()->mi();
for (Set<MultipleDeriv>::const_iterator i=W1.begin(); i!=W1.end(); i++)
{
const MultipleDeriv& md = *i;
TEUCHOS_TEST_FOR_EXCEPT(md.order() != 1);
Deriv lambda = *(md.begin());
if (lambda.isCoordDeriv()) continue;
TEUCHOS_TEST_FOR_EXCEPT(!lambda.isFunctionalDeriv());
FunctionIdentifier lambda_fid = lambda.fid();
const MultiIndex& lambda_mi = lambda.opOnFunc().mi();
for (MultipleDeriv::const_iterator j=mu.begin(); j!=mu.end(); j++)
{
const Deriv& d = *j;
if (d.isCoordDeriv()) continue;
FunctionIdentifier d_fid = d.fid();
const MultiIndex& d_mi = d.opOnFunc().mi();
if (d_fid != lambda_fid) continue;
if (!(alpha + lambda_mi == d_mi)) continue;
MultipleDeriv z = mu.factorOutDeriv(d);
z.put(lambda);
rtn.put(z);
}
}
}
SUNDANCE_MSG3(verb, tabs << "backed-out derivs = " << rtn);
return rtn;
}
