// ========================================================================
void process_nodesets(Ioss::Region ®ion, stk::mesh::MetaData &meta)
{
const Ioss::NodeSetContainer& node_sets = region.get_nodesets();
stk::io::default_part_processing(node_sets, meta, stk::topology::NODE_RANK);
/** \todo REFACTOR should "distribution_factor" be a default field
* that is automatically declared on all objects that it exists
* on as is done in current framework?
*/
stk::mesh::Field<double> & distribution_factors_field =
meta.declare_field<stk::mesh::Field<double> >(stk::topology::NODE_RANK, "distribution_factors");
/** \todo REFACTOR How to associate distribution_factors field
* with the nodeset part if a node is a member of multiple
* nodesets
*/
for(Ioss::NodeSetContainer::const_iterator it = node_sets.begin();
it != node_sets.end(); ++it) {
Ioss::NodeSet *entity = *it;
if (stk::io::include_entity(entity)) {
stk::mesh::Part* const part = meta.get_part(entity->name());
STKIORequire(part != NULL);
STKIORequire(entity->field_exists("distribution_factors"));
stk::mesh::put_field(distribution_factors_field, *part);
/** \todo IMPLEMENT truly handle fields... For this case we
* are just defining a field for each transient field that is
* present in the mesh...
*/
stk::io::define_io_fields(entity, Ioss::Field::TRANSIENT,
*part,
part->primary_entity_rank() ) ;
}
}
}
