fealpy
fealpy copied to clipboard
Bugs in ScaledMonomialSpace2d
在 ScaledMonomialSpace2d.py
中, 如上图中标注的,
第 194 行应该为 ldof = self.number_of_local_dofs(p=p)
;
第 205 行, np.eisum
应该为 np.einsum
.
好的,这个程序我都没有测试过,谢谢。
其实这个关于基函数的实现还可以改进的,我找到时间会重新实现。
已经修改了。
魏老师, 在 ScaledMonomialSpace2d.py
中
-
def mass_matrix(self):
可以改为def mass_matrix(self, p=None):
的形式吗, - 还需要添加一个计算刚度矩阵的函数
stiff_matrix(self, p=None)
, 需要组装不同多项式阶数的刚度矩阵.
好的,我今天就加上去
我已经加上去了。 我想了一下,我还需要把这个缩放单项式空间重构一下,最近会找时间去做。