Thanks, Daniel,
I used FESystem<dim> to wrap the scalar FiniteElement and the FE_Nothing
together. Here are the hp-relevant variables:
FESystem<dim> fe_continuous;
FESystem<dim> fe_discontinuous;
hp::FECollection<dim> fe_c;
hp::QCollection<dim> q_c;
hp::QCollection<dim-1> q_face_c;
This is how I do initialization
Equation<dim>::Equation (ParameterHandler &prm)
:
ProblemDefinition<dim>(prm),
fe_order(prm.get_integer("finite element polynomial degree")),
fe_continuous(FE_Q<dim>(fe_order),1,
FE_Nothing<dim>(),1),
fe_discontinuous(FE_Nothing<dim>(),1,
FE_DGQ<dim>(fe_order),1),
fe_c(hp::FECollection<dim>(fe_continuous,fe_discontinuous)),
dof_handler (triangulation)
{
q_c.push_back(QGauss<dim>(fe_order+1));
q_c.push_back(QGauss<dim>(fe_order+1));
q_face_c.push_back(QGauss<dim>(fe_order+1));
q_face_c.push_back(QGauss<dim>(fe_order+1));
}
As the problem is a scalar problem, fv.shape_value and shape_grad work as
they are supposed to, but not the get_function_values.
在 2017年6月13日星期二 UTC-7上午2:58:09，Daniel Arndt写道：
>
> Weixiong,
>
> [...]
>>
> Then I got confused:
>> 1. I used hp class a year ago or so, fv.get_function_values works
>> directly. Now I have to use the above change. Did I do something wrong?
>> 2. Then why fv.shape_value still works as it did? Isn't fv still a
>> FESystem object?
>>
> How did you initialize hp::DoFHandler? What is the hp::FECollectionyou
> used?
> If you are solving a scalar problem using the same equation in different
> subdomains with different finite elements,
> I would assume that the Finite Elemnts in the hp::FECollection are also
> scalar valued and have just one component.
> However, from what you write above it seems that you initialized the
> hp::FECollection with vector valued FiniteElements.
> Can you confirm either?
>
> Best,
> Daniel
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
