add growth factor on DX1

Man-Yat Chu · Man-Yat Chu · commit e29bbe42806d · 2017-09-27T16:32:07.000-07:00
diff --git a/fastpm/multi.py b/fastpm/multi.py
@@ -73,7 +73,7 @@ def apply_density(k, v):
             source1[...] += sp.Omega(a) * source
             Omega_tot += sp.Omega(a)
 
-            DX1 = lpt1(source, Q)
+            DX1 = lpt1(source, Q)*pt.D1(a)
             sp.S[...] = DX1
 
             def apply_velocity(k, v):
