1

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」を更新するたびに「a_tt」の値が更新されるため、最後にすべて同じ要素を持つ配列があります。これを回避するには、代わりに何をすべきですか?

4

1 に答える 1