/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c

Bug Summary

File:	dm/impls/plex/plexgmsh.c
Warning:	line 853, column 23 The right operand of '<' is a garbage value

Annotated Source Code

[?] Use j/k keys for keyboard navigation

#define PETSCDM_DLL

#include <petsc/private/dmpleximpl.h> /*I "petscdmplex.h" I*/

#include <petsc/private/hashmapi.h>

typedef struct {

PetscViewer viewer;

int fileFormat;

int dataSize;

PetscBool binary;

PetscBool byteSwap;

size_t wlen;

void *wbuf;

size_t slen;

void *sbuf;

} GmshFile;

static PetscErrorCode GmshBufferGet(GmshFile *gmsh, size_t count, size_t eltsize, void *buf)

{

size_t size = count * eltsize;

PetscErrorCode ierr;

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 22; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

if (gmsh->wlen < size) {

ierr = PetscFree(gmsh->wbuf)((*PetscTrFree)((void*)(gmsh->wbuf),24,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((gmsh->wbuf) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),24,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

ierr = PetscMalloc(size, &gmsh->wbuf)((*PetscTrMalloc)((size),PETSC_FALSE,25,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(void**)(&gmsh->wbuf)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),25,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

gmsh->wlen = size;

}

*(void**)buf = size ? gmsh->wbuf : NULL((void*)0);

PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);

}

static PetscErrorCode GmshBufferSizeGet(GmshFile *gmsh, size_t count, void *buf)

{

size_t dataSize = (size_t)gmsh->dataSize;

size_t size = count * dataSize;

PetscErrorCode ierr;

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 38; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

if (gmsh->slen < size) {

ierr = PetscFree(gmsh->sbuf)((*PetscTrFree)((void*)(gmsh->sbuf),40,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((gmsh->sbuf) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),40,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

ierr = PetscMalloc(size, &gmsh->sbuf)((*PetscTrMalloc)((size),PETSC_FALSE,41,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(void**)(&gmsh->sbuf)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),41,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

gmsh->slen = size;

}

*(void**)buf = size ? gmsh->sbuf : NULL((void*)0);

}

static PetscErrorCode GmshRead(GmshFile *gmsh, void *buf, PetscInt count, PetscDataType dtype)

{

PetscErrorCode ierr;

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 51; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

ierr = PetscViewerRead(gmsh->viewer, buf, count, NULL((void*)0), dtype);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),52,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

if (gmsh->byteSwap) {ierr = PetscByteSwap(buf, dtype, count);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),53,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

}

static PetscErrorCode GmshReadString(GmshFile *gmsh, char *buf, PetscInt count)

{

PetscErrorCode ierr;

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 60; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

ierr = PetscViewerRead(gmsh->viewer, buf, count, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),61,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

}

static PetscErrorCode GmshMatch(PETSC_UNUSED__attribute((unused)) GmshFile *gmsh, const char Section[], char line[PETSC_MAX_PATH_LEN4096], PetscBool *match)

{

PetscErrorCode ierr;

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 68; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

ierr = PetscStrcmp(line, Section, match);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),69,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

}

static PetscErrorCode GmshExpect(GmshFile *gmsh, const char Section[], char line[PETSC_MAX_PATH_LEN4096])

{

PetscBool match;

PetscErrorCode ierr;

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 78; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

ierr = GmshMatch(gmsh, Section, line, &match);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),79,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

if (!match) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file, expecting %s",Section)return PetscError(((MPI_Comm)0x44000001),80,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file, expecting %s"
,Section);

}

static PetscErrorCode GmshReadSection(GmshFile *gmsh, char line[PETSC_MAX_PATH_LEN4096])

{

PetscBool match;

PetscErrorCode ierr;

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 89; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

while (PETSC_TRUE) {

ierr = GmshReadString(gmsh, line, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),91,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

ierr = GmshMatch(gmsh, "$Comments", line, &match);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),92,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

if (!match) break;

while (PETSC_TRUE) {

ierr = GmshReadString(gmsh, line, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),95,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

ierr = GmshMatch(gmsh, "$EndComments", line, &match);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),96,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

if (match) break;

}

100

101

}

102

103

static PetscErrorCode GmshReadEndSection(GmshFile *gmsh, const char EndSection[], char line[PETSC_MAX_PATH_LEN4096])

104

{

105

PetscErrorCode ierr;

106

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 106; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

107

ierr = GmshReadString(gmsh, line, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),107,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

108

ierr = GmshExpect(gmsh, EndSection, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),108,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

109

110

}

111

112

static PetscErrorCode GmshReadSize(GmshFile *gmsh, PetscInt *buf, PetscInt count)

113

{

114

PetscInt i;

115

size_t dataSize = (size_t)gmsh->dataSize;

116

PetscErrorCode ierr;

117

118

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 118; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

←

Within the expansion of the macro 'PetscFunctionBegin':

→

119

if (dataSize == sizeof(PetscInt)) {

←

Assuming the condition is false

→

←

Taking false branch

→

120

ierr = GmshRead(gmsh, buf, count, PETSC_INT);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),120,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

121

} else if (dataSize == sizeof(int)) {

←

Taking false branch

→

122

int *ibuf = NULL((void*)0);

123

ierr = GmshBufferSizeGet(gmsh, count, &ibuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),123,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

124

ierr = GmshRead(gmsh, ibuf, count, PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),124,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

125

for (i = 0; i < count; ++i) buf[i] = (PetscInt)ibuf[i];

126

} else if (dataSize == sizeof(long)) {

←

Assuming the condition is false

→

←

Taking false branch

→

127

long *ibuf = NULL((void*)0);

128

ierr = GmshBufferSizeGet(gmsh, count, &ibuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),128,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

129

ierr = GmshRead(gmsh, ibuf, count, PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),129,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

130

for (i = 0; i < count; ++i) buf[i] = (PetscInt)ibuf[i];

131

} else if (dataSize == sizeof(PetscInt64)) {

←

Taking false branch

→

132

PetscInt64 *ibuf = NULL((void*)0);

133

ierr = GmshBufferSizeGet(gmsh, count, &ibuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),133,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

134

ierr = GmshRead(gmsh, ibuf, count, PETSC_INT64);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),134,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

135

for (i = 0; i < count; ++i) buf[i] = (PetscInt)ibuf[i];

136

}

137

←

Within the expansion of the macro 'PetscFunctionReturn':

→

138

}

139

140

static PetscErrorCode GmshReadInt(GmshFile *gmsh, int *buf, PetscInt count)

141

{

142

PetscErrorCode ierr;

143

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 143; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

144

ierr = GmshRead(gmsh, buf, count, PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),144,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

145

146

}

147

148

static PetscErrorCode GmshReadDouble(GmshFile *gmsh, double *buf, PetscInt count)

149

{

150

PetscErrorCode ierr;

151

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 151; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

152

ierr = GmshRead(gmsh, buf, count, PETSC_DOUBLE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),152,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

153

154

}

155

156

typedef struct {

157

PetscInt id; /* Entity number */

158

PetscInt dim; /* Entity dimension */

159

double bbox[6]; /* Bounding box */

160

PetscInt numTags; /* Size of tag array */

161

int tags[4]; /* Tag array */

162

} GmshEntity;

163

164

typedef struct {

165

GmshEntity *entity[4];

166

PetscHMapI entityMap[4];

167

} GmshEntities;

168

169

static PetscErrorCode GmshEntitiesCreate(PetscInt count[4], GmshEntities **entities)

170

{

171

PetscInt dim;

172

PetscErrorCode ierr;

173

174

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 174; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

175

ierr = PetscNew(entities)PetscMallocA(1,PETSC_TRUE,175,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(1)*sizeof(**((entities))),((entities)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),175,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

176

for (dim = 0; dim < 4; ++dim) {

177

ierr = PetscCalloc1(count[dim], &(*entities)->entity[dim])PetscMallocA(1,PETSC_TRUE,177,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(count[dim])*sizeof(**(&(*entities)->entity[dim
])),(&(*entities)->entity[dim]));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),177,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

178

ierr = PetscHMapICreate(&(*entities)->entityMap[dim]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),178,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

179

}

180

181

}

182

183

static PetscErrorCode GmshEntitiesAdd(GmshEntities *entities, PetscInt index, PetscInt dim, PetscInt eid, GmshEntity** entity)

184

{

185

PetscErrorCode ierr;

186

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 186; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

187

ierr = PetscHMapISet(entities->entityMap[dim], eid, index);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),187,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

188

entities->entity[dim][index].dim = dim;

189

entities->entity[dim][index].id = eid;

190

if (entity) *entity = &entities->entity[dim][index];

191

192

}

193

194

static PetscErrorCode GmshEntitiesGet(GmshEntities *entities, PetscInt dim, PetscInt eid, GmshEntity** entity)

195

{

196

PetscInt index;

197

PetscErrorCode ierr;

198

199

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 199; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

200

ierr = PetscHMapIGet(entities->entityMap[dim], eid, &index);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),200,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

201

*entity = &entities->entity[dim][index];

202

203

}

204

205

static PetscErrorCode GmshEntitiesDestroy(GmshEntities **entities)

206

{

207

PetscInt dim;

208

PetscErrorCode ierr;

209

210

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 210; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

211

if (!*entities) PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);

212

