Dear all,
I am writing code using hp class. When accessing the solution as:
/*fe_c and q_c are fe and q collections*/
hp::FEValue hp_fv(fe_c,q_c,update_values | update_quadrature_points);
for(cells)
{
const FEValues &fv = hp_fv.get_present_fe_values();
std::vector local_solution(q_c[0].size())
Okay, I found sth interesting: I have to use sort of component mask to make
the error gone:
for (cells)
{
hp_fv.reinit(cell);
const FEValues &fv = hp_fv.get_present_fe_values();
FEValuesExtractors::Scalar tmp(cell->active_fe_index());
fv[tmp].get_function_values(global_solution,local_solut
All right, my understanding is that fv now is a FESystem object and have to
be used correspondingly. So I tried and it worked.
hp::FEValue hp_fv(fe_c,q_c,update_values | update_quadrature_points);
for(cells)
{
const FEValues &fv = hp_fv.get_present_fe_values();
FEValuesExtractors::Scalar tmp(