ImplDetails.cpp

/*------------------------------------------------------------------------*/
/*                 Copyright 2010 Sandia Corporation.                     */
/*  Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive   */
/*  license for use of this work by or on behalf of the U.S. Government.  */
/*  Export of this program may require a license from the                 */
/*  United States Government.                                             */
/*------------------------------------------------------------------------*/
 
 
#include <limits>
#include <sstream>
#include <stdexcept>
 
#include <stk_linsys/ImplDetails.hpp>
 
namespace stk_classic {
namespace linsys {
namespace impl {
 
int
map_field_to_int(FieldIdMap& field_id_map,
                 const stk_classic::mesh::FieldBase& field)
{
  FieldIdMap::iterator iter = field_id_map.find(&field);
 
  if (iter == field_id_map.end()) {
    iter = field_id_map.insert(iter, std::make_pair(&field,field_id_map.size()));
  }
 
  return iter->second;
}
 
int
query_field_to_int_mapping(const FieldIdMap& field_id_map,
                           const stk_classic::mesh::FieldBase& field)
{
  FieldIdMap::const_iterator iter = field_id_map.find(&field);
 
  if (iter == field_id_map.end()) {
    std::ostringstream msg;
    msg << "stk_classic::linsys::query_field_to_int_mapping ERROR: "
        << " field with name '"<<field.name()<<"' not found in field-to-int map.";
    std::string str = msg.str();
    throw std::runtime_error(str);
  }
 
  return iter->second;
}
 
const stk_classic::mesh::FieldBase*
get_field(const FieldIdMap& field_id_map,
          int field_id)
{
  FieldIdMap::const_iterator
    iter = field_id_map.begin(), iter_end = field_id_map.end();
 
  while(iter!=iter_end && iter->second != field_id) ++iter;
 
  if (iter == iter_end) {
    std::ostringstream msg;
    msg << "stk_classic::linsys::get_dof ERROR: "
     << "field_id ("<<field_id<<") returned from fei query is not mapped to a stk_classic::mesh::Field.";
    std::string str = msg.str();
    throw std::runtime_error(str);
  }
 
  return iter->first;
}
 
int entitytype_to_int(stk_classic::mesh::EntityRank entity_rank)
{
  int int_type = static_cast<int>(entity_rank);
 
  verify_convertible_to_int(entity_rank, "stk_classic::linsys::entitytype_to_int");
 
  return int_type;
}
 
int entityid_to_int(stk_classic::mesh::EntityId id)
{
  int int_id = static_cast<int>(id);
 
  verify_convertible_to_int(id, "stk_classic::linsys::entityid_to_int");
 
  return int_id;
}
 
}//namespace impl
}//namespace linsys
}//namespace stk_classic