Lines Matching refs:genCtx
27 static PetscInt DMMoab_SetTensorElementConnectivity_Private(DMMoabMeshGeneratorCtx &genCtx, PetscIn… in DMMoab_SetTensorElementConnectivity_Private() argument
29 switch (genCtx.dim) { in DMMoab_SetTensorElementConnectivity_Private()
41 connectivity[offset + subent_conn[2]] = corner + 1 + genCtx.ystride; in DMMoab_SetTensorElementConnectivity_Private()
42 connectivity[offset + subent_conn[3]] = corner + genCtx.ystride; in DMMoab_SetTensorElementConnectivity_Private()
50 connectivity[offset + subent_conn[2]] = corner + 1 + genCtx.ystride; in DMMoab_SetTensorElementConnectivity_Private()
51 connectivity[offset + subent_conn[3]] = corner + genCtx.ystride; in DMMoab_SetTensorElementConnectivity_Private()
52 connectivity[offset + subent_conn[4]] = corner + genCtx.zstride; in DMMoab_SetTensorElementConnectivity_Private()
53 connectivity[offset + subent_conn[5]] = corner + 1 + genCtx.zstride; in DMMoab_SetTensorElementConnectivity_Private()
54 connectivity[offset + subent_conn[6]] = corner + 1 + genCtx.ystride + genCtx.zstride; in DMMoab_SetTensorElementConnectivity_Private()
55 connectivity[offset + subent_conn[7]] = corner + genCtx.ystride + genCtx.zstride; in DMMoab_SetTensorElementConnectivity_Private()
61 static PetscInt DMMoab_SetSimplexElementConnectivity_Private(DMMoabMeshGeneratorCtx &genCtx, PetscI… in DMMoab_SetSimplexElementConnectivity_Private() argument
67 switch (genCtx.dim) { in DMMoab_SetSimplexElementConnectivity_Private()
75 C = corner + 1 + genCtx.ystride; in DMMoab_SetSimplexElementConnectivity_Private()
76 D = corner + genCtx.ystride; in DMMoab_SetSimplexElementConnectivity_Private()
127 C = corner + 1 + genCtx.ystride; in DMMoab_SetSimplexElementConnectivity_Private()
128 D = corner + genCtx.ystride; in DMMoab_SetSimplexElementConnectivity_Private()
129 E = corner + genCtx.zstride; in DMMoab_SetSimplexElementConnectivity_Private()
130 F = corner + 1 + genCtx.zstride; in DMMoab_SetSimplexElementConnectivity_Private()
131 G = corner + 1 + genCtx.ystride + genCtx.zstride; in DMMoab_SetSimplexElementConnectivity_Private()
132 H = corner + genCtx.ystride + genCtx.zstride; in DMMoab_SetSimplexElementConnectivity_Private()
178 … PetscInt> DMMoab_SetElementConnectivity_Private(DMMoabMeshGeneratorCtx &genCtx, PetscInt offset, … in DMMoab_SetElementConnectivity_Private() argument
185 if (genCtx.simplex) { in DMMoab_SetElementConnectivity_Private()
186 subelem = simplices_per_tensor[genCtx.dim]; in DMMoab_SetElementConnectivity_Private()
188 …vcount = DMMoab_SetSimplexElementConnectivity_Private(genCtx, m, offset, corner, subent_conn, conn… in DMMoab_SetElementConnectivity_Private()
193 …vcount = DMMoab_SetTensorElementConnectivity_Private(genCtx, offset, corner, subent_conn, connect… in DMMoab_SetElementConnectivity_Private()
198 …::Interface *mbImpl, moab::ReadUtilIface *iface, DMMoabMeshGeneratorCtx &genCtx, PetscInt m, Petsc… in DMMoab_GenerateVertices_Private() argument
211 …nnodes = genCtx.blockSizeVertexXYZ[0] * (genCtx.dim > 1 ? genCtx.blockSizeVertexXYZ[1] * (genCtx.d… in DMMoab_GenerateVertices_Private()
222 x = (m * genCtx.A + a) * genCtx.q; in DMMoab_GenerateVertices_Private()
223 y = (n * genCtx.B + b) * genCtx.q; in DMMoab_GenerateVertices_Private()
224 z = (k * genCtx.C + c) * genCtx.q; in DMMoab_GenerateVertices_Private()
228 for (kk = 0; kk < (genCtx.dim > 2 ? genCtx.blockSizeVertexXYZ[2] : 1); kk++) { in DMMoab_GenerateVertices_Private()
229 for (jj = 0; jj < (genCtx.dim > 1 ? genCtx.blockSizeVertexXYZ[1] : 1); jj++) { in DMMoab_GenerateVertices_Private()
230 for (ii = 0; ii < genCtx.blockSizeVertexXYZ[0]; ii++, ix++) { in DMMoab_GenerateVertices_Private()
232 arrays[0][ix] = (x + ii) * genCtx.dx + genCtx.xyzbounds[0]; in DMMoab_GenerateVertices_Private()
233 arrays[1][ix] = (y + jj) * genCtx.dy + genCtx.xyzbounds[2]; in DMMoab_GenerateVertices_Private()
234 arrays[2][ix] = (z + kk) * genCtx.dz + genCtx.xyzbounds[4]; in DMMoab_GenerateVertices_Private()
241 gids[ix] = 1 + (x + ii) + (y + jj) * genCtx.NX + (z + kk) * (genCtx.NX * genCtx.NY); in DMMoab_GenerateVertices_Private()
252 …::Interface *mbImpl, moab::ReadUtilIface *iface, DMMoabMeshGeneratorCtx &genCtx, PetscInt m, Petsc… in DMMoab_GenerateElements_Private() argument
258 …nsorelems = genCtx.blockSizeElementXYZ[0] * (genCtx.dim > 1 ? genCtx.blockSizeElementXY… in DMMoab_GenerateElements_Private()
264 switch (genCtx.dim) { in DMMoab_GenerateElements_Private()
271 if (genCtx.simplex) { in DMMoab_GenerateElements_Private()
273 nelems = ntensorelems * simplices_per_tensor[genCtx.dim]; in DMMoab_GenerateElements_Private()
284 if (genCtx.simplex) { in DMMoab_GenerateElements_Private()
286 nelems = ntensorelems * simplices_per_tensor[genCtx.dim]; in DMMoab_GenerateElements_Private()
304 xe = m * genCtx.A * genCtx.blockSizeElementXYZ[0] + a * genCtx.blockSizeElementXYZ[0]; in DMMoab_GenerateElements_Private()
305 …ye = (genCtx.dim > 1 ? n * genCtx.B * genCtx.blockSizeElementXYZ[1] + b * genCtx.blockSizeElementX… in DMMoab_GenerateElements_Private()
306 …ze = (genCtx.dim > 2 ? k * genCtx.C * genCtx.blockSizeElementXYZ[2] + c * genCtx.blockSizeElementX… in DMMoab_GenerateElements_Private()
309 for (kk = 0; kk < (genCtx.dim > 2 ? genCtx.blockSizeElementXYZ[2] : 1); kk++) { in DMMoab_GenerateElements_Private()
310 for (jj = 0; jj < (genCtx.dim > 1 ? genCtx.blockSizeElementXYZ[1] : 1); jj++) { in DMMoab_GenerateElements_Private()
311 for (ii = 0; ii < genCtx.blockSizeElementXYZ[0]; ii++) { in DMMoab_GenerateElements_Private()
312 …moab::EntityHandle corner = startv + genCtx.q * ii + genCtx.q * jj * genCtx.ystride + genCtx.q * k… in DMMoab_GenerateElements_Private()
314 …std::pair<PetscInt, PetscInt> entoffset = DMMoab_SetElementConnectivity_Private(genCtx, ix, corner… in DMMoab_GenerateElements_Private()
320 … gids[ie + j] = 1 + ((xe + ii) + (ye + jj) * genCtx.nex + (ze + kk) * (genCtx.nex * genCtx.ney)); in DMMoab_GenerateElements_Private()
331 if (genCtx.adjEnts) { /* we need to update adjacencies now, because some elements are new */ in DMMoab_GenerateElements_Private()
341 static PetscErrorCode DMMBUtil_InitializeOptions(DMMoabMeshGeneratorCtx &genCtx, PetscInt dim, Pets… in DMMBUtil_InitializeOptions() argument
345 genCtx.dim = dim; in DMMBUtil_InitializeOptions()
346 genCtx.simplex = simplex; in DMMBUtil_InitializeOptions()
347 genCtx.newMergeMethod = genCtx.keep_skins = genCtx.adjEnts = true; in DMMBUtil_InitializeOptions()
349 genCtx.q = 1; in DMMBUtil_InitializeOptions()
350 genCtx.fraction = genCtx.remainder = genCtx.cumfraction = 0; in DMMBUtil_InitializeOptions()
352 …if (!genCtx.usrxyzgrid) { /* not overridden by genCtx - assume nele equally and that genCtx wants … in DMMBUtil_InitializeOptions()
354 genCtx.fraction = nelems / nprocs; /* partition only by the largest dimension */ in DMMBUtil_InitializeOptions()
355 …genCtx.remainder = nelems % nprocs; /* remainder after partition which gets evenly distributed b… in DMMBUtil_InitializeOptions()
356 …genCtx.cumfraction = (rank > 0 ? (genCtx.fraction) * (rank) + (rank - 1 < genCtx.remainder ? rank … in DMMBUtil_InitializeOptions()
357 if (rank < genCtx.remainder) /* This process gets "fraction+1" elements */ in DMMBUtil_InitializeOptions()
358 genCtx.fraction++; in DMMBUtil_InitializeOptions()
360 …cInt_FMT ", Cumulative fraction = %" PetscInt_FMT "\n", genCtx.fraction, genCtx.remainder, genCtx.… in DMMBUtil_InitializeOptions()
361 switch (genCtx.dim) { in DMMBUtil_InitializeOptions()
363 genCtx.blockSizeElementXYZ[0] = genCtx.fraction; in DMMBUtil_InitializeOptions()
364 genCtx.blockSizeElementXYZ[1] = 1; in DMMBUtil_InitializeOptions()
365 genCtx.blockSizeElementXYZ[2] = 1; in DMMBUtil_InitializeOptions()
368 genCtx.blockSizeElementXYZ[0] = nelems; in DMMBUtil_InitializeOptions()
369 genCtx.blockSizeElementXYZ[1] = genCtx.fraction; in DMMBUtil_InitializeOptions()
370 genCtx.blockSizeElementXYZ[2] = 1; in DMMBUtil_InitializeOptions()
374 genCtx.blockSizeElementXYZ[0] = nelems; in DMMBUtil_InitializeOptions()
375 genCtx.blockSizeElementXYZ[1] = nelems; in DMMBUtil_InitializeOptions()
376 genCtx.blockSizeElementXYZ[2] = genCtx.fraction; in DMMBUtil_InitializeOptions()
384 for (PetscInt i = 0; i < 6; i++) genCtx.xyzbounds[i] = bounds[i]; in DMMBUtil_InitializeOptions()
386 genCtx.xyzbounds[0] = genCtx.xyzbounds[2] = genCtx.xyzbounds[4] = 0.0; in DMMBUtil_InitializeOptions()
387 genCtx.xyzbounds[1] = genCtx.xyzbounds[3] = genCtx.xyzbounds[5] = 1.0; in DMMBUtil_InitializeOptions()
390 if (!genCtx.usrprocgrid) { in DMMBUtil_InitializeOptions()
391 switch (genCtx.dim) { in DMMBUtil_InitializeOptions()
393 genCtx.M = nprocs; in DMMBUtil_InitializeOptions()
394 genCtx.N = genCtx.K = 1; in DMMBUtil_InitializeOptions()
397 genCtx.N = nprocs; in DMMBUtil_InitializeOptions()
398 genCtx.M = genCtx.K = 1; in DMMBUtil_InitializeOptions()
401 genCtx.K = nprocs; in DMMBUtil_InitializeOptions()
402 genCtx.M = genCtx.N = 1; in DMMBUtil_InitializeOptions()
407 if (!genCtx.usrrefgrid) genCtx.A = genCtx.B = genCtx.C = 1; in DMMBUtil_InitializeOptions()
410 genCtx.nex = genCtx.ney = genCtx.nez = 0; in DMMBUtil_InitializeOptions()
411 genCtx.xstride = genCtx.ystride = genCtx.zstride = 0; in DMMBUtil_InitializeOptions()
412 genCtx.NX = genCtx.NY = genCtx.NZ = 0; in DMMBUtil_InitializeOptions()
413 genCtx.nex = genCtx.ney = genCtx.nez = 0; in DMMBUtil_InitializeOptions()
414 genCtx.blockSizeVertexXYZ[0] = genCtx.blockSizeVertexXYZ[1] = genCtx.blockSizeVertexXYZ[2] = 1; in DMMBUtil_InitializeOptions()
416 switch (genCtx.dim) { in DMMBUtil_InitializeOptions()
418 genCtx.blockSizeVertexXYZ[0] = genCtx.q * genCtx.blockSizeElementXYZ[0] + 1; in DMMBUtil_InitializeOptions()
419 genCtx.blockSizeVertexXYZ[1] = genCtx.q * genCtx.blockSizeElementXYZ[1] + 1; in DMMBUtil_InitializeOptions()
420 genCtx.blockSizeVertexXYZ[2] = genCtx.q * genCtx.blockSizeElementXYZ[2] + 1; in DMMBUtil_InitializeOptions()
422 …genCtx.nex = genCtx.M * genCtx.A * genCtx.blockSizeElementXYZ[0]; /* number of e… in DMMBUtil_InitializeOptions()
423 …genCtx.dx = (genCtx.xyzbounds[1] - genCtx.xyzbounds[0]) / (nelems * genCtx.q); /* distance be… in DMMBUtil_InitializeOptions()
424 genCtx.NX = (genCtx.q * genCtx.nex + 1); in DMMBUtil_InitializeOptions()
425 genCtx.xstride = 1; in DMMBUtil_InitializeOptions()
426 …genCtx.ney = genCtx.N * genCtx.B * genCtx.blockSizeElementXYZ[1]; /* number of e… in DMMBUtil_InitializeOptions()
427 …genCtx.dy = (genCtx.xyzbounds[3] - genCtx.xyzbounds[2]) / (nelems * genCtx.q); /* distance be… in DMMBUtil_InitializeOptions()
428 genCtx.NY = (genCtx.q * genCtx.ney + 1); in DMMBUtil_InitializeOptions()
429 genCtx.ystride = genCtx.blockSizeVertexXYZ[0]; in DMMBUtil_InitializeOptions()
430 …genCtx.nez = genCtx.K * genCtx.C * genCtx.blockSizeElementXYZ[2]; /* number of e… in DMMBUtil_InitializeOptions()
431 …genCtx.dz = (genCtx.xyzbounds[5] - genCtx.xyzbounds[4]) / (nelems * genCtx.q); /* distance be… in DMMBUtil_InitializeOptions()
432 genCtx.NZ = (genCtx.q * genCtx.nez + 1); in DMMBUtil_InitializeOptions()
433 genCtx.zstride = genCtx.blockSizeVertexXYZ[0] * genCtx.blockSizeVertexXYZ[1]; in DMMBUtil_InitializeOptions()
436 genCtx.blockSizeVertexXYZ[0] = genCtx.q * genCtx.blockSizeElementXYZ[0] + 1; in DMMBUtil_InitializeOptions()
437 genCtx.blockSizeVertexXYZ[1] = genCtx.q * genCtx.blockSizeElementXYZ[1] + 1; in DMMBUtil_InitializeOptions()
438 genCtx.blockSizeVertexXYZ[2] = 0; in DMMBUtil_InitializeOptions()
440 …genCtx.nex = genCtx.M * genCtx.A * genCtx.blockSizeElementXYZ[0]; /* number … in DMMBUtil_InitializeOptions()
441 …genCtx.dx = (genCtx.xyzbounds[1] - genCtx.xyzbounds[0]) / (genCtx.nex * genCtx.q); /* distanc… in DMMBUtil_InitializeOptions()
442 genCtx.NX = (genCtx.q * genCtx.nex + 1); in DMMBUtil_InitializeOptions()
443 genCtx.xstride = 1; in DMMBUtil_InitializeOptions()
444 …genCtx.ney = genCtx.N * genCtx.B * genCtx.blockSizeElementXYZ[1]; /* number of e… in DMMBUtil_InitializeOptions()
445 …genCtx.dy = (genCtx.xyzbounds[3] - genCtx.xyzbounds[2]) / (nelems * genCtx.q); /* distance be… in DMMBUtil_InitializeOptions()
446 genCtx.NY = (genCtx.q * genCtx.ney + 1); in DMMBUtil_InitializeOptions()
447 genCtx.ystride = genCtx.blockSizeVertexXYZ[0]; in DMMBUtil_InitializeOptions()
450 genCtx.blockSizeVertexXYZ[1] = genCtx.blockSizeVertexXYZ[2] = 0; in DMMBUtil_InitializeOptions()
451 …genCtx.blockSizeVertexXYZ[0] = genCtx.q * genCtx.blockSizeElementXY… in DMMBUtil_InitializeOptions()
453 …genCtx.nex = genCtx.M * genCtx.A * genCtx.blockSizeElementXYZ[0]; /* number of e… in DMMBUtil_InitializeOptions()
454 …genCtx.dx = (genCtx.xyzbounds[1] - genCtx.xyzbounds[0]) / (nelems * genCtx.q); /* distance be… in DMMBUtil_InitializeOptions()
455 genCtx.NX = (genCtx.q * genCtx.nex + 1); in DMMBUtil_InitializeOptions()
456 genCtx.xstride = 1; in DMMBUtil_InitializeOptions()
461 …k(genCtx.dim >= 1 && genCtx.dim <= 3, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Invalid topolog… in DMMBUtil_InitializeOptions()
462 …genCtx.M * genCtx.N * genCtx.K == nprocs, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Invalid [m,… in DMMBUtil_InitializeOptions()
463 genCtx.N, genCtx.K, nprocs); in DMMBUtil_InitializeOptions()
465 …genCtx.xyzbounds[0] < genCtx.xyzbounds[1], PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "X-dim: Lef… in DMMBUtil_InitializeOptions()
466 …genCtx.dim <= 1 || genCtx.xyzbounds[2] < genCtx.xyzbounds[3], PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTO… in DMMBUtil_InitializeOptions()
467 …genCtx.dim <= 2 || genCtx.xyzbounds[4] < genCtx.xyzbounds[5], PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTO… in DMMBUtil_InitializeOptions()
469 … %" PetscInt_FMT ", %" PetscInt_FMT "\n", genCtx.blockSizeElementXYZ[0], genCtx.blockSizeElementXY… in DMMBUtil_InitializeOptions()
470 …, %" PetscInt_FMT ", %" PetscInt_FMT "\n", genCtx.blockSizeVertexXYZ[0], genCtx.blockSizeVertexXYZ… in DMMBUtil_InitializeOptions()
471 …sors := %" PetscInt_FMT ", %" PetscInt_FMT ", %" PetscInt_FMT "\n", genCtx.A, genCtx.B, genCtx.C)); in DMMBUtil_InitializeOptions()
472 …sors := %" PetscInt_FMT ", %" PetscInt_FMT ", %" PetscInt_FMT "\n", genCtx.M, genCtx.N, genCtx.K)); in DMMBUtil_InitializeOptions()
473 …= %" PetscInt_FMT ", %" PetscInt_FMT ", %" PetscInt_FMT "\n", genCtx.nex, genCtx.ney, genCtx.nez)); in DMMBUtil_InitializeOptions()
474 PetscCall(PetscInfo(NULL, "Local delxyz:= %g, %g, %g\n", genCtx.dx, genCtx.dy, genCtx.dz)); in DMMBUtil_InitializeOptions()
475 …PetscInt_FMT ", %" PetscInt_FMT ", %" PetscInt_FMT "\n", genCtx.xstride, genCtx.ystride, genCtx.zs… in DMMBUtil_InitializeOptions()
510 DMMoabMeshGeneratorCtx genCtx; in DMMoabCreateBoxMesh() local
521 PetscCall(PetscLogEventRegister("GenerateMesh", DM_CLASSID, &genCtx.generateMesh)); in DMMoabCreateBoxMesh()
522 PetscCall(PetscLogEventRegister("AddVertices", DM_CLASSID, &genCtx.generateVertices)); in DMMoabCreateBoxMesh()
523 PetscCall(PetscLogEventRegister("AddElements", DM_CLASSID, &genCtx.generateElements)); in DMMoabCreateBoxMesh()
524 PetscCall(PetscLogEventRegister("ParResolve", DM_CLASSID, &genCtx.parResolve)); in DMMoabCreateBoxMesh()
525 PetscCall(PetscLogEventBegin(genCtx.generateMesh, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
561 genCtx.M = genCtx.N = genCtx.K = 1; in DMMoabCreateBoxMesh()
562 genCtx.A = genCtx.B = genCtx.C = 1; in DMMoabCreateBoxMesh()
563 genCtx.blockSizeElementXYZ[0] = 0; in DMMoabCreateBoxMesh()
564 genCtx.blockSizeElementXYZ[1] = 0; in DMMoabCreateBoxMesh()
565 genCtx.blockSizeElementXYZ[2] = 0; in DMMoabCreateBoxMesh()
569 …oints in x direction", "DMMoabSetSizes", genCtx.blockSizeElementXYZ[0], &genCtx.blockSizeElementXY… in DMMoabCreateBoxMesh()
570 …oints in y direction", "DMMoabSetSizes", genCtx.blockSizeElementXYZ[1], &genCtx.blockSizeElementXY… in DMMoabCreateBoxMesh()
571 …oints in z direction", "DMMoabSetSizes", genCtx.blockSizeElementXYZ[2], &genCtx.blockSizeElementXY… in DMMoabCreateBoxMesh()
574 …, "Number of processors in x direction", "DMMoabSetNumProcs", genCtx.M, &genCtx.M, &genCtx.usrproc… in DMMoabCreateBoxMesh()
575 …, "Number of processors in y direction", "DMMoabSetNumProcs", genCtx.N, &genCtx.N, &genCtx.usrproc… in DMMoabCreateBoxMesh()
576 …, "Number of processors in z direction", "DMMoabSetNumProcs", genCtx.K, &genCtx.K, &genCtx.usrproc… in DMMoabCreateBoxMesh()
579 … of refinement blocks in x direction", "DMMoabSetRefinement", genCtx.A, &genCtx.A, &genCtx.usrrefg… in DMMoabCreateBoxMesh()
580 … of refinement blocks in y direction", "DMMoabSetRefinement", genCtx.B, &genCtx.B, &genCtx.usrrefg… in DMMoabCreateBoxMesh()
581 … of refinement blocks in z direction", "DMMoabSetRefinement", genCtx.C, &genCtx.C, &genCtx.usrrefg… in DMMoabCreateBoxMesh()
584 …PetscCall(DMMBUtil_InitializeOptions(genCtx, dim, useSimplex, global_rank, global_size, bounds, ne… in DMMoabCreateBoxMesh()
588 if (genCtx.adjEnts) genCtx.keep_skins = true; /* do not delete anything - consumes more memory */ in DMMoabCreateBoxMesh()
592 switch (genCtx.dim) { in DMMoabCreateBoxMesh()
594 ml = (genCtx.cumfraction); in DMMoabCreateBoxMesh()
597 nl = (genCtx.cumfraction); in DMMoabCreateBoxMesh()
600 kl = (genCtx.cumfraction) / genCtx.q / genCtx.blockSizeElementXYZ[2] / genCtx.C; //genCtx.K in DMMoabCreateBoxMesh()
634 PetscCall(PetscLogEventEnd(genCtx.generateMesh, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
636 for (a = 0; a < (genCtx.dim > 0 ? genCtx.A : genCtx.A); a++) { in DMMoabCreateBoxMesh()
637 for (b = 0; b < (genCtx.dim > 1 ? genCtx.B : 1); b++) { in DMMoabCreateBoxMesh()
638 for (c = 0; c < (genCtx.dim > 2 ? genCtx.C : 1); c++) { in DMMoabCreateBoxMesh()
641 PetscCall(PetscLogEventBegin(genCtx.generateVertices, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
642 …PetscCall(DMMoab_GenerateVertices_Private(mbImpl, readMeshIface, genCtx, ml, nl, kl, a, b, c, glob… in DMMoabCreateBoxMesh()
643 PetscCall(PetscLogEventEnd(genCtx.generateVertices, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
645 PetscCall(PetscLogEventBegin(genCtx.generateElements, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
646 …PetscCall(DMMoab_GenerateElements_Private(mbImpl, readMeshIface, genCtx, ml, nl, kl, a, b, c, glob… in DMMoabCreateBoxMesh()
647 PetscCall(PetscLogEventEnd(genCtx.generateElements, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
650 switch (genCtx.dim) { in DMMoabCreateBoxMesh()
652 part_num += (c + kl * genCtx.C) * (genCtx.M * genCtx.A * genCtx.N * genCtx.B); in DMMoabCreateBoxMesh()
654 part_num += (b + nl * genCtx.B) * (genCtx.M * genCtx.A); in DMMoabCreateBoxMesh()
656 part_num += (a + ml * genCtx.A); in DMMoabCreateBoxMesh()
672 if (genCtx.adjEnts) { in DMMoabCreateBoxMesh()
673 if (genCtx.dim > 1) { in DMMoabCreateBoxMesh()
679 if (genCtx.dim > 2) { in DMMoabCreateBoxMesh()
710 PetscCall(PetscLogEventBegin(genCtx.parResolve, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
712 merr = mbImpl->get_entities_by_dimension(dmmoab->fileset, genCtx.dim, cells); in DMMoabCreateBoxMesh()
717 if (genCtx.A * genCtx.B * genCtx.C != 1) { // merge needed in DMMoabCreateBoxMesh()
719 if (genCtx.newMergeMethod) { in DMMoabCreateBoxMesh()
749 PetscCall(PetscLogEventEnd(genCtx.parResolve, 0, 0, 0, 0)); in DMMoabCreateBoxMesh()
752 if (!genCtx.keep_skins) { // default is to delete the 1- and 2-dimensional entities in DMMoabCreateBoxMesh()
755 if (genCtx.dim > 1) { in DMMoabCreateBoxMesh()
760 if (genCtx.dim > 2) { in DMMoabCreateBoxMesh()