On 12/20/20 2:58 AM, Alberto Salvadori wrote:
Point<dim> centroid = cell->face(face_number)->barycenter() ;
provides always centroid = {0,0,0}, whereas the following
Point<dim> centroid = cell->face(face_number)->center() ;
returns the correct amount.
Am I misunderstanding something or may it be a bug?
I suspect that you forgot to try running your code in debug mode :-) In debug
mode, I'm pretty sure you would run into this assertion:
template <int structdim, int dim, int spacedim>
Point<spacedim>
barycenter(const TriaAccessor<structdim, dim, spacedim> &)
{
// this function catches all the cases not
// explicitly handled above
Assert(false, ExcNotImplemented());
return Point<spacedim>();
}
In other words, the function is simply not implemented -- though I suspect
that it would not be very difficult to do so.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/c6c8111e-a53a-da7a-4b5e-40b870546283%40colostate.edu.
