Loop only finite faces in the volume computation

For outward wrapping
This commit is contained in:
Mael Rouxel-Labbé 2023-10-16 17:02:05 +02:00
parent 116bb2beb9
commit 18154a758b
1 changed files with 1 additions and 1 deletions

View File

@ -1632,7 +1632,7 @@ public:
auto wrap_volume = [&]() auto wrap_volume = [&]()
{ {
FT vol = 0; FT vol = 0;
for(Cell_handle ch : m_tr.all_cell_handles()) for(Cell_handle ch : m_tr.finite_cell_handles())
if(!ch->is_outside()) if(!ch->is_outside())
vol += volume(m_tr.point(ch, 0), m_tr.point(ch, 1), m_tr.point(ch, 2), m_tr.point(ch, 3)); vol += volume(m_tr.point(ch, 0), m_tr.point(ch, 1), m_tr.point(ch, 2), m_tr.point(ch, 3));