Huayi Wei

已经合并进来了

不传 p 的时候，返回的是 QuadrangleMesh 对象。传 p 进去的时候是 LagrangeQuadrangleMesh 对象，这是四边形的高阶网格对象，其单个单元顶点的编号规则和 QuadrangleMesh 是不一样的。

四边形的高阶网格对象是为参数有限元服务的，因为其边可能是曲边的，其画图需要特别的处理方法。我会在首先加入对 p = 1 四边形网格的支持，p > 1 的情形需要再想一下。

``` import numpy as np import matplotlib.pyplot as plt from fealpy.mesh import MeshFactory as MF box = [0, 1, 0, 1] #mesh = MF.boxmesh2d(box, nx=2, ny=2, meshtype='tri', p=3) mesh =...

请不要用 oldexample 中的例子，那些代码已经不用了。

很抱歉回复的有点晚，最近太忙了，请见谅。

现在在 functionspace 模块中已经没有 tools 这个文件了

因为内部程序更新后，有一些例子程序没有及时更新，所以会出问题。你可以看下 doc/fealpycopyright.pdf, 还有 doc/fealpy_ppt.pdf， 先了解一下 fealpy 的设计思想和基本结构。

mesh 对象有一个 cell_data 的字典数据，你可以建立一个数组 比如 domain_flag, 形状是 （NC，），假设你有 5 种材料，分别用 1, 2, 3, 4, 5 五个整数来标示 domain_flag[i] == 3, 则标记第 i 个单元为 3 号材料， 这个数据可以放在 mesh 对象的 cell_data 字典中 mesh.cell_data['domain_flag']...

很抱歉, 刚看到你的这个消息. 这个问题解决了吗? 你可以试着更新到最新的 FEALPy.