确认一下HexMeshDataStructure的local排序

[scaomath]

@weihuayi 似乎现在FEALPy里面face的顺序不再像魏老师在ifem里面实现的那样了

这个是ifem

elem = [s s+ni s+ni+nij s+nij s+1 s+ni+1 s+ni+nij+1 s+nij+1]; 
totalFace = uint32([elem(:,[1 4 3 2]); elem(:,[1 2 6 5]);...
                               elem(:,[5 6 7 8]); elem(:,[8 7 3 4]);...
                               elem(:,[4 1 5 8]); elem(:,[2 3 7 6])]);

这个是现在的FEALPy

fealpy/fealpy/mesh/StructureHexMesh.py

Lines 74 to 77 in ea12ef7

	localFace = np.array([
	(0, 3, 2, 1), (4, 5, 6, 7), # bottom and top faces
	(0, 4, 7, 3), (1, 2, 6, 5), # left and right faces
	(0, 1, 5, 4), (2, 3, 7, 6)])# front and back faces

如果老ifem里面的排序建立这个mesh类，是不是那个类里面的face2cell或者face_to_cell的函数不能直接用？

[weihuayi]

是的，不能直接用。StructureHexMesh.py 中实现 face_to_cell 是按照前面的 localFace 顺序来实现的，但是比较手工化的实现。在非结构的六面体网格中，你给个不同的 localFace，后面的 face_to_cell 仍然可以正确的工作。