Return complex fc3_sum instead of fc3_sum squared. This is for implementing frequency shift of instantaneous phonon.
Preparing implementation of instantaneous phonon
Minor change