YcZhang

Two questions in the demo MixedPoisson.py

7

Dear Mikael, In the [MixedPoisson.py](https://github.com/spectralDNS/shenfun/blob/4faa1fbd364380dd52ff0bbce88f54e8d42bce08/demo/MixedPoisson.py#L1), I have two questions: 1 . How to modify it to periodic boundary conditions in both x- and y-directions. I tried to change `SD =...

TriangleMesh 中 coarsen() 对特定网格的 bug ?

1

如果生成 `米` 字型网格, `coarsen()` 后出现 bug: 1. 原始网格  2. 只标记第 `42` 个单元, 进行 `coarsen()` 后的网格  --- 下面是测试代码 ```python from fealpy.mesh import MeshFactory as MF import numpy as np...

LagrangeFiniteElementSpace 中 grad_basis 不支持 TD-1 维的积分点

11

需要求 `grad_basis()` 在边(面)上的值. 能否将 `grad_basis()` 中的 `TD=self.TD` 改为类似于 `basis()` 中的 `TD = bc.shape[-1] - 1 `?

ScaledMonomialSpace2d 中 value, grad_value 加入 p=None 参数

1

对应 `ScaledMonomialSpace2d` 中 `basis`, `grad_basis` 函数, `value`, `grad_value` 能否加入 `p=None` 参数?

the `integral` function in `class FEMeshIntegralAlg()`

11

从 line-621 --- line-640 ```python # old api def integral(self, u, celltype=False, barycenter=True): """ """ qf = self.integrator bcs = qf.quadpts # 积分点 (NQ, 3) ws = qf.weights # 积分点对应的权重...

the FEMeshIntegralAlg bug in ScaledMonomialSpace2d

11

魏老师, 可以先把 `FEMeshIntegralAlg` 从 `ScaledMonomialSpace2d` 中去掉吗? 对于 tri mesh, 做积分时 `FEMeshIntegralAlg' 仍然存在 bug, 魏老师可以尝试运行一下下面的测试程序 ```python from fealpy.functionspace.ScaledMonomialSpace2d import ScaledMonomialSpace2d import numpy as np from fealpy.mesh import TriangleMesh # --- mesh...

```python from fealpy.mesh import StructureQuadMesh box = [0, 1, 0, 1] nx = 2 ny = 2 mesh = StructureQuadMesh(box, nx, ny) cellarea = mesh.entity_measure('cell') ``` then, the error is...

add the `number_of_edges_of_cells()` function to the class StructureQuadMeshDataStructure?

13

魏老师, 可以在 `class StructureQuadMeshDataStructure` 中添加 `number_of_edges_of_cells()` 函数么? 添加的时, `number_of_edges_of_cells()` 可以返回 `(NC,)` 个值, 而不是单值 `4` 么?

The bug in ScaledMonomialSpace2d's stiff_matrix()??

1

魏老师, `ScaledMonomialSpace2d()` 中 `stiff_matrix()` ```python def f(x, index): gphi = self.grad_basis(x, index=index, p=p) return np.einsum('jkm, jpm->jkp', gphi, gphi) ``` 是否应该为 ```python def f(x, index): gphi = self.grad_basis(x, index=index, p=p) return...

mesh.ds.cell_to_edge() 中 keyword 没统一

1

`mesh.ds.cell_to_edge()` 中 keyword 没统一 --- 例如 `fealpy/functionspace/NonConformingVirtualElementSpace2d.py -- line 50, cell2edge = mesh.ds.cell_to_edge(sparse=False)` 如果用 polygonmesh 的话, 会提示 keyword 错误