The Higher Education and Research forge

Home My Page Projects Code Snippets Project Openings MMD
Summary Activity Forums Tracker Lists Tasks Docs Surveys News SCM Files

SCM Repository

authoralainbastide <alain.bastide@univ-reunion.fr>
Sat, 18 Jul 2020 02:37:45 +0000 (06:37 +0400)
committeralainbastide <alain.bastide@univ-reunion.fr>
Sat, 18 Jul 2020 02:37:45 +0000 (06:37 +0400)
src/primal.c

index 9b92a62..4aefbf0 100644 (file)
@@ -42,7 +42,29 @@ DEF_hexaedron_localNodeListNumbers
 DEF_hexaedron_localSegmentList
 DEF_hexaedron_localSegmentsListNumbers
 
+
+
+
+
 void setHexahedreVertex(struct primalMesh * myPrimalMesh)
+/*
+ * void setHexahedreVertex(struct primalMesh * myPrimalMesh)
+ * 
+ * Build vertex for each point without overlap
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use : L= myPrimalMesh->L
+ *          myPrimalMesh->M
+ *          myPrimalMesh->l
+ *          myPrimalMesh->N
+ *          myPrimalMesh->H
+ *          myPrimalMesh->P
+ * 
+ * Allocate : myPrimalMesh->vertex
+ * Produce : myPrimalMesh->vertex
+ * 
+ */
 {
   
   //    clock_t t;
@@ -80,6 +102,21 @@ void setHexahedreVertex(struct primalMesh * myPrimalMesh)
 }
 
 void setHexahedreCellToVertexNumber(struct primalMesh * myPrimalMesh)
+/*
+ * void setHexahedreCellToVertexNumber(struct primalMesh * myPrimalMesh)
+ * 
+ * This function builds a vector with the number of vertices of each cell; 
+ * here, we only use the hexahedron (HEXAHEDRON) but the code could be 
+ * updated to use several cell types
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use : 
+ * 
+ * Allocate : myPrimalMesh->cellToVertexNumber
+ * Produce : myPrimalMesh->cellToVertexNumber
+ * 
+ */
 {
   clock_t t=startFunction(__FUNCTION__);
   
@@ -104,54 +141,23 @@ void setHexahedreCellToVertexNumber(struct primalMesh * myPrimalMesh)
 }
 
 
-//void setHexahedreFaceToVertexNumber(struct primalMesh * myPrimalMesh)
-//{
-//  // vertex number on faces
-//  clock_t t=startFunction(__FUNCTION__);
-
-//  myPrimalMesh->faceToVertexNumber = (connectivity_short *) malloc(myPrimalMesh->faceNumber * sizeof(connectivity_short));
-
-//  for(connectivity_int i=0;i<myPrimalMesh->faceNumber;i++)
-//    {
-
-
-//      myPrimalMesh->faceToVertexNumber[i] = QUAD;
-
-//      debug_print( "%-40s : %ld (%d)\n",
-//                   "myPrimalMesh.faceToVertexNumber",
-//                   i,
-//                   myPrimalMesh->faceToVertexNumber[i]
-//                   );
-
-
-//    }
-//  endFunction(__FUNCTION__, t);
-//}
-
-//void setHexahedreCellToFacesNumber(struct primalMesh * myPrimalMesh)
-//{
-
-//  clock_t t=startFunction(__FUNCTION__);
-
-//  myPrimalMesh->cellToFacesNumber = (connectivity_short*) malloc(myPrimalMesh->cellNumber * sizeof(connectivity_short));
-
-
-//  for (connectivity_int i=0;i<myPrimalMesh->cellNumber;i++)
-//    {
-//      if(myPrimalMesh->cellToVertexNumber[i] == HEXAHEDRON)
-//        myPrimalMesh->cellToFacesNumber[i] = HEXAHEDRON_FACES;
-
-//    }
-
-//  for(connectivity_int ii=0;ii<myPrimalMesh->cellNumber;ii++)
-//    {
-//      debug_print( "%-40s : %ld (%d) \n" ,"myPrimalMesh->cellToFacesNumber", ii, myPrimalMesh->cellToFacesNumber[ii]);
-//    }
-//  endFunction(__FUNCTION__, t);
-//}
-
 
 void setHexahedreCellToVertex(struct primalMesh * myPrimalMesh)
