Dear Deepak,
Have you run your code in debug mode to confirm that no errors are thrown?
I have used the Kelly error estimator in the context of an hp problem with
no issues. Have you tried increasing the number of face quadrature points
(I notice that qrule starts at 1, but without knowing which type and degree
FE you are using, its hard to judge if the way you choose the quadrature
order is correct).
Regards,
J-P
On Friday, August 26, 2016 at 1:23:03 PM UTC+2, Deepak Gupta wrote:
>
> Dear All,
>
> I am trying to use KellyErrorEstimator of deal.II and I saw various
> tutorials which use it. I ran those examples and they work well.
> Recently, I tried it with my work and when I print the
> estimated_error_per_cell vector, all of the values in that vector are
> zero. I have been trying to find the bug for past two days but all in
> vain. May be someone can comment (if they also encountered ever) on why
> this error can be zero, based on the way it is implemented in deal.II.
> The piece of code is as follows: (this is an hp-adaptive framework):
>
> //Quadrature collection for FE
> hp::QCollection<dim-1> quad_collection;
> for (unsigned int qrule = 1; qrule <= (fem->fe_collection).size();
> ++qrule){
> quad_collection.push_back(QGauss<dim-1>(qrule));
> }
>
> std::cout<<"No. of dofs : "<<fem->dof_handler->n_dofs()<<std::endl;
> Vector<float> estimated_error_per_cell
> ((*(fem->triangulation)).n_active_cells());
> std::cout<<"Size of error vector :
> "<<estimated_error_per_cell.size()<<std::endl;
> KellyErrorEstimator<dim>::estimate (*(fem->dof_handler),
> quad_collection,
> typename
> FunctionMap<dim>::type(),
> fem->solution,
> estimated_error_per_cell);
>
>
> dof_hander->n_dofs() gives correct output. Size of
> estimated_error_per_cell is correct. fem-> solution vector is correct
> (cross-checked). A minimal working example for my problem would be
> similar to step-27 but then step-27 works fine.
>
>
> I hope someone has seen such an issue and can help.
>
> Best
>
> Deepak
>
> --
> Deepak K. Gupta
>
> PhD Candidate/Researcher
> Structural Optimization & Mechanics
> TU Delft / Precision & Microsystems Engineering
> Faculty of Mechanical, Maritime & Materials Engineering (3mE)
> Mekelweg 2, 2628 CD Delft, The Netherlands
>
> T +31 (0)15 27 86818
> [email protected]
>
> Room 3mE.G.1.150
>
>
--
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.
