decl_curve_rnn_trainingdata.py

isample = 0
for time_stage_change in np.linspace(*bounds_stage_change_time, num_discrete_stage_changes):
    for _ in range(num_realizations_per_stage_change):
        _, production, stage = calc_two_stage_decline(p0, exp_stage_zero, exp_stage_one, time_max, t
          time_stage_change, num=num_timesteps) 
        for num_sample_points in range(1, num_timesteps):
            features[isample, :, ifeature_stage] = stage[:] 
            features[isample, :num_sample_points, ifeature_production] = production[:num_sample_points]
            targets[isample, 1:, itarget_production] = production[1:] 
            isample += 1