Panzer_BCStrategy_Dirichlet_Utilities.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//           Panzer: A partial differential equation assembly
//       engine for strongly coupled complex multiphysics systems
//
// Copyright 2011 NTESS and the Panzer contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef PANZER_BCSTRATEGY_DIRCHLET_UTILITIES_HPP
#define PANZER_BCSTRATEGY_DIRCHLET_UTILITIES_HPP
 
#include "Phalanx_FieldManager.hpp"
#include "Panzer_PhysicsBlock.hpp"
#include "Panzer_String_Utilities.hpp"
#include "Teuchos_ParameterList.hpp"
#include <vector>
#include <list>
#include <string>
#include <algorithm>
 
namespace panzer {
 
 
template<typename EvalT>
void 
buildAndRegisterSubsetOfClosureModelEvaluatorsForType(PHX::FieldManager<panzer::Traits>& fm,
                  const panzer::PhysicsBlock& pb,
                  const panzer::ClosureModelFactory_TemplateManager<panzer::Traits>& factory,
                  const std::string comma_sep_closure_model_list,
                  const Teuchos::ParameterList& models,
                  const Teuchos::ParameterList& user_data)
{
  
  std::vector<std::string> closure_model_vector;
  panzer::StringTokenizer(closure_model_vector,comma_sep_closure_model_list,",");
  
  // copy into list
  std::list<std::string> closure_model_list;
  for (std::vector<std::string>::iterator i=closure_model_vector.begin(); i != closure_model_vector.end(); ++i)
    closure_model_list.push_back(*i);
  
  Teuchos::ParameterList models_to_build;
  
  for (Teuchos::ParameterList::ConstIterator model = models.begin(); model != models.end(); ++model) {
    
    std::list<std::string>::iterator search = 
      std::find(closure_model_list.begin(), closure_model_list.end(), model->first);
    
    if (search != closure_model_list.end()) {
      closure_model_list.erase(search);
      models_to_build.sublist(model->first) = models.sublist(model->first);
    }
    else
      models_to_build.sublist(model->first);
  }
  
  TEUCHOS_TEST_FOR_EXCEPTION(closure_model_list.size() != 0, std::logic_error,
           "Error - the list of closure models \"" << comma_sep_closure_model_list << "\" contains an invalid model.");
 
  pb.buildAndRegisterClosureModelEvaluatorsForType<EvalT>(fm,factory,models_to_build,user_data);
 
}
  
} // namespace panzer
 
#endif