Huayi Wei

integral 这个 api 是没有引入 coordtype 以前的 API， 目前为兼容以前的例子，就留在这里。单元上的积分， 建议使用 cell_integral. 或者，integral 保留， barycenter 这个参数也去掉 ``` def integral(self, u, etype='cell', q=None): """ Notes ----- 计算函数 u 在指定网格实体上的整体积分。 """ if etype ==...

或者命名为 `mesh_integral(...)`?

指明是在整个网格上的积分

原始的 `integral`, 仍然保留， 只是不再判断 coordtype, 这样还和一些老的程序兼容。 而 PolygonMesh 上的积分算法也会考虑和 FEM 网格积分算法的接口调整为一致的。

好的， 我再测试一下。     ------------------ Original ------------------ From:  "YcZhang"

你定义函数的时候， index 不能默认为 None， 应该是 `np.s_[:]` ``` @cartesian def f1(x, index=np.s_[:]): return uh.value(x, index) ```

quad 的积分和三角形的模式是一样的。     ------------------ Original ------------------ From:  "YcZhang"

是的， 我需要把逻辑修正一下，就可以了。

这个程序还没写好, 更新的版本还在修改中. 改完就会上传上来.

我刚上传一个新的版本, 但运行的结果还不是太好, 还需要进一步的优化.