for (dim = 0; dim < 4; ++dim) {

213

ierr = PetscFree((*entities)->entity[dim])((*PetscTrFree)((void*)((*entities)->entity[dim]),213,__func__
,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c") ||
(((*entities)->entity[dim]) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),213,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

214

ierr = PetscHMapIDestroy(&(*entities)->entityMap[dim]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),214,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

215

}

216

ierr = PetscFree((*entities))((*PetscTrFree)((void*)((*entities)),216,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || (((*entities)) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),216,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

217

218

}

219

220

typedef struct {

221

PetscInt id; /* Entity number */

222

PetscInt dim; /* Entity dimension */

223

PetscInt cellType; /* Cell type */

224

PetscInt numNodes; /* Size of node array */

225

PetscInt nodes[8]; /* Node array */

226

PetscInt numTags; /* Size of tag array */

227

int tags[4]; /* Tag array */

228

} GmshElement;

229

230

static PetscErrorCode DMPlexCreateGmsh_ReadNodes_v22(GmshFile *gmsh, int shift, PetscInt *numVertices, double **gmsh_nodes)

231

{

232

PetscViewer viewer = gmsh->viewer;

233

PetscBool byteSwap = gmsh->byteSwap;

234

char line[PETSC_MAX_PATH_LEN4096];

235

int v, num, nid, snum;

236

double *coordinates;

237

PetscErrorCode ierr;

238

239

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 239; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

240

ierr = PetscViewerRead(viewer, line, 1, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),240,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

241

snum = sscanf(line, "%d", &num);

242

if (snum != 1) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),242,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

243

ierr = PetscMalloc1(num*3, &coordinates)PetscMallocA(1,PETSC_FALSE,243,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(num*3)*sizeof(**(&coordinates)),(&coordinates
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),243,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

244

*numVertices = num;

245

*gmsh_nodes = coordinates;

246

for (v = 0; v < num; ++v) {

247

double *xyz = coordinates + v*3;

248

ierr = PetscViewerRead(viewer, &nid, 1, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),248,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

249

if (byteSwap) {ierr = PetscByteSwap(&nid, PETSC_ENUM, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),249,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

250

if (nid != v+shift) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unexpected node number %d should be %d", nid, v+shift)return PetscError(((MPI_Comm)0x44000001),250,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unexpected node number %d should be %d"
,nid,v+shift);

251

ierr = PetscViewerRead(viewer, xyz, 3, NULL((void*)0), PETSC_DOUBLE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),251,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

252

if (byteSwap) {ierr = PetscByteSwap(xyz, PETSC_DOUBLE, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),252,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

253

}

254

255

}

256

257

/* Gmsh elements can be of any dimension/co-dimension, so we need to traverse the

258

file contents multiple times to figure out the true number of cells and facets

259

in the given mesh. To make this more efficient we read the file contents only

260

once and store them in memory, while determining the true number of cells. */

261

static PetscErrorCode DMPlexCreateGmsh_ReadElements_v22(GmshFile* gmsh, PETSC_UNUSED__attribute((unused)) int shift, PetscInt *numCells, GmshElement **gmsh_elems)

262

{

263

PetscViewer viewer = gmsh->viewer;

264

PetscBool binary = gmsh->binary;

265

PetscBool byteSwap = gmsh->byteSwap;

266

char line[PETSC_MAX_PATH_LEN4096];

267

GmshElement *elements;

268

int i, c, p, num, ibuf[1+4+512], snum;

269

int cellType, dim, numNodes, numNodesIgnore, numElem, numTags;

270

PetscErrorCode ierr;

271

272

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 272; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

273

ierr = PetscViewerRead(viewer, line, 1, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),273,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);;

274

snum = sscanf(line, "%d", &num);

275

if (snum != 1) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),275,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

276

ierr = PetscMalloc1(num, &elements)PetscMallocA(1,PETSC_FALSE,276,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(num)*sizeof(**(&elements)),(&elements));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),276,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

277

*numCells = num;

278

*gmsh_elems = elements;

279

for (c = 0; c < num;) {

280

ierr = PetscViewerRead(viewer, ibuf, 3, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),280,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

281

if (byteSwap) {ierr = PetscByteSwap(ibuf, PETSC_ENUM, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),281,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

282

if (binary) {

283

cellType = ibuf[0];

284

numElem = ibuf[1];

285

numTags = ibuf[2];

286

} else {

287

elements[c].id = ibuf[0];

288

cellType = ibuf[1];

289

numTags = ibuf[2];

290

numElem = 1;

291

}

292

/* http://gmsh.info/doc/texinfo/gmsh.html#MSH-ASCII-file-format */

293

numNodesIgnore = 0;

294

switch (cellType) {

295

case 1: /* 2-node line */

296

dim = 1;

297

numNodes = 2;

298

break;

299

case 2: /* 3-node triangle */

300

dim = 2;

301

numNodes = 3;

302

break;

303

case 3: /* 4-node quadrangle */

304

dim = 2;

305

numNodes = 4;

306

break;

307

case 4: /* 4-node tetrahedron */

308

dim = 3;

309

numNodes = 4;

310

break;

311

case 5: /* 8-node hexahedron */

312

dim = 3;

313

numNodes = 8;

314

break;

315

case 6: /* 6-node wedge */

316

dim = 3;

317

numNodes = 6;

318

break;

319

case 8: /* 3-node 2nd order line */

320

dim = 1;

321

numNodes = 2;

322

numNodesIgnore = 1;

323

break;

324

case 9: /* 6-node 2nd order triangle */

325

dim = 2;

326

numNodes = 3;

327

numNodesIgnore = 3;

328

break;

329

case 10: /* 9-node 2nd order quadrangle */

330

dim = 2;

331

numNodes = 4;

332

numNodesIgnore = 5;

333

break;

334

case 11: /* 10-node 2nd order tetrahedron */

335

dim = 3;

336

numNodes = 4;

337

numNodesIgnore = 6;

338

break;

339

case 12: /* 27-node 2nd order hexhedron */

340

dim = 3;

341

numNodes = 8;

342

numNodesIgnore = 19;

343

break;

344

case 13: /* 18-node 2nd wedge */

345

dim = 3;

346

numNodes = 6;

347

numNodesIgnore = 12;

348

break;

349

case 15: /* 1-node vertex */

350

dim = 0;

351

numNodes = 1;

352

break;

353

case 7: /* 5-node pyramid */

354

case 14: /* 14-node 2nd order pyramid */

355

default:

356

SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unsupported Gmsh element type %d", cellType)return PetscError(((MPI_Comm)0x44000001),356,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unsupported Gmsh element type %d",cellType
);

357

}

358

if (binary) {

359

const int nint = 1 + numTags + numNodes + numNodesIgnore;

360

/* Loop over element blocks */

361

for (i = 0; i < numElem; ++i, ++c) {

362

ierr = PetscViewerRead(viewer, ibuf, nint, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),362,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

363

if (byteSwap) {ierr = PetscByteSwap(ibuf, PETSC_ENUM, nint);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),363,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

364

elements[c].dim = dim;

365

elements[c].numNodes = numNodes;

366

elements[c].numTags = numTags;

367

elements[c].id = ibuf[0];

368

elements[c].cellType = cellType;

369

for (p = 0; p < numTags; p++) elements[c].tags[p] = ibuf[1 + p];

370

for (p = 0; p < numNodes; p++) elements[c].nodes[p] = ibuf[1 + numTags + p];

371

}

372

} else {

373

const int nint = numTags + numNodes + numNodesIgnore;

374

elements[c].dim = dim;

375

elements[c].numNodes = numNodes;

376

elements[c].numTags = numTags;

377

elements[c].cellType = cellType;

378

ierr = PetscViewerRead(viewer, ibuf, nint, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),378,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

379

for (p = 0; p < numTags; p++) elements[c].tags[p] = ibuf[p];

380

for (p = 0; p < numNodes; p++) elements[c].nodes[p] = ibuf[numTags + p];

381

c++;

382

}

383

}

384

385

}

386

387

388

$Entities

389

numPoints(unsigned long) numCurves(unsigned long)

390

numSurfaces(unsigned long) numVolumes(unsigned long)

391

// points

392

tag(int) boxMinX(double) boxMinY(double) boxMinZ(double)

393

boxMaxX(double) boxMaxY(double) boxMaxZ(double)

394

numPhysicals(unsigned long) phyisicalTag[...](int)

395

...

396

// curves

397

tag(int) boxMinX(double) boxMinY(double) boxMinZ(double)

398

boxMaxX(double) boxMaxY(double) boxMaxZ(double)

399

numPhysicals(unsigned long) physicalTag[...](int)

400

numBREPVert(unsigned long) tagBREPVert[...](int)

401

...

402

// surfaces

403

tag(int) boxMinX(double) boxMinY(double) boxMinZ(double)

404

boxMaxX(double) boxMaxY(double) boxMaxZ(double)

405

numPhysicals(unsigned long) physicalTag[...](int)

406

numBREPCurve(unsigned long) tagBREPCurve[...](int)

407

...

408

// volumes

409

tag(int) boxMinX(double) boxMinY(double) boxMinZ(double)

410

boxMaxX(double) boxMaxY(double) boxMaxZ(double)

411

numPhysicals(unsigned long) physicalTag[...](int)

412

numBREPSurfaces(unsigned long) tagBREPSurfaces[...](int)

413

...

414

$EndEntities

415

416

static PetscErrorCode DMPlexCreateGmsh_ReadEntities_v40(GmshFile *gmsh, GmshEntities **entities)

417

{

418

PetscViewer viewer = gmsh->viewer;

419

PetscBool byteSwap = gmsh->byteSwap;

420

long index, num, lbuf[4];

421

int dim, eid, numTags, *ibuf, t;

422

PetscInt count[4], i;

423

GmshEntity *entity = NULL((void*)0);

424

PetscErrorCode ierr;

425

426

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 426; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

427

ierr = PetscViewerRead(viewer, lbuf, 4, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),427,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

428

if (byteSwap) {ierr = PetscByteSwap(lbuf, PETSC_LONG, 4);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),428,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

429

for (i = 0; i < 4; ++i) count[i] = lbuf[i];

430

ierr = GmshEntitiesCreate(count, entities);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),430,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

431

for (dim = 0; dim < 4; ++dim) {

432

for (index = 0; index < count[dim]; ++index) {

433

ierr = PetscViewerRead(viewer, &eid, 1, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),433,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

434

if (byteSwap) {ierr = PetscByteSwap(&eid, PETSC_ENUM, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),434,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

435

ierr = GmshEntitiesAdd(*entities, (PetscInt)index, dim, eid, &entity);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),435,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

436

ierr = PetscViewerRead(viewer, entity->bbox, 6, NULL((void*)0), PETSC_DOUBLE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),436,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

437

if (byteSwap) {ierr = PetscByteSwap(entity->bbox, PETSC_DOUBLE, 6);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),437,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

438

ierr = PetscViewerRead(viewer, &num, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),438,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

439

if (byteSwap) {ierr = PetscByteSwap(&num, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),439,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

440

ierr = GmshBufferGet(gmsh, num, sizeof(int), &ibuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),440,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

441

ierr = PetscViewerRead(viewer, ibuf, num, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),441,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

442

if (byteSwap) {ierr = PetscByteSwap(ibuf, PETSC_ENUM, num);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),442,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

443

entity->numTags = numTags = (int) PetscMin(num, 4)(((num)<(4)) ? (num) : (4));

444

for (t = 0; t < numTags; ++t) entity->tags[t] = ibuf[t];

445

if (dim == 0) continue;

446

ierr = PetscViewerRead(viewer, &num, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),446,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

447

if (byteSwap) {ierr = PetscByteSwap(&num, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),447,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

448

ierr = GmshBufferGet(gmsh, num, sizeof(int), &ibuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),448,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

449

ierr = PetscViewerRead(viewer, ibuf, num, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),449,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

450

if (byteSwap) {ierr = PetscByteSwap(ibuf, PETSC_ENUM, num);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),450,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

451

}

452

}

453

454

}

455

456

457

$Nodes

458

numEntityBlocks(unsigned long) numNodes(unsigned long)

459

tagEntity(int) dimEntity(int) typeNode(int) numNodes(unsigned long)

460

tag(int) x(double) y(double) z(double)

461

...

462

...

463

$EndNodes

464

465

static PetscErrorCode DMPlexCreateGmsh_ReadNodes_v40(GmshFile *gmsh, int shift, PetscInt *numVertices, double **gmsh_nodes)

466

{

467

PetscViewer viewer = gmsh->viewer;

468

PetscBool byteSwap = gmsh->byteSwap;

469

long block, node, v, numEntityBlocks, numTotalNodes, numNodes;

470

int info[3], nid;

471

double *coordinates;

472

char *cbuf;

473

PetscErrorCode ierr;

474

475

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 475; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

476

ierr = PetscViewerRead(viewer, &numEntityBlocks, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),476,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

477

if (byteSwap) {ierr = PetscByteSwap(&numEntityBlocks, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),477,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

478

ierr = PetscViewerRead(viewer, &numTotalNodes, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),478,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

479

if (byteSwap) {ierr = PetscByteSwap(&numTotalNodes, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),479,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

480

ierr = PetscMalloc1(numTotalNodes*3, &coordinates)PetscMallocA(1,PETSC_FALSE,480,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numTotalNodes*3)*sizeof(**(&coordinates)),(&
coordinates));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),480,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

481

*numVertices = numTotalNodes;

482

*gmsh_nodes = coordinates;

483

for (v = 0, block = 0; block < numEntityBlocks; ++block) {

484

ierr = PetscViewerRead(viewer, info, 3, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),484,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

485

ierr = PetscViewerRead(viewer, &numNodes, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),485,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

486

if (byteSwap) {ierr = PetscByteSwap(&numNodes, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),486,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

487

if (gmsh->binary) {

488

int nbytes = sizeof(int) + 3*sizeof(double);

489

ierr = GmshBufferGet(gmsh, numNodes, nbytes, &cbuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),489,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

490

ierr = PetscViewerRead(viewer, cbuf, numNodes*nbytes, NULL((void*)0), PETSC_CHAR);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),490,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

491

for (node = 0; node < numNodes; ++node, ++v) {

492

char *cnid = cbuf + node*nbytes, *cxyz = cnid + sizeof(int);

493

double *xyz = coordinates + v*3;

494

if (!PetscBinaryBigEndian()) {ierr = PetscByteSwap(cnid, PETSC_ENUM, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),494,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

495

if (!PetscBinaryBigEndian()) {ierr = PetscByteSwap(cxyz, PETSC_DOUBLE, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),495,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

496

ierr = PetscMemcpy(&nid, cnid, sizeof(int));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),496,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

497

ierr = PetscMemcpy(xyz, cxyz, 3*sizeof(double));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),497,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

498

if (byteSwap) {ierr = PetscByteSwap(&nid, PETSC_ENUM, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),498,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

499

if (byteSwap) {ierr = PetscByteSwap(xyz, PETSC_DOUBLE, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),499,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

500

if (nid != v+shift) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unexpected node number %d should be %d", nid, v+shift)return PetscError(((MPI_Comm)0x44000001),500,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unexpected node number %d should be %d"
,nid,v+shift);

501

}

502

} else {

503

for (node = 0; node < numNodes; ++node, ++v) {

504

double *xyz = coordinates + v*3;

505

ierr = PetscViewerRead(viewer, &nid, 1, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),505,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

506

if (byteSwap) {ierr = PetscByteSwap(&nid, PETSC_ENUM, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),506,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

507

if (nid != v+shift) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unexpected node number %d should be %d", nid, v+shift)return PetscError(((MPI_Comm)0x44000001),507,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unexpected node number %d should be %d"
,nid,v+shift);

508

ierr = PetscViewerRead(viewer, xyz, 3, NULL((void*)0), PETSC_DOUBLE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),508,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

509

if (byteSwap) {ierr = PetscByteSwap(xyz, PETSC_DOUBLE, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),509,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

510

}

511

}

512

}

513

514

}

515

516

517

$Elements

518

numEntityBlocks(unsigned long) numElements(unsigned long)

519

tagEntity(int) dimEntity(int) typeEle(int) numElements(unsigned long)

520

tag(int) numVert[...](int)

521

...

522

...

523

$EndElements

524

525

static PetscErrorCode DMPlexCreateGmsh_ReadElements_v40(GmshFile *gmsh, PETSC_UNUSED__attribute((unused)) int shift, GmshEntities *entities, PetscInt *numCells, GmshElement **gmsh_elems)

526

{

527

PetscViewer viewer = gmsh->viewer;

528

PetscBool byteSwap = gmsh->byteSwap;

529

long c, block, numEntityBlocks, numTotalElements, elem, numElements;

530

int p, info[3], *ibuf = NULL((void*)0);

531

int eid, dim, numTags, *tags, cellType, numNodes;

532

GmshEntity *entity = NULL((void*)0);

533

GmshElement *elements;

534

PetscErrorCode ierr;

535

536

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 536; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

537

ierr = PetscViewerRead(viewer, &numEntityBlocks, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),537,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

538

if (byteSwap) {ierr = PetscByteSwap(&numEntityBlocks, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),538,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

539

ierr = PetscViewerRead(viewer, &numTotalElements, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),539,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

540

if (byteSwap) {ierr = PetscByteSwap(&numTotalElements, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),540,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

541

ierr = PetscCalloc1(numTotalElements, &elements)PetscMallocA(1,PETSC_TRUE,541,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numTotalElements)*sizeof(**(&elements)),(&elements
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),541,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

542

*numCells = numTotalElements;

543

*gmsh_elems = elements;

544

for (c = 0, block = 0; block < numEntityBlocks; ++block) {

545

ierr = PetscViewerRead(viewer, info, 3, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),545,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

546

if (byteSwap) {ierr = PetscByteSwap(info, PETSC_ENUM, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),546,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

547

eid = info[0]; dim = info[1]; cellType = info[2];

548

ierr = GmshEntitiesGet(entities, dim, eid, &entity);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),548,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

549

numTags = entity->numTags;

550

tags = entity->tags;

551

switch (cellType) {

552

case 1: /* 2-node line */

553

numNodes = 2;

554

break;

555

case 2: /* 3-node triangle */

556

numNodes = 3;

557

break;

558

case 3: /* 4-node quadrangle */

559

numNodes = 4;

560

break;

561

case 4: /* 4-node tetrahedron */

562

numNodes = 4;

563

break;

564

case 5: /* 8-node hexahedron */

565

numNodes = 8;

566

break;

567

case 6: /* 6-node wedge */

568

numNodes = 6;

569

break;

570

case 15: /* 1-node vertex */

571

numNodes = 1;

572

break;

573

default:

574

SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unsupported Gmsh element type %d", cellType)return PetscError(((MPI_Comm)0x44000001),574,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unsupported Gmsh element type %d",cellType
);

575

}

576

ierr = PetscViewerRead(viewer, &numElements, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),576,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

577

if (byteSwap) {ierr = PetscByteSwap(&numElements, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),577,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

578

ierr = GmshBufferGet(gmsh, (1+numNodes)*numElements, sizeof(int), &ibuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),578,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

579

ierr = PetscViewerRead(viewer, ibuf, (1+numNodes)*numElements, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),579,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

580

if (byteSwap) {ierr = PetscByteSwap(ibuf, PETSC_ENUM, (1+numNodes)*numElements);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),580,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

581

for (elem = 0; elem < numElements; ++elem, ++c) {

582

int *id = ibuf + elem*(1+numNodes), *nodes = id + 1;

583

GmshElement *element = elements + c;

584

element->dim = dim;

585

element->cellType = cellType;

586

element->numNodes = numNodes;

587

element->numTags = numTags;

588

element->id = id[0];

589

for (p = 0; p < numNodes; p++) element->nodes[p] = nodes[p];

590

for (p = 0; p < numTags; p++) element->tags[p] = tags[p];

591

}

592

}

593

594

}

595

596

static PetscErrorCode DMPlexCreateGmsh_ReadPeriodic_v40(GmshFile *gmsh, int shift, PetscInt slaveMap[], PetscBT bt)

597

{

598

PetscViewer viewer = gmsh->viewer;

599

int fileFormat = gmsh->fileFormat;

600

PetscBool binary = gmsh->binary;

601

PetscBool byteSwap = gmsh->byteSwap;

602

int numPeriodic, snum, i;

603

char line[PETSC_MAX_PATH_LEN4096];

604

PetscErrorCode ierr;

605

606

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 606; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

607

if (fileFormat == 22 || !binary) {

608

ierr = PetscViewerRead(viewer, line, 1, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),608,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

609

snum = sscanf(line, "%d", &numPeriodic);

610

if (snum != 1) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),610,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

611

} else {

612

ierr = PetscViewerRead(viewer, &numPeriodic, 1, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),612,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

613

if (byteSwap) {ierr = PetscByteSwap(&numPeriodic, PETSC_ENUM, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),613,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

614

}

615

for (i = 0; i < numPeriodic; i++) {

616

int ibuf[3], slaveDim = -1, slaveTag = -1, masterTag = -1, slaveNode, masterNode;

617

long j, nNodes;

618

double affine[16];

619

620

if (fileFormat == 22 || !binary) {

621

ierr = PetscViewerRead(viewer, line, 3, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),621,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

622

snum = sscanf(line, "%d %d %d", &slaveDim, &slaveTag, &masterTag);

623

if (snum != 3) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),623,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

624

} else {

625

ierr = PetscViewerRead(viewer, ibuf, 3, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),625,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

626

if (byteSwap) {ierr = PetscByteSwap(ibuf, PETSC_ENUM, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),626,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

627

slaveDim = ibuf[0]; slaveTag = ibuf[1]; masterTag = ibuf[2];

628

}

629

(void)slaveDim; (void)slaveTag; (void)masterTag; /* unused */

630

631

if (fileFormat == 22 || !binary) {

632

ierr = PetscViewerRead(viewer, line, 1, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),632,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

633

snum = sscanf(line, "%ld", &nNodes);

634

if (snum != 1) { /* discard transformation and try again */

635

ierr = PetscViewerRead(viewer, line, -PETSC_MAX_PATH_LEN4096, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),635,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

636

ierr = PetscViewerRead(viewer, line, 1, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),636,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

637

snum = sscanf(line, "%ld", &nNodes);

638

if (snum != 1) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),638,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

639

}

640

} else {

641

ierr = PetscViewerRead(viewer, &nNodes, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),641,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

642

if (byteSwap) {ierr = PetscByteSwap(&nNodes, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),642,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

643

if (nNodes == -1) { /* discard transformation and try again */

644

ierr = PetscViewerRead(viewer, affine, 16, NULL((void*)0), PETSC_DOUBLE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),644,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

645

ierr = PetscViewerRead(viewer, &nNodes, 1, NULL((void*)0), PETSC_LONG);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),645,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

646

if (byteSwap) {ierr = PetscByteSwap(&nNodes, PETSC_LONG, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),646,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

647

}

648

}

649

650

for (j = 0; j < nNodes; j++) {

651

if (fileFormat == 22 || !binary) {

652

ierr = PetscViewerRead(viewer, line, 2, NULL((void*)0), PETSC_STRING);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),652,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

653

snum = sscanf(line, "%d %d", &slaveNode, &masterNode);

654

if (snum != 2) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),654,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

655

} else {

656

ierr = PetscViewerRead(viewer, ibuf, 2, NULL((void*)0), PETSC_ENUM);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),656,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

657

if (byteSwap) {ierr = PetscByteSwap(ibuf, PETSC_ENUM, 2);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),657,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

658

slaveNode = ibuf[0]; masterNode = ibuf[1];

659

}

660

slaveMap[slaveNode - shift] = masterNode - shift;

661

ierr = PetscBTSet(bt, slaveNode - shift);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),661,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

662

ierr = PetscBTSet(bt, masterNode - shift);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),662,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

663

}

664

}

665

666

}

667

668

669

$Entities

670

numPoints(size_t) numCurves(size_t)

671

numSurfaces(size_t) numVolumes(size_t)

672

pointTag(int) X(double) Y(double) Z(double)

673

numPhysicalTags(size_t) physicalTag(int) ...

674

...

675

curveTag(int) minX(double) minY(double) minZ(double)

676

maxX(double) maxY(double) maxZ(double)

677

numPhysicalTags(size_t) physicalTag(int) ...

678

numBoundingPoints(size_t) pointTag(int) ...

679

...

680

surfaceTag(int) minX(double) minY(double) minZ(double)

681

maxX(double) maxY(double) maxZ(double)

682

numPhysicalTags(size_t) physicalTag(int) ...

683

numBoundingCurves(size_t) curveTag(int) ...

684

...

685

volumeTag(int) minX(double) minY(double) minZ(double)

686

maxX(double) maxY(double) maxZ(double)

687

numPhysicalTags(size_t) physicalTag(int) ...

688

numBoundngSurfaces(size_t) surfaceTag(int) ...

689

...

690

$EndEntities

691

692

static PetscErrorCode DMPlexCreateGmsh_ReadEntities_v41(GmshFile *gmsh, GmshEntities **entities)

693

{

694

PetscInt count[4], index, numTags, i;

695

int dim, eid, *tags = NULL((void*)0);

696

GmshEntity *entity = NULL((void*)0);

697

PetscErrorCode ierr;

698

699

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 699; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

700

ierr = GmshReadSize(gmsh, count, 4);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),700,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

701

ierr = GmshEntitiesCreate(count, entities);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),701,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

