ElemAssembly& RBAssemblyExpansion::get_F_assembly(unsigned int q) { if(q >= get_n_F_terms()) libmesh_error_msg("Error: We must have q < get_n_F_terms in get_F_assembly."); return *_F_assembly_vector[q]; }
void RBAssemblyExpansion::perform_F_boundary_assembly(unsigned int q, FEMContext& context) { if(q >= get_n_F_terms()) libmesh_error_msg("Error: We must have q < get_n_F_terms in perform_F_interior_assembly."); libmesh_assert(_A_assembly_vector[q]); return _F_assembly_vector[q]->boundary_assembly( context ); }
Number RBThetaExpansion::eval_F_theta(unsigned int q, const RBParameters & mu) { if (q >= get_n_F_terms()) libmesh_error_msg("Error: We must have q < get_n_F_terms in eval_F_theta."); libmesh_assert(_F_theta_vector[q]); return _F_theta_vector[q]->evaluate( mu ); }
Number RBThetaExpansion::eval_F_theta(unsigned int q, const RBParameters& mu) { if(q >= get_n_F_terms()) { libMesh::err << "Error: We must have q < get_n_F_terms in eval_F_theta." << std::endl; libmesh_error(); } libmesh_assert(_F_theta_vector[q]); return _F_theta_vector[q]->evaluate( mu ); }