@@ -42,13 +42,15 @@ namespace
4242//
4343// 1. Available at <https://pi.math.cornell.edu/~hatcher/AT/ATpage.html>.
4444constexpr std::tuple<std::int64_t , std::int64_t , std::int64_t , std::int64_t >
45- num_entities (std::int64_t i, std::int64_t j, std::int64_t k, int nrefine) {
45+ num_entities (std::int64_t i, std::int64_t j, std::int64_t k, int nrefine)
46+ {
4647 i <<= nrefine;
4748 j <<= nrefine;
4849 k <<= nrefine;
4950 std::int64_t vertices = (i + 1 ) * (j + 1 ) * (k + 1 );
50- std::int64_t edges = 7 *i*j*k + 3 *(i*j + i*k + j*k) + (i + j + k);
51- std::int64_t faces = 12 *i*j*k + 2 *(i*j + i*k + j*k);
51+ std::int64_t edges
52+ = 7 * i * j * k + 3 * (i * j + i * k + j * k) + (i + j + k);
53+ std::int64_t faces = 12 * i * j * k + 2 * (i * j + i * k + j * k);
5254 std::int64_t cells = 6 * (i * j * k);
5355 return {vertices, edges, faces, cells};
5456}
@@ -179,7 +181,8 @@ create_cube_mesh(MPI_Comm comm, std::size_t target_dofs, bool target_dofs_total,
179181 else
180182 MPI_Comm_dup (comm, &sub_comm);
181183
182- auto cell_part = dolfinx::mesh::create_cell_partitioner (dolfinx::mesh::GhostMode::none, graph_part, 2 );
184+ auto cell_part = dolfinx::mesh::create_cell_partitioner (
185+ dolfinx::mesh::GhostMode::none, graph_part, 2 );
183186 auto mesh = dolfinx::mesh::create_box (
184187 comm, sub_comm, {{{0.0 , 0.0 , 0.0 }, {1.0 , 1.0 , 1.0 }}}, {Nx, Ny, Nz},
185188 dolfinx::mesh::CellType::tetrahedron, cell_part);
@@ -196,9 +199,10 @@ create_cube_mesh(MPI_Comm comm, std::size_t target_dofs, bool target_dofs_total,
196199 {
197200 mesh.topology_mutable ()->create_connectivity (3 , 1 );
198201 auto [new_mesh, _parent_edges, _parent_facet] = dolfinx::refinement::refine (
199- mesh, std::nullopt ,
200- dolfinx::mesh::create_cell_partitioner (dolfinx::mesh::GhostMode::shared_facet, 2 ),
201- dolfinx::refinement::Option::parent_cell_and_facet);
202+ mesh, std::nullopt ,
203+ dolfinx::mesh::create_cell_partitioner (
204+ dolfinx::mesh::GhostMode::shared_facet, 2 ),
205+ dolfinx::refinement::Option::parent_cell_and_facet);
202206 mesh = std::move (new_mesh);
203207 }
204208
@@ -373,9 +377,10 @@ create_spoke_mesh(MPI_Comm comm, std::size_t target_dofs,
373377 < target)
374378 {
375379 auto [new_mesh, _parent_edges, _parent_facet] = dolfinx::refinement::refine (
376- *mesh, std::nullopt ,
377- dolfinx::mesh::create_cell_partitioner (dolfinx::mesh::GhostMode::shared_facet, 2 ),
378- dolfinx::refinement::Option::parent_cell_and_facet);
380+ *mesh, std::nullopt ,
381+ dolfinx::mesh::create_cell_partitioner (
382+ dolfinx::mesh::GhostMode::shared_facet, 2 ),
383+ dolfinx::refinement::Option::parent_cell_and_facet);
379384 mesh = std::make_shared<dolfinx::mesh::Mesh<double >>(new_mesh);
380385 mesh->topology_mutable ()->create_entities (1 );
381386 }
@@ -411,9 +416,10 @@ create_spoke_mesh(MPI_Comm comm, std::size_t target_dofs,
411416 marked_edges.push_back (i);
412417
413418 auto [new_mesh, _parent_edges, _parent_facet] = dolfinx::refinement::refine (
414- *mesh, marked_edges,
415- dolfinx::mesh::create_cell_partitioner (dolfinx::mesh::GhostMode::shared_facet, 2 ),
416- dolfinx::refinement::Option::parent_cell_and_facet);
419+ *mesh, marked_edges,
420+ dolfinx::mesh::create_cell_partitioner (
421+ dolfinx::mesh::GhostMode::shared_facet, 2 ),
422+ dolfinx::refinement::Option::parent_cell_and_facet);
417423 meshi = std::make_shared<dolfinx::mesh::Mesh<double >>(new_mesh);
418424
419425 double actual_fraction
0 commit comments