Mapped restart fails to initialize multiblock mesh
When generating a mapped restart (write_restart.py -m
), if the mesh has multiple blocks I see the failure:
PANIC[1]: Assertion failed at /home/zjibben/codes/telluride/truchas/src/truchas/distributed_mesh/face_neighbor_table_type.F90:0111
Here is a small reproducer: mesh.jou mesh.gen run test.inp. This generates a simple 2-block hex mesh, runs a couple steps of a trivial heat transfer problem, then attempts to generate a restart mapped onto the same mesh. I've tried Cubit 15.5, 15.7, and 15.8, and Truchas 21.05 - 21.08 with GCC, Intel, and NAG, all with the same result.
With some print statements, it seems the problem is xcnode
somehow ends up with a cell containing a negative number of nodes.