ElemAssembly& RBAssemblyExpansion::get_A_assembly(unsigned int q) { if(q >= get_n_A_terms()) libmesh_error_msg("Error: We must have q < get_n_A_terms in get_A_assembly."); return *_A_assembly_vector[q]; }
void RBAssemblyExpansion::perform_A_boundary_assembly(unsigned int q, FEMContext& context) { if(q >= get_n_A_terms()) libmesh_error_msg("Error: We must have q < get_n_A_terms in perform_A_boundary_assembly."); libmesh_assert(_A_assembly_vector[q]); return _A_assembly_vector[q]->boundary_assembly( context ); }
Number RBThetaExpansion::eval_A_theta(unsigned int q, const RBParameters & mu) { if (q >= get_n_A_terms()) libmesh_error_msg("Error: We must have q < get_n_A_terms in eval_A_theta."); libmesh_assert(_A_theta_vector[q]); return _A_theta_vector[q]->evaluate( mu ); }
Number RBThetaExpansion::eval_A_theta(unsigned int q, const RBParameters& mu) { if(q >= get_n_A_terms()) { libMesh::err << "Error: We must have q < get_n_A_terms in eval_A_theta." << std::endl; libmesh_error(); } libmesh_assert(_A_theta_vector[q]); return _A_theta_vector[q]->evaluate( mu ); }