On 23 May 2016, at 10:54, David Ham <David.Ham@imperial.ac.uk> wrote:
For the special case of a periodic box, doing it using mesh._from_cell_list and constructing the connectivity and coordinates yourself is not that hard.
Note that if you do this you will need to create a 3D cube embedded in 4D and then remap to a 3D periodic cube. The reason for this is that from_cell_list assumes a P1 coordinate field. If you want to directly make the 3D periodic cube you need to interact with PETSc's DMPlex at a slightly lower level to first make the topology you want, then add the geometry and then build the firedrake mesh. This will probably also run into a problem with us because I think we currently assume that the coordinate field that comes out of DMPlex is in P1. But changing that is an easy fix. Cheers, Lawrence