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 "SundanceCurveNormExpr.hpp"
00032 #include "SundanceEvalManager.hpp"
00033 #include "SundanceSparsitySuperset.hpp"
00034 #include "SundanceOut.hpp"
00035 #include "SundanceObjectWithVerbosity.hpp"
00036
00037 using namespace Sundance;
00038 using namespace Teuchos;
00039
00040
00041 CurveNormExpr::CurveNormExpr(int dir, const std::string& name)
00042 : EvaluatableExpr(),
00043 dir_(dir),
00044 name_(coordName(dir, name))
00045 {}
00046
00047 bool CurveNormExpr::lessThan(const ScalarExpr* other) const
00048 {
00049 const CurveNormExpr* c = dynamic_cast<const CurveNormExpr*>(other);
00050 TEUCHOS_TEST_FOR_EXCEPTION(c==0, std::logic_error, "cast should never fail at this point");
00051 return dir() < c->dir();
00052 }
00053
00054 XMLObject CurveNormExpr::toXML() const
00055 {
00056 XMLObject rtn("CurveNormExpr");
00057 rtn.addAttribute("dir", Teuchos::toString(dir_));
00058 rtn.addAttribute("name", name());
00059 return rtn;
00060 }
00061
00062 string CurveNormExpr::coordName(int dir, const std::string& name)
00063 {
00064 if (name.length() > 0) return name;
00065 switch(dir)
00066 {
00067 case 0:
00068 return "nx";
00069 case 1:
00070 return "ny";
00071 case 2:
00072 return "nz";
00073 default:
00074 TEUCHOS_TEST_FOR_EXCEPTION(true, std::runtime_error,
00075 "CurveNormExpr::coordName direction out of range [0,2]");
00076 return "error";
00077 }
00078 }
00079
00080
00081 Set<MultipleDeriv>
00082 CurveNormExpr::internalFindW(int order, const EvalContext& context) const
00083 {
00084 Tabs tab0;
00085 SUNDANCE_VERB_HIGH(tab0 << "CurveNormExpr::internalFindW() for " << toString());
00086 Set<MultipleDeriv> rtn;
00087
00088 if (order==0) rtn.put(MultipleDeriv());
00089
00090 if (order==1)
00091 {
00092
00093
00094
00095
00096
00097 }
00098
00099 SUNDANCE_VERB_HIGH(tab0 << "W[" << order << "]=" << rtn);
00100 return rtn;
00101 }
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141