message: Simplify the algorithm to delete adjacenceies for facets on internal borders
diff: @@ -1215,29 +1215,21 @@ namespace RINGMesh {
  const std::vector< Box3d >& surface_boxes )
  {
  const Surface& S = geomodel_.surface( surface_id );
  - std::vector< index_t > adjacet_facets_id( 3 * S.nb_mesh_elements() );
  for( index_t f = 0; f < S.nb_mesh_elements(); ++f ) {
  + std::vector< index_t > adjacent_facets_id( 3 );
  for( index_t e = 0; e < 3; ++e ) {
  - adjacet_facets_id[3 * f + e] = S.facet_adjacent_index( f, e );
  + adjacent_facets_id[e] = S.facet_adjacent_index( f, e );
  if( !S.is_on_border( f, e ) ) {
  bool internal_border = is_edge_in_several_surfaces( geomodel_,
  surface_id, f, e, surface_nns, surface_boxes );
  if( internal_border ) {
  - adjacet_facets_id[3 * f + e] = GEO::NO_FACET;
  + adjacent_facets_id[e] = NO_ID;
  }
  }
  }
  - }
  -
  - for( index_t f = 0; f < S.nb_mesh_elements(); f++ ) {
  - std::vector< index_t >::const_iterator first = adjacet_facets_id.begin()
  - + 3 * f;
  - std::vector< index_t >::const_iterator last = adjacet_facets_id.begin()
  - + 3 * f + 3;
  - std::vector< index_t > cur_facet_adjacencies( first, last );
  geometry.set_surface_element_adjacency( surface_id, f,
  - cur_facet_adjacencies );
  + adjacent_facets_id );
  }
  }
  