# A_n: amplitude of spiral modes Â
  n_modes = 20  # Number of radial modes Â
  map_Q = np.zeros(12 L2) Â
  map_U = np.zeros(12  L2) Â
  for n in range(1, n_modes + 1): Â
    for m in [-2, 0, 2]:  # Even-parity modes for EB correlation Â
      psi_nlm = spiral_harmonic(n, m, beta, L) Â
      # Add to Stokes parameters (Q iU ~ E iB) Â
      map_Q += A_n  np.real(psi_nlm) Â
      map_U += A_n  np.imag(psi_nlm) Â
  return map_Q, map_U Â
Beri Komentar
Belum ada komentar. Jadilah yang pertama untuk memberikan komentar!