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

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");