add calculate_G and calculate_Sigma
change some multiplications of complex numbers
create function d_times_cd; change exp to cdexp
slightly changes in self_consistent.f90