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 "SundanceExpr.hpp"
00032 #include "SundanceScalarExpr.hpp"
00033 #include "SundanceSymbolicTransformation.hpp"
00034
00035
00036 using namespace Sundance;
00037 using namespace Sundance;
00038
00039 using namespace Teuchos;
00040 using namespace Sundance;
00041
00042 SymbolicTransformation::SymbolicTransformation()
00043 {}
00044
00045 RCP<ScalarExpr> SymbolicTransformation::chooseSign(int sign,
00046 const RCP<ScalarExpr>& expr)
00047 {
00048
00049
00050 switch(sign)
00051 {
00052 case 1:
00053 return expr;
00054 case -1:
00055 {
00056 Expr e = -Expr::handle(expr);
00057 RCP<ScalarExpr> rtn = rcp_dynamic_cast<ScalarExpr>(e.ptr());
00058 TEUCHOS_TEST_FOR_EXCEPTION(rtn.get() == NULL, std::logic_error,
00059 "Non-scalar expr "
00060 << e.toString()
00061 << " detected in SymbolicTransformation::chooseSign");
00062 return rtn;
00063 }
00064 default:
00065 TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error,
00066 "sign != +/- 1 in Expr::transformSign()");
00067 }
00068 return expr;
00069 }
00070
00071 Expr SymbolicTransformation::chooseSign(int sign,
00072 const Expr& expr)
00073 {
00074
00075
00076 switch(sign)
00077 {
00078 case 1:
00079 return expr;
00080 case -1:
00081 return -expr;
00082 default:
00083 TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error,
00084 "sign != +/- 1 in Expr::transformSign()");
00085 }
00086 return expr;
00087 }
00088
00089 RCP<ScalarExpr> SymbolicTransformation::getScalar(const Expr& expr)
00090 {
00091 RCP<ScalarExpr> s = rcp_dynamic_cast<ScalarExpr>(expr.ptr());
00092
00093 TEUCHOS_TEST_FOR_EXCEPTION(s.get()==NULL, std::logic_error,
00094 "non-scalar detected in SymbolicTransformation::getScalar");
00095
00096 return s;
00097 }