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

authorAlain <alain.bastide@univ-reunion.fr>
Thu, 2 Jul 2020 01:54:29 +0000 (05:54 +0400)
committerAlain <alain.bastide@univ-reunion.fr>
Sat, 4 Jul 2020 02:55:14 +0000 (06:55 +0400)
.gitignore
src/main.c
src/primal.h

index 1ef9424..780d182 100644 (file)
@@ -44,3 +44,4 @@ Makefile
 *.cmake
 test/CMakeFiles
 CMakeFiles
+*.autosave
index ddaabde..81bc536 100644 (file)
@@ -41,7 +41,7 @@ int main (void)
 #else
   myPrimalMesh.M = 200; // cell number on X
   myPrimalMesh.N = 200; // cell number on Y
-  myPrimalMesh.P = 50; // cell number on Z
+  myPrimalMesh.P = 200; // cell number on Z
   
 #endif
   
@@ -149,6 +149,19 @@ int main (void)
   
   setDualFacesCentres(&myPrimalMesh, &myDualMesh);
   
+  /* To do dual segments to generate global dual surface
+   * Faces = myPrimalMesh->segmentToFaceOwner union myPrimalMesh->segmentToFaceNeighbour
+   * myPrimalMesh->faceToCells
+   * Faces2 = myPrimalMesh->cellToFacesOwner union myPrimalMesh->cellToFacesNeighbour
+   * Faces3 = Faces intersect Face2
+   * 
+   * TODO : 
+   * Dual faces on primal faces and domain boundary 
+   * 
+   * Volumes of dual mesh
+   * 
+   * test Volume dual = volume primal
+  */
   
   //  FILE *fp;
   //  fp = fopen("one.dat", "wb");
index d472e1a..0a24252 100644 (file)
@@ -4,7 +4,7 @@
 #include "mmd.h"
 
 #define DEF_hexaedron_localNodeList \
-connectivity_int  hexaedron_localNodeList[HEXAHEDRON_FACES][QUAD]={ \
+static const connectivity_int  hexaedron_localNodeList[HEXAHEDRON_FACES][QUAD]={ \
 {0,3,2,1}, \
 {0,1,5,4}, \
 {0,4,7,3}, \
@@ -15,7 +15,7 @@ connectivity_int  hexaedron_localNodeList[HEXAHEDRON_FACES][QUAD]={ \
 
 
 #define DEF_hexaedron_localNodeListNumbers \
-connectivity_int  hexaedron_localNodeListNumbers[HEXAHEDRON_FACES]={ \
+static const connectivity_int  hexaedron_localNodeListNumbers[HEXAHEDRON_FACES]={ \
   QUAD, \
   QUAD, \
   QUAD, \
@@ -26,7 +26,7 @@ connectivity_int  hexaedron_localNodeListNumbers[HEXAHEDRON_FACES]={ \
 
 
 #define DEF_hexaedron_localSegmentList \
-connectivity_int  hexaedron_localSegmentList[HEXAHEDRON_SEGMENTS][SEGMENTVERTEX]={ \
+static const connectivity_int  hexaedron_localSegmentList[HEXAHEDRON_SEGMENTS][SEGMENTVERTEX]={ \
 {0,1}, \
 {1,5}, \
 {5,4}, \
@@ -44,7 +44,7 @@ connectivity_int  hexaedron_localSegmentList[HEXAHEDRON_SEGMENTS][SEGMENTVERTEX]
 
 
 #define DEF_hexaedron_localSegmentsListNumbers \
-connectivity_int  hexaedron_localSegmentsListNumbers[HEXAHEDRON_SEGMENTS]={ \
+static const connectivity_int  hexaedron_localSegmentsListNumbers[HEXAHEDRON_SEGMENTS]={ \
   SEGMENTVERTEX, \
   SEGMENTVERTEX, \
   SEGMENTVERTEX, \