Skip to content

Instantly share code, notes, and snippets.

@BrianMartell
Created December 5, 2025 20:45
Show Gist options
  • Select an option

  • Save BrianMartell/c86b0598c5d43a849adc3d53da543f98 to your computer and use it in GitHub Desktop.

Select an option

Save BrianMartell/c86b0598c5d43a849adc3d53da543f98 to your computer and use it in GitHub Desktop.
PUH-BrianMartell PUH_v6_DESI_Euclid_Sim.py-Updated New Py Code
import numpy as np
import matplotlib.pyplot as plt
theta = np.linspace(0, np.pi, 500)
A_H = 0.09
H_ratio = 1 + A_H * np.cos(theta)
plt.figure(figsize=(8,5))
plt.plot(np.rad2deg(theta), H_ratio, 'b-')
plt.axhline(1.09, color='r', ls='--', label='Fast pole +9%')
plt.axhline(0.91, color='k', ls='--', label='Slow pole –9%')
plt.xlabel('Angle from fast pole (deg)')
plt.ylabel('Normalized H(z)')
plt.title('PUH v6: Dipolar Hubble Flow (DESI/Euclid test)')
plt.legend()
plt.grid(True)
plt.savefig('puh_desi_euclid.png', dpi=300)
plt.show()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment