void Albany::PoissonsEquationProblem:: buildProblem( Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct> > meshSpecs, Albany::StateManager& stateMgr) { int physSets = meshSpecs.size(); *out << "Num MeshSpecs: " << physSets << '\n'; fm.resize(physSets); bool haveSidesets = false; *out << "Calling PoissonsEquationProblem::buildEvaluators" << '\n'; for (int ps = 0; ps < physSets; ++ps) { fm[ps] = Teuchos::rcp(new PHX::FieldManager<PHAL::AlbanyTraits>); buildEvaluators(*fm[ps], *meshSpecs[ps], stateMgr, BUILD_RESID_FM, Teuchos::null); if (meshSpecs[ps]->ssNames.size() > 0) haveSidesets = true; } constructDirichletEvaluators(*meshSpecs[0]); if( haveSidesets ) constructNeumannEvaluators(meshSpecs[0]); if( params->isType<Teuchos::RCP<ATO::Topology> >("Topology") ) setupTopOpt(meshSpecs,stateMgr); }
//------------------------------------------------------------------------------ void Albany::MechanicsProblem:: buildProblem( Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct>> meshSpecs, Albany::StateManager& stateMgr) { // Construct All Phalanx Evaluators int physSets = meshSpecs.size(); *out << "Num MeshSpecs: " << physSets << '\n'; fm.resize(physSets); bool haveSidesets = false; *out << "Calling MechanicsProblem::buildEvaluators" << '\n'; for (int ps = 0; ps < physSets; ++ps) { fm[ps] = Teuchos::rcp(new PHX::FieldManager<PHAL::AlbanyTraits>); buildEvaluators(*fm[ps], *meshSpecs[ps], stateMgr, BUILD_RESID_FM, Teuchos::null); if (meshSpecs[ps]->ssNames.size() > 0) haveSidesets = true; } constructDirichletEvaluators(*meshSpecs[0]); if (haveSidesets) constructNeumannEvaluators(meshSpecs[0]); }
void Albany::HeatProblem:: buildProblem( Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct> > meshSpecs, Albany::StateManager& stateMgr) { /* Construct All Phalanx Evaluators */ int physSets = meshSpecs.size(); std::cout << "Heat Problem Num MeshSpecs: " << physSets << std::endl; fm.resize(physSets); for (int ps=0; ps<physSets; ps++) { fm[ps] = Teuchos::rcp(new PHX::FieldManager<PHAL::AlbanyTraits>); buildEvaluators(*fm[ps], *meshSpecs[ps], stateMgr, BUILD_RESID_FM, Teuchos::null); } if(meshSpecs[0]->nsNames.size() > 0) // Build a nodeset evaluator if nodesets are present constructDirichletEvaluators(meshSpecs[0]->nsNames); if(meshSpecs[0]->ssNames.size() > 0) // Build a sideset evaluator if sidesets are present constructNeumannEvaluators(meshSpecs[0]); }
void Hydrology::buildProblem (Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct> > meshSpecs, Albany::StateManager& stateMgr) { // Building cell basis and cubature const CellTopologyData * const cell_top = &meshSpecs[0]->ctd; intrepidBasis = Albany::getIntrepid2Basis(*cell_top); cellType = Teuchos::rcp(new shards::CellTopology (cell_top)); Intrepid2::DefaultCubatureFactory cubFactory; cubature = cubFactory.create<PHX::Device, RealType, RealType>(*cellType, meshSpecs[0]->cubatureDegree); elementBlockName = meshSpecs[0]->ebName; const int worksetSize = meshSpecs[0]->worksetSize; const int numCellVertices = cellType->getNodeCount(); const int numCellNodes = intrepidBasis->getCardinality(); const int numCellQPs = cubature->getNumPoints(); dl = Teuchos::rcp(new Albany::Layouts(worksetSize,numCellVertices,numCellNodes,numCellQPs,numDim)); /* Construct All Phalanx Evaluators */ TEUCHOS_TEST_FOR_EXCEPTION(meshSpecs.size()!=1,std::logic_error,"Problem supports one Material Block"); fm.resize(1); fm[0] = Teuchos::rcp(new PHX::FieldManager<PHAL::AlbanyTraits>); buildEvaluators(*fm[0], *meshSpecs[0], stateMgr, Albany::BUILD_RESID_FM,Teuchos::null); if(meshSpecs[0]->nsNames.size() > 0) { // Build a nodeset evaluator if nodesets are present constructDirichletEvaluators(*meshSpecs[0]); } if(meshSpecs[0]->ssNames.size() > 0) { // Build a sideset evaluator if sidesets are present constructNeumannEvaluators(meshSpecs[0]); } }
void Albany::HMCProblem:: buildProblem( Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct>> meshSpecs, Albany::StateManager& stateMgr) { /* Construct All Phalanx Evaluators */ TEUCHOS_TEST_FOR_EXCEPTION(meshSpecs.size()!=1,std::logic_error,"Problem supports one Material Block"); fm.resize(1); fm[0] = Teuchos::rcp(new PHX::FieldManager<PHAL::AlbanyTraits>); buildEvaluators(*fm[0], *meshSpecs[0], stateMgr, BUILD_RESID_FM, Teuchos::null); if(meshSpecs[0]->nsNames.size() > 0) // Build a nodeset evaluator if nodesets are present constructDirichletEvaluators(*meshSpecs[0]); if(meshSpecs[0]->ssNames.size() > 0) // Build a sideset evaluator if sidesets are present constructNeumannEvaluators(meshSpecs[0]); #ifdef ALBANY_ATO if( params->isType<Teuchos::RCP<ATO::Topology>>("Topology") ) setupTopOpt(meshSpecs,stateMgr); #endif }
void FELIX::Hydrology::buildProblem (Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct> > meshSpecs, Albany::StateManager& stateMgr) { TEUCHOS_TEST_FOR_EXCEPTION (meshSpecs.size()!=1,std::logic_error,"Problem supports one Material Block"); fm.resize(1); fm[0] = Teuchos::rcp(new PHX::FieldManager<PHAL::AlbanyTraits>); // Build evaluators buildEvaluators(*fm[0], *meshSpecs[0], stateMgr, Albany::BUILD_RESID_FM, Teuchos::null); // Build bc evaluators if(meshSpecs[0]->nsNames.size() > 0) // Build a nodeset evaluator if nodesets are present constructDirichletEvaluators(*meshSpecs[0]); if(meshSpecs[0]->ssNames.size() > 0) // Build a sideset evaluator if sidesets are present constructNeumannEvaluators(meshSpecs[0]); }