Skip to content

Instantly share code, notes, and snippets.

Show Gist options
  • Select an option

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

Select an option

Save BrianMartell/539dbf9987cd7ec37f72722a1be60559 to your computer and use it in GitHub Desktop.
PUH-BrianMartell PUH_T242_Cyclic_Rebound_Trigger_Sim.py-Updated New Py Code for Conceptual Simulation
import numpy as np
import matplotlib.pyplot as plt
# Two processing rates: external vs core-frame near saturation
folding_fraction = np.linspace(0, 0.999, 400)
external_rate = np.ones_like(folding_fraction) # normal external processing
core_rate = (1 - folding_fraction)**4 # collapses near saturation (T238)
plt.figure(figsize=(12, 7))
plt.plot(folding_fraction, external_rate, 'b-', label='External Substrate Processing Rate')
plt.plot(folding_fraction, core_rate, 'r--', label='Core-Frame Processing Rate (→0 at saturation)')
plt.axvline(0.99, color='orange', linestyle=':', label='T173 Snap Threshold')
plt.title('PUH T242: Two Timescales in Cyclic Rebound')
plt.xlabel('Folding Fraction (approaching 1)')
plt.ylabel('Relative Processing Rate (Time)')
plt.legend(); plt.grid(True); plt.show()
print("Shows external 'trillions of years' vs core 'near-instant' as same event.")
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment