fealpy icon indicating copy to clipboard operation
fealpy copied to clipboard

Bugs in ScaledMonomialSpace2d

Open yoczhang opened this issue 5 years ago • 6 comments

image

ScaledMonomialSpace2d.py 中, 如上图中标注的, 第 194 行应该为 ldof = self.number_of_local_dofs(p=p); 第 205 行, np.eisum 应该为 np.einsum.

yoczhang avatar Feb 14 '20 13:02 yoczhang

好的,这个程序我都没有测试过,谢谢。

weihuayi avatar Feb 15 '20 09:02 weihuayi

其实这个关于基函数的实现还可以改进的,我找到时间会重新实现。

weihuayi avatar Feb 15 '20 09:02 weihuayi

已经修改了。

weihuayi avatar Feb 15 '20 09:02 weihuayi

魏老师, 在 ScaledMonomialSpace2d.py

  1. def mass_matrix(self): 可以改为 def mass_matrix(self, p=None): 的形式吗,
  2. 还需要添加一个计算刚度矩阵的函数 stiff_matrix(self, p=None), 需要组装不同多项式阶数的刚度矩阵.

yoczhang avatar Feb 16 '20 12:02 yoczhang

好的,我今天就加上去

weihuayi avatar Feb 17 '20 08:02 weihuayi

我已经加上去了。 我想了一下,我还需要把这个缩放单项式空间重构一下,最近会找时间去做。

weihuayi avatar Feb 17 '20 12:02 weihuayi