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.
