00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #include "SundanceDerivOfSymbFunc.hpp"
00032
00033 #include "SundanceTestFuncElement.hpp"
00034
00035 #include "SundanceCoordExpr.hpp"
00036 #include "PlayaTabs.hpp"
00037 #include "SundanceOut.hpp"
00038
00039 using namespace Sundance;
00040 using namespace Sundance;
00041
00042 using namespace Sundance;
00043 using namespace Teuchos;
00044
00045
00046 DerivOfSymbFunc::DerivOfSymbFunc(const MultiIndex& op,
00047 const RCP<ScalarExpr>& arg)
00048 : DiffOp(op, arg), argFid_(FunctionIdentifier())
00049 {
00050 const SymbolicFuncElement* f
00051 = dynamic_cast<const SymbolicFuncElement*>(evaluatableArg());
00052 TEUCHOS_TEST_FOR_EXCEPTION(f==0, std::logic_error, "argument to DerivOfSymbFunc ctor "
00053 "is not a symbolic function");
00054 argFid_ = f->fid();
00055 }
00056
00057
00058 Deriv DerivOfSymbFunc::representMeAsFunctionalDeriv() const
00059 {
00060 const SymbolicFuncElement* f
00061 = dynamic_cast<const SymbolicFuncElement*>(evaluatableArg());
00062 TEUCHOS_TEST_FOR_EXCEPTION(f==0, std::logic_error, "DerivOfSymbFunc::"
00063 "representMeAsFunctionalDeriv(), 'this' pointer "
00064 "is not a symbolic function");
00065 return Deriv(f, SpatialDerivSpecifier(mi()));
00066 }
00067
00068 Evaluator* DerivOfSymbFunc::createEvaluator(const EvaluatableExpr* expr,
00069 const EvalContext& context) const
00070 {
00071 return new DerivOfSymbFuncEvaluator(dynamic_cast<const DerivOfSymbFunc*>(expr), context);
00072 }
00073
00074
00075 bool DerivOfSymbFunc::lessThan(const ScalarExpr* other) const
00076 {
00077 const DerivOfSymbFunc* d = dynamic_cast<const DerivOfSymbFunc*>(other);
00078 TEUCHOS_TEST_FOR_EXCEPTION(d==0, std::logic_error, "cast should never fail at this point");
00079
00080 if (argFid_ < d->argFid_) return true;
00081 if (d->argFid_ < argFid_) return false;
00082
00083 return DiffOp::lessThan(other);
00084 }