Huayi Wei

我统一了一下, 已经更新到仓库了

正确的用法是 ``` from fealpy.mesh import MeshFactory as MF box = [0, 1, 0, 1] mesh = MF.boxmesh2d(box, nx=10, ny=10, meshtype='tri') ```

应该不可以，关键你说的 grad_basis 在边（面）上的值的具体数学含义，是单元的梯度值在边界上的限制呢？还是边界上切身梯度值？ 或者你需求背景是什么？

所以你的 \nabla q 是单元值在 \Gamma 上的限制，对吗? 你会给出边上的积分点，然后计算 \nabla q 在些积分点上的值？

@yoczhang 你看一下 这个代码 [edge_grad_basis ](https://github.com/weihuayi/fealpy/blob/cd1be360c6cafa7aa8f8c025882841712869fa1e/fealpy/functionspace/LagrangeFiniteElementSpace.py#L312)

另外，你的数值离散必须要在三角形网格上做吗？

是的

index 和 lidx 长度应该一样， 你要计算一个单元的基函数在某条边上积分点处的函数值，你需要知道单元的编号 index 和 这条边在这个单元上的局部编号。 我问 “另外，你的数值离散必须要在三角形网格上做吗？" 的意思是： 如果你的的离散方法适用于多边形网格的，最好不用 Lagrange 有限元空间。

好的，这个程序我都没有测试过，谢谢。

其实这个关于基函数的实现还可以改进的，我找到时间会重新实现。