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 "SundanceSymbolicFunc.hpp"
00032 #include "SundanceSymbolicFuncElement.hpp"
00033 #include "SundanceZeroExpr.hpp"
00034
00035 #include "SundanceDerivSet.hpp"
00036 #include "PlayaTabs.hpp"
00037
00038 using namespace Sundance;
00039 using namespace Sundance;
00040
00041 using namespace Sundance;
00042 using namespace Teuchos;
00043
00044 SymbolicFunc::SymbolicFunc(const FunctionWithID& fid,
00045 const RCP<const CommonFuncDataStub>& data)
00046 : ListExpr(), FunctionWithID(fid), commonData_(data)
00047 {}
00048
00049
00050 void SymbolicFunc::substituteZero() const
00051 {
00052 for (int i=0; i<this->size(); i++)
00053 {
00054 const SymbolicFuncElement* u
00055 = dynamic_cast<const SymbolicFuncElement*>(element(i).ptr().get());
00056 TEUCHOS_TEST_FOR_EXCEPTION(u==0, std::logic_error,
00057 "Non-symbolic function "
00058 << element(i).toString()
00059 << " detected in SymbolicFunc::substituteZero()");
00060 u->substituteZero();
00061 }
00062 }
00063
00064 void SymbolicFunc
00065 ::substituteFunction(const RCP<DiscreteFunctionStub>& u0) const
00066 {
00067 TEUCHOS_TEST_FOR_EXCEPTION(this->size() != u0->size(), std::logic_error,
00068 "Mismatch between sizes of symbolic " << toString()
00069 << " and discrete func " << u0->toString()
00070 << " in substituteFunction()");
00071
00072 for (int i=0; i<this->size(); i++)
00073 {
00074 const SymbolicFuncElement* u
00075 = dynamic_cast<const SymbolicFuncElement*>(element(i).ptr().get());
00076 TEUCHOS_TEST_FOR_EXCEPTION(u==0, std::logic_error,
00077 "Non-symbolic function "
00078 << element(i).toString()
00079 << " detected in SymbolicFunc::substituteFunction()");
00080
00081 RCP<DiscreteFuncElement> df
00082 = rcp_dynamic_cast<DiscreteFuncElement>(u0->element(i).ptr());
00083 u->substituteFunction(df);
00084 }
00085 }
00086