code: def rotated_Gamma_ml(m,l,phi1,phi2,theta1,theta2,gamma_ml):
      
      
      rotated_gamma = 0
      
      for ii in range(2*l+1):
          rotated_gamma += Dlmk(l,m,ii-l,phi1,phi2,theta1,theta2).conjugate()*gamma_ml[ii]
  
      return rotated_gamma
docs: This function takes any gamma in the computational frame and rotates it to the
      cosmic frame.