Refs #11. fixed a bug in mp.forward that causes weird I(E=0)

neutrons · May 11, 2016 · 99f16d9 · 99f16d9
1 parent aaec434
commit 99f16d9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/multiphonon/forward.py b/multiphonon/forward.py
@@ -201,7 +201,7 @@ def computeA1E(E,g, beta, dE):
     t = g/(E*g0)*t
     z = zero_ind
     # remove NaN
-    t[z] = 2.0*( t[z+1] + ( t[z+1] - t[z+2] ) ) # XXX: why 2.0* ?
+    t[z] = t[z+1] + ( t[z+1] - t[z+2] )
     # XXX: normalize?
     # t /= t.sum()
     return E, t