50 std::string qoi_list = input(
"QoI/enabled_qois",
"none" );
52 std::vector<std::string> qoi_names;
54 if( qoi_list != std::string(
"none") )
59 std::tr1::shared_ptr<CompositeQoI> qois(
new CompositeQoI );
61 if( !qoi_names.empty() )
63 for( std::vector<std::string>::const_iterator name = qoi_names.begin();
64 name != qoi_names.end(); ++name )
66 this->
add_qoi( input, *name, qois );
71 if( input(
"screen-options/echo_qoi",
false ) )
80 void QoIFactory::add_qoi(
const GetPot& ,
const std::string& qoi_name, std::tr1::shared_ptr<CompositeQoI>& qois )
101 libMesh::err <<
"Error: Invalid QoI name " << qoi_name << std::endl;
107 qois->add_qoi( *qoi );
113 const std::string& qoi_name )
115 int num_physics = input.vector_variable_size(
"Physics/enabled_physics");
118 libmesh_assert(num_physics > 0);
120 std::set<std::string> requested_physics;
121 std::set<std::string> required_physics;
124 for(
int i = 0; i < num_physics; i++ )
126 requested_physics.insert( input(
"Physics/enabled_physics",
"NULL", i ) );
144 std::cout <<
"==========================================================" << std::endl
145 <<
"List of Enabled QoIs:" << std::endl;
147 for(
unsigned int q = 0; q < qois->n_qois(); q++ )
149 std::cout << qois->get_qoi(q).name() << std::endl;
152 std::cout <<
"==========================================================" << std::endl;
158 const std::set<std::string>& required_physics,
159 const std::string& qoi_name )
161 bool physics_found =
false;
162 for( std::set<std::string>::const_iterator name = required_physics.begin();
163 name != required_physics.end();
166 if( requested_physics.find( (*name) ) != requested_physics.end() )
167 physics_found =
true;
177 const std::set<std::string>& required_physics )
179 libMesh::err <<
"================================================================" << std::endl
180 <<
"QoI " << qoi_name << std::endl
181 <<
"requires one of the following physics which were not found:" <<std::endl;
183 for( std::set<std::string>::const_iterator name = required_physics.begin();
184 name != required_physics.end();
187 libMesh::err << *name << std::endl;
190 libMesh::err <<
"================================================================" << std::endl;
virtual void check_qoi_physics_consistency(const GetPot &input, const std::string &qoi_name)
void consistency_error_msg(const std::string &qoi_name, const std::set< std::string > &required_physics)
const std::string avg_nusselt
virtual void add_qoi(const GetPot &input, const std::string &qoi_name, std::tr1::shared_ptr< CompositeQoI > &qois)
const std::string vorticity
virtual std::tr1::shared_ptr< CompositeQoI > build(const GetPot &input)
virtual void echo_qoi_list(std::tr1::shared_ptr< CompositeQoI > &qois)
const PhysicsName low_mach_navier_stokes
const PhysicsName heat_transfer
void split_string(const std::string &input, const std::string &delimiter, std::vector< std::string > &results)
const std::string parsed_interior
void consistency_helper(const std::set< std::string > &requested_physics, const std::set< std::string > &required_physics, const std::string &qoi_name)