Creating a mesh with UnitIntervalMesh causes crash

[funsim]

Hi,

this simple script results in a *** FLUIDITY ERROR ***:

  from firedrake import *
  UnitIntervalMesh(2)

The error message I get is:

*** FLUIDITY ERROR ***
Source location: (Fields_Allocates.F90, 1467)
Error message: Surface element in sndgln is not on the surface.
Backtrace will follow if it is available:
/home/firedrake/src/firedrake/lib/libfluidity.so.1(fprint_backtrace_+0x1f) [0x7fa07dec9a59]
/home/firedrake/src/firedrake/lib/libfluidity.so.1(__fldebug_MOD_flabort_pinpoint+0x277) [0x7fa07dee6ef7]
/home/firedrake/src/firedrake/lib/libfluidity.so.1(__fields_allocates_MOD_add_faces_face_list+0xd55) [0x7fa07df53c15]
/home/firedrake/src/firedrake/lib/libfluidity.so.1(__fields_allocates_MOD_add_faces+0x2056) [0x7fa07df5afd6]
/home/firedrake/src/firedrake/lib/libfluidity.so.1(__read_triangle_MOD_read_triangle_files_to_field+0x1f33) [0x7fa07e02dc43]
/home/firedrake/src/firedrake/lib/libfluidity.so.1(__read_triangle_MOD_read_triangle_simple+0x29b) [0x7fa07e030afb]
/home/firedrake/src/firedrake/lib/libfluidity.so.1(__mesh_files_MOD_read_mesh_simple+0x49e) [0x7fa07dfde81e]
/home/firedrake/src/firedrake/lib/libfluidity.so.1(read_mesh_f+0x447) [0x7fa07e003077]
/home/firedrake/src/firedrake/python/firedrake/core_types.so(+0x24cc5) [0x7fa07e5d3cc5]
python() [0x4ce25e]
python(PyObject_Call+0x36) [0x4abf96]
/home/firedrake/src/firedrake/python/firedrake/core_types.so(+0x3005d) [0x7fa07e5df05d]
python(PyEval_EvalFrameEx+0xa5b) [0x47bf2b]
python(PyEval_EvalCodeEx+0x19a) [0x4e048a]
python() [0x4e0eea]
python() [0x4ce25e]
python() [0x4a0534]
python() [0x4ae24e]
python(PyEval_EvalFrameEx+0xa5b) [0x47bf2b]
python(PyEval_EvalCodeEx+0x19a) [0x4e048a]
python(PyEval_EvalCode+0x32) [0x53ff12]
python() [0x54038b]
python(PyRun_FileExFlags+0x92) [0x4658f4]
python(PyRun_SimpleFileExFlags+0x2d8) [0x465e24]
python(Py_Main+0xb4e) [0x466b99]
/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7fa07ef14ea5]
python() [0x4e1115]
Use addr2line -e <binary> <address> to decipher.
Error is terminal.
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 16.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.

I installed the entire PyOP2 + Firedrake system today, so I should be on the most up to date version.

[kynan]

I can reproduce the issue.