@@ -25,16 +25,14 @@ class recomb_model(eqx.Module):
2525
2626 integration_spacing : jnp .float64
2727 lna_axis_full : jnp .array
28- concrete_axis_size : jnp .array
29- concrete_axis_size_postSahaHe : jnp .array
3028
3129 z1 : jnp .float64
3230
3331 twog_redshift : jnp .float64
3432 He4equil_redshift : jnp .float64
3533 idx_4He_equil : jnp .array
3634
37- def __init__ (self ,integration_spacing = 5.0e-4 , Nsteps = 800 , Nsteps_postSahaHe = 4000 , z0 = 8000. , z1 = 0. ):
35+ def __init__ (self ,integration_spacing = 5.0e-4 , z0 = 8000. , z1 = 0. ):
3836 """
3937 Initialize complete recombination model.
4038
@@ -45,10 +43,6 @@ def __init__(self,integration_spacing = 5.0e-4, Nsteps=800, Nsteps_postSahaHe=40
4543 -----------
4644 integration_spacing : float, optional
4745 Step size for integration (default: 5.0e-4)
48- Nsteps : int, optional
49- Maximum integration steps (default: 800)
50- Nsteps_postSahaHe : int, optional
51- Maximum steps for post-Saha helium phase (default: 4000)
5246 z0 : float, optional
5347 Initial redshift (default: 8000.)
5448 z1 : float, optional
@@ -59,8 +53,6 @@ def __init__(self,integration_spacing = 5.0e-4, Nsteps=800, Nsteps_postSahaHe=40
5953
6054 # Define time axes
6155 self .lna_axis_full = jnp .arange (- jnp .log (1 + z0 ), - jnp .log (1 + z1 ), self .integration_spacing )
62- self .concrete_axis_size = jnp .zeros (Nsteps )
63- self .concrete_axis_size_postSahaHe = jnp .zeros (Nsteps_postSahaHe )
6456
6557 self .twog_redshift = 701.
6658 self .He4equil_redshift = 3601. # generous
0 commit comments