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 #ifndef SUNDANCE_DISCRETEFUNCELEMENT_H
00032 #define SUNDANCE_DISCRETEFUNCELEMENT_H
00033
00034
00035 #include "SundanceDefs.hpp"
00036 #include "SundanceFuncElementBase.hpp"
00037 #include "SundanceDiscreteFuncEvaluator.hpp"
00038 #include "SundanceDiscreteFuncDataStub.hpp"
00039
00040 namespace Sundance
00041 {
00042 using namespace Sundance;
00043
00044
00045 using namespace Teuchos;
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059 class DiscreteFuncElement : public virtual EvaluatableExpr,
00060 public FuncElementBase,
00061 public virtual GenericEvaluatorFactory<DiscreteFuncElement, DiscreteFuncElementEvaluator>
00062 {
00063 public:
00064
00065 DiscreteFuncElement(const RCP<DiscreteFuncDataStub>& data,
00066 const std::string& name,
00067 const std::string& suffix,
00068 const FunctionIdentifier& fid,
00069 int myIndexIntoVector);
00070
00071
00072 virtual ~DiscreteFuncElement() {;}
00073
00074
00075
00076
00077 RCP<const DiscreteFuncDataStub> commonData() const {return commonData_;}
00078
00079
00080
00081 DiscreteFuncDataStub* commonData() {return commonData_.get();}
00082
00083
00084 int myIndex() const {return myIndex_;}
00085
00086
00087
00088 void addMultiIndex(const MultiIndex& newMi) const ;
00089
00090
00091
00092
00093
00094 int maxDiffOrderOnDiscreteFunctions() const {return 0;}
00095
00096
00097
00098
00099
00100 virtual bool hasDiscreteFunctions() const {return true;}
00101
00102
00103
00104
00105
00106 virtual bool hasTestFunctions() const {return false;}
00107
00108
00109 virtual Set<MultipleDeriv>
00110 internalFindW(int order, const EvalContext& context) const ;
00111
00112 virtual Set<MultipleDeriv>
00113 internalFindV(int order, const EvalContext& context) const ;
00114
00115 virtual Set<MultipleDeriv>
00116 internalFindC(int order, const EvalContext& context) const ;
00117
00118
00119 virtual RCP<Array<Set<MultipleDeriv> > >
00120 internalDetermineR(const EvalContext& context,
00121 const Array<Set<MultipleDeriv> >& RInput) const ;
00122
00123
00124 virtual XMLObject toXML() const ;
00125
00126
00127 const Set<MultiIndex>& multiIndexSet() const {return miSet_;}
00128
00129
00130 virtual RCP<ExprBase> getRcp() {return rcp(this);}
00131
00132
00133 bool lessThan(const ScalarExpr* other) const ;
00134
00135 private:
00136
00137 RCP<DiscreteFuncDataStub> commonData_;
00138
00139 mutable Set<MultiIndex> miSet_;
00140
00141 int myIndex_;
00142
00143
00144 };
00145 }
00146
00147 #endif