702

for (dim = 0; dim < 4; ++dim) {

703

for (index = 0; index < count[dim]; ++index) {

704

ierr = GmshReadInt(gmsh, &eid, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),704,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

705

ierr = GmshEntitiesAdd(*entities, (PetscInt)index, dim, eid, &entity);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),705,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

706

ierr = GmshReadDouble(gmsh, entity->bbox, (dim == 0) ? 3 : 6);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),706,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

707

ierr = GmshReadSize(gmsh, &numTags, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),707,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

708

ierr = GmshBufferGet(gmsh, numTags, sizeof(int), &tags);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),708,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

709

ierr = GmshReadInt(gmsh, tags, numTags);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),709,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

710

entity->numTags = PetscMin(numTags, 4)(((numTags)<(4)) ? (numTags) : (4));

711

for (i = 0; i < entity->numTags; ++i) entity->tags[i] = tags[i];

712

if (dim == 0) continue;

713

ierr = GmshReadSize(gmsh, &numTags, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),713,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

714

ierr = GmshBufferGet(gmsh, numTags, sizeof(int), &tags);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),714,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

715

ierr = GmshReadInt(gmsh, tags, numTags);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),715,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

716

}

717

}

718

719

}

720

721

722

$Nodes

723

numEntityBlocks(size_t) numNodes(size_t)

724

minNodeTag(size_t) maxNodeTag(size_t)

725

entityDim(int) entityTag(int) parametric(int; 0 or 1) numNodesBlock(size_t)

726

nodeTag(size_t)

727

...

728

x(double) y(double) z(double)

729

< u(double; if parametric and entityDim = 1 or entityDim = 2) >

730

< v(double; if parametric and entityDim = 2) >

731

...

732

...

733

$EndNodes

734

735

static PetscErrorCode DMPlexCreateGmsh_ReadNodes_v41(GmshFile *gmsh, int shift, PetscInt *numVertices, double **gmsh_nodes)

736

{

737

int info[3];

738

PetscInt sizes[4], numEntityBlocks, numNodes, numNodesBlock = 0, *nodeTag = NULL((void*)0), block, node, i;

739

double *coordinates;

740

PetscErrorCode ierr;

741

742

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 742; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

743

ierr = GmshReadSize(gmsh, sizes, 4);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),743,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

744

numEntityBlocks = sizes[0]; numNodes = sizes[1];

745

ierr = PetscMalloc1(numNodes*3, &coordinates)PetscMallocA(1,PETSC_FALSE,745,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numNodes*3)*sizeof(**(&coordinates)),(&coordinates
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),745,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

746

*numVertices = numNodes;

747

*gmsh_nodes = coordinates;

748

for (block = 0, node = 0; block < numEntityBlocks; ++block, node += numNodesBlock) {

749

ierr = GmshReadInt(gmsh, info, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),749,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

750

ierr = GmshReadSize(gmsh, &numNodesBlock, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),750,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

751

ierr = GmshBufferGet(gmsh, numNodesBlock, sizeof(PetscInt), &nodeTag);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),751,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

752

ierr = GmshReadSize(gmsh, nodeTag, numNodesBlock);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),752,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

753

for (i = 0; i < numNodesBlock; ++i) if (nodeTag[i] != node+i+shift) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unexpected node number %d should be %d", nodeTag[i], node+i+shift)return PetscError(((MPI_Comm)0x44000001),753,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unexpected node number %d should be %d"
,nodeTag[i],node+i+shift);

754

if (info[2] != 0) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Parametric coordinates not supported")return PetscError(((MPI_Comm)0x44000001),754,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Parametric coordinates not supported"
);

755

ierr = GmshReadDouble(gmsh, coordinates+node*3, numNodesBlock*3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),755,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

756

}

757

758

}

759

760

761

$Elements

762

numEntityBlocks(size_t) numElements(size_t)

763

minElementTag(size_t) maxElementTag(size_t)

764

entityDim(int) entityTag(int) elementType(int; see below) numElementsBlock(size_t)

765

elementTag(size_t) nodeTag(size_t) ...

766

...

767

...

768

$EndElements

769

770

static PetscErrorCode DMPlexCreateGmsh_ReadElements_v41(GmshFile *gmsh, PETSC_UNUSED__attribute((unused)) int shift, GmshEntities *entities, PetscInt *numCells, GmshElement **gmsh_elems)

771

{

772

int info[3], eid, dim, cellType, *tags;

773

PetscInt sizes[4], *ibuf = NULL((void*)0), numEntityBlocks, numElements, numBlockElements, numNodes, numTags, block, elem, c, p;

774

GmshEntity *entity = NULL((void*)0);

775

GmshElement *elements;

776

PetscErrorCode ierr;

777

778

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 778; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

779

ierr = GmshReadSize(gmsh, sizes, 4);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),779,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

780

numEntityBlocks = sizes[0]; numElements = sizes[1];

781

ierr = PetscCalloc1(numElements, &elements)PetscMallocA(1,PETSC_TRUE,781,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numElements)*sizeof(**(&elements)),(&elements
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),781,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

782

*numCells = numElements;

783

*gmsh_elems = elements;

784

for (c = 0, block = 0; block < numEntityBlocks; ++block) {

785

ierr = GmshReadInt(gmsh, info, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),785,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

786

dim = info[0]; eid = info[1]; cellType = info[2];

787

ierr = GmshEntitiesGet(entities, dim, eid, &entity);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),787,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

788

numTags = entity->numTags;

789

tags = entity->tags;

790

switch (cellType) {

791

case 1: /* 2-node line */

792

numNodes = 2;

793

break;

794

case 2: /* 3-node triangle */

795

numNodes = 3;

796

break;

797

case 3: /* 4-node quadrangle */

798

numNodes = 4;

799

break;

800

case 4: /* 4-node tetrahedron */

801

numNodes = 4;

802

break;

803

case 5: /* 8-node hexahedron */

804

numNodes = 8;

805

break;

806

case 6: /* 6-node wedge */

807

numNodes = 6;

808

break;

809

case 15: /* 1-node vertex */

810

numNodes = 1;

811

break;

812

default:

813

SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unsupported Gmsh element type %d", cellType)return PetscError(((MPI_Comm)0x44000001),813,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unsupported Gmsh element type %d",cellType
);

814

}

815

ierr = GmshReadSize(gmsh, &numBlockElements, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),815,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

816

ierr = GmshBufferGet(gmsh, (1+numNodes)*numBlockElements, sizeof(PetscInt), &ibuf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),816,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

817

ierr = GmshReadSize(gmsh, ibuf, (1+numNodes)*numBlockElements);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),817,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

818

for (elem = 0; elem < numBlockElements; ++elem, ++c) {

819

GmshElement *element = elements + c;

820

PetscInt *id = ibuf + elem*(1+numNodes), *nodes = id + 1;

821

element->id = id[0];

822

element->dim = dim;

823

element->cellType = cellType;

824

element->numNodes = numNodes;

825

element->numTags = numTags;

826

for (p = 0; p < numNodes; p++) element->nodes[p] = nodes[p];

827

for (p = 0; p < numTags; p++) element->tags[p] = tags[p];

828

}

829

}

830

831

}

832

833

834

$Periodic

835

numPeriodicLinks(size_t)

836

entityDim(int) entityTag(int) entityTagMaster(int)

837

numAffine(size_t) value(double) ...

838

numCorrespondingNodes(size_t)

839

nodeTag(size_t) nodeTagMaster(size_t)

840

...

841

...

842

$EndPeriodic

843

844

static PetscErrorCode DMPlexCreateGmsh_ReadPeriodic_v41(GmshFile *gmsh, int shift, PetscInt slaveMap[], PetscBT bt)

845

{

846

int info[3];

847

PetscInt numPeriodicLinks, numAffine, numCorrespondingNodes, *nodeTags = NULL((void*)0), link, node;

'numPeriodicLinks' declared without an initial value

→

848

double dbuf[16];

849

PetscErrorCode ierr;

850

851

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 851; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

852

ierr = GmshReadSize(gmsh, &numPeriodicLinks, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),852,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

←

Calling 'GmshReadSize'

→

←

Returning from 'GmshReadSize'

→

853

for (link = 0; link < numPeriodicLinks; ++link) {

←

The right operand of '<' is a garbage value

854

ierr = GmshReadInt(gmsh, info, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),854,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

855

ierr = GmshReadSize(gmsh, &numAffine, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),855,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

856

ierr = GmshReadDouble(gmsh, dbuf, numAffine);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),856,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

857

ierr = GmshReadSize(gmsh, &numCorrespondingNodes, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),857,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

858

ierr = GmshBufferGet(gmsh, numCorrespondingNodes, sizeof(PetscInt), &nodeTags);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),858,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

859

ierr = GmshReadSize(gmsh, nodeTags, numCorrespondingNodes*2);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),859,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

860

for (node = 0; node < numCorrespondingNodes; ++node) {

861

PetscInt slaveNode = nodeTags[node*2+0] - shift;

862

PetscInt masterNode = nodeTags[node*2+1] - shift;

863

slaveMap[slaveNode] = masterNode;

864

ierr = PetscBTSet(bt, slaveNode);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),864,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

865

ierr = PetscBTSet(bt, masterNode);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),865,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

866

}

867

}

868

869

}

870

871

872

$MeshFormat // same as MSH version 2

873

version(ASCII double; currently 4.1)

874

file-type(ASCII int; 0 for ASCII mode, 1 for binary mode)

875

data-size(ASCII int; sizeof(size_t))

876

