хранить старые значения FiPy

Question

Я пытаюсь решить некоторые дифференциальные уравнения, используя FiPy в Python и как новичок, у меня все еще есть некоторые проблемы. Я делаю следующее: я определяю переменную ячейки, я решаю уравнение для этой переменной, и я ее обновляю. Я хочу сохранить его значения после каждой итерации. Вот пример:

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» обновляется каждый раз, когда я обновляю «а», так что я в конце концов массив со всеми теми же элементами. Что я могу сделать, чтобы избежать этого?

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))