Huayi Wei

它引入了半边数据结构, 可以统一表示 tri, quad, polygon 数据结构. 核心数据结构采用了动态数组的技术, 更适合自适应算法.

原始的数据结构中, 有些接口返回的数据是有区别的, 比如单元边的个数, 三角形返回一个整数 3, 四边形返回整数 4, 多边形返回一个长度为 (NC, ) 的一维数组. 这样返回有效率方面的考虑, 这样也更切合用户对不同网格类型的预期.

但也许还是可以改进一下的, 还是上面的单元边个数的例子, 三角形返回一个长度为 1 的一维数组 np.array([3]) 即可. 这样可以以数值离散中利用数组的广播功能, 简化数值离散程序对网格进行区分编程的繁琐.

> 是的, 魏老师, 并不影响计算, 我改一下自己的程序就好了, 这个空间是 hybrid high-order 方法的空间, 类似于 WG 的空间, 不过重构算子时差别比较大 正常情况下, 三角形, 四边形是特殊的多边形, 所以我感觉重构算子时不应该差别很大的, 应该是一样的, 你再考虑一下.

可以实现，直接可以用半边数据结构，也可以重新写一个复合网格类（这个需要一些工作）

@Geniuswangwei 你的这个网格看起来很怪

这个完全可以做到，我们要讨论一下，如何更好支持这种带 PML 层的应用

我把时间取消了，你更新一下 FEALPy

我不太理解你的问题，能说的理详细一点吗?     ------------------ Original ------------------ From: ***@***.***>; Date:  Fri, Mar 19, 2021 09:57 AM To: ***@***.***>; Cc: ***@***.***>; Subject:  [weihuayi/fealpy] 插值函数 (#33)   魏老师，通过有限元数值解插值全域函数，只要给出一个点的位置，就可以知道插值全域函数值，您开发的包里面的这个封装函数是如何实现的，能仔细讲解一下吗？ — You are receiving this because you...

p 属于什么空间？     ------------------ Original ------------------ From: ***@***.***>; Date:  Sat, Mar 20, 2021 12:23 PM To: ***@***.***>; Cc: ***@***.***>; ***@***.***>; Subject:  Re: [weihuayi/fealpy] 插值函数 (#33)   魏老师，我在求一个非线性方程需要迭代求解，首通过有限元求解方程得到数值离散解uh，然后非线性项的变量是关于p，而p=2|u|^2，因此我需要把离散解插值回计算区域的函数，这样只需要给出函数点，就可以直接得到p的值，好像您开发包里面有关于数值解构造函数的，所以想问问您实现的具体细节？ — You are receiving...