+/*
+ * void setHexahedreCellToVertex(struct primalMesh * myPrimalMesh)
+ * 
+ * This function builds a vector by assigning for each cell the 
+ * number of vertices. The order of the vertices is ordered with 
+ * a positive sense of membership.
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use : 
+ * 
+ * Allocate : myPrimalMesh->cellToVertex
+ * Produce : myPrimalMesh->cellToVertex
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -237,122 +243,25 @@ void setHexahedreCellToVertex(struct primalMesh * myPrimalMesh)
 }
 
 
-//void setHexahedreFaceToVertex(struct primalMesh * myPrimalMesh)
-//{
-
-//  clock_t t=startFunction(__FUNCTION__);
-//  affiche("\n");
-
-//  connectivity_int it_cell=0;
-
-//  // vertex number on faces
-//  connectivity_int cN1 =  myPrimalMesh->P*myPrimalMesh->N;
-
-//  connectivity_int it_surface=0;
-//  connectivity_int it_face=0;
-
-
-//  allocMemoryForSpecialMatrixConnectivity(&(myPrimalMesh->faceToVertex),myPrimalMesh->faceNumber, myPrimalMesh->faceToVertexNumber);
-
-
-
-//  for(connectivity_int i=0;i<myPrimalMesh->M;i++)
-//    for(connectivity_int j=0;j<myPrimalMesh->N;j++)
-//      for(connectivity_int k=0;k<myPrimalMesh->P;k++)
-//        {
-//          //it_surface = k + j*(myPrimalMesh->P+1) + i*(myPrimalMesh->P+1)*(myPrimalMesh->N+1);
-//          it_cell     = k + j*myPrimalMesh->P + i*cN1; //cN1 = myPrimalMesh->P*myPrimalMesh->N;
-//          it_face     = 0;
-
-//          if(myPrimalMesh->cellToVertexNumber[it_cell] == HEXAHEDRON)
-//            {
-//              connectivity_int vertexList[HEXAHEDRON];
-//              // FACE 1
-//              myPrimalMesh->faceToVertex[it_surface][0] = myPrimalMesh->cellToVertex[it_cell][0];
-//              myPrimalMesh->faceToVertex[it_surface][1] = myPrimalMesh->cellToVertex[it_cell][3];
-//              myPrimalMesh->faceToVertex[it_surface][2] = myPrimalMesh->cellToVertex[it_cell][2];
-//              myPrimalMesh->faceToVertex[it_surface][3] = myPrimalMesh->cellToVertex[it_cell][1];
-
-//              //                    vertexList[0] = myPrimalMesh->cellToVertex[it_cell][0];
-//              //                    vertexList[1] = myPrimalMesh->cellToVertex[it_cell][3];
-//              //                    vertexList[2] = myPrimalMesh->cellToVertex[it_cell][2];
-//              //                    vertexList[3] = myPrimalMesh->cellToVertex[it_cell][1];
-
-//              //                    // FAIRE TRIER LISTE
-//              //                    // INVERT LISTE et TRIER LISTE
-
-//              //                    myPrimalMesh->faceToVertex[it_surface][0] = vertexList[0];
-//              //                    myPrimalMesh->faceToVertex[it_surface][1] = vertexList[3];
-//              //                    myPrimalMesh->faceToVertex[it_surface][2] = vertexList[2];
-//              //                    myPrimalMesh->faceToVertex[it_surface][3] = vertexList[1];
-
-//              it_surface+=1;
-
-//              // FACE 2
-//              myPrimalMesh->faceToVertex[it_surface][0] = myPrimalMesh->cellToVertex[it_cell][0];
-//              myPrimalMesh->faceToVertex[it_surface][1] = myPrimalMesh->cellToVertex[it_cell][1];
-//              myPrimalMesh->faceToVertex[it_surface][2] = myPrimalMesh->cellToVertex[it_cell][5];
-//              myPrimalMesh->faceToVertex[it_surface][3] = myPrimalMesh->cellToVertex[it_cell][4];
-
-//              it_surface+=1;
-
-
-//              // FACE 3
-//              myPrimalMesh->faceToVertex[it_surface][0] = myPrimalMesh->cellToVertex[it_cell][0];
-//              myPrimalMesh->faceToVertex[it_surface][1] = myPrimalMesh->cellToVertex[it_cell][4];
-//              myPrimalMesh->faceToVertex[it_surface][2] = myPrimalMesh->cellToVertex[it_cell][7];
-//              myPrimalMesh->faceToVertex[it_surface][3] = myPrimalMesh->cellToVertex[it_cell][3];
-
-//              it_surface+=1;
-
-
-//              // FACE 4
-//              myPrimalMesh->faceToVertex[it_surface][0] = myPrimalMesh->cellToVertex[it_cell][1];
-//              myPrimalMesh->faceToVertex[it_surface][1] = myPrimalMesh->cellToVertex[it_cell][2];
-//              myPrimalMesh->faceToVertex[it_surface][2] = myPrimalMesh->cellToVertex[it_cell][6];
-//              myPrimalMesh->faceToVertex[it_surface][3] = myPrimalMesh->cellToVertex[it_cell][5];
-
-//              it_surface+=1;
-
-
-//              // FACE 5
-//              myPrimalMesh->faceToVertex[it_surface][0] = myPrimalMesh->cellToVertex[it_cell][4];
-//              myPrimalMesh->faceToVertex[it_surface][1] = myPrimalMesh->cellToVertex[it_cell][5];
-//              myPrimalMesh->faceToVertex[it_surface][2] = myPrimalMesh->cellToVertex[it_cell][6];
-//              myPrimalMesh->faceToVertex[it_surface][3] = myPrimalMesh->cellToVertex[it_cell][7];
-
-//              it_surface+=1;
-
-
-//              // FACE 6
-//              myPrimalMesh->faceToVertex[it_surface][0] = myPrimalMesh->cellToVertex[it_cell][3];
-//              myPrimalMesh->faceToVertex[it_surface][1] = myPrimalMesh->cellToVertex[it_cell][7];
-//              myPrimalMesh->faceToVertex[it_surface][2] = myPrimalMesh->cellToVertex[it_cell][6];
-//              myPrimalMesh->faceToVertex[it_surface][3] = myPrimalMesh->cellToVertex[it_cell][2];
-
-//              it_surface+=1;
-
-//            }
-
-//        }
-
-//  for(connectivity_int k=0;k<myPrimalMesh->faceNumber;k++)
-//    {
-//      debug_print( "%-40s : %ld ( ","myPrimalMesh->faceToVertex", k);
-
-//      for(connectivity_int ii=0;ii<myPrimalMesh->faceToVertexNumber[k];ii++)
-//        {
-//          debug_print( "%ld " , myPrimalMesh->faceToVertex[k][ii]);
-//        }
-//      debug_print(")\n");
-//    }
-
-
-//  endFunction(__FUNCTION__, t);
-
-//}
 
 void setHexahedreVertexToCellNumbers(struct primalMesh * myPrimalMesh)
+/*
+ * void setHexahedreVertexToCellNumbers(struct primalMesh * myPrimalMesh)
+ * 
+ * This function constructs a vector by assigning for each vertex 
+ * the number of cells attached to it. 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  myPrimalMesh->cellNumber
+ *        myPrimalMesh->vertexNumber 
+ *        myPrimalMesh->cellToVertex
+ *        myPrimalMesh->cellToVertexNumber
+ * 
+ * Allocate : myPrimalMesh->vertexToCellNumber
+ * Produce : myPrimalMesh->vertexToCellNumber
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -389,6 +298,21 @@ void setHexahedreVertexToCellNumbers(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreVertexToCells(struct primalMesh * myPrimalMesh)
+/*
+ * void setHexahedreVertexToCells(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  myPrimalMesh->vertexNumber
+ *        myPrimalMesh->vertexToCellNumber
+ *        myPrimalMesh->cellToVertex
+ *        
+ * 
+ * Allocate : myPrimalMesh->vertexToCells
+ * Produce : myPrimalMesh->vertexToCells
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -441,6 +365,28 @@ void setHexahedreVertexToCells(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreCellToCells(struct primalMesh * myPrimalMesh)
+/*
+ * void setHexahedreCellToCells(struct primalMesh * myPrimalMesh)
+ * 
+ * This function searches the cells around a cell using vertices 
+ * to find them.
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  myPrimalMesh->M
+ *        myPrimalMesh->N
+ *        myPrimalMesh->P
+ *        myPrimalMesh->cellNumber
+ *        myPrimalMesh->cellToVertexNumber
+ *        myPrimalMesh->cellToVertex
+ *        
+ * 
+ * Allocate : myPrimalMesh->cellToCells
+ *            myPrimalMesh->cellToCellsNumbers
+ * Produce : myPrimalMesh->cellToCells
+ *           myPrimalMesh->cellToCellsNumbers
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -451,10 +397,6 @@ void setHexahedreCellToCells(struct primalMesh * myPrimalMesh)
   // vertex number on faces
   connectivity_int cN1 =  myPrimalMesh->P*myPrimalMesh->N;
   
-  //  connectivity_int it_face=0;
-  //    connectivity_int it_segment=0;
-  
-  //  allocMemoryForMatrixConnectivity(&(myPrimalMesh->segmentToVertex),myPrimalMesh->segmentNumber,SEGMENTVERTEX);
   connectivity_int * vertexList= (connectivity_int *) malloc(SEGMENTVERTEX * sizeof(connectivity_int));
   
   myPrimalMesh->cellToCells = (connectivity_int **) malloc( (myPrimalMesh->cellNumber) * sizeof(connectivity_int *));
@@ -466,9 +408,8 @@ void setHexahedreCellToCells(struct primalMesh * myPrimalMesh)
     for(connectivity_int j=0;j<myPrimalMesh->N;j++)
       for(connectivity_int k=0;k<myPrimalMesh->P;k++)
         {
-          //it_surface = k + j*(myPrimalMesh->P+1) + i*(myPrimalMesh->P+1)*(myPrimalMesh->N+1);
+          
           it_cell     = k + j*myPrimalMesh->P + i*cN1; //cN1 = myPrimalMesh->P*myPrimalMesh->N;
-          //          it_face     = 0;
           
           if(myPrimalMesh->cellToVertexNumber[it_cell] == HEXAHEDRON)
             {
@@ -495,19 +436,7 @@ void setHexahedreCellToCells(struct primalMesh * myPrimalMesh)
                           
                         }
                     }
-                  
-                  
-                  /*
-                        debug_print( "%-40s : Celli (%ld)   -> VERTEX (%7ld) -> CELL LIST ( ","cellList", it_cell, vertexSelect);
-                        
-                        for(connectivity_int ii=0;ii<cellNumbers;ii++)
-                        {
-                            debug_print("%ld ",cellList[ii]);
-                        }
-                        
-                        debug_print(")\n");
-*/
-                  
+           
                 }                        //                        cellNumbers+=jj;
               
               sortListConnectivity(cellList, cellNumbers);
