YcZhang
YcZhang
The same problem.
Yes, it works for me.
@mikaem Hi Mikael and other question about all periodic bcs, I don't quite understand the choice of solver, if I want to solve the [biharmonic2D](https://github.com/spectralDNS/shenfun/blob/f2a9078b04ffa813b68498566da168fdf4302cab/demo/biharmonic2D.py#L1) problem with all periodic bcs...
OK, thank you!
>这里不是程序的 Bug, 是少输入了一个参数, FEMeshIntegralAlg 中默认处理的是 三角形, 四边形, 四面体, 六面体这样规则的单元, 它的输入点, 都是重心坐标的形式, 所以 barycenter 默认设置为 True, 你要处理输入是笛卡尔坐标点时, 就要设置为 False. 这里确实不是 `FEMeshIntegralAlg` 本身的 bug, 而是放在 `ScaledMonomialSpace2d` 中会对统一的程序有影响. 我猜想 `ScaledMonomialSpace2d` 的设计之初是为任意多边形分片间断 (DG) 空间考虑的...
>我正在尝试一种新坐标系 Decorator, 用户要积分的函数, 加上一个坐标系的 Decorator, 就会自动在你的函数对象里加一个属性, 这样数值积分程序, 可以自动判断你输入的函数需要什么样的坐标点了. 但还在尝试中, 但还有一些问题没有完全解决掉. 期待.
嗯嗯 好的~. 不过我还是觉得在 `ScaledMonomialSpace2d` 中用 meshtype 来区分不同的空间不太合理 😂, 并且 fealpy 下 tri mesh 是可以直接用 `PolygonMeshIntegralAlg` 做积分的. 也就是程序本来不用考虑 meshtype, 现在需要区别一下, 影响也不大.
不过 `value`, `grad_value` 加入 `p=None` 后, 相应的 `cell2dof` 也需要修改, 貌似不是一个好的建议 ... ? -_-!