之前采用混合元求解非定常偏微分方程时,更新数值解都没哟问题,最近算了一个类似方程,更新一直报错。
solve(a == L, w) (c, m) = w.split() # # update the numerical solution c_n.assign(c)
程序源代码如上,找了好久没发现错误,最后更改为一下写法可算
solve(a == L, w) (c, m) = w.split(deepcopy=True) # # update the numerical solution c_n.assign(c)