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 "SundanceConstantExpr.hpp"
00032
00033 using namespace Sundance;
00034 using namespace Sundance;
00035
00036 using namespace Sundance;
00037 using namespace Teuchos;
00038
00039 ConstantExpr::ConstantExpr(const double& value)
00040 : SpatiallyConstantExpr(), value_(value)
00041 {}
00042
00043
00044
00045 Set<MultipleDeriv>
00046 ConstantExpr::internalFindW(int order, const EvalContext& context) const
00047 {
00048 Set<MultipleDeriv> rtn;
00049
00050 if (order==0) rtn.put(MultipleDeriv());
00051 Tabs tab;
00052 SUNDANCE_MSG2(context.setupVerbosity(),
00053 tab << "ConstantExpr::internalFindW found" << rtn << " for order="
00054 << order);
00055
00056 return rtn;
00057 }
00058
00059 Set<MultipleDeriv>
00060 ConstantExpr::internalFindV(int order, const EvalContext& context) const
00061 {
00062 Set<MultipleDeriv> rtn;
00063 Tabs tab;
00064 SUNDANCE_MSG2(context.setupVerbosity(),
00065 tab << "ConstantExpr::internalFindV is a no-op");
00066 return rtn;
00067 }
00068
00069
00070 Set<MultipleDeriv>
00071 ConstantExpr::internalFindC(int order, const EvalContext& context) const
00072 {
00073 Tabs tab;
00074 SUNDANCE_MSG2(context.setupVerbosity(),
00075 tab << "ConstantExpr::internalFindC is forwarding to findR()");
00076 return findR(order, context);
00077 }
00078
00079
00080 bool ConstantExpr::lessThan(const ScalarExpr* other) const
00081 {
00082 const ConstantExpr* c = dynamic_cast<const ConstantExpr*>(other);
00083 TEUCHOS_TEST_FOR_EXCEPTION(c==0, std::logic_error, "cast should never fail at this point");
00084 return value() < c->value();
00085 }
00086
00087
00088 std::ostream& ConstantExpr::toText(std::ostream& os, bool ) const
00089 {
00090 os << value();
00091 return os;
00092 }
00093
00094
00095 XMLObject ConstantExpr::toXML() const
00096 {
00097 XMLObject rtn("Constant");
00098 rtn.addAttribute("value", Teuchos::toString(value()));
00099 return rtn;
00100 }
00101
00102