Correction of a bug causing duplicated edges in opnbdy mode

Author: dapogny

src/mmg2d/hash_2d.c

@@ -600,14 +600,15 @@ int MMG2D_pack(MMG5_pMesh mesh,MMG5_pSol sol,MMG5_pSol met) {
       else if ( (pt->ref==pt1->ref) && MG_SIN(pt->tag[i]) ) {
         ++mesh->na;
       }
-      else if ( mesh->info.opnbdy ) {
+      else if ( mesh->info.opnbdy && pt->ref == pt1->ref && k < iel ) {
         if ( (pt->tag[i] & MG_REF) || pt->tag[i] & MG_BDY ) {
           assert ( pt->tag[i] & (MG_REF+MG_BDY) );
           ++mesh->na;
         }
       }
     }
   }
+    
   /** Count edges stored in quadrangles */
   for (k=1; k<=mesh->nquad; k++) {
     pq = &mesh->quadra[k];
@@ -676,7 +677,7 @@ int MMG2D_pack(MMG5_pMesh mesh,MMG5_pSol sol,MMG5_pSol met) {
           pt1 = &mesh->tria[iel];
           if ( !iel || (pt->ref > pt1->ref) ||
                ((pt->ref==pt1->ref) && MG_SIN(pt->tag[i])) ||
-               (mesh->info.opnbdy && ((pt->tag[i] & MG_REF) || (pt->tag[i] & MG_BDY)))) {
+               (mesh->info.opnbdy && (pt->ref==pt1->ref) && (k<iel) && ((pt->tag[i] & MG_REF) || (pt->tag[i] & MG_BDY)))) {
             ++ned;
             ped = &mesh->edge[ned];
             ped->a = pt->v[i1];

src/mmg2d/inout_2d.c

@@ -1304,7 +1304,7 @@ int MMG2D_saveMesh(MMG5_pMesh mesh,const char *filename) {
     }
   }
 
-  /* edges */
+  /* Edges */
   nedreq = 0;
   if ( mesh->na ) {
     if(!bin) {

src/mmg3d/hash_3d.c

@@ -1216,6 +1216,7 @@ int MMG5_hGeom(MMG5_pMesh mesh) {
         i2  = MMG5_iprv2[i];
         kk  = adja[i] / 3;
         if ( (!kk) || pt->tag[i] & MG_NOM ) {
+          if ( !kk ) pt->tag[i] |= (MG_BDY + MG_GEO); // Dunno what to do here
           if ( pt->tag[i] & MG_NOM ) {
             if ( mesh->info.iso )
               pt->edg[i] = ( pt->edg[i] != 0 ) ?  -MMG5_abs(pt->edg[i]) : mesh->info.isoref;