@@ -537,18 +466,7 @@ void setHexahedreCellToCells(struct primalMesh * myPrimalMesh)
     }
 #endif
   
-  //  for(connectivity_int k=0;k<myPrimalMesh->segmentNumber;k++)
-  //    {
-  //      debug_print( "%-40s : %ld ( ","myPrimalMesh->segmentToVertex", k);
-  
-  //      for(connectivity_int ii=0;ii<SEGMENTVERTEX;ii++)
-  //        {
-  //          debug_print( "%ld " , myPrimalMesh->segmentToVertex[k][ii]);
-  //        }
-  
-  //      debug_print(")\n");
-  //    }
-  
+
   free(cellList);
   free(vertexList);
   
@@ -560,6 +478,33 @@ void setHexahedreCellToCells(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreCellToFacesOwnerNeighbour(struct primalMesh * myPrimalMesh)
+/*
+ * void setHexahedreCellToFacesOwnerNeighbour(struct primalMesh * myPrimalMesh)
+ * 
+ * XXX
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  myPrimalMesh->cellNumber
+ *        
+ * 
+ * Allocate : myPrimalMesh->cellToFacesOwner
+ *            myPrimalMesh->cellToFacesNeighbour
+ *            myPrimalMesh->cellToFacesOwnerNumber
+ *            myPrimalMesh->cellToFacesNeighbourNumber
+ *            
+ *            
+ *            
+ * Produce : myPrimalMesh->faces 
+ *           myPrimalMesh->faceNumber
+ *           myPrimalMesh->cellToFacesOwner
+ *           myPrimalMesh->cellToFacesOwnerNumber
+ *           myPrimalMesh->cellToFacesNeighbour
+ *           myPrimalMesh->cellToFacesNeighbourNumber
+ *          
+ *           
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -773,6 +718,19 @@ void setHexahedreCellToFacesOwnerNeighbour(struct primalMesh * myPrimalMesh)
 }
 
 void setHexahedreFaceToCells(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -853,6 +811,19 @@ void setHexahedreFaceToCells(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreSegments(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1043,6 +1014,19 @@ void setHexahedreSegments(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreSegmentsCentres(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1099,6 +1083,19 @@ void setHexahedreSegmentsCentres(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreSegments2(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1205,6 +1202,19 @@ void setHexahedreSegments2(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreVertexToSegments(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1296,6 +1306,19 @@ void setHexahedreVertexToSegments(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreSegmentToFaces(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1529,6 +1552,19 @@ void setHexahedreSegmentToFaces(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreSegmentToFaces2(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1658,6 +1694,19 @@ void setHexahedreSegmentToFaces2(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreSortFaceToVertex(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   //////////////////// TRIE LES NOEUDS
   
@@ -1686,29 +1735,7 @@ void setHexahedreSortFaceToVertex(struct primalMesh * myPrimalMesh)
       
     }
   
-  
-  
-  
-  //        for(conctivity_int k=0;k<myPrimalMesh->faceNumber;k++)
-  //        {
-  //            conctivity_int it_surface=0;
-  //            int             bool=0;
-  //            do
-  //            {
-  //                for(conctivity_int ii=0;ii<myPrimalMesh->faceToVertexNumber[k];ii++)
-  //                {
-  //                    //     ((ii + 1) % myPrimalMesh->faceToVertexNumber[k]) ;
-  
-  //                    if(myPrimalMesh->faceToVertex[k][ii] == myPrimalMesh->faceToVertex[it_surface][0])
-  //                    {
-  //                        bool = 1;
-  //                        debug_print("TROUVE : (%ld) myPrimalMesh->faceToVertex[%ld][%ld] = %ld\n", it_surface,k,ii,myPrimalMesh->faceToVertex[k][ii])
-  //                    }
-  //                }
-  //                it_surface++;
-  //            }
-  //            while( (it_surface<myPrimalMesh->faceNumber) && (bool == 0) );
-  //        }
   endFunction(__FUNCTION__, t);
   
 }
@@ -1745,6 +1772,19 @@ void setHexahedreAllocVertexToSegment(struct primalMesh * myPrimalMesh)
 }
 
 void setHexahedreFaceCentersAreas(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1898,6 +1938,19 @@ void setHexahedreFaceCentersAreas(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreEstimateVolumeCentroid(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   
   clock_t t=startFunction(__FUNCTION__);
@@ -1972,6 +2025,19 @@ void setHexahedreEstimateVolumeCentroid(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreVolumeCentroid(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   clock_t t=startFunction(__FUNCTION__);
   
@@ -2082,6 +2148,19 @@ void setHexahedreVolumeCentroid(struct primalMesh * myPrimalMesh)
 
 
 void setHexahedreBoundaryFaces(struct primalMesh * myPrimalMesh)
+/*
+ * void XXX(struct primalMesh * myPrimalMesh)
+ * 
+ * 
+ * Input : myPrimalMesh struct
+ * 
+ * Use :  
+ *        
+ * 
+ * Allocate : 
+ * Produce : 
+ * 
+ */
 {
   clock_t t=startFunction(__FUNCTION__);