Huayi Wei

@lalala753 我已经安排人来实现这个功能

@jibingquanm 第一个问题，在 TriangleMesh 中，有一个 `point_to_bc(point)` 可以用，该方法可以找到 point 所在的单元，并计算出它对应的重心坐标。 `index, bc = mesh.point_to_bc(point)`

@jibingquanm 你用的有限元是线性元，还是高次的元。

要用到 LagrangeTriangleMesh（https://github.com/weihuayi/fealpy/blob/master/fealpy/mesh/lagrange_triangle_mesh.py）。 ``` from fealpy.mesh import TriangleMesh, LagrangeTriangleMesh mesh = TriangleMesh.from_box(nx=2, ny=2) node = mesh.entity('node') cell = mesh.entity('cell') lmesh = LagrangeTriangleMesh(node, cell, p=3) # 3 次 Lagrange 网格 lmesh.nodedata['uh'] =...

@jibingquanm 上面的代码比较旧了，我安排了一个学生基于最新的接口来写上面的代码，这一周应该可以完成。非常抱歉，这一段事情很多，没有及时回复您。

@zhf-0 单元节点的顺序要保证满足右手法则。另外，你的加密方式是什么？

给定一个四面体单元 $T = (x_0, x_1, x_2, x_3)$, 满足右手法则，向量 $x_0x_1$, $x_0x_2$ 和 $x_0x_3$ 的混合积是一个正数

@zhf-0 我已经安排同学来检查测试了。

@zhf-0 程序已经更新了，你测试一下吧

@zhf-0 你是说 bisect 自适应吗？