// Start with the shell bounded by
// two nested cubes
- if (n <= 6)
+ if (n == 6)
{
for (unsigned int i=0;i<8;++i)
vertices.push_back(p+hexahedron[i]*irad);
cells[i].vertices[j] = cell_vertices[i][j];
cells[i].material_id = 0;
}
+
+ tria.create_triangulation (vertices, cells, SubCellData());
}
// A more regular subdivision can
// be obtained by two nested
// rhombic dodecahedra
- else if (n <= 12)
+ else if (n == 12)
{
for (unsigned int i=0;i<8;++i)
vertices.push_back(p+hexahedron[i]*irad);
}
cells[i].material_id = 0;
}
+
+ tria.create_triangulation (vertices, cells, SubCellData());
}
else
{
- Assert(false, ExcIndexRange(n, 1, 12));
+ Assert(false, ExcMessage ("Invalid number of coarse mesh cells."));
}
- tria.create_triangulation (vertices, cells,
- SubCellData()); // no boundary
- // information
-
if (colorize)
colorize_hyper_shell(tria, p, inner_radius, outer_radius);
}