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 "PlayaTabs.hpp"
00032 #include "SundanceOut.hpp"
00033 #include "SundanceRegionQuadCombo.hpp"
00034
00035 using namespace Sundance;
00036 using namespace Teuchos;
00037
00038
00039 RegionQuadCombo::RegionQuadCombo()
00040 : id_(-1), domain_(), quad_(), paramCurve_() , watch_()
00041 {;}
00042
00043 RegionQuadCombo::RegionQuadCombo(const RCP<CellFilterStub>& domain,
00044 const RCP<QuadratureFamilyStub>& quad,
00045 const ParametrizedCurve& paramCurve ,
00046 const WatchFlag& watch)
00047 : id_(getID(domain, quad,watch)), domain_(domain), quad_(quad),
00048 paramCurve_(paramCurve) , watch_(watch)
00049 {;}
00050
00051
00052 int RegionQuadCombo::getID(const RCP<CellFilterStub>& domain,
00053 const RCP<QuadratureFamilyStub>& quad,
00054 const WatchFlag& watch)
00055 {
00056 RegTriple p(domain, quad, watch);
00057
00058 if (!domainAndQuadToIDMap().containsKey(p))
00059 {
00060 int id = topID();
00061 domainAndQuadToIDMap().put(p, id);
00062 }
00063 return domainAndQuadToIDMap().get(p);
00064 }
00065
00066 string RegionQuadCombo::toString() const
00067 {
00068 TeuchosOStringStream oss;
00069 Tabs tabs;
00070
00071 oss << "Integration Region" << std::endl;
00072 {
00073 oss << tabs << "cell filter=" << domain_->description() << std::endl;
00074 oss << tabs << "quadrature rule=" << quad_->description() << std::endl;
00075 oss << tabs << "watchpoint=[" << watch().name() << "]" << std::endl;
00076 }
00077 return oss.str();
00078 }
00079
00080 Sundance::Map<RegTriple, int>& RegionQuadCombo::domainAndQuadToIDMap()
00081 {
00082 static Sundance::Map<RegTriple, int> rtn = Sundance::Map<RegTriple, int>();
00083 return rtn;
00084 }
00085
00086
00087
00088
00089