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 "SundanceEssentialBC.hpp"
00032 #include "SundanceSumOfBCs.hpp"
00033 #include "SundanceZeroExpr.hpp"
00034
00035 using namespace Sundance;
00036 using namespace Sundance;
00037
00038 using namespace Teuchos;
00039 using namespace Sundance;
00040
00041 using namespace Sundance;
00042
00043 Expr Sundance::EssentialBC(const Handle<CellFilterStub>& domain,
00044 const Expr& integrand,
00045 const WatchFlag& watch)
00046 {
00047 const ZeroExpr* z = dynamic_cast<const ZeroExpr*>(integrand.ptr().get());
00048 const ConstantExpr* c = dynamic_cast<const ConstantExpr*>(integrand.ptr().get());
00049 if (z != 0 || (c != 0 && c->value()==0.0))
00050 {
00051 return integrand;
00052 }
00053 RCP<QuadratureFamilyStub> quad = QuadratureFamilyStub::defaultQuadrature();
00054 return new SumOfBCs(domain.ptr(), Re(integrand), quad,
00055 ParametrizedCurve::returnDummyCurve() , watch);
00056 }
00057
00058 Expr Sundance::EssentialBC(const Handle<CellFilterStub>& domain,
00059 const Expr& integrand,
00060 const Handle<QuadratureFamilyStub>& quad,
00061 const WatchFlag& watch)
00062 {
00063 const ZeroExpr* z = dynamic_cast<const ZeroExpr*>(integrand.ptr().get());
00064 const ConstantExpr* c = dynamic_cast<const ConstantExpr*>(integrand.ptr().get());
00065 if (z != 0 || (c != 0 && c->value()==0.0))
00066 {
00067 return integrand;
00068 }
00069 return new SumOfBCs(domain.ptr(), Re(integrand), quad.ptr(),
00070 ParametrizedCurve::returnDummyCurve(), watch);
00071 }
00072
00073
00074 Expr Sundance::EssentialBC(const Handle<CellFilterStub>& domain,
00075 const Expr& integrand,
00076 const Handle<QuadratureFamilyStub>& quad,
00077 const ParametrizedCurve& curve,
00078 const WatchFlag& watch)
00079 {
00080 const ZeroExpr* z = dynamic_cast<const ZeroExpr*>(integrand.ptr().get());
00081 const ConstantExpr* c = dynamic_cast<const ConstantExpr*>(integrand.ptr().get());
00082 if (z != 0 || (c != 0 && c->value()==0.0))
00083 {
00084 return integrand;
00085 }
00086 return new SumOfBCs(domain.ptr(), Re(integrand), quad.ptr(), curve, watch);
00087 }