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 [email protected].
For more options, visit https://groups.google.com/d/optout.
