YcZhang
YcZhang
[edge_grad_basis ](https://github.com/weihuayi/fealpy/blob/cd1be360c6cafa7aa8f8c025882841712869fa1e/fealpy/functionspace/LagrangeFiniteElementSpace.py#L312) 中 `lidx` 表示什么? 我目前测试中出现错误: `gphi = np.einsum('k...ij, kjm->k...im', R, Dlambda[index, :, :])`, `index` 中的指标超出了 `Dlambda` 的维数
嗯嗯, 目前就是用最简单的有限元来实现一下
魏老师, 在 `ScaledMonomialSpace2d.py` 中 1. `def mass_matrix(self):` 可以改为 `def mass_matrix(self, p=None):` 的形式吗, 2. 还需要添加一个计算刚度矩阵的函数 `stiff_matrix(self, p=None)`, 需要组装不同多项式阶数的刚度矩阵.
? 我这里可以正常运行 😂
好的. 魏老师这几天可不可也更新一下其他算例, 先有一个较稳定的版本, 后面更新时可以以不同的版本号给出. 另外, 我 python 升级到了 3.8, 给出了一堆 SyntaxWarning, 全部都是关于 `is` 的使用, 如 `if spacetype is 'C':` 中 `is` 更准确的应该为 `==`, 是否可以考虑换成 `==`?
3.8 版本的警告如下 ```python # python3.8 >>> spacetype='C' >>> spacetype is 'C' :1: SyntaxWarning: "is" with a literal. Did you mean "=="? True >>> ``` 而在 3.7 版本中没有此类警告 ```python # python3.7...
`sys.argv[]` 中其实就是一个列表，里面是存放的是程序外部输入的参数，例如 `example/PoissonFEMWithDirichletBC_example.py`，运行时可以用 `>python PoissonFEMWithDirichletBC_example.py 1 2 3 2`，其中 `1 2 3 2` 分别对应 ```python p = int(sys.argv[1]) n = int(sys.argv[2]) maxit = int(sys.argv[3]) d = int(sys.argv[4]) ``` --- 如果你直接在...
如果直接在终端下命令行运行的话，要看你调用的是哪个 python 了，比如，你是通过 Anaconda (虚拟环境) 安装的 fealpy，那么你需要先激活 fealpy 这个虚拟环境，然后在此环境下 `python fealpyexample.py`。你可以查看你的 IDE 调用的是哪个路径下的 python，然后可以在终端下输入绝对路径下的 python 也是可以的。