< int with value one; only in binary mode, to detect endianness >

877

$EndMeshFormat

878

879

static PetscErrorCode DMPlexCreateGmsh_ReadMeshFormat(GmshFile *gmsh)

880

{

881

char line[PETSC_MAX_PATH_LEN4096];

882

int snum, fileType, fileFormat, dataSize, checkEndian;

883

float version;

884

PetscErrorCode ierr;

885

886

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 886; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

887

ierr = GmshReadString(gmsh, line, 3);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),887,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

888

snum = sscanf(line, "%f %d %d", &version, &fileType, &dataSize);

889

if (snum != 3) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unable to parse Gmsh file header: %s", line)return PetscError(((MPI_Comm)0x44000001),889,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unable to parse Gmsh file header: %s"
,line);

890

if (version < 2.2) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Gmsh file version %3.1f must be at least 2.2", (double)version)return PetscError(((MPI_Comm)0x44000001),890,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Gmsh file version %3.1f must be at least 2.2"
,(double)version);

891

if ((int)version == 3) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Gmsh file version %3.1f not supported", (double)version)return PetscError(((MPI_Comm)0x44000001),891,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Gmsh file version %3.1f not supported"
,(double)version);

892

if (version > 4.1) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Gmsh file version %3.1f must be at most 4.1", (double)version)return PetscError(((MPI_Comm)0x44000001),892,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Gmsh file version %3.1f must be at most 4.1"
,(double)version);

893

if (gmsh->binary && !fileType) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Viewer is binary but Gmsh file is ASCII")return PetscError(((MPI_Comm)0x44000001),893,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Viewer is binary but Gmsh file is ASCII"
);

894

if (!gmsh->binary && fileType) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Viewer is ASCII but Gmsh file is binary")return PetscError(((MPI_Comm)0x44000001),894,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Viewer is ASCII but Gmsh file is binary"
);

895

fileFormat = (int)(version*10); /* XXX Should use (int)roundf(version*10) ? */

896

if (fileFormat <= 40 && dataSize != sizeof(double)) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Data size %d is not valid for a Gmsh file", dataSize)return PetscError(((MPI_Comm)0x44000001),896,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Data size %d is not valid for a Gmsh file"
,dataSize);

897

if (fileFormat >= 41 && dataSize != sizeof(int) && dataSize != sizeof(PetscInt64)) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Data size %d is not valid for a Gmsh file", dataSize)return PetscError(((MPI_Comm)0x44000001),897,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Data size %d is not valid for a Gmsh file"
,dataSize);

898

gmsh->fileFormat = fileFormat;

899

gmsh->dataSize = dataSize;

900

gmsh->byteSwap = PETSC_FALSE;

901

if (gmsh->binary) {

902

ierr = GmshReadInt(gmsh, &checkEndian, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),902,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

903

if (checkEndian != 1) {

904

ierr = PetscByteSwap(&checkEndian, PETSC_ENUM, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),904,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

905

if (checkEndian != 1) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unable to detect endianness in Gmsh file header: %s", line)return PetscError(((MPI_Comm)0x44000001),905,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unable to detect endianness in Gmsh file header: %s"
,line);

906

gmsh->byteSwap = PETSC_TRUE;

907

}

908

}

909

910

}

911

912

913

PhysicalNames

914

numPhysicalNames(ASCII int)

915

dimension(ASCII int) physicalTag(ASCII int) "name"(127 characters max)

916

...

917

$EndPhysicalNames

918

919

static PetscErrorCode DMPlexCreateGmsh_ReadPhysicalNames(GmshFile *gmsh)

920

{

921

char line[PETSC_MAX_PATH_LEN4096], name[128+2], *p, *q;

922

int snum, numRegions, region, dim, tag;

923

PetscErrorCode ierr;

924

925

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 925; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

926

ierr = GmshReadString(gmsh, line, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),926,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

927

snum = sscanf(line, "%d", &numRegions);

928

if (snum != 1) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),928,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

929

for (region = 0; region < numRegions; ++region) {

930

ierr = GmshReadString(gmsh, line, 2);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),930,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

931

snum = sscanf(line, "%d %d", &dim, &tag);

932

if (snum != 2) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),932,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

933

ierr = GmshReadString(gmsh, line, -(PetscInt)sizeof(line));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),933,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

934

ierr = PetscStrchr(line, '"', &p);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),934,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

935

if (!p) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),935,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

936

ierr = PetscStrrchr(line, '"', &q);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),936,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

937

if (q == p) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "File is not a valid Gmsh file")return PetscError(((MPI_Comm)0x44000001),937,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"File is not a valid Gmsh file");

938

ierr = PetscStrncpy(name, p+1, (size_t)(q-p-1));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),938,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

939

}

940

941

}

942

943

/*@C

944

DMPlexCreateGmshFromFile - Create a DMPlex mesh from a Gmsh file

945

946

+ comm - The MPI communicator

947

. filename - Name of the Gmsh file

948

- interpolate - Create faces and edges in the mesh

949

950

Output Parameter:

951

. dm - The DM object representing the mesh

952

953

Level: beginner

954

955

.seealso: DMPlexCreateFromFile(), DMPlexCreateGmsh(), DMPlexCreate()

956

@*/

957

PetscErrorCode DMPlexCreateGmshFromFile(MPI_Comm comm, const char filename[], PetscBool interpolate, DM *dm)

958

{

959

PetscViewer viewer;

960

PetscMPIInt rank;

961

int fileType;

962

PetscViewerType vtype;

963

PetscErrorCode ierr;

964

965

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 965; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

966

ierr = MPI_Comm_rank(comm, &rank);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),966,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

967

968

/* Determine Gmsh file type (ASCII or binary) from file header */

969

if (!rank) {

970

GmshFile gmsh_, *gmsh = &gmsh_;

971

char line[PETSC_MAX_PATH_LEN4096];

972

int snum;

973

float version;

974

975

ierr = PetscArrayzero(gmsh,1)PetscMemzero(gmsh,(1)*sizeof(*(gmsh)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),975,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

976

ierr = PetscViewerCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001), &gmsh->viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),976,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

977

ierr = PetscViewerSetType(gmsh->viewer, PETSCVIEWERASCII"ascii");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),977,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

978

ierr = PetscViewerFileSetMode(gmsh->viewer, FILE_MODE_READ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),978,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

979

ierr = PetscViewerFileSetName(gmsh->viewer, filename);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),979,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

980

/* Read only the first two lines of the Gmsh file */

981

ierr = GmshReadSection(gmsh, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),981,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

982

ierr = GmshExpect(gmsh, "$MeshFormat", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),982,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

983

ierr = GmshReadString(gmsh, line, 2);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),983,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

984

snum = sscanf(line, "%f %d", &version, &fileType);

985

if (snum != 2) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unable to parse Gmsh file header: %s", line)return PetscError(((MPI_Comm)0x44000001),985,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Unable to parse Gmsh file header: %s"
,line);

986

if (version < 2.2) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Gmsh file version %3.1f must be at least 2.2", (double)version)return PetscError(((MPI_Comm)0x44000001),986,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Gmsh file version %3.1f must be at least 2.2"
,(double)version);

987

if ((int)version == 3) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Gmsh file version %3.1f not supported", (double)version)return PetscError(((MPI_Comm)0x44000001),987,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Gmsh file version %3.1f not supported"
,(double)version);

988

if (version > 4.1) SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Gmsh file version %3.1f must be at most 4.1", (double)version)return PetscError(((MPI_Comm)0x44000001),988,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Gmsh file version %3.1f must be at most 4.1"
,(double)version);

989

ierr = PetscViewerDestroy(&gmsh->viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),989,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

990

}

991

ierr = MPI_Bcast(&fileType, 1, MPI_INT, 0, comm)((petsc_allreduce_ct += PetscMPIParallelComm((comm)),0) || MPI_Bcast
((&fileType),(1),(((MPI_Datatype)0x4c000405)),(0),(comm))
);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),991,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

992

vtype = (fileType == 0) ? PETSCVIEWERASCII"ascii" : PETSCVIEWERBINARY"binary";

993

994

/* Create appropriate viewer and build plex */

995

ierr = PetscViewerCreate(comm, &viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),995,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

996

ierr = PetscViewerSetType(viewer, vtype);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),996,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

997

ierr = PetscViewerFileSetMode(viewer, FILE_MODE_READ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),997,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

998

ierr = PetscViewerFileSetName(viewer, filename);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),998,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

999

ierr = DMPlexCreateGmsh(comm, viewer, interpolate, dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),999,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1000

ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1000,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1001

1002

}

1003

1004

/*@

1005

DMPlexCreateGmsh - Create a DMPlex mesh from a Gmsh file viewer

1006

1007

Collective

1008

1009

Input Parameters:

1010

+ comm - The MPI communicator

1011

. viewer - The Viewer associated with a Gmsh file

1012

- interpolate - Create faces and edges in the mesh

1013

1014

Output Parameter:

1015

. dm - The DM object representing the mesh

1016

1017

Note: http://gmsh.info/doc/texinfo/gmsh.html#MSH-file-format

1018

1019

Level: beginner

1020

1021

.seealso: DMPLEX, DMCreate()

1022

@*/

1023

PetscErrorCode DMPlexCreateGmsh(MPI_Comm comm, PetscViewer viewer, PetscBool interpolate, DM *dm)

1024

