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