udpate

fealpy/mesh/HexahedronMesh.py

@@ -23,6 +23,7 @@ class HexahedronMeshDataStructure(Mesh3dDataStructure):
         [4, 0], [3, 0], [5, 0], [0, 2], 
         [2, 4], [4, 3], [3, 5], [5, 2], 
         [1, 4], [1, 3], [1, 5], [2, 1]])
+    ccw = np.array([0, 1, 2, 3])
     NVC = 8
     NEC = 12
     NFC = 6

fealpy/mesh/Mesh2d.py

@@ -200,11 +200,6 @@ def find_node(self, axes, node=None,
         if node is None:
             node = self.node
 
-
-        if (index is None) and (showindex == True):
-            NN = self.number_of_nodes() 
-            index = np.array(range(NN))
-
         find_node(axes, node,
                 index=index, showindex=showindex,
                 color=color, markersize=markersize,

fealpy/mesh/Mesh3d.py

@@ -102,8 +102,7 @@ def add_plot(
         if isinstance(plot, ModuleType):
             from mpl_toolkits.mplot3d import Axes3D
             fig = plot.figure()
-            fig.set_facecolor('white')
-            axes = fig.gca(projection='3d')
+            axes = plot.subplot(projection='3d')
         else:
             axes = plot
 

fealpy/mesh/TetrahedronMesh.py

@@ -10,6 +10,7 @@
 class TetrahedronMeshDataStructure(Mesh3dDataStructure):
     OFace = np.array([(1, 2, 3),  (0, 3, 2), (0, 1, 3), (0, 2, 1)])
     SFace = np.array([(1, 2, 3),  (0, 2, 3), (0, 1, 3), (0, 1, 2)])  
+    ccw = np.array([0, 1, 2])
 
     localFace = np.array([(1, 2, 3),  (0, 3, 2), (0, 1, 3), (0, 2, 1)])
     localEdge = np.array([(0, 1), (0, 2), (0, 3), (1, 2), (1, 3), (2, 3)])

fealpy/mesh/mesh_tools.py

@@ -30,9 +30,9 @@ def find_node(
         index = range(node.shape[0])
     elif (type(index) is np.int_):
         index = np.array([index], dtype=np.int_)
-    elif (type(index) is np.ndarray) and (index.dtype == np.bool):
+    elif (type(index) is np.ndarray) and (index.dtype == np.bool_):
         index, = np.nonzero(index)
-    elif (type(index) is list) and (type(index[0]) is np.bool):
+    elif (type(index) is list) and (type(index[0]) is np.bool_):
         index, = np.nonzero(index)
     else:
         pass
@@ -120,7 +120,7 @@ def find_entity(
         else:
             pass #TODO: raise a error
     elif (type(index) is np.ndarray) :
-        if index.dtype == np.bool:
+        if index.dtype == np.bool_:
             index, = np.nonzero(index)
     elif (type(index) is list) & (type(index[0]) is np.bool_):
         index, = np.nonzero(index)
@@ -438,8 +438,6 @@ def show_mesh_3d(
                color=edgecolor)
         return axes.add_collection3d(edges)
 
-    face = mesh.boundary_face(threshold=threshold)
-
     face = mesh.entity('face')
     isBdFace = mesh.ds.boundary_face_flag()
     if threshold is None:

tutorial/lagrange_space.py

@@ -3,7 +3,7 @@
 import matplotlib.pyplot as plt
 from mpl_toolkits.mplot3d import Axes3D
 
-from fealpy.decorator import baricentric
+from fealpy.decorator import barycentric
 from fealpy.mesh import MeshFactory as MF
 from fealpy.functionspace import LagrangeFiniteElementSpace