{

1025

PetscViewer parentviewer = NULL((void*)0);

1026

double *coordsIn = NULL((void*)0);

1027

GmshEntities *entities = NULL((void*)0);

1028

GmshElement *gmsh_elem = NULL((void*)0);

1029

PetscSection coordSection;

1030

Vec coordinates;

1031

PetscBT periodicV = NULL((void*)0), periodicC = NULL((void*)0);

1032

PetscScalar *coords;

1033

PetscInt dim = 0, embedDim = -1, coordSize, c, v, d, cell, *periodicMap = NULL((void*)0), *periodicMapI = NULL((void*)0), *hybridMap = NULL((void*)0), cMax = PETSC_DETERMINE-1;

1034

PetscInt numVertices = 0, numCells = 0, trueNumCells = 0;

1035

int i, shift = 1;

1036

PetscMPIInt rank;

1037

PetscBool binary, zerobase = PETSC_FALSE, usemarker = PETSC_FALSE;

1038

PetscBool enable_hybrid = interpolate, periodic = PETSC_TRUE;

1039

PetscErrorCode ierr;

1040

1041

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
; petscstack->line[petscstack->currentsize] = 1041; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1042

ierr = MPI_Comm_rank(comm, &rank);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1042,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1043

ierr = PetscObjectOptionsBegin((PetscObject)viewer)0; do { PetscOptionItems PetscOptionsObjectBase; PetscOptionItems
*PetscOptionsObject = &PetscOptionsObjectBase; PetscOptionsObject
->options = ((PetscObject)(PetscObject)viewer)->options
; for (PetscOptionsObject->count=(PetscOptionsPublish?-1:1
); PetscOptionsObject->count<2; PetscOptionsObject->
count++) { PetscErrorCode _5_ierr = PetscObjectOptionsBegin_Private
(PetscOptionsObject,(PetscObject)viewer);do {if (__builtin_expect
(!!(_5_ierr),0)) return PetscError(((MPI_Comm)0x44000001),1043
,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,_5_ierr,PETSC_ERROR_REPEAT," ");} while (0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1043,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1044

ierr = PetscOptionsHead(PetscOptionsObject,"DMPlex Gmsh options");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1044,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1045

ierr = PetscOptionsBool("-dm_plex_gmsh_hybrid", "Generate hybrid cell bounds", "DMPlexCreateGmsh", enable_hybrid, &enable_hybrid, NULL)PetscOptionsBool_Private(PetscOptionsObject,"-dm_plex_gmsh_hybrid"
,"Generate hybrid cell bounds","DMPlexCreateGmsh",enable_hybrid
,&enable_hybrid,((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1045,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1046

ierr = PetscOptionsBool("-dm_plex_gmsh_periodic","Read Gmsh periodic section", "DMPlexCreateGmsh", periodic, &periodic, NULL)PetscOptionsBool_Private(PetscOptionsObject,"-dm_plex_gmsh_periodic"
,"Read Gmsh periodic section","DMPlexCreateGmsh",periodic,&
periodic,((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1046,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1047

ierr = PetscOptionsBool("-dm_plex_gmsh_use_marker", "Generate marker label", "DMPlexCreateGmsh", usemarker, &usemarker, NULL)PetscOptionsBool_Private(PetscOptionsObject,"-dm_plex_gmsh_use_marker"
,"Generate marker label","DMPlexCreateGmsh",usemarker,&usemarker
,((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1047,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1048

ierr = PetscOptionsBool("-dm_plex_gmsh_zero_base", "Read Gmsh file with zero base indices", "DMPlexCreateGmsh", zerobase, &zerobase, NULL)PetscOptionsBool_Private(PetscOptionsObject,"-dm_plex_gmsh_zero_base"
,"Read Gmsh file with zero base indices","DMPlexCreateGmsh",zerobase
,&zerobase,((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1048,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1049

ierr = PetscOptionsBoundedInt("-dm_plex_gmsh_spacedim", "Embedding space dimension", "DMPlexCreateGmsh", embedDim, &embedDim, NULL,PETSC_DECIDE)PetscOptionsInt_Private(PetscOptionsObject,"-dm_plex_gmsh_spacedim"
,"Embedding space dimension","DMPlexCreateGmsh",embedDim,&
embedDim,((void*)0),-1,2147483647);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1049,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1050

ierr = PetscOptionsTail()0; do {if (PetscOptionsObject->count != 1) do { do { ; if (
petscstack && petscstack->currentsize > 0) { petscstack
->currentsize--; petscstack->function[petscstack->currentsize
] = 0; petscstack->file[petscstack->currentsize] = 0; petscstack
->line[petscstack->currentsize] = 0; petscstack->petscroutine
[petscstack->currentsize] = PETSC_FALSE; } if (petscstack)
{ petscstack->hotdepth = (((petscstack->hotdepth-1)<
(0)) ? (0) : (petscstack->hotdepth-1)); } ; } while (0); return
(0);} while (0);} while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1050,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1051

ierr = PetscOptionsEnd()_5_ierr = PetscOptionsEnd_Private(PetscOptionsObject);do {if (
__builtin_expect(!!(_5_ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1051,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,_5_ierr,PETSC_ERROR_REPEAT," ");} while (0);}} while (0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1051,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1052

if (zerobase) shift = 0;

1053

1054

ierr = DMCreate(comm, dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1054,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1055

ierr = DMSetType(*dm, DMPLEX"plex");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1055,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1056

ierr = PetscLogEventBegin(DMPLEX_CreateGmsh,*dm,0,0,0)(((PetscLogPLB && petsc_stageLog->stageInfo[petsc_stageLog
->curStage].perfInfo.active && petsc_stageLog->
stageInfo[petsc_stageLog->curStage].eventLog->eventInfo
[DMPLEX_CreateGmsh].active) ? (*PetscLogPLB)((DMPLEX_CreateGmsh
),0,(PetscObject)(*dm),(PetscObject)(0),(PetscObject)(0),(PetscObject
)(0)) : 0 ));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1056,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1057

1058

ierr = PetscObjectTypeCompare((PetscObject)viewer, PETSCVIEWERBINARY"binary", &binary);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1058,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1059

1060

/* Binary viewers read on all ranks, get subviewer to read only in rank 0 */

1061

if (binary) {

1062

parentviewer = viewer;

1063

ierr = PetscViewerGetSubViewer(parentviewer, PETSC_COMM_SELF((MPI_Comm)0x44000001), &viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1063,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1064

}

1065

1066

if (!rank) {

1067

GmshFile gmsh_, *gmsh = &gmsh_;

1068

char line[PETSC_MAX_PATH_LEN4096];

1069

PetscBool match;

1070

1071

ierr = PetscArrayzero(gmsh,1)PetscMemzero(gmsh,(1)*sizeof(*(gmsh)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1071,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1072

gmsh->viewer = viewer;

1073

gmsh->binary = binary;

1074

1075

/* Read mesh format */

1076

ierr = GmshReadSection(gmsh, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1076,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1077

ierr = GmshExpect(gmsh, "$MeshFormat", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1077,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1078

ierr = DMPlexCreateGmsh_ReadMeshFormat(gmsh);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1078,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1079

ierr = GmshReadEndSection(gmsh, "$EndMeshFormat", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1079,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1080

1081

/* OPTIONAL Read physical names */

1082

ierr = GmshReadSection(gmsh, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1082,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1083

ierr = GmshMatch(gmsh,"$PhysicalNames", line, &match);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1083,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1084

if (match) {

1085

ierr = DMPlexCreateGmsh_ReadPhysicalNames(gmsh);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1085,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1086

ierr = GmshReadEndSection(gmsh, "$EndPhysicalNames", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1086,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1087

/* Initial read for entity section */

1088

ierr = GmshReadSection(gmsh, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1088,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1089

}

1090

1091

/* Read entities */

1092

if (gmsh->fileFormat >= 40) {

1093

ierr = GmshExpect(gmsh, "$Entities", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1093,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1094

switch (gmsh->fileFormat) {

1095

case 41: ierr = DMPlexCreateGmsh_ReadEntities_v41(gmsh, &entities);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1095,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1096

default: ierr = DMPlexCreateGmsh_ReadEntities_v40(gmsh, &entities);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1096,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1097

}

1098

ierr = GmshReadEndSection(gmsh, "$EndEntities", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1098,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1099

/* Initial read for nodes section */

1100

ierr = GmshReadSection(gmsh, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1100,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1101

}

1102

1103

/* Read nodes */

1104

ierr = GmshExpect(gmsh, "$Nodes", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1104,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1105

switch (gmsh->fileFormat) {

1106

case 41: ierr = DMPlexCreateGmsh_ReadNodes_v41(gmsh, shift, &numVertices, &coordsIn);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1106,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1107

case 40: ierr = DMPlexCreateGmsh_ReadNodes_v40(gmsh, shift, &numVertices, &coordsIn);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1107,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1108

default: ierr = DMPlexCreateGmsh_ReadNodes_v22(gmsh, shift, &numVertices, &coordsIn);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1108,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1109

}

1110

ierr = GmshReadEndSection(gmsh, "$EndNodes", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1110,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1111

1112

/* Read elements */

1113

ierr = GmshReadSection(gmsh, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1113,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);;

1114

ierr = GmshExpect(gmsh, "$Elements", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1114,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1115

switch (gmsh->fileFormat) {

1116

case 41: ierr = DMPlexCreateGmsh_ReadElements_v41(gmsh, shift, entities, &numCells, &gmsh_elem);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1116,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1117

case 40: ierr = DMPlexCreateGmsh_ReadElements_v40(gmsh, shift, entities, &numCells, &gmsh_elem);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1117,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1118

default: ierr = DMPlexCreateGmsh_ReadElements_v22(gmsh, shift, &numCells, &gmsh_elem);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1118,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1119

}

1120

ierr = GmshReadEndSection(gmsh, "$EndElements", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1120,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1121

1122

/* OPTIONAL Read periodic section */

1123

if (periodic) {

1124

ierr = GmshReadSection(gmsh, line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1124,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1125

ierr = GmshMatch(gmsh, "$Periodic", line, &periodic);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1125,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1126

}

1127

if (periodic) {

1128

PetscInt pVert, *periodicMapT, *aux;

1129

1130

ierr = PetscMalloc1(numVertices, &periodicMapT)PetscMallocA(1,PETSC_FALSE,1130,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numVertices)*sizeof(**(&periodicMapT)),(&periodicMapT
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1130,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1131

ierr = PetscBTCreate(numVertices, &periodicV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1131,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1132

for (i = 0; i < numVertices; i++) periodicMapT[i] = i;

1133

1134

ierr = GmshExpect(gmsh, "$Periodic", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1134,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1135

switch (gmsh->fileFormat) {

1136

case 41: ierr = DMPlexCreateGmsh_ReadPeriodic_v41(gmsh, shift, periodicMapT, periodicV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1136,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1137

default: ierr = DMPlexCreateGmsh_ReadPeriodic_v40(gmsh, shift, periodicMapT, periodicV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1137,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); break;

1138

}

1139

ierr = GmshReadEndSection(gmsh, "$EndPeriodic", line);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1139,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1140

1141

/* we may have slaves of slaves */

1142

for (i = 0; i < numVertices; i++) {

1143

while (periodicMapT[periodicMapT[i]] != periodicMapT[i]) {

1144

periodicMapT[i] = periodicMapT[periodicMapT[i]];

1145

}

1146

}

1147

/* periodicMap : from old to new numbering (periodic vertices excluded)

1148

periodicMapI: from new to old numbering */

1149

ierr = PetscMalloc1(numVertices, &periodicMap)PetscMallocA(1,PETSC_FALSE,1149,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numVertices)*sizeof(**(&periodicMap)),(&periodicMap
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1149,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1150

ierr = PetscMalloc1(numVertices, &periodicMapI)PetscMallocA(1,PETSC_FALSE,1150,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numVertices)*sizeof(**(&periodicMapI)),(&periodicMapI
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1150,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1151

ierr = PetscMalloc1(numVertices, &aux)PetscMallocA(1,PETSC_FALSE,1151,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(numVertices)*sizeof(**(&aux)),(&aux));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1151,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1152

for (i = 0, pVert = 0; i < numVertices; i++) {

1153

if (periodicMapT[i] != i) {

1154

pVert++;

1155

} else {

1156

aux[i] = i - pVert;

1157

periodicMapI[i - pVert] = i;

1158

}

1159

}

1160

for (i = 0 ; i < numVertices; i++) {

1161

periodicMap[i] = aux[periodicMapT[i]];

1162

}

1163

ierr = PetscFree(periodicMapT)((*PetscTrFree)((void*)(periodicMapT),1163,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((periodicMapT) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1163,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1164

ierr = PetscFree(aux)((*PetscTrFree)((void*)(aux),1164,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((aux) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1164,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1165

/* remove periodic vertices */

1166

numVertices = numVertices - pVert;

1167

}

1168

1169

ierr = GmshEntitiesDestroy(&entities);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1169,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1170

ierr = PetscFree(gmsh->wbuf)((*PetscTrFree)((void*)(gmsh->wbuf),1170,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((gmsh->wbuf) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1170,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1171

ierr = PetscFree(gmsh->sbuf)((*PetscTrFree)((void*)(gmsh->sbuf),1171,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((gmsh->sbuf) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1171,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1172

}

1173

1174

if (parentviewer) {

1175

ierr = PetscViewerRestoreSubViewer(parentviewer, PETSC_COMM_SELF((MPI_Comm)0x44000001), &viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1175,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1176

}

1177

1178

if (!rank) {

1179

PetscBool hybrid = PETSC_FALSE;

1180

PetscInt cellType = -1;

1181

1182

for (trueNumCells = 0, c = 0; c < numCells; ++c) {

1183

if (gmsh_elem[c].dim > dim) {dim = gmsh_elem[c].dim; trueNumCells = 0; hybrid = PETSC_FALSE; cellType = -1;}

1184

if (gmsh_elem[c].dim < dim) continue;

1185

if (cellType == -1) cellType = gmsh_elem[c].cellType;

1186

/* different cell type indicate an hybrid mesh in PLEX */

1187

if (cellType != gmsh_elem[c].cellType) hybrid = PETSC_TRUE;

1188

/* wedges always indicate an hybrid mesh in PLEX */

1189

if (cellType == 6 || cellType == 13) hybrid = PETSC_TRUE;

1190

trueNumCells++;

1191

}

1192

/* Renumber cells for hybrid grids */

1193

if (hybrid && enable_hybrid) {

1194

PetscInt hc1 = 0, hc2 = 0, *hybridCells1 = NULL((void*)0), *hybridCells2 = NULL((void*)0);

1195

PetscInt cell, tn, *tp;

1196

int n1 = 0,n2 = 0;

1197

1198

ierr = PetscMalloc1(trueNumCells, &hybridCells1)PetscMallocA(1,PETSC_FALSE,1198,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(trueNumCells)*sizeof(**(&hybridCells1)),(&hybridCells1
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1198,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1199

ierr = PetscMalloc1(trueNumCells, &hybridCells2)PetscMallocA(1,PETSC_FALSE,1199,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(trueNumCells)*sizeof(**(&hybridCells2)),(&hybridCells2
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1199,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1200

for (cell = 0, c = 0; c < numCells; ++c) {

1201

if (gmsh_elem[c].dim == dim) {

1202

if (!n1) n1 = gmsh_elem[c].cellType;

1203

else if (!n2 && n1 != gmsh_elem[c].cellType) n2 = gmsh_elem[c].cellType;

1204

1205

if (gmsh_elem[c].cellType == n1) hybridCells1[hc1++] = cell;

1206

else if (gmsh_elem[c].cellType == n2) hybridCells2[hc2++] = cell;

1207

else SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Cannot handle more than 2 cell types")return PetscError(((MPI_Comm)0x44000001),1207,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Cannot handle more than 2 cell types"
);

1208

cell++;

1209

}

1210

}

1211

1212

switch (n1) {

1213

case 2: /* triangles */

1214

case 9:

1215

switch (n2) {

1216

case 0: /* single type mesh */

1217

case 3: /* quads */

1218

case 10:

1219

break;

1220

default:

1221

SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Unsupported cell types %d and %d",n1, n2)return PetscError(((MPI_Comm)0x44000001),1221,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Unsupported cell types %d and %d",n1
,n2);

1222

}

1223

break;

1224

case 3: /* quadrilateral */

1225

case 10:

1226

switch (n2) {

1227

case 0: /* single type mesh */

1228

case 2: /* swap since we list simplices first */

1229

case 9:

1230

tn = hc1;

1231

hc1 = hc2;

1232

hc2 = tn;

1233

1234

tp = hybridCells1;

1235

hybridCells1 = hybridCells2;

1236

hybridCells2 = tp;

1237

break;

1238

default:

1239

SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Unsupported cell types %d and %d",n1, n2)return PetscError(((MPI_Comm)0x44000001),1239,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Unsupported cell types %d and %d",n1
,n2);

1240

}

1241

break;

1242

case 4: /* tetrahedra */

1243

case 11:

1244

switch (n2) {

1245

case 0: /* single type mesh */

1246

case 6: /* wedges */

1247

case 13:

1248

break;

1249

default:

1250

SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Unsupported cell types %d and %d",n1, n2)return PetscError(((MPI_Comm)0x44000001),1250,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Unsupported cell types %d and %d",n1
,n2);

1251

}

1252

break;

1253

case 5: /* hexahedra */

1254

case 12:

1255

switch (n2) {

1256

case 0: /* single type mesh */

1257

case 6: /* wedges */

1258

case 13:

1259

break;

1260

default:

1261

SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Unsupported cell types %d and %d",n1, n2)return PetscError(((MPI_Comm)0x44000001),1261,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Unsupported cell types %d and %d",n1
,n2);

1262

}

1263

break;

1264

case 6: /* wedge */

1265

case 13:

1266

switch (n2) {

1267

case 0: /* single type mesh */

1268

case 4: /* tetrahedra: swap since we list simplices first */

1269

case 11:

1270

case 5: /* hexahedra */

1271

case 12:

1272

tn = hc1;

1273

hc1 = hc2;

1274

hc2 = tn;

1275

1276

tp = hybridCells1;

1277

hybridCells1 = hybridCells2;

1278

hybridCells2 = tp;

1279

break;

1280

default:

1281

SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Unsupported cell types %d and %d",n1, n2)return PetscError(((MPI_Comm)0x44000001),1281,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Unsupported cell types %d and %d",n1
,n2);

1282

}

1283

break;

1284

default:

1285

SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Unsupported cell types %d and %d",n1, n2)return PetscError(((MPI_Comm)0x44000001),1285,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Unsupported cell types %d and %d",n1
,n2);

1286

}

1287

cMax = hc1;

1288

ierr = PetscMalloc1(trueNumCells, &hybridMap)PetscMallocA(1,PETSC_FALSE,1288,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,(size_t)(trueNumCells)*sizeof(**(&hybridMap)),(&hybridMap
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1288,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1289

for (cell = 0; cell < hc1; cell++) hybridMap[hybridCells1[cell]] = cell;

1290

for (cell = 0; cell < hc2; cell++) hybridMap[hybridCells2[cell]] = cell + hc1;

1291

ierr = PetscFree(hybridCells1)((*PetscTrFree)((void*)(hybridCells1),1291,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((hybridCells1) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1291,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1292

ierr = PetscFree(hybridCells2)((*PetscTrFree)((void*)(hybridCells2),1292,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((hybridCells2) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1292,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1293

}

1294

1295

}

1296

1297

/* Allocate the cell-vertex mesh */

1298

ierr = DMPlexSetChart(*dm, 0, trueNumCells+numVertices);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1298,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1299

for (cell = 0, c = 0; c < numCells; ++c) {

1300

if (gmsh_elem[c].dim == dim) {

1301

ierr = DMPlexSetConeSize(*dm, hybridMap ? hybridMap[cell] : cell, gmsh_elem[c].numNodes);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1301,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1302

cell++;

1303

}

1304

}

1305

ierr = DMSetUp(*dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1305,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1306

/* Add cell-vertex connections */

1307

for (cell = 0, c = 0; c < numCells; ++c) {

1308

if (gmsh_elem[c].dim == dim) {

1309

PetscInt pcone[8], corner;

1310

for (corner = 0; corner < gmsh_elem[c].numNodes; ++corner) {

1311

const PetscInt cc = gmsh_elem[c].nodes[corner] - shift;

1312

pcone[corner] = (periodicMap ? periodicMap[cc] : cc) + trueNumCells;

1313

}

1314

if (dim == 3) {

1315

/* Tetrahedra are inverted */

1316

if (gmsh_elem[c].cellType == 4 || gmsh_elem[c].cellType == 11) {

1317

PetscInt tmp = pcone[0];

1318

pcone[0] = pcone[1];

1319

pcone[1] = tmp;

1320

}

1321

/* Hexahedra are inverted */

1322

if (gmsh_elem[c].cellType == 5 || gmsh_elem[c].cellType == 12) {

1323

PetscInt tmp = pcone[1];

1324

pcone[1] = pcone[3];

1325

pcone[3] = tmp;

1326

}

1327

/* Prisms are inverted */

1328

if (gmsh_elem[c].cellType == 6 || gmsh_elem[c].cellType == 13) {

1329

PetscInt tmp;

1330

1331

tmp = pcone[1];

1332

pcone[1] = pcone[2];

1333

pcone[2] = tmp;

1334

tmp = pcone[4];

1335

pcone[4] = pcone[5];

1336

pcone[5] = tmp;

1337

}

1338

} else if (dim == 2 && hybridMap && hybridMap[cell] >= cMax) { /* hybrid cells */

1339

/* quads are input to PLEX as prisms */

1340

if (gmsh_elem[c].cellType == 3 || gmsh_elem[c].cellType == 10) {

1341

PetscInt tmp = pcone[2];

1342

pcone[2] = pcone[3];

1343

pcone[3] = tmp;

1344

}

1345

}

1346

ierr = DMPlexSetCone(*dm, hybridMap ? hybridMap[cell] : cell, pcone);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1346,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1347

cell++;

1348

}

1349

}

1350

ierr = MPI_Bcast(&dim, 1, MPIU_INT, 0, comm)((petsc_allreduce_ct += PetscMPIParallelComm((comm)),0) || MPI_Bcast
((&dim),(1),(((MPI_Datatype)0x4c000405)),(0),(comm)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1350,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1351

ierr = DMSetDimension(*dm, dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1351,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1352

ierr = DMPlexSetHybridBounds(*dm, cMax, PETSC_DETERMINE-1, PETSC_DETERMINE-1, PETSC_DETERMINE-1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1352,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1353

ierr = DMPlexSymmetrize(*dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1353,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1354

ierr = DMPlexStratify(*dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1354,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1355

if (interpolate) {

1356

DM idm;

1357

1358

ierr = DMPlexInterpolate(*dm, &idm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1358,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1359

ierr = DMDestroy(dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1359,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1360

*dm = idm;

1361

}

1362

1363

if (usemarker && !interpolate && dim > 1) SETERRQ(comm,PETSC_ERR_SUP,"Cannot create marker label without interpolation")return PetscError(comm,1363,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,56,PETSC_ERROR_INITIAL,"Cannot create marker label without interpolation"
);

1364

if (!rank && usemarker) {

1365

PetscInt f, fStart, fEnd;

1366

1367

ierr = DMCreateLabel(*dm, "marker");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1367,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1368

ierr = DMPlexGetHeightStratum(*dm, 1, &fStart, &fEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1368,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1369

for (f = fStart; f < fEnd; ++f) {

1370

PetscInt suppSize;

1371

1372

ierr = DMPlexGetSupportSize(*dm, f, &suppSize);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1372,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1373

if (suppSize == 1) {

1374

PetscInt *cone = NULL((void*)0), coneSize, p;

1375

1376

ierr = DMPlexGetTransitiveClosure(*dm, f, PETSC_TRUE, &coneSize, &cone);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1376,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1377

for (p = 0; p < coneSize; p += 2) {

1378

ierr = DMSetLabelValue(*dm, "marker", cone[p], 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1378,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1379

}

1380

ierr = DMPlexRestoreTransitiveClosure(*dm, f, PETSC_TRUE, &coneSize, &cone);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1380,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1381

}

1382

}

1383

}

1384

1385

if (!rank) {

1386

PetscInt vStart, vEnd;

1387

1388

ierr = DMPlexGetDepthStratum(*dm, 0, &vStart, &vEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1388,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1389

for (cell = 0, c = 0; c < numCells; ++c) {

1390

1391

/* Create face sets */

1392

if (interpolate && gmsh_elem[c].dim == dim-1) {

1393

const PetscInt *join;

1394

PetscInt joinSize, pcone[8], corner;

1395

/* Find the relevant facet with vertex joins */

1396

for (corner = 0; corner < gmsh_elem[c].numNodes; ++corner) {

1397

const PetscInt cc = gmsh_elem[c].nodes[corner] - shift;

1398

pcone[corner] = (periodicMap ? periodicMap[cc] : cc) + vStart;

1399

}

1400

ierr = DMPlexGetFullJoin(*dm, gmsh_elem[c].numNodes, pcone, &joinSize, &join);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1400,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1401

if (joinSize != 1) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Could not determine Plex facet for GMsh element %d (Plex cell %D)", gmsh_elem[c].id, c)return PetscError(((MPI_Comm)0x44000001),1401,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,62,PETSC_ERROR_INITIAL,"Could not determine Plex facet for GMsh element %d (Plex cell %D)"
,gmsh_elem[c].id,c);

1402

ierr = DMSetLabelValue(*dm, "Face Sets", join[0], gmsh_elem[c].tags[0]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1402,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1403

ierr = DMPlexRestoreJoin(*dm, gmsh_elem[c].numNodes, (const PetscInt *) pcone, &joinSize, &join);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1403,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1404

}

1405

1406

/* Create cell sets */

1407

if (gmsh_elem[c].dim == dim) {

1408

if (gmsh_elem[c].numTags > 0) {

1409

ierr = DMSetLabelValue(*dm, "Cell Sets", hybridMap ? hybridMap[cell] : cell, gmsh_elem[c].tags[0]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1409,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1410

}

1411

cell++;

1412

}

1413

1414

/* Create vertex sets */

1415

if (gmsh_elem[c].dim == 0) {

1416

if (gmsh_elem[c].numTags > 0) {

1417

const PetscInt cc = gmsh_elem[c].nodes[0] - shift;

1418

const PetscInt vid = (periodicMap ? periodicMap[cc] : cc) + vStart;

1419

ierr = DMSetLabelValue(*dm, "Vertex Sets", vid, gmsh_elem[c].tags[0]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1419,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1420

}

1421

}

1422

}

1423

}

1424

1425

/* Create coordinates */

1426

if (embedDim < 0) embedDim = dim;

1427

ierr = DMSetCoordinateDim(*dm, embedDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1427,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1428

ierr = DMGetCoordinateSection(*dm, &coordSection);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1428,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1429

ierr = PetscSectionSetNumFields(coordSection, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1429,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1430

ierr = PetscSectionSetFieldComponents(coordSection, 0, embedDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1430,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1431

if (periodicMap) { /* we need to localize coordinates on cells */

1432

ierr = PetscSectionSetChart(coordSection, 0, trueNumCells + numVertices);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1432,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1433

} else {

1434

ierr = PetscSectionSetChart(coordSection, trueNumCells, trueNumCells + numVertices);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1434,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1435

}

1436

for (v = trueNumCells; v < trueNumCells+numVertices; ++v) {

1437

ierr = PetscSectionSetDof(coordSection, v, embedDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1437,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1438

ierr = PetscSectionSetFieldDof(coordSection, v, 0, embedDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1438,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1439

}

1440

if (periodicMap) {

1441

ierr = PetscBTCreate(trueNumCells, &periodicC);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1441,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1442

for (cell = 0, c = 0; c < numCells; ++c) {

1443

if (gmsh_elem[c].dim == dim) {

1444

PetscInt corner;

1445

PetscBool pc = PETSC_FALSE;

1446

for (corner = 0; corner < gmsh_elem[c].numNodes; ++corner) {

1447

pc = (PetscBool)(pc || PetscBTLookup(periodicV, gmsh_elem[c].nodes[corner] - shift));

1448

}

1449

if (pc) {

1450

PetscInt dof = gmsh_elem[c].numNodes*embedDim;

1451

PetscInt ucell = hybridMap ? hybridMap[cell] : cell;

1452

ierr = PetscBTSet(periodicC, ucell);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1452,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1453

ierr = PetscSectionSetDof(coordSection, ucell, dof);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1453,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1454

ierr = PetscSectionSetFieldDof(coordSection, ucell, 0, dof);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1454,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1455

}

1456

cell++;

1457

}

1458

}

1459

}

1460

ierr = PetscSectionSetUp(coordSection);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1460,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1461

ierr = PetscSectionGetStorageSize(coordSection, &coordSize);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1461,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1462

ierr = VecCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001), &coordinates);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1462,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1463

ierr = PetscObjectSetName((PetscObject) coordinates, "coordinates");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1463,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1464

ierr = VecSetSizes(coordinates, coordSize, PETSC_DETERMINE-1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1464,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1465

ierr = VecSetBlockSize(coordinates, embedDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1465,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1466

ierr = VecSetType(coordinates, VECSTANDARD"standard");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1466,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1467

ierr = VecGetArray(coordinates, &coords);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1467,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1468

if (periodicMap) {

1469

PetscInt off;

1470

1471

for (cell = 0, c = 0; c < numCells; ++c) {

1472

PetscInt pcone[8], corner;

1473

if (gmsh_elem[c].dim == dim) {

1474

PetscInt ucell = hybridMap ? hybridMap[cell] : cell;

1475

if (PetscUnlikely(PetscBTLookup(periodicC, ucell))__builtin_expect(!!(PetscBTLookup(periodicC, ucell)),0)) {

1476

for (corner = 0; corner < gmsh_elem[c].numNodes; ++corner) {

1477

pcone[corner] = gmsh_elem[c].nodes[corner] - shift;

1478

}

1479

if (dim == 3) {

1480

/* Tetrahedra are inverted */

1481

if (gmsh_elem[c].cellType == 4 || gmsh_elem[c].cellType == 11) {

1482

PetscInt tmp = pcone[0];

1483

pcone[0] = pcone[1];

1484

pcone[1] = tmp;

1485

}

1486

/* Hexahedra are inverted */

1487

if (gmsh_elem[c].cellType == 5 || gmsh_elem[c].cellType == 12) {

1488

PetscInt tmp = pcone[1];

1489

pcone[1] = pcone[3];

1490

pcone[3] = tmp;

1491

}

1492

/* Prisms are inverted */

1493

if (gmsh_elem[c].cellType == 6 || gmsh_elem[c].cellType == 13) {

1494

PetscInt tmp;

1495

1496

tmp = pcone[1];

1497

pcone[1] = pcone[2];

1498

pcone[2] = tmp;

1499

tmp = pcone[4];

1500

pcone[4] = pcone[5];

1501

pcone[5] = tmp;

1502

}

1503

} else if (dim == 2 && hybridMap && hybridMap[cell] >= cMax) { /* hybrid cells */

1504

/* quads are input to PLEX as prisms */

1505

if (gmsh_elem[c].cellType == 3 || gmsh_elem[c].cellType == 10) {

1506

PetscInt tmp = pcone[2];

1507

pcone[2] = pcone[3];

1508

pcone[3] = tmp;

1509

}

1510

}

1511

ierr = PetscSectionGetOffset(coordSection, ucell, &off);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1511,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1512

for (corner = 0; corner < gmsh_elem[c].numNodes; ++corner) {

1513

v = pcone[corner];

1514

for (d = 0; d < embedDim; ++d) {

1515

coords[off++] = (PetscReal) coordsIn[v*3+d];

1516

}

1517

}

1518

}

1519

cell++;

1520

}

1521

}

1522

for (v = 0; v < numVertices; ++v) {

1523

ierr = PetscSectionGetOffset(coordSection, v + trueNumCells, &off);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1523,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1524

for (d = 0; d < embedDim; ++d) {

1525

coords[off+d] = (PetscReal) coordsIn[periodicMapI[v]*3+d];

1526

}

1527

}

1528

} else {

1529

for (v = 0; v < numVertices; ++v) {

1530

for (d = 0; d < embedDim; ++d) {

1531

coords[v*embedDim+d] = (PetscReal) coordsIn[v*3+d];

1532

}

1533

}

1534

}

1535

ierr = VecRestoreArray(coordinates, &coords);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1535,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1536

ierr = DMSetCoordinatesLocal(*dm, coordinates);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1536,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1537

1538

periodic = periodicMap ? PETSC_TRUE : PETSC_FALSE;

1539

ierr = MPI_Bcast(&periodic, 1, MPIU_BOOL, 0, comm)((petsc_allreduce_ct += PetscMPIParallelComm((comm)),0) || MPI_Bcast
((&periodic),(1),(MPIU_BOOL),(0),(comm)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1539,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1540

ierr = DMSetPeriodicity(*dm, periodic, NULL((void*)0), NULL((void*)0), NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1540,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1541

1542

ierr = PetscFree(coordsIn)((*PetscTrFree)((void*)(coordsIn),1542,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((coordsIn) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1542,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1543

ierr = PetscFree(gmsh_elem)((*PetscTrFree)((void*)(gmsh_elem),1543,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((gmsh_elem) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1543,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1544

ierr = PetscFree(hybridMap)((*PetscTrFree)((void*)(hybridMap),1544,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((hybridMap) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1544,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1545

ierr = PetscFree(periodicMap)((*PetscTrFree)((void*)(periodicMap),1545,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((periodicMap) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1545,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1546

ierr = PetscFree(periodicMapI)((*PetscTrFree)((void*)(periodicMapI),1546,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
) || ((periodicMapI) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1546,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1547

ierr = PetscBTDestroy(&periodicV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1547,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1548

ierr = PetscBTDestroy(&periodicC);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1548,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1549

ierr = VecDestroy(&coordinates);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1549,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1550

1551

ierr = PetscLogEventEnd(DMPLEX_CreateGmsh,*dm,0,0,0)(((PetscLogPLE && petsc_stageLog->stageInfo[petsc_stageLog
->curStage].perfInfo.active && petsc_stageLog->
stageInfo[petsc_stageLog->curStage].eventLog->eventInfo
[DMPLEX_CreateGmsh].active) ? (*PetscLogPLE)((DMPLEX_CreateGmsh
),0,(PetscObject)(*dm),(PetscObject)(0),(PetscObject)(0),(PetscObject
)(0)) : 0 ));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1551,__func__,"/sandbox/petsc/petsc.master/src/dm/impls/plex/plexgmsh.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1552

1553

}