Schematic representation of a single timestep in our full network code. { Y₀(Z₁,A₁),Y₀(Z₂,A₂),... } denotes the set of initial nuclei abundances, ρ₀ is the initial mass density, ρ(t) defines the density evolution, T₀ is the initial temperature, Y_e,0 denotes the initial electron fraction, and n_n,0 is the initial neutron density. First the neutron decay is computed before the reaction network is solved using the Crank-Nicholson algorithm. Next the fission contribution is calculated along with the new physical parameters. If the changes in abundance or n_n are too large, the timestep is reattempted with dt = 0.1dt. Adaptive timesteps are used to maximize dt.