存儲舊值FiPy

Question

我想解決一些在Python中使用FiPy的微分方程，並且作爲新手，我仍然有一些問題。我所做的是以下幾點：我定義了一個單元變量，我爲這個變量求解了一個方程，並且我更新了它。我想在每次迭代後存儲它的值。這裏有一個例子：

a = CellVariable(mesh,name='a', value=0., hasOld=True) 
# eq is an equation involving 'a' 
# define an array to store the values of 'a' after solving 'eq' 
a_tt = [] 
for t in range(10): 
     eq.sweep(dt=0.01) 
     a.updateOld() 
     a_tt.append(a) 

我意識到我的錯誤 - 在「a_tt」的價值觀我每次更新時間更新「A」，所以我必須在結束所有相同元素的數組。 爲了避免這種情況，我還需要做些什麼？

Answer 1

我認爲a_tt.append(a.copy())可能工作。

否則，在http://www.ctcms.nist.gov/fipy/examples/diffusion/generated/examples.diffusion.mesh1D.html掃描部分中使用的方法應該工作。例如：

a_tt.append(CellVariable(mesh=m, value=a.value))