Esempio n. 1
0
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];
}
Esempio n. 2
0
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 );
}
Esempio n. 3
0
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 );
}
Esempio n. 4
0
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 );
}