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_SYMBOLICFUNCELEMENT_H
00032 #define SUNDANCE_SYMBOLICFUNCELEMENT_H
00033
00034
00035 #include "SundanceDefs.hpp"
00036 #include "SundanceFuncElementBase.hpp"
00037 #include "SundanceEvaluatableExpr.hpp"
00038 #include "SundanceDiscreteFuncElement.hpp"
00039 #include "SundanceDiscreteFunctionStub.hpp"
00040 #include "SundanceSymbolicFuncEvaluator.hpp"
00041 #include "SundanceSymbolicFuncDescriptor.hpp"
00042 #include "SundanceCommonFuncDataStub.hpp"
00043
00044 namespace Sundance
00045 {
00046 using namespace Sundance;
00047 class DiscreteFuncElement;
00048 using namespace Teuchos;
00049
00050
00051
00052
00053
00054
00055
00056
00057 class SymbolicFuncElement : public FuncElementBase,
00058 public SymbolicFuncDescriptor,
00059 virtual public EvaluatableExpr,
00060 public GenericEvaluatorFactory<SymbolicFuncElement, SymbolicFuncElementEvaluator>
00061 {
00062 public:
00063
00064 SymbolicFuncElement(const std::string& name,
00065 const std::string& suffix,
00066 const FunctionIdentifier& fid,
00067 const RCP<const CommonFuncDataStub>& data);
00068
00069
00070 virtual ~SymbolicFuncElement() {;}
00071
00072
00073 void accumulateFuncSet(Set<int>& funcDofIDs,
00074 const Set<int>& activeSet) const ;
00075
00076
00077 virtual bool hasTestFunctions() const {return false;}
00078
00079
00080
00081
00082
00083
00084
00085
00086 void substituteZero() const ;
00087
00088
00089
00090
00091
00092
00093
00094 void substituteFunction(const RCP<DiscreteFuncElement>& u0) const ;
00095
00096
00097
00098 const EvaluatableExpr* evalPt() const {return evalPt_.get();}
00099
00100
00101
00102 EvaluatableExpr* evalPt() {return evalPt_.get();}
00103
00104
00105
00106 bool evalPtIsZero() const ;
00107
00108
00109 const RCP<const CommonFuncDataStub>& commonData() const {return commonData_;}
00110
00111
00112
00113
00114
00115 virtual Set<MultipleDeriv>
00116 internalFindW(int order, const EvalContext& context) const ;
00117
00118
00119 virtual Set<MultipleDeriv>
00120 internalFindV(int order, const EvalContext& context) const ;
00121
00122
00123 virtual Set<MultipleDeriv>
00124 internalFindC(int order, const EvalContext& context) const ;
00125
00126
00127 virtual RCP<Array<Set<MultipleDeriv> > >
00128 internalDetermineR(const EvalContext& context,
00129 const Array<Set<MultipleDeriv> >& RInput) const ;
00130
00131 virtual void registerSpatialDerivs(const EvalContext& context,
00132 const Set<MultiIndex>& miSet) const ;
00133
00134
00135
00136
00137
00138 virtual bool isIndependentOf(const Expr& u) const ;
00139
00140
00141
00142
00143 virtual bool isLinearForm(const Expr& u) const ;
00144
00145
00146 virtual RCP<ExprBase> getRcp() {return rcp(this);}
00147
00148 private:
00149 RCP<const CommonFuncDataStub> commonData_;
00150
00151 mutable RCP<EvaluatableExpr> evalPt_;
00152
00153 mutable Array<int> evalPtDerivSetIndices_;
00154 };
00155 }
00156
00157 #endif