Commit 847d159f authored by Roland Haas's avatar Roland Haas

POWER: bugfix correct complex phase of strain

parent 98decfb7
......@@ -126,7 +126,9 @@ def psi4ToStrain(mp_psi4, f0):
hf = ffi(freq, dhf, f0)
time, h = myFourierTransformInverse(freq, hf, t0[0])
hTable = np.column_stack((time, h))
# we have d^2(h+ -i hx)/dt^2 = psi4 ie the sign of the imaginary part
# needs to be flipped
hTable = np.column_stack((time, h.conj()))
return hTable
#Fixed frequency integration
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment