在fipy中相邻/连接顶点的单元格的方法？

from fipy import Grid2D

# 创建一个2D网格
mesh = Grid2D(nx=3, ny=3)

# 获取单元格的面的ID
cell_face_ids = mesh.getCellFaceIDs()

# 遍历每个单元格
for cell in mesh.cells:
    # 获取当前单元格的索引
    cell_index = cell.mesh.cellFaceIDs[cell]
    
    # 获取相邻单元格的面的ID
    neighbor_face_ids = cell_face_ids[cell_index]
    
    # 遍历相邻单元格的面的ID
    for face_id in neighbor_face_ids:
        # 获取相邻单元格的索引
        neighbor_index = mesh.faceCellIDs[face_id]
        
        # 根据相邻单元格的索引进行相应操作
        # ...