Rebased onto master and merged.
DirichletBC was not working correctly for quad and hex cells, because there was a mismatch in facet ordering between DOLFIN and FIAT. Now edges and faces numbering matches the numbering from entity_closure_dofs
DirichletBC
entity_closure_dofs
entity_closure_dofs for the quad element:
{0: {0: [0], 1: [1], 2: [2], 3: [3]}, 1: {0: [0, 1], 1: [2, 3], 2: [0, 2], 3: [1, 3]}, 2: {0: [0, 1, 2, 3]}}
entity_closure_dofs for the hex element:
{0: {0: [0], 1: [1], 2: [2], 3: [3], 4: [4], 5: [5], 6: [6], 7: [7]}, 1: {0: [0, 1], 1: [2, 3], 2: [4, 5], 3: [6, 7], 4: [0, 2], 5: [1, 3], 6: [4, 6], 7: [5, 7], 8: [0, 4], 9: [1, 5], 10: [2, 6], 11: [3, 7]}, 2: {0: [0, 1, 2, 3], 1: [4, 5, 6, 7], 2: [0, 1, 4, 5], 3: [2, 3, 6, 7], 4: [0, 2, 4, 6], 5: [1, 3, 5, 7]}, 3: {0: [0, 1, 2, 3, 4, 5, 6, 7]}}
DirichletBC
was not working correctly for quad and hex cells, because there was a mismatch in facet ordering between DOLFIN and FIAT. Now edges and faces numbering matches the numbering fromentity_closure_dofs
entity_closure_dofs
for the quad element:{0: {0: [0], 1: [1], 2: [2], 3: [3]}, 1: {0: [0, 1], 1: [2, 3], 2: [0, 2], 3: [1, 3]}, 2: {0: [0, 1, 2, 3]}}
entity_closure_dofs
for the hex element:{0: {0: [0], 1: [1], 2: [2], 3: [3], 4: [4], 5: [5], 6: [6], 7: [7]}, 1: {0: [0, 1], 1: [2, 3], 2: [4, 5], 3: [6, 7], 4: [0, 2], 5: [1, 3], 6: [4, 6], 7: [5, 7], 8: [0, 4], 9: [1, 5], 10: [2, 6], 11: [3, 7]}, 2: {0: [0, 1, 2, 3], 1: [4, 5, 6, 7], 2: [0, 1, 4, 5], 3: [2, 3, 6, 7], 4: [0, 2, 4, 6], 5: [1, 3, 5, 7]}, 3: {0: [0, 1, 2, 3, 4, 5, 6, 7]}}