1 static char help[] = "Tests DMPlex Gmsh reader.\n\n"; 2 3 #include <petscdmplex.h> 4 5 #if !defined(PETSC_GMSH_EXE) 6 #define PETSC_GMSH_EXE "gmsh" 7 #endif 8 9 int main(int argc, char **argv) 10 { 11 DM dm; 12 const char *const mshlist[] = {"seg", "tri", "qua", "tet", "wed", "hex"}; 13 const char *const fmtlist[] = {"msh22", "msh40", "msh41"}; 14 PetscInt msh = 5; 15 PetscInt fmt = 2; 16 PetscBool bin = PETSC_FALSE; 17 PetscInt order = 2; 18 19 const char cmdtemplate[] = "%s -3 -format %s %s -order %d %s -o %s"; 20 char gmsh[PETSC_MAX_PATH_LEN] = PETSC_GMSH_EXE; 21 char tag[PETSC_MAX_PATH_LEN], path[PETSC_MAX_PATH_LEN]; 22 char geo[PETSC_MAX_PATH_LEN], geodir[PETSC_MAX_PATH_LEN] = "."; 23 char out[PETSC_MAX_PATH_LEN], outdir[PETSC_MAX_PATH_LEN] = "."; 24 char cmd[PETSC_MAX_PATH_LEN*4]; 25 PetscBool set,flg; 26 FILE *fp; 27 PetscErrorCode ierr; 28 29 ierr = PetscInitialize(&argc, &argv, NULL, help);if (ierr) return ierr; 30 31 ierr = PetscStrncpy(geodir, "${PETSC_DIR}/share/petsc/datafiles/meshes", sizeof(geodir));CHKERRQ(ierr); 32 ierr = PetscOptionsGetenv(PETSC_COMM_SELF, "GMSH", path, sizeof(path), &set);CHKERRQ(ierr); 33 if (set) {ierr = PetscStrncpy(gmsh, path, sizeof(gmsh));CHKERRQ(ierr);} 34 ierr = PetscOptionsGetString(NULL, NULL, "-gmsh", gmsh, sizeof(gmsh), NULL);CHKERRQ(ierr); 35 ierr = PetscOptionsGetString(NULL, NULL, "-dir", geodir, sizeof(geodir), NULL);CHKERRQ(ierr); 36 ierr = PetscOptionsGetString(NULL, NULL, "-out", outdir, sizeof(outdir), NULL);CHKERRQ(ierr); 37 ierr = PetscOptionsGetEList(NULL, NULL, "-msh", mshlist, 6, &msh, NULL);CHKERRQ(ierr); 38 ierr = PetscOptionsGetEList(NULL, NULL, "-fmt", fmtlist, 3, &fmt, NULL);CHKERRQ(ierr); 39 ierr = PetscOptionsGetBool(NULL, NULL, "-bin", &bin, NULL);CHKERRQ(ierr); 40 ierr = PetscOptionsGetInt(NULL, NULL, "-order", &order, NULL);CHKERRQ(ierr); 41 if (fmt == 1) bin = PETSC_FALSE; /* Recent Gmsh releases cannot generate msh40+binary format*/ 42 43 { /* This test requires Gmsh >= 4.2.0 */ 44 int inum = 0, major = 0, minor = 0, micro = 0; 45 ierr = PetscSNPrintf(cmd, sizeof(cmd), "%s -info", gmsh);CHKERRQ(ierr); 46 ierr = PetscPOpen(PETSC_COMM_SELF, NULL, cmd, "r", &fp);CHKERRQ(ierr); 47 if (fp) {inum = fscanf(fp, "Version : %d.%d.%d", &major, &minor, µ);} 48 ierr = PetscPClose(PETSC_COMM_SELF, fp);CHKERRQ(ierr); 49 if (inum != 3 || major < 4 || (major == 4 && minor < 2)) { 50 ierr = PetscPrintf(PETSC_COMM_SELF, "Gmsh>=4.2.0 not available\n");CHKERRQ(ierr); goto finish; 51 } 52 } 53 54 ierr = PetscSNPrintf(tag, sizeof(tag), "%s-%d-%s%s", mshlist[msh], order, fmtlist[fmt], bin?"-bin":"");CHKERRQ(ierr); 55 ierr = PetscSNPrintf(geo, sizeof(geo), "%s/gmsh-%s.geo", geodir, mshlist[msh]);CHKERRQ(ierr); 56 ierr = PetscSNPrintf(out, sizeof(out), "%s/mesh-%s.msh", outdir, tag);CHKERRQ(ierr); 57 ierr = PetscStrreplace(PETSC_COMM_SELF, geo, path, sizeof(path));CHKERRQ(ierr); 58 ierr = PetscFixFilename(path, geo);CHKERRQ(ierr); 59 ierr = PetscStrreplace(PETSC_COMM_SELF, out, path, sizeof(path));CHKERRQ(ierr); 60 ierr = PetscFixFilename(path, out);CHKERRQ(ierr); 61 ierr = PetscTestFile(geo, 'r', &flg);CHKERRQ(ierr); 62 if (!flg) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_USER_INPUT, "File not found: %s", geo); 63 64 ierr = PetscSNPrintf(cmd, sizeof(cmd), cmdtemplate, gmsh, fmtlist[fmt], bin?"-bin":"", order, geo, out);CHKERRQ(ierr); 65 ierr = PetscPOpen(PETSC_COMM_SELF, NULL, cmd, "r", &fp);CHKERRQ(ierr); 66 ierr = PetscPClose(PETSC_COMM_SELF, fp);CHKERRQ(ierr); 67 68 ierr = DMPlexCreateFromFile(PETSC_COMM_SELF, out, PETSC_TRUE, &dm);CHKERRQ(ierr); 69 ierr = PetscSNPrintf(tag, sizeof(tag), "mesh-%s", mshlist[msh]);CHKERRQ(ierr); 70 ierr = PetscObjectSetName((PetscObject)dm, tag);CHKERRQ(ierr); 71 ierr = DMViewFromOptions(dm, NULL, "-dm_view");CHKERRQ(ierr); 72 ierr = DMSetFromOptions(dm);CHKERRQ(ierr); 73 ierr = DMDestroy(&dm);CHKERRQ(ierr); 74 75 finish: 76 ierr = PetscFinalize(); 77 return ierr; 78 } 79 80 /*TEST 81 82 build: 83 requires: define(PETSC_HAVE_POPEN) 84 85 test: 86 requires: define(PETSC_GMSH_EXE) 87 args: -dir ${wPETSC_DIR}/share/petsc/datafiles/meshes 88 args: -msh {{seg tri qua tet wed hex}separate_output} 89 args: -order {{1 2 3 7}} 90 args: -fmt {{msh22 msh40 msh41}} -bin {{0 1}} 91 args: -dm_view ::ascii_info_detail 92 args: -dm_plex_check_all 93 94 TEST*/ 95