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<dim> &fv = hp_fv.get_present_fe_values();
FEValuesExtractors::Scalar tmp(cell->active_fe_index());
fv[tmp].get_function_values(global_solution,local_solution);
/*do something to the local solution*/
}
When using hp class previously (a year ago), I directly used
fv.get_function_values, which worked.
在 2017年6月12日星期一 UTC-7下午3:25:13，Weixiong Zheng写道：
>
> 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<dim> hp_fv(fe_c,q_c,update_values | update_quadrature_points);
> for(cells)
> {
> const FEValues<dim> &fv = hp_fv.get_present_fe_values();
> std::vector<double> local_solution(q_c[0].size());//I am using the same
> quadrature rule for all
> fv.get_function_values(/*Vector<double>*/ global_solution,
> local_solution);
> /*do something to local_solution*/
> }
> The problem is a scalar problem. It is about solving the same equation in
> different subdomains with different finite elements.
>
> I then got the following error message. Any ideas?
>
> Thanks.
> Weixiong
>
> An error occurred in line <2717> of file <../source/fe/fe_values.cc> in
> function
>
> void dealii::FEValuesBase<2, 2>::get_function_values(const InputVector
> &, std::vector<typename InputVector::value_type> &) const [dim = 2,
> spacedim = 2, InputVector = dealii::Vector<double>]
>
> The violated condition was:
>
> (fe->n_components()) == (1)
>
> The name and call sequence of the exception was:
>
> dealii::ExcDimensionMismatch((fe->n_components()),(1))
>
> Additional Information:
>
> Dimension 2 not equal to 1
>
--
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.
