; LLVM IR from CLASS-SZ/class_sz :: class-sz/hyrec (history)

; Unoptimized IR (clang -disable-llvm-passes):
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-i128:128-f80:128-n8:16:32:64-S128"
target triple = "x86_64-unknown-linux-gnu"

%struct.REC_COSMOPARAMS = type { double, double, double, double, double, double, double, double, double, double, double, double, double, double, i64, double, i16, double, double, double, double, double, double, double }

; Function Attrs: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.start.p0(i64 immarg, ptr nocapture) #0

; Function Attrs: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.end.p0(i64 immarg, ptr nocapture) #0

; Function Attrs: nounwind uwtable
declare dso_local double @rec_HubbleConstant(ptr noundef, double noundef) #1

; Function Attrs: nounwind
declare double @exp(double noundef) #2

; Function Attrs: nounwind uwtable
declare dso_local double @rec_Tmss(double noundef, double noundef, double noundef, double noundef, double noundef, double noundef) #1

; Function Attrs: nounwind uwtable
declare dso_local double @rec_dTmdlna(double noundef, double noundef, double noundef, double noundef, double noundef, double noundef, double noundef) #1

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.fmuladd.f64(double, double, double) #3

; Function Attrs: nounwind uwtable
declare dso_local void @rec_get_xe_next1(ptr noundef, double noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, i32 noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef) #1

; Function Attrs: nounwind uwtable
declare dso_local double @energy_injection_rate(ptr noundef, double noundef) #1

; Function Attrs: nounwind uwtable
declare dso_local void @rec_get_xe_next2(ptr noundef, double noundef, double noundef, double noundef, ptr noundef, ptr noundef, ptr noundef, i32 noundef, i32 noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef) #1

; Function Attrs: nounwind uwtable
define dso_local void @rec_build_history(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, ptr noundef %4) #1 {
  %6 = alloca ptr, align 8
  %7 = alloca ptr, align 8
  %8 = alloca ptr, align 8
  %9 = alloca ptr, align 8
  %10 = alloca ptr, align 8
  %11 = alloca i64, align 8
  %12 = alloca ptr, align 8
  %13 = alloca [3 x ptr], align 16
  %14 = alloca double, align 8
  %15 = alloca double, align 8
  %16 = alloca double, align 8
  %17 = alloca double, align 8
  %18 = alloca double, align 8
  %19 = alloca double, align 8
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca double, align 8
  store ptr %0, ptr %6, align 8, !tbaa !5
  store ptr %1, ptr %7, align 8, !tbaa !5
  store ptr %2, ptr %8, align 8, !tbaa !5
  store ptr %3, ptr %9, align 8, !tbaa !5
  store ptr %4, ptr %10, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %11) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #5
  call void @llvm.lifetime.start.p0(i64 24, ptr %13) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %17) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %18) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %19) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %20) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %22) #5
  %23 = load ptr, ptr %6, align 8, !tbaa !5
  %24 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %23, i32 0, i32 14
  %25 = load i64, ptr %24, align 8, !tbaa !9
  %26 = trunc i64 %25 to i32
  %27 = call ptr @create_2D_array(i32 noundef 311, i32 noundef %26)
  store ptr %27, ptr %12, align 8, !tbaa !5
  %28 = load ptr, ptr %6, align 8, !tbaa !5
  %29 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %28, i32 0, i32 14
  %30 = load i64, ptr %29, align 8, !tbaa !9
  %31 = trunc i64 %30 to i32
  %32 = call ptr @create_1D_array(i32 noundef %31)
  %33 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  store ptr %32, ptr %33, align 16, !tbaa !5
  %34 = load ptr, ptr %6, align 8, !tbaa !5
  %35 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %34, i32 0, i32 14
  %36 = load i64, ptr %35, align 8, !tbaa !9
  %37 = trunc i64 %36 to i32
  %38 = call ptr @create_1D_array(i32 noundef %37)
  %39 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 1
  store ptr %38, ptr %39, align 8, !tbaa !5
  %40 = load ptr, ptr %6, align 8, !tbaa !5
  %41 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %40, i32 0, i32 14
  %42 = load i64, ptr %41, align 8, !tbaa !9
  %43 = trunc i64 %42 to i32
  %44 = call ptr @create_1D_array(i32 noundef %43)
  %45 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 2
  store ptr %44, ptr %45, align 16, !tbaa !5
  %46 = load ptr, ptr %6, align 8, !tbaa !5
  %47 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %46, i32 0, i32 11
  %48 = load double, ptr %47, align 8, !tbaa !14
  store double %48, ptr %15, align 8, !tbaa !15
  store double 1.000000e+00, ptr %22, align 8, !tbaa !15
  store i64 0, ptr %11, align 8, !tbaa !16
  br label %49

49:                                               ; preds = %97, %5
  %50 = load i64, ptr %11, align 8, !tbaa !16
  %51 = load ptr, ptr %6, align 8, !tbaa !5
  %52 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %51, i32 0, i32 14
  %53 = load i64, ptr %52, align 8, !tbaa !9
  %54 = icmp slt i64 %50, %53
  br i1 %54, label %55, label %58

55:                                               ; preds = %49
  %56 = load double, ptr %22, align 8, !tbaa !15
  %57 = fcmp ogt double %56, 1.000000e-09
  br label %58

58:                                               ; preds = %55, %49
  %59 = phi i1 [ false, %49 ], [ %57, %55 ]
  br i1 %59, label %60, label %100

60:                                               ; preds = %58
  %61 = load ptr, ptr %6, align 8, !tbaa !5
  %62 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %61, i32 0, i32 11
  %63 = load double, ptr %62, align 8, !tbaa !14
  %64 = fadd double 1.000000e+00, %63
  %65 = load ptr, ptr %6, align 8, !tbaa !5
  %66 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %65, i32 0, i32 13
  %67 = load double, ptr %66, align 8, !tbaa !17
  %68 = fneg double %67
  %69 = load i64, ptr %11, align 8, !tbaa !16
  %70 = sitofp i64 %69 to double
  %71 = fmul double %68, %70
  %72 = call double @exp(double noundef %71) #5
  %73 = call double @llvm.fmuladd.f64(double %64, double %72, double -1.000000e+00)
  store double %73, ptr %15, align 8, !tbaa !15
  %74 = load ptr, ptr %6, align 8, !tbaa !5
  %75 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %74, i32 0, i32 9
  %76 = load double, ptr %75, align 8, !tbaa !18
  %77 = load ptr, ptr %6, align 8, !tbaa !5
  %78 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %77, i32 0, i32 0
  %79 = load double, ptr %78, align 8, !tbaa !19
  %80 = load ptr, ptr %6, align 8, !tbaa !5
  %81 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %80, i32 0, i32 10
  %82 = load double, ptr %81, align 8, !tbaa !20
  %83 = load double, ptr %15, align 8, !tbaa !15
  %84 = call double @rec_sahaHeII(double noundef %76, double noundef %79, double noundef %82, double noundef %83, ptr noundef %22)
  %85 = load ptr, ptr %9, align 8, !tbaa !5
  %86 = load i64, ptr %11, align 8, !tbaa !16
  %87 = getelementptr inbounds double, ptr %85, i64 %86
  store double %84, ptr %87, align 8, !tbaa !15
  %88 = load ptr, ptr %6, align 8, !tbaa !5
  %89 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %88, i32 0, i32 0
  %90 = load double, ptr %89, align 8, !tbaa !19
  %91 = load double, ptr %15, align 8, !tbaa !15
  %92 = fadd double 1.000000e+00, %91
  %93 = fmul double %90, %92
  %94 = load ptr, ptr %10, align 8, !tbaa !5
  %95 = load i64, ptr %11, align 8, !tbaa !16
  %96 = getelementptr inbounds double, ptr %94, i64 %95
  store double %93, ptr %96, align 8, !tbaa !15
  br label %97

97:                                               ; preds = %60
  %98 = load i64, ptr %11, align 8, !tbaa !16
  %99 = add nsw i64 %98, 1
  store i64 %99, ptr %11, align 8, !tbaa !16
  br label %49, !llvm.loop !21

100:                                              ; preds = %58
  store double 0.000000e+00, ptr %22, align 8, !tbaa !15
  br label %101

101:                                              ; preds = %152, %100
  %102 = load i64, ptr %11, align 8, !tbaa !16
  %103 = load ptr, ptr %6, align 8, !tbaa !5
  %104 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %103, i32 0, i32 14
  %105 = load i64, ptr %104, align 8, !tbaa !9
  %106 = icmp slt i64 %102, %105
  br i1 %106, label %107, label %110

107:                                              ; preds = %101
  %108 = load double, ptr %22, align 8, !tbaa !15
  %109 = fcmp olt double %108, 5.000000e-04
  br label %110

110:                                              ; preds = %107, %101
  %111 = phi i1 [ false, %101 ], [ %109, %107 ]
  br i1 %111, label %112, label %155

112:                                              ; preds = %110
  %113 = load ptr, ptr %6, align 8, !tbaa !5
  %114 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %113, i32 0, i32 11
  %115 = load double, ptr %114, align 8, !tbaa !14
  %116 = fadd double 1.000000e+00, %115
  %117 = load ptr, ptr %6, align 8, !tbaa !5
  %118 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %117, i32 0, i32 13
  %119 = load double, ptr %118, align 8, !tbaa !17
  %120 = fneg double %119
  %121 = load i64, ptr %11, align 8, !tbaa !16
  %122 = sitofp i64 %121 to double
  %123 = fmul double %120, %122
  %124 = call double @exp(double noundef %123) #5
  %125 = call double @llvm.fmuladd.f64(double %116, double %124, double -1.000000e+00)
  store double %125, ptr %15, align 8, !tbaa !15
  %126 = load ptr, ptr %6, align 8, !tbaa !5
  %127 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %126, i32 0, i32 9
  %128 = load double, ptr %127, align 8, !tbaa !18
  %129 = load ptr, ptr %6, align 8, !tbaa !5
  %130 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %129, i32 0, i32 0
  %131 = load double, ptr %130, align 8, !tbaa !19
  %132 = load ptr, ptr %6, align 8, !tbaa !5
  %133 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %132, i32 0, i32 10
  %134 = load double, ptr %133, align 8, !tbaa !20
  %135 = load ptr, ptr %6, align 8, !tbaa !5
  %136 = load double, ptr %15, align 8, !tbaa !15
  %137 = call double @rec_HubbleConstant(ptr noundef %135, double noundef %136)
  %138 = load double, ptr %15, align 8, !tbaa !15
  %139 = call double @xe_PostSahaHe(double noundef %128, double noundef %131, double noundef %134, double noundef %137, double noundef %138, ptr noundef %22)
  %140 = load ptr, ptr %9, align 8, !tbaa !5
  %141 = load i64, ptr %11, align 8, !tbaa !16
  %142 = getelementptr inbounds double, ptr %140, i64 %141
  store double %139, ptr %142, align 8, !tbaa !15
  %143 = load ptr, ptr %6, align 8, !tbaa !5
  %144 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %143, i32 0, i32 0
  %145 = load double, ptr %144, align 8, !tbaa !19
  %146 = load double, ptr %15, align 8, !tbaa !15
  %147 = fadd double 1.000000e+00, %146
  %148 = fmul double %145, %147
  %149 = load ptr, ptr %10, align 8, !tbaa !5
  %150 = load i64, ptr %11, align 8, !tbaa !16
  %151 = getelementptr inbounds double, ptr %149, i64 %150
  store double %148, ptr %151, align 8, !tbaa !15
  br label %152

152:                                              ; preds = %112
  %153 = load i64, ptr %11, align 8, !tbaa !16
  %154 = add nsw i64 %153, 1
  store i64 %154, ptr %11, align 8, !tbaa !16
  br label %101, !llvm.loop !23

155:                                              ; preds = %110
  %156 = load ptr, ptr %6, align 8, !tbaa !5
  %157 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %156, i32 0, i32 11
  %158 = load double, ptr %157, align 8, !tbaa !14
  %159 = fadd double 1.000000e+00, %158
  %160 = load ptr, ptr %6, align 8, !tbaa !5
  %161 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %160, i32 0, i32 13
  %162 = load double, ptr %161, align 8, !tbaa !17
  %163 = fneg double %162
  %164 = load i64, ptr %11, align 8, !tbaa !16
  %165 = sub nsw i64 %164, 3
  %166 = sitofp i64 %165 to double
  %167 = fmul double %163, %166
  %168 = call double @exp(double noundef %167) #5
  %169 = call double @llvm.fmuladd.f64(double %159, double %168, double -1.000000e+00)
  store double %169, ptr %18, align 8, !tbaa !15
  %170 = load ptr, ptr %9, align 8, !tbaa !5
  %171 = load i64, ptr %11, align 8, !tbaa !16
  %172 = sub nsw i64 %171, 2
  %173 = getelementptr inbounds double, ptr %170, i64 %172
  %174 = load double, ptr %173, align 8, !tbaa !15
  %175 = load ptr, ptr %9, align 8, !tbaa !5
  %176 = load i64, ptr %11, align 8, !tbaa !16
  %177 = sub nsw i64 %176, 4
  %178 = getelementptr inbounds double, ptr %175, i64 %177
  %179 = load double, ptr %178, align 8, !tbaa !15
  %180 = fsub double %174, %179
  %181 = fdiv double %180, 2.000000e+00
  %182 = load ptr, ptr %6, align 8, !tbaa !5
  %183 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %182, i32 0, i32 13
  %184 = load double, ptr %183, align 8, !tbaa !17
  %185 = fdiv double %181, %184
  store double %185, ptr %19, align 8, !tbaa !15
  %186 = load ptr, ptr %6, align 8, !tbaa !5
  %187 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %186, i32 0, i32 11
  %188 = load double, ptr %187, align 8, !tbaa !14
  %189 = fadd double 1.000000e+00, %188
  %190 = load ptr, ptr %6, align 8, !tbaa !5
  %191 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %190, i32 0, i32 13
  %192 = load double, ptr %191, align 8, !tbaa !17
  %193 = fneg double %192
  %194 = load i64, ptr %11, align 8, !tbaa !16
  %195 = sub nsw i64 %194, 2
  %196 = sitofp i64 %195 to double
  %197 = fmul double %193, %196
  %198 = call double @exp(double noundef %197) #5
  %199 = call double @llvm.fmuladd.f64(double %189, double %198, double -1.000000e+00)
  store double %199, ptr %16, align 8, !tbaa !15
  %200 = load ptr, ptr %9, align 8, !tbaa !5
  %201 = load i64, ptr %11, align 8, !tbaa !16
  %202 = sub nsw i64 %201, 1
  %203 = getelementptr inbounds double, ptr %200, i64 %202
  %204 = load double, ptr %203, align 8, !tbaa !15
  %205 = load ptr, ptr %9, align 8, !tbaa !5
  %206 = load i64, ptr %11, align 8, !tbaa !16
  %207 = sub nsw i64 %206, 3
  %208 = getelementptr inbounds double, ptr %205, i64 %207
  %209 = load double, ptr %208, align 8, !tbaa !15
  %210 = fsub double %204, %209
  %211 = fdiv double %210, 2.000000e+00
  %212 = load ptr, ptr %6, align 8, !tbaa !5
  %213 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %212, i32 0, i32 13
  %214 = load double, ptr %213, align 8, !tbaa !17
  %215 = fdiv double %211, %214
  store double %215, ptr %17, align 8, !tbaa !15
  store double 1.000000e+00, ptr %22, align 8, !tbaa !15
  br label %216

216:                                              ; preds = %338, %155
  %217 = load i64, ptr %11, align 8, !tbaa !16
  %218 = load ptr, ptr %6, align 8, !tbaa !5
  %219 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %218, i32 0, i32 14
  %220 = load i64, ptr %219, align 8, !tbaa !9
  %221 = icmp slt i64 %217, %220
  br i1 %221, label %222, label %230

222:                                              ; preds = %216
  %223 = load double, ptr %22, align 8, !tbaa !15
  %224 = fcmp ogt double %223, 1.000000e-04
  br i1 %224, label %228, label %225

225:                                              ; preds = %222
  %226 = load double, ptr %15, align 8, !tbaa !15
  %227 = fcmp ogt double %226, 1.650000e+03
  br label %228

228:                                              ; preds = %225, %222
  %229 = phi i1 [ true, %222 ], [ %227, %225 ]
  br label %230

230:                                              ; preds = %228, %216
  %231 = phi i1 [ false, %216 ], [ %229, %228 ]
  br i1 %231, label %232, label %341

232:                                              ; preds = %230
  %233 = load ptr, ptr %6, align 8, !tbaa !5
  %234 = load double, ptr %15, align 8, !tbaa !15
  %235 = load ptr, ptr %9, align 8, !tbaa !5
  %236 = load i64, ptr %11, align 8, !tbaa !16
  %237 = sub nsw i64 %236, 1
  %238 = getelementptr inbounds double, ptr %235, i64 %237
  %239 = load double, ptr %238, align 8, !tbaa !15
  %240 = load ptr, ptr %9, align 8, !tbaa !5
  %241 = load i64, ptr %11, align 8, !tbaa !16
  %242 = getelementptr inbounds double, ptr %240, i64 %241
  %243 = load ptr, ptr %7, align 8, !tbaa !5
  %244 = load i64, ptr %11, align 8, !tbaa !16
  %245 = sub nsw i64 %244, 1
  %246 = trunc i64 %245 to i32
  %247 = load ptr, ptr %8, align 8, !tbaa !5
  %248 = load ptr, ptr %12, align 8, !tbaa !5
  %249 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  call void @rec_get_xe_next1(ptr noundef %233, double noundef %234, double noundef %239, ptr noundef %242, ptr noundef %243, i32 noundef 1, i32 noundef %246, ptr noundef %247, ptr noundef %248, ptr noundef %249, ptr noundef %16, ptr noundef %17, ptr noundef %18, ptr noundef %19)
  %250 = load ptr, ptr %6, align 8, !tbaa !5
  %251 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %250, i32 0, i32 11
  %252 = load double, ptr %251, align 8, !tbaa !14
  %253 = fadd double 1.000000e+00, %252
  %254 = load ptr, ptr %6, align 8, !tbaa !5
  %255 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %254, i32 0, i32 13
  %256 = load double, ptr %255, align 8, !tbaa !17
  %257 = fneg double %256
  %258 = load i64, ptr %11, align 8, !tbaa !16
  %259 = sitofp i64 %258 to double
  %260 = fmul double %257, %259
  %261 = call double @exp(double noundef %260) #5
  %262 = call double @llvm.fmuladd.f64(double %253, double %261, double -1.000000e+00)
  store double %262, ptr %15, align 8, !tbaa !15
  %263 = load ptr, ptr %9, align 8, !tbaa !5
  %264 = load i64, ptr %11, align 8, !tbaa !16
  %265 = getelementptr inbounds double, ptr %263, i64 %264
  %266 = load double, ptr %265, align 8, !tbaa !15
  %267 = load ptr, ptr %6, align 8, !tbaa !5
  %268 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %267, i32 0, i32 0
  %269 = load double, ptr %268, align 8, !tbaa !19
  %270 = load double, ptr %15, align 8, !tbaa !15
  %271 = fadd double 1.000000e+00, %270
  %272 = fmul double %269, %271
  %273 = load ptr, ptr %6, align 8, !tbaa !5
  %274 = load double, ptr %15, align 8, !tbaa !15
  %275 = call double @rec_HubbleConstant(ptr noundef %273, double noundef %274)
  %276 = load ptr, ptr %6, align 8, !tbaa !5
  %277 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %276, i32 0, i32 10
  %278 = load double, ptr %277, align 8, !tbaa !20
  %279 = load ptr, ptr %6, align 8, !tbaa !5
  %280 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %279, i32 0, i32 9
  %281 = load double, ptr %280, align 8, !tbaa !18
  %282 = load double, ptr %15, align 8, !tbaa !15
  %283 = fadd double 1.000000e+00, %282
  %284 = call double @cube(double noundef %283)
  %285 = fmul double %281, %284
  %286 = load ptr, ptr %6, align 8, !tbaa !5
  %287 = load double, ptr %15, align 8, !tbaa !15
  %288 = call double @energy_injection_rate(ptr noundef %286, double noundef %287)
  %289 = call double @rec_Tmss(double noundef %266, double noundef %272, double noundef %275, double noundef %278, double noundef %285, double noundef %288)
  %290 = load ptr, ptr %10, align 8, !tbaa !5
  %291 = load i64, ptr %11, align 8, !tbaa !16
  %292 = getelementptr inbounds double, ptr %290, i64 %291
  store double %289, ptr %292, align 8, !tbaa !15
  %293 = load ptr, ptr %12, align 8, !tbaa !5
  %294 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  %295 = load ptr, ptr %9, align 8, !tbaa !5
  %296 = load i64, ptr %11, align 8, !tbaa !16
  %297 = getelementptr inbounds double, ptr %295, i64 %296
  %298 = load double, ptr %297, align 8, !tbaa !15
  %299 = load ptr, ptr %6, align 8, !tbaa !5
  %300 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %299, i32 0, i32 0
  %301 = load double, ptr %300, align 8, !tbaa !19
  %302 = load double, ptr %15, align 8, !tbaa !15
  %303 = fadd double 1.000000e+00, %302
  %304 = fmul double %301, %303
  %305 = fmul double %304, 0x3F16970041294F34
  %306 = load ptr, ptr %6, align 8, !tbaa !5
  %307 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %306, i32 0, i32 9
  %308 = load double, ptr %307, align 8, !tbaa !18
  %309 = load double, ptr %15, align 8, !tbaa !15
  %310 = fadd double 1.000000e+00, %309
  %311 = call double @cube(double noundef %310)
  %312 = fmul double %308, %311
  %313 = fmul double %312, 0x3EB0C6F7A0B5ED8D
  %314 = load ptr, ptr %8, align 8, !tbaa !5
  %315 = load ptr, ptr %6, align 8, !tbaa !5
  %316 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %315, i32 0, i32 11
  %317 = load double, ptr %316, align 8, !tbaa !14
  %318 = load ptr, ptr %6, align 8, !tbaa !5
  %319 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %318, i32 0, i32 13
  %320 = load double, ptr %319, align 8, !tbaa !17
  %321 = load i64, ptr %11, align 8, !tbaa !16
  %322 = trunc i64 %321 to i32
  %323 = load double, ptr %15, align 8, !tbaa !15
  call void @update_fminus_Saha(ptr noundef %293, ptr noundef %294, double noundef %298, double noundef %305, double noundef %313, ptr noundef %314, double noundef %317, double noundef %320, i32 noundef %322, double noundef %323, i32 noundef 0)
  %324 = load ptr, ptr %9, align 8, !tbaa !5
  %325 = load i64, ptr %11, align 8, !tbaa !16
  %326 = getelementptr inbounds double, ptr %324, i64 %325
  %327 = load double, ptr %326, align 8, !tbaa !15
  %328 = load ptr, ptr %6, align 8, !tbaa !5
  %329 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %328, i32 0, i32 9
  %330 = load double, ptr %329, align 8, !tbaa !18
  %331 = load ptr, ptr %6, align 8, !tbaa !5
  %332 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %331, i32 0, i32 0
  %333 = load double, ptr %332, align 8, !tbaa !19
  %334 = load double, ptr %15, align 8, !tbaa !15
  %335 = call double @rec_saha_xe_H(double noundef %330, double noundef %333, double noundef %334)
  %336 = fsub double %327, %335
  %337 = call double @llvm.fabs.f64(double %336)
  store double %337, ptr %22, align 8, !tbaa !15
  br label %338

338:                                              ; preds = %232
  %339 = load i64, ptr %11, align 8, !tbaa !16
  %340 = add nsw i64 %339, 1
  store i64 %340, ptr %11, align 8, !tbaa !16
  br label %216, !llvm.loop !24

341:                                              ; preds = %230
  store double 0.000000e+00, ptr %22, align 8, !tbaa !15
  br label %342

342:                                              ; preds = %434, %341
  %343 = load i64, ptr %11, align 8, !tbaa !16
  %344 = load ptr, ptr %6, align 8, !tbaa !5
  %345 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %344, i32 0, i32 14
  %346 = load i64, ptr %345, align 8, !tbaa !9
  %347 = icmp slt i64 %343, %346
  br i1 %347, label %348, label %351

348:                                              ; preds = %342
  %349 = load double, ptr %22, align 8, !tbaa !15
  %350 = fcmp olt double %349, 5.000000e-05
  br label %351

351:                                              ; preds = %348, %342
  %352 = phi i1 [ false, %342 ], [ %350, %348 ]
  br i1 %352, label %353, label %437

353:                                              ; preds = %351
  %354 = load ptr, ptr %6, align 8, !tbaa !5
  %355 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %354, i32 0, i32 11
  %356 = load double, ptr %355, align 8, !tbaa !14
  %357 = fadd double 1.000000e+00, %356
  %358 = load ptr, ptr %6, align 8, !tbaa !5
  %359 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %358, i32 0, i32 13
  %360 = load double, ptr %359, align 8, !tbaa !17
  %361 = fneg double %360
  %362 = load i64, ptr %11, align 8, !tbaa !16
  %363 = sitofp i64 %362 to double
  %364 = fmul double %361, %363
  %365 = call double @exp(double noundef %364) #5
  %366 = call double @llvm.fmuladd.f64(double %357, double %365, double -1.000000e+00)
  store double %366, ptr %15, align 8, !tbaa !15
  %367 = load ptr, ptr %6, align 8, !tbaa !5
  %368 = load double, ptr %15, align 8, !tbaa !15
  %369 = call double @rec_HubbleConstant(ptr noundef %367, double noundef %368)
  store double %369, ptr %14, align 8, !tbaa !15
  %370 = load ptr, ptr %6, align 8, !tbaa !5
  %371 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %370, i32 0, i32 9
  %372 = load double, ptr %371, align 8, !tbaa !18
  %373 = load double, ptr %15, align 8, !tbaa !15
  %374 = fadd double 1.000000e+00, %373
  %375 = call double @cube(double noundef %374)
  %376 = fmul double %372, %375
  %377 = fmul double %376, 0x3EB0C6F7A0B5ED8D
  %378 = load double, ptr %14, align 8, !tbaa !15
  %379 = load ptr, ptr %6, align 8, !tbaa !5
  %380 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %379, i32 0, i32 0
  %381 = load double, ptr %380, align 8, !tbaa !19
  %382 = fmul double 0x3F16970041294F34, %381
  %383 = load double, ptr %15, align 8, !tbaa !15
  %384 = fadd double 1.000000e+00, %383
  %385 = fmul double %382, %384
  %386 = load ptr, ptr %7, align 8, !tbaa !5
  %387 = load ptr, ptr %8, align 8, !tbaa !5
  %388 = load ptr, ptr %6, align 8, !tbaa !5
  %389 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %388, i32 0, i32 11
  %390 = load double, ptr %389, align 8, !tbaa !14
  %391 = load ptr, ptr %6, align 8, !tbaa !5
  %392 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %391, i32 0, i32 13
  %393 = load double, ptr %392, align 8, !tbaa !17
  %394 = load ptr, ptr %12, align 8, !tbaa !5
  %395 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  %396 = load i64, ptr %11, align 8, !tbaa !16
  %397 = trunc i64 %396 to i32
  %398 = load double, ptr %15, align 8, !tbaa !15
  %399 = load ptr, ptr %6, align 8, !tbaa !5
  %400 = load double, ptr %15, align 8, !tbaa !15
  %401 = call double @energy_injection_rate(ptr noundef %399, double noundef %400)
  %402 = call double @xe_PostSahaH(double noundef %377, double noundef %378, double noundef %385, ptr noundef %386, ptr noundef %387, double noundef %390, double noundef %393, ptr noundef %394, ptr noundef %395, i32 noundef %397, double noundef %398, ptr noundef %22, i32 noundef 1, double noundef %401)
  %403 = load ptr, ptr %9, align 8, !tbaa !5
  %404 = load i64, ptr %11, align 8, !tbaa !16
  %405 = getelementptr inbounds double, ptr %403, i64 %404
  store double %402, ptr %405, align 8, !tbaa !15
  %406 = load ptr, ptr %9, align 8, !tbaa !5
  %407 = load i64, ptr %11, align 8, !tbaa !16
  %408 = getelementptr inbounds double, ptr %406, i64 %407
  %409 = load double, ptr %408, align 8, !tbaa !15
  %410 = load ptr, ptr %6, align 8, !tbaa !5
  %411 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %410, i32 0, i32 0
  %412 = load double, ptr %411, align 8, !tbaa !19
  %413 = load double, ptr %15, align 8, !tbaa !15
  %414 = fadd double 1.000000e+00, %413
  %415 = fmul double %412, %414
  %416 = load double, ptr %14, align 8, !tbaa !15
  %417 = load ptr, ptr %6, align 8, !tbaa !5
  %418 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %417, i32 0, i32 10
  %419 = load double, ptr %418, align 8, !tbaa !20
  %420 = load ptr, ptr %6, align 8, !tbaa !5
  %421 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %420, i32 0, i32 9
  %422 = load double, ptr %421, align 8, !tbaa !18
  %423 = load double, ptr %15, align 8, !tbaa !15
  %424 = fadd double 1.000000e+00, %423
  %425 = call double @cube(double noundef %424)
  %426 = fmul double %422, %425
  %427 = load ptr, ptr %6, align 8, !tbaa !5
  %428 = load double, ptr %15, align 8, !tbaa !15
  %429 = call double @energy_injection_rate(ptr noundef %427, double noundef %428)
  %430 = call double @rec_Tmss(double noundef %409, double noundef %415, double noundef %416, double noundef %419, double noundef %426, double noundef %429)
  %431 = load ptr, ptr %10, align 8, !tbaa !5
  %432 = load i64, ptr %11, align 8, !tbaa !16
  %433 = getelementptr inbounds double, ptr %431, i64 %432
  store double %430, ptr %433, align 8, !tbaa !15
  br label %434

434:                                              ; preds = %353
  %435 = load i64, ptr %11, align 8, !tbaa !16
  %436 = add nsw i64 %435, 1
  store i64 %436, ptr %11, align 8, !tbaa !16
  br label %342, !llvm.loop !25

437:                                              ; preds = %351
  %438 = load ptr, ptr %6, align 8, !tbaa !5
  %439 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %438, i32 0, i32 11
  %440 = load double, ptr %439, align 8, !tbaa !14
  %441 = fadd double 1.000000e+00, %440
  %442 = load ptr, ptr %6, align 8, !tbaa !5
  %443 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %442, i32 0, i32 13
  %444 = load double, ptr %443, align 8, !tbaa !17
  %445 = fneg double %444
  %446 = load i64, ptr %11, align 8, !tbaa !16
  %447 = sub nsw i64 %446, 3
  %448 = sitofp i64 %447 to double
  %449 = fmul double %445, %448
  %450 = call double @exp(double noundef %449) #5
  %451 = call double @llvm.fmuladd.f64(double %441, double %450, double -1.000000e+00)
  store double %451, ptr %18, align 8, !tbaa !15
  %452 = load ptr, ptr %9, align 8, !tbaa !5
  %453 = load i64, ptr %11, align 8, !tbaa !16
  %454 = sub nsw i64 %453, 2
  %455 = getelementptr inbounds double, ptr %452, i64 %454
  %456 = load double, ptr %455, align 8, !tbaa !15
  %457 = load ptr, ptr %9, align 8, !tbaa !5
  %458 = load i64, ptr %11, align 8, !tbaa !16
  %459 = sub nsw i64 %458, 4
  %460 = getelementptr inbounds double, ptr %457, i64 %459
  %461 = load double, ptr %460, align 8, !tbaa !15
  %462 = fsub double %456, %461
  %463 = fdiv double %462, 2.000000e+00
  %464 = load ptr, ptr %6, align 8, !tbaa !5
  %465 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %464, i32 0, i32 13
  %466 = load double, ptr %465, align 8, !tbaa !17
  %467 = fdiv double %463, %466
  store double %467, ptr %19, align 8, !tbaa !15
  %468 = load ptr, ptr %6, align 8, !tbaa !5
  %469 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %468, i32 0, i32 11
  %470 = load double, ptr %469, align 8, !tbaa !14
  %471 = fadd double 1.000000e+00, %470
  %472 = load ptr, ptr %6, align 8, !tbaa !5
  %473 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %472, i32 0, i32 13
  %474 = load double, ptr %473, align 8, !tbaa !17
  %475 = fneg double %474
  %476 = load i64, ptr %11, align 8, !tbaa !16
  %477 = sub nsw i64 %476, 2
  %478 = sitofp i64 %477 to double
  %479 = fmul double %475, %478
  %480 = call double @exp(double noundef %479) #5
  %481 = call double @llvm.fmuladd.f64(double %471, double %480, double -1.000000e+00)
  store double %481, ptr %16, align 8, !tbaa !15
  %482 = load ptr, ptr %9, align 8, !tbaa !5
  %483 = load i64, ptr %11, align 8, !tbaa !16
  %484 = sub nsw i64 %483, 1
  %485 = getelementptr inbounds double, ptr %482, i64 %484
  %486 = load double, ptr %485, align 8, !tbaa !15
  %487 = load ptr, ptr %9, align 8, !tbaa !5
  %488 = load i64, ptr %11, align 8, !tbaa !16
  %489 = sub nsw i64 %488, 3
  %490 = getelementptr inbounds double, ptr %487, i64 %489
  %491 = load double, ptr %490, align 8, !tbaa !15
  %492 = fsub double %486, %491
  %493 = fdiv double %492, 2.000000e+00
  %494 = load ptr, ptr %6, align 8, !tbaa !5
  %495 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %494, i32 0, i32 13
  %496 = load double, ptr %495, align 8, !tbaa !17
  %497 = fdiv double %493, %496
  store double %497, ptr %17, align 8, !tbaa !15
  br label %498

498:                                              ; preds = %586, %437
  %499 = load i64, ptr %11, align 8, !tbaa !16
  %500 = load ptr, ptr %6, align 8, !tbaa !5
  %501 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %500, i32 0, i32 14
  %502 = load i64, ptr %501, align 8, !tbaa !9
  %503 = icmp slt i64 %499, %502
  br i1 %503, label %504, label %522

504:                                              ; preds = %498
  %505 = load ptr, ptr %10, align 8, !tbaa !5
  %506 = load i64, ptr %11, align 8, !tbaa !16
  %507 = sub nsw i64 %506, 1
  %508 = getelementptr inbounds double, ptr %505, i64 %507
  %509 = load double, ptr %508, align 8, !tbaa !15
  %510 = load ptr, ptr %6, align 8, !tbaa !5
  %511 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %510, i32 0, i32 0
  %512 = load double, ptr %511, align 8, !tbaa !19
  %513 = fdiv double %509, %512
  %514 = load double, ptr %15, align 8, !tbaa !15
  %515 = fadd double 1.000000e+00, %514
  %516 = fdiv double %513, %515
  %517 = fsub double 1.000000e+00, %516
  %518 = fcmp olt double %517, 5.000000e-04
  br i1 %518, label %519, label %522

519:                                              ; preds = %504
  %520 = load double, ptr %15, align 8, !tbaa !15
  %521 = fcmp ogt double %520, 7.000000e+02
  br label %522

522:                                              ; preds = %519, %504, %498
  %523 = phi i1 [ false, %504 ], [ false, %498 ], [ %521, %519 ]
  br i1 %523, label %524, label %589

524:                                              ; preds = %522
  %525 = load ptr, ptr %6, align 8, !tbaa !5
  %526 = load double, ptr %15, align 8, !tbaa !15
  %527 = load ptr, ptr %9, align 8, !tbaa !5
  %528 = load i64, ptr %11, align 8, !tbaa !16
  %529 = sub nsw i64 %528, 1
  %530 = getelementptr inbounds double, ptr %527, i64 %529
  %531 = load double, ptr %530, align 8, !tbaa !15
  %532 = load ptr, ptr %9, align 8, !tbaa !5
  %533 = load i64, ptr %11, align 8, !tbaa !16
  %534 = getelementptr inbounds double, ptr %532, i64 %533
  %535 = load ptr, ptr %7, align 8, !tbaa !5
  %536 = load i64, ptr %11, align 8, !tbaa !16
  %537 = sub nsw i64 %536, 1
  %538 = trunc i64 %537 to i32
  %539 = load ptr, ptr %8, align 8, !tbaa !5
  %540 = load ptr, ptr %12, align 8, !tbaa !5
  %541 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  call void @rec_get_xe_next1(ptr noundef %525, double noundef %526, double noundef %531, ptr noundef %534, ptr noundef %535, i32 noundef 2, i32 noundef %538, ptr noundef %539, ptr noundef %540, ptr noundef %541, ptr noundef %16, ptr noundef %17, ptr noundef %18, ptr noundef %19)
  %542 = load ptr, ptr %6, align 8, !tbaa !5
  %543 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %542, i32 0, i32 11
  %544 = load double, ptr %543, align 8, !tbaa !14
  %545 = fadd double 1.000000e+00, %544
  %546 = load ptr, ptr %6, align 8, !tbaa !5
  %547 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %546, i32 0, i32 13
  %548 = load double, ptr %547, align 8, !tbaa !17
  %549 = fneg double %548
  %550 = load i64, ptr %11, align 8, !tbaa !16
  %551 = sitofp i64 %550 to double
  %552 = fmul double %549, %551
  %553 = call double @exp(double noundef %552) #5
  %554 = call double @llvm.fmuladd.f64(double %545, double %553, double -1.000000e+00)
  store double %554, ptr %15, align 8, !tbaa !15
  %555 = load ptr, ptr %9, align 8, !tbaa !5
  %556 = load i64, ptr %11, align 8, !tbaa !16
  %557 = getelementptr inbounds double, ptr %555, i64 %556
  %558 = load double, ptr %557, align 8, !tbaa !15
  %559 = load ptr, ptr %6, align 8, !tbaa !5
  %560 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %559, i32 0, i32 0
  %561 = load double, ptr %560, align 8, !tbaa !19
  %562 = load double, ptr %15, align 8, !tbaa !15
  %563 = fadd double 1.000000e+00, %562
  %564 = fmul double %561, %563
  %565 = load ptr, ptr %6, align 8, !tbaa !5
  %566 = load double, ptr %15, align 8, !tbaa !15
  %567 = call double @rec_HubbleConstant(ptr noundef %565, double noundef %566)
  %568 = load ptr, ptr %6, align 8, !tbaa !5
  %569 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %568, i32 0, i32 10
  %570 = load double, ptr %569, align 8, !tbaa !20
  %571 = load ptr, ptr %6, align 8, !tbaa !5
  %572 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %571, i32 0, i32 9
  %573 = load double, ptr %572, align 8, !tbaa !18
  %574 = load double, ptr %15, align 8, !tbaa !15
  %575 = fadd double 1.000000e+00, %574
  %576 = call double @cube(double noundef %575)
  %577 = fmul double %573, %576
  %578 = fmul double %577, 0x3EB0C6F7A0B5ED8D
  %579 = load ptr, ptr %6, align 8, !tbaa !5
  %580 = load double, ptr %15, align 8, !tbaa !15
  %581 = call double @energy_injection_rate(ptr noundef %579, double noundef %580)
  %582 = call double @rec_Tmss(double noundef %558, double noundef %564, double noundef %567, double noundef %570, double noundef %578, double noundef %581)
  %583 = load ptr, ptr %10, align 8, !tbaa !5
  %584 = load i64, ptr %11, align 8, !tbaa !16
  %585 = getelementptr inbounds double, ptr %583, i64 %584
  store double %582, ptr %585, align 8, !tbaa !15
  br label %586

586:                                              ; preds = %524
  %587 = load i64, ptr %11, align 8, !tbaa !16
  %588 = add nsw i64 %587, 1
  store i64 %588, ptr %11, align 8, !tbaa !16
  br label %498, !llvm.loop !26

589:                                              ; preds = %522
  %590 = load ptr, ptr %9, align 8, !tbaa !5
  %591 = load i64, ptr %11, align 8, !tbaa !16
  %592 = sub nsw i64 %591, 3
  %593 = getelementptr inbounds double, ptr %590, i64 %592
  %594 = load double, ptr %593, align 8, !tbaa !15
  %595 = load ptr, ptr %10, align 8, !tbaa !5
  %596 = load i64, ptr %11, align 8, !tbaa !16
  %597 = sub nsw i64 %596, 3
  %598 = getelementptr inbounds double, ptr %595, i64 %597
  %599 = load double, ptr %598, align 8, !tbaa !15
  %600 = load ptr, ptr %6, align 8, !tbaa !5
  %601 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %600, i32 0, i32 0
  %602 = load double, ptr %601, align 8, !tbaa !19
  %603 = load double, ptr %18, align 8, !tbaa !15
  %604 = fadd double 1.000000e+00, %603
  %605 = fmul double %602, %604
  %606 = load ptr, ptr %6, align 8, !tbaa !5
  %607 = load double, ptr %18, align 8, !tbaa !15
  %608 = call double @rec_HubbleConstant(ptr noundef %606, double noundef %607)
  %609 = load ptr, ptr %6, align 8, !tbaa !5
  %610 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %609, i32 0, i32 10
  %611 = load double, ptr %610, align 8, !tbaa !20
  %612 = load ptr, ptr %6, align 8, !tbaa !5
  %613 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %612, i32 0, i32 9
  %614 = load double, ptr %613, align 8, !tbaa !18
  %615 = load double, ptr %18, align 8, !tbaa !15
  %616 = fadd double 1.000000e+00, %615
  %617 = call double @cube(double noundef %616)
  %618 = fmul double %614, %617
  %619 = load ptr, ptr %6, align 8, !tbaa !5
  %620 = load double, ptr %18, align 8, !tbaa !15
  %621 = call double @energy_injection_rate(ptr noundef %619, double noundef %620)
  %622 = call double @rec_dTmdlna(double noundef %594, double noundef %599, double noundef %605, double noundef %608, double noundef %611, double noundef %618, double noundef %621)
  store double %622, ptr %21, align 8, !tbaa !15
  %623 = load ptr, ptr %9, align 8, !tbaa !5
  %624 = load i64, ptr %11, align 8, !tbaa !16
  %625 = sub nsw i64 %624, 2
  %626 = getelementptr inbounds double, ptr %623, i64 %625
  %627 = load double, ptr %626, align 8, !tbaa !15
  %628 = load ptr, ptr %10, align 8, !tbaa !5
  %629 = load i64, ptr %11, align 8, !tbaa !16
  %630 = sub nsw i64 %629, 2
  %631 = getelementptr inbounds double, ptr %628, i64 %630
  %632 = load double, ptr %631, align 8, !tbaa !15
  %633 = load ptr, ptr %6, align 8, !tbaa !5
  %634 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %633, i32 0, i32 0
  %635 = load double, ptr %634, align 8, !tbaa !19
  %636 = load double, ptr %16, align 8, !tbaa !15
  %637 = fadd double 1.000000e+00, %636
  %638 = fmul double %635, %637
  %639 = load ptr, ptr %6, align 8, !tbaa !5
  %640 = load double, ptr %16, align 8, !tbaa !15
  %641 = call double @rec_HubbleConstant(ptr noundef %639, double noundef %640)
  %642 = load ptr, ptr %6, align 8, !tbaa !5
  %643 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %642, i32 0, i32 10
  %644 = load double, ptr %643, align 8, !tbaa !20
  %645 = load ptr, ptr %6, align 8, !tbaa !5
  %646 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %645, i32 0, i32 9
  %647 = load double, ptr %646, align 8, !tbaa !18
  %648 = load double, ptr %16, align 8, !tbaa !15
  %649 = fadd double 1.000000e+00, %648
  %650 = call double @cube(double noundef %649)
  %651 = fmul double %647, %650
  %652 = load ptr, ptr %6, align 8, !tbaa !5
  %653 = load double, ptr %16, align 8, !tbaa !15
  %654 = call double @energy_injection_rate(ptr noundef %652, double noundef %653)
  %655 = call double @rec_dTmdlna(double noundef %627, double noundef %632, double noundef %638, double noundef %641, double noundef %644, double noundef %651, double noundef %654)
  store double %655, ptr %20, align 8, !tbaa !15
  br label %656

656:                                              ; preds = %706, %589
  %657 = load i64, ptr %11, align 8, !tbaa !16
  %658 = load ptr, ptr %6, align 8, !tbaa !5
  %659 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %658, i32 0, i32 14
  %660 = load i64, ptr %659, align 8, !tbaa !9
  %661 = icmp slt i64 %657, %660
  br i1 %661, label %662, label %665

662:                                              ; preds = %656
  %663 = load double, ptr %15, align 8, !tbaa !15
  %664 = fcmp ogt double %663, 7.000000e+02
  br label %665

665:                                              ; preds = %662, %656
  %666 = phi i1 [ false, %656 ], [ %664, %662 ]
  br i1 %666, label %667, label %709

667:                                              ; preds = %665
  %668 = load ptr, ptr %6, align 8, !tbaa !5
  %669 = load double, ptr %15, align 8, !tbaa !15
  %670 = load ptr, ptr %9, align 8, !tbaa !5
  %671 = load i64, ptr %11, align 8, !tbaa !16
  %672 = sub nsw i64 %671, 1
  %673 = getelementptr inbounds double, ptr %670, i64 %672
  %674 = load double, ptr %673, align 8, !tbaa !15
  %675 = load ptr, ptr %10, align 8, !tbaa !5
  %676 = load i64, ptr %11, align 8, !tbaa !16
  %677 = sub nsw i64 %676, 1
  %678 = getelementptr inbounds double, ptr %675, i64 %677
  %679 = load double, ptr %678, align 8, !tbaa !15
  %680 = load ptr, ptr %9, align 8, !tbaa !5
  %681 = load i64, ptr %11, align 8, !tbaa !16
  %682 = getelementptr inbounds double, ptr %680, i64 %681
  %683 = load ptr, ptr %10, align 8, !tbaa !5
  %684 = load i64, ptr %11, align 8, !tbaa !16
  %685 = getelementptr inbounds double, ptr %683, i64 %684
  %686 = load ptr, ptr %7, align 8, !tbaa !5
  %687 = load i64, ptr %11, align 8, !tbaa !16
  %688 = sub nsw i64 %687, 1
  %689 = trunc i64 %688 to i32
  %690 = load ptr, ptr %8, align 8, !tbaa !5
  %691 = load ptr, ptr %12, align 8, !tbaa !5
  %692 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  call void @rec_get_xe_next2(ptr noundef %668, double noundef %669, double noundef %674, double noundef %679, ptr noundef %682, ptr noundef %685, ptr noundef %686, i32 noundef 2, i32 noundef %689, ptr noundef %690, ptr noundef %691, ptr noundef %692, ptr noundef %16, ptr noundef %17, ptr noundef %20, ptr noundef %18, ptr noundef %19, ptr noundef %21)
  %693 = load ptr, ptr %6, align 8, !tbaa !5
  %694 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %693, i32 0, i32 11
  %695 = load double, ptr %694, align 8, !tbaa !14
  %696 = fadd double 1.000000e+00, %695
  %697 = load ptr, ptr %6, align 8, !tbaa !5
  %698 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %697, i32 0, i32 13
  %699 = load double, ptr %698, align 8, !tbaa !17
  %700 = fneg double %699
  %701 = load i64, ptr %11, align 8, !tbaa !16
  %702 = sitofp i64 %701 to double
  %703 = fmul double %700, %702
  %704 = call double @exp(double noundef %703) #5
  %705 = call double @llvm.fmuladd.f64(double %696, double %704, double -1.000000e+00)
  store double %705, ptr %15, align 8, !tbaa !15
  br label %706

706:                                              ; preds = %667
  %707 = load i64, ptr %11, align 8, !tbaa !16
  %708 = add nsw i64 %707, 1
  store i64 %708, ptr %11, align 8, !tbaa !16
  br label %656, !llvm.loop !27

709:                                              ; preds = %665
  br label %710

710:                                              ; preds = %760, %709
  %711 = load i64, ptr %11, align 8, !tbaa !16
  %712 = load ptr, ptr %6, align 8, !tbaa !5
  %713 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %712, i32 0, i32 14
  %714 = load i64, ptr %713, align 8, !tbaa !9
  %715 = icmp slt i64 %711, %714
  br i1 %715, label %716, label %719

716:                                              ; preds = %710
  %717 = load double, ptr %15, align 8, !tbaa !15
  %718 = fcmp ogt double %717, 2.000000e+01
  br label %719

719:                                              ; preds = %716, %710
  %720 = phi i1 [ false, %710 ], [ %718, %716 ]
  br i1 %720, label %721, label %763

721:                                              ; preds = %719
  %722 = load ptr, ptr %6, align 8, !tbaa !5
  %723 = load double, ptr %15, align 8, !tbaa !15
  %724 = load ptr, ptr %9, align 8, !tbaa !5
  %725 = load i64, ptr %11, align 8, !tbaa !16
  %726 = sub nsw i64 %725, 1
  %727 = getelementptr inbounds double, ptr %724, i64 %726
  %728 = load double, ptr %727, align 8, !tbaa !15
  %729 = load ptr, ptr %10, align 8, !tbaa !5
  %730 = load i64, ptr %11, align 8, !tbaa !16
  %731 = sub nsw i64 %730, 1
  %732 = getelementptr inbounds double, ptr %729, i64 %731
  %733 = load double, ptr %732, align 8, !tbaa !15
  %734 = load ptr, ptr %9, align 8, !tbaa !5
  %735 = load i64, ptr %11, align 8, !tbaa !16
  %736 = getelementptr inbounds double, ptr %734, i64 %735
  %737 = load ptr, ptr %10, align 8, !tbaa !5
  %738 = load i64, ptr %11, align 8, !tbaa !16
  %739 = getelementptr inbounds double, ptr %737, i64 %738
  %740 = load ptr, ptr %7, align 8, !tbaa !5
  %741 = load i64, ptr %11, align 8, !tbaa !16
  %742 = sub nsw i64 %741, 1
  %743 = trunc i64 %742 to i32
  %744 = load ptr, ptr %8, align 8, !tbaa !5
  %745 = load ptr, ptr %12, align 8, !tbaa !5
  %746 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  call void @rec_get_xe_next2(ptr noundef %722, double noundef %723, double noundef %728, double noundef %733, ptr noundef %736, ptr noundef %739, ptr noundef %740, i32 noundef 3, i32 noundef %743, ptr noundef %744, ptr noundef %745, ptr noundef %746, ptr noundef %16, ptr noundef %17, ptr noundef %20, ptr noundef %18, ptr noundef %19, ptr noundef %21)
  %747 = load ptr, ptr %6, align 8, !tbaa !5
  %748 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %747, i32 0, i32 11
  %749 = load double, ptr %748, align 8, !tbaa !14
  %750 = fadd double 1.000000e+00, %749
  %751 = load ptr, ptr %6, align 8, !tbaa !5
  %752 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %751, i32 0, i32 13
  %753 = load double, ptr %752, align 8, !tbaa !17
  %754 = fneg double %753
  %755 = load i64, ptr %11, align 8, !tbaa !16
  %756 = sitofp i64 %755 to double
  %757 = fmul double %754, %756
  %758 = call double @exp(double noundef %757) #5
  %759 = call double @llvm.fmuladd.f64(double %750, double %758, double -1.000000e+00)
  store double %759, ptr %15, align 8, !tbaa !15
  br label %760

760:                                              ; preds = %721
  %761 = load i64, ptr %11, align 8, !tbaa !16
  %762 = add nsw i64 %761, 1
  store i64 %762, ptr %11, align 8, !tbaa !16
  br label %710, !llvm.loop !28

763:                                              ; preds = %719
  br label %764

764:                                              ; preds = %809, %763
  %765 = load i64, ptr %11, align 8, !tbaa !16
  %766 = load ptr, ptr %6, align 8, !tbaa !5
  %767 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %766, i32 0, i32 14
  %768 = load i64, ptr %767, align 8, !tbaa !9
  %769 = icmp slt i64 %765, %768
  br i1 %769, label %770, label %812

770:                                              ; preds = %764
  %771 = load ptr, ptr %6, align 8, !tbaa !5
  %772 = load double, ptr %15, align 8, !tbaa !15
  %773 = load ptr, ptr %9, align 8, !tbaa !5
  %774 = load i64, ptr %11, align 8, !tbaa !16
  %775 = sub nsw i64 %774, 1
  %776 = getelementptr inbounds double, ptr %773, i64 %775
  %777 = load double, ptr %776, align 8, !tbaa !15
  %778 = load ptr, ptr %10, align 8, !tbaa !5
  %779 = load i64, ptr %11, align 8, !tbaa !16
  %780 = sub nsw i64 %779, 1
  %781 = getelementptr inbounds double, ptr %778, i64 %780
  %782 = load double, ptr %781, align 8, !tbaa !15
  %783 = load ptr, ptr %9, align 8, !tbaa !5
  %784 = load i64, ptr %11, align 8, !tbaa !16
  %785 = getelementptr inbounds double, ptr %783, i64 %784
  %786 = load ptr, ptr %10, align 8, !tbaa !5
  %787 = load i64, ptr %11, align 8, !tbaa !16
  %788 = getelementptr inbounds double, ptr %786, i64 %787
  %789 = load ptr, ptr %7, align 8, !tbaa !5
  %790 = load i64, ptr %11, align 8, !tbaa !16
  %791 = sub nsw i64 %790, 1
  %792 = trunc i64 %791 to i32
  %793 = load ptr, ptr %8, align 8, !tbaa !5
  %794 = load ptr, ptr %12, align 8, !tbaa !5
  %795 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  call void @rec_get_xe_next2(ptr noundef %771, double noundef %772, double noundef %777, double noundef %782, ptr noundef %785, ptr noundef %788, ptr noundef %789, i32 noundef 4, i32 noundef %792, ptr noundef %793, ptr noundef %794, ptr noundef %795, ptr noundef %16, ptr noundef %17, ptr noundef %20, ptr noundef %18, ptr noundef %19, ptr noundef %21)
  %796 = load ptr, ptr %6, align 8, !tbaa !5
  %797 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %796, i32 0, i32 11
  %798 = load double, ptr %797, align 8, !tbaa !14
  %799 = fadd double 1.000000e+00, %798
  %800 = load ptr, ptr %6, align 8, !tbaa !5
  %801 = getelementptr inbounds %struct.REC_COSMOPARAMS, ptr %800, i32 0, i32 13
  %802 = load double, ptr %801, align 8, !tbaa !17
  %803 = fneg double %802
  %804 = load i64, ptr %11, align 8, !tbaa !16
  %805 = sitofp i64 %804 to double
  %806 = fmul double %803, %805
  %807 = call double @exp(double noundef %806) #5
  %808 = call double @llvm.fmuladd.f64(double %799, double %807, double -1.000000e+00)
  store double %808, ptr %15, align 8, !tbaa !15
  br label %809

809:                                              ; preds = %770
  %810 = load i64, ptr %11, align 8, !tbaa !16
  %811 = add nsw i64 %810, 1
  store i64 %811, ptr %11, align 8, !tbaa !16
  br label %764, !llvm.loop !29

812:                                              ; preds = %764
  %813 = load ptr, ptr %12, align 8, !tbaa !5
  call void @free_2D_array(ptr noundef %813, i32 noundef 311)
  %814 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 0
  %815 = load ptr, ptr %814, align 16, !tbaa !5
  call void @free(ptr noundef %815) #5
  %816 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 1
  %817 = load ptr, ptr %816, align 8, !tbaa !5
  call void @free(ptr noundef %817) #5
  %818 = getelementptr inbounds [3 x ptr], ptr %13, i64 0, i64 2
  %819 = load ptr, ptr %818, align 16, !tbaa !5
  call void @free(ptr noundef %819) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %22) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %20) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %19) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %18) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %17) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %14) #5
  call void @llvm.lifetime.end.p0(i64 24, ptr %13) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %11) #5
  ret void
}

declare ptr @create_2D_array(i32 noundef, i32 noundef) #4

declare ptr @create_1D_array(i32 noundef) #4

declare double @rec_sahaHeII(double noundef, double noundef, double noundef, double noundef, ptr noundef) #4

declare double @xe_PostSahaHe(double noundef, double noundef, double noundef, double noundef, double noundef, ptr noundef) #4

declare double @cube(double noundef) #4

declare void @update_fminus_Saha(ptr noundef, ptr noundef, double noundef, double noundef, double noundef, ptr noundef, double noundef, double noundef, i32 noundef, double noundef, i32 noundef) #4

declare double @rec_saha_xe_H(double noundef, double noundef, double noundef) #4

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.fabs.f64(double) #3

declare double @xe_PostSahaH(double noundef, double noundef, double noundef, ptr noundef, ptr noundef, double noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, i32 noundef, double noundef) #4

declare void @free_2D_array(ptr noundef, i32 noundef) #4

; Function Attrs: nounwind
declare void @free(ptr noundef) #2

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { nounwind "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #4 = { "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #5 = { nounwind }

!llvm.module.flags = !{!0, !1, !2, !3}
!llvm.ident = !{!4}

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"any pointer", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !12, i64 112}
!10 = !{!"", !11, i64 0, !11, i64 8, !11, i64 16, !11, i64 24, !11, i64 32, !11, i64 40, !11, i64 48, !11, i64 56, !11, i64 64, !11, i64 72, !11, i64 80, !11, i64 88, !11, i64 96, !11, i64 104, !12, i64 112, !11, i64 120, !13, i64 128, !11, i64 136, !11, i64 144, !11, i64 152, !11, i64 160, !11, i64 168, !11, i64 176, !11, i64 184}
!11 = !{!"double", !7, i64 0}
!12 = !{!"long", !7, i64 0}
!13 = !{!"short", !7, i64 0}
!14 = !{!10, !11, i64 88}
!15 = !{!11, !11, i64 0}
!16 = !{!12, !12, i64 0}
!17 = !{!10, !11, i64 104}
!18 = !{!10, !11, i64 72}
!19 = !{!10, !11, i64 0}
!20 = !{!10, !11, i64 80}
!21 = distinct !{!21, !22}
!22 = !{!"llvm.loop.mustprogress"}
!23 = distinct !{!23, !22}
!24 = distinct !{!24, !22}
!25 = distinct !{!25, !22}
!26 = distinct !{!26, !22}
!27 = distinct !{!27, !22}
!28 = distinct !{!28, !22}
!29 = distinct !{!29, !22}


; Optimized IR (opt -passes='default<O3>'):
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-i128:128-f80:128-n8:16:32:64-S128"
target triple = "x86_64-unknown-linux-gnu"

; Function Attrs: mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.start.p0(i64 immarg, ptr nocapture) #0

; Function Attrs: mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.end.p0(i64 immarg, ptr nocapture) #0

; Function Attrs: nounwind uwtable
declare dso_local double @rec_HubbleConstant(ptr noundef, double noundef) local_unnamed_addr #1

; Function Attrs: mustprogress nofree nounwind willreturn memory(write)
declare double @exp(double noundef) local_unnamed_addr #2

; Function Attrs: nounwind uwtable
declare dso_local double @rec_Tmss(double noundef, double noundef, double noundef, double noundef, double noundef, double noundef) local_unnamed_addr #1

; Function Attrs: nounwind uwtable
declare dso_local double @rec_dTmdlna(double noundef, double noundef, double noundef, double noundef, double noundef, double noundef, double noundef) local_unnamed_addr #1

; Function Attrs: mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.fmuladd.f64(double, double, double) #3

; Function Attrs: nounwind uwtable
declare dso_local void @rec_get_xe_next1(ptr noundef, double noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, i32 noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef) local_unnamed_addr #1

; Function Attrs: nounwind uwtable
declare dso_local double @energy_injection_rate(ptr noundef, double noundef) local_unnamed_addr #1

; Function Attrs: nounwind uwtable
declare dso_local void @rec_get_xe_next2(ptr noundef, double noundef, double noundef, double noundef, ptr noundef, ptr noundef, ptr noundef, i32 noundef, i32 noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef, ptr noundef) local_unnamed_addr #1

; Function Attrs: nounwind uwtable
define dso_local void @rec_build_history(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, ptr noundef %4) local_unnamed_addr #1 {
  %6 = alloca [3 x ptr], align 16
  %7 = alloca double, align 8
  %8 = alloca double, align 8
  %9 = alloca double, align 8
  %10 = alloca double, align 8
  %11 = alloca double, align 8
  %12 = alloca double, align 8
  %13 = alloca double, align 8
  call void @llvm.lifetime.start.p0(i64 24, ptr nonnull %6) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %7) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %8) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %9) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %10) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %11) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %12) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %13) #6
  %14 = getelementptr inbounds i8, ptr %0, i64 112
  %15 = load i64, ptr %14, align 8, !tbaa !5
  %16 = trunc i64 %15 to i32
  %17 = tail call ptr @create_2D_array(i32 noundef 311, i32 noundef %16) #6
  %18 = load i64, ptr %14, align 8, !tbaa !5
  %19 = trunc i64 %18 to i32
  %20 = tail call ptr @create_1D_array(i32 noundef %19) #6
  store ptr %20, ptr %6, align 16, !tbaa !12
  %21 = load i64, ptr %14, align 8, !tbaa !5
  %22 = trunc i64 %21 to i32
  %23 = tail call ptr @create_1D_array(i32 noundef %22) #6
  %24 = getelementptr inbounds i8, ptr %6, i64 8
  store ptr %23, ptr %24, align 8, !tbaa !12
  %25 = load i64, ptr %14, align 8, !tbaa !5
  %26 = trunc i64 %25 to i32
  %27 = tail call ptr @create_1D_array(i32 noundef %26) #6
  %28 = getelementptr inbounds i8, ptr %6, i64 16
  store ptr %27, ptr %28, align 16, !tbaa !12
  %29 = getelementptr inbounds i8, ptr %0, i64 88
  %30 = load double, ptr %29, align 8, !tbaa !14
  store double 1.000000e+00, ptr %13, align 8, !tbaa !15
  %31 = load i64, ptr %14, align 8, !tbaa !5
  %32 = icmp sgt i64 %31, 0
  br i1 %32, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %5
  %33 = getelementptr inbounds i8, ptr %0, i64 104
  %34 = getelementptr inbounds i8, ptr %0, i64 72
  %35 = getelementptr inbounds i8, ptr %0, i64 80
  br label %36

36:                                               ; preds = %.lr.ph, %36
  %.0271280 = phi i64 [ 0, %.lr.ph ], [ %54, %36 ]
  %37 = load double, ptr %29, align 8, !tbaa !14
  %38 = fadd double %37, 1.000000e+00
  %39 = load double, ptr %33, align 8, !tbaa !16
  %40 = fneg double %39
  %41 = uitofp nneg i64 %.0271280 to double
  %42 = fmul double %41, %40
  %43 = call double @exp(double noundef %42) #6
  %44 = call double @llvm.fmuladd.f64(double %38, double %43, double -1.000000e+00)
  %45 = load double, ptr %34, align 8, !tbaa !17
  %46 = load double, ptr %0, align 8, !tbaa !18
  %47 = load double, ptr %35, align 8, !tbaa !19
  %48 = call double @rec_sahaHeII(double noundef %45, double noundef %46, double noundef %47, double noundef %44, ptr noundef nonnull %13) #6
  %49 = getelementptr inbounds double, ptr %3, i64 %.0271280
  store double %48, ptr %49, align 8, !tbaa !15
  %50 = load double, ptr %0, align 8, !tbaa !18
  %51 = fadd double %44, 1.000000e+00
  %52 = fmul double %51, %50
  %53 = getelementptr inbounds double, ptr %4, i64 %.0271280
  store double %52, ptr %53, align 8, !tbaa !15
  %54 = add nuw nsw i64 %.0271280, 1
  %55 = load i64, ptr %14, align 8, !tbaa !5
  %56 = icmp slt i64 %54, %55
  %57 = load double, ptr %13, align 8
  %58 = fcmp ogt double %57, 1.000000e-09
  %59 = select i1 %56, i1 %58, i1 false
  br i1 %59, label %36, label %._crit_edge, !llvm.loop !20

._crit_edge:                                      ; preds = %36, %5
  %60 = phi i64 [ %31, %5 ], [ %55, %36 ]
  %.0271.lcssa = phi i64 [ 0, %5 ], [ %54, %36 ]
  %.0.lcssa = phi double [ %30, %5 ], [ %44, %36 ]
  store double 0.000000e+00, ptr %13, align 8, !tbaa !15
  %61 = icmp slt i64 %.0271.lcssa, %60
  br i1 %61, label %.lr.ph284, label %._crit_edge285

.lr.ph284:                                        ; preds = %._crit_edge
  %62 = getelementptr inbounds i8, ptr %0, i64 104
  %63 = getelementptr inbounds i8, ptr %0, i64 72
  %64 = getelementptr inbounds i8, ptr %0, i64 80
  br label %65

65:                                               ; preds = %.lr.ph284, %65
  %.1272282 = phi i64 [ %.0271.lcssa, %.lr.ph284 ], [ %84, %65 ]
  %66 = load double, ptr %29, align 8, !tbaa !14
  %67 = fadd double %66, 1.000000e+00
  %68 = load double, ptr %62, align 8, !tbaa !16
  %69 = fneg double %68
  %70 = uitofp nneg i64 %.1272282 to double
  %71 = fmul double %70, %69
  %72 = call double @exp(double noundef %71) #6
  %73 = call double @llvm.fmuladd.f64(double %67, double %72, double -1.000000e+00)
  %74 = load double, ptr %63, align 8, !tbaa !17
  %75 = load double, ptr %0, align 8, !tbaa !18
  %76 = load double, ptr %64, align 8, !tbaa !19
  %77 = call double @rec_HubbleConstant(ptr noundef nonnull %0, double noundef %73)
  %78 = call double @xe_PostSahaHe(double noundef %74, double noundef %75, double noundef %76, double noundef %77, double noundef %73, ptr noundef nonnull %13) #6
  %79 = getelementptr inbounds double, ptr %3, i64 %.1272282
  store double %78, ptr %79, align 8, !tbaa !15
  %80 = load double, ptr %0, align 8, !tbaa !18
  %81 = fadd double %73, 1.000000e+00
  %82 = fmul double %81, %80
  %83 = getelementptr inbounds double, ptr %4, i64 %.1272282
  store double %82, ptr %83, align 8, !tbaa !15
  %84 = add nuw nsw i64 %.1272282, 1
  %85 = load i64, ptr %14, align 8, !tbaa !5
  %86 = icmp slt i64 %84, %85
  %87 = load double, ptr %13, align 8
  %88 = fcmp olt double %87, 5.000000e-04
  %89 = select i1 %86, i1 %88, i1 false
  br i1 %89, label %65, label %._crit_edge285, !llvm.loop !22

._crit_edge285:                                   ; preds = %65, %._crit_edge
  %.1272.lcssa = phi i64 [ %.0271.lcssa, %._crit_edge ], [ %84, %65 ]
  %.1.lcssa = phi double [ %.0.lcssa, %._crit_edge ], [ %73, %65 ]
  %90 = load double, ptr %29, align 8, !tbaa !14
  %91 = fadd double %90, 1.000000e+00
  %92 = getelementptr inbounds i8, ptr %0, i64 104
  %93 = load double, ptr %92, align 8, !tbaa !16
  %94 = fneg double %93
  %95 = add nsw i64 %.1272.lcssa, -3
  %96 = sitofp i64 %95 to double
  %97 = fmul double %96, %94
  %98 = call double @exp(double noundef %97) #6
  %99 = call double @llvm.fmuladd.f64(double %91, double %98, double -1.000000e+00)
  store double %99, ptr %9, align 8, !tbaa !15
  %100 = add nsw i64 %.1272.lcssa, -2
  %101 = getelementptr inbounds double, ptr %3, i64 %100
  %102 = load double, ptr %101, align 8, !tbaa !15
  %103 = getelementptr double, ptr %3, i64 %.1272.lcssa
  %104 = getelementptr i8, ptr %103, i64 -32
  %105 = load double, ptr %104, align 8, !tbaa !15
  %106 = fsub double %102, %105
  %107 = fmul double %106, 5.000000e-01
  %108 = load double, ptr %92, align 8, !tbaa !16
  %109 = fdiv double %107, %108
  store double %109, ptr %10, align 8, !tbaa !15
  %110 = load double, ptr %29, align 8, !tbaa !14
  %111 = fadd double %110, 1.000000e+00
  %112 = fneg double %108
  %113 = sitofp i64 %100 to double
  %114 = fmul double %113, %112
  %115 = call double @exp(double noundef %114) #6
  %116 = call double @llvm.fmuladd.f64(double %111, double %115, double -1.000000e+00)
  store double %116, ptr %7, align 8, !tbaa !15
  %117 = getelementptr i8, ptr %103, i64 -8
  %118 = load double, ptr %117, align 8, !tbaa !15
  %119 = getelementptr inbounds double, ptr %3, i64 %95
  %120 = load double, ptr %119, align 8, !tbaa !15
  %121 = fsub double %118, %120
  %122 = fmul double %121, 5.000000e-01
  %123 = load double, ptr %92, align 8, !tbaa !16
  %124 = fdiv double %122, %123
  store double %124, ptr %8, align 8, !tbaa !15
  store double 1.000000e+00, ptr %13, align 8, !tbaa !15
  %125 = load i64, ptr %14, align 8, !tbaa !5
  %126 = icmp slt i64 %.1272.lcssa, %125
  br i1 %126, label %.lr.ph291, label %.critedge

.lr.ph291:                                        ; preds = %._crit_edge285
  %127 = getelementptr inbounds i8, ptr %0, i64 80
  %128 = getelementptr inbounds i8, ptr %0, i64 72
  br label %133

129:                                              ; preds = %133
  %130 = fcmp ogt double %175, 1.000000e-04
  %131 = fcmp ogt double %146, 1.650000e+03
  %132 = or i1 %130, %131
  br i1 %132, label %133, label %.critedge, !llvm.loop !23

133:                                              ; preds = %.lr.ph291, %129
  %.2273288362 = phi i64 [ %.1272.lcssa, %.lr.ph291 ], [ %176, %129 ]
  %.2289361 = phi double [ %.1.lcssa, %.lr.ph291 ], [ %146, %129 ]
  %134 = add nsw i64 %.2273288362, -1
  %135 = getelementptr inbounds double, ptr %3, i64 %134
  %136 = load double, ptr %135, align 8, !tbaa !15
  %137 = getelementptr inbounds double, ptr %3, i64 %.2273288362
  %138 = trunc i64 %134 to i32
  call void @rec_get_xe_next1(ptr noundef nonnull %0, double noundef %.2289361, double noundef %136, ptr noundef %137, ptr noundef %1, i32 noundef 1, i32 noundef %138, ptr noundef %2, ptr noundef %17, ptr noundef nonnull %6, ptr noundef nonnull %7, ptr noundef nonnull %8, ptr noundef nonnull %9, ptr noundef nonnull %10)
  %139 = load double, ptr %29, align 8, !tbaa !14
  %140 = fadd double %139, 1.000000e+00
  %141 = load double, ptr %92, align 8, !tbaa !16
  %142 = fneg double %141
  %143 = uitofp nneg i64 %.2273288362 to double
  %144 = fmul double %143, %142
  %145 = call double @exp(double noundef %144) #6
  %146 = call double @llvm.fmuladd.f64(double %140, double %145, double -1.000000e+00)
  %147 = load double, ptr %137, align 8, !tbaa !15
  %148 = load double, ptr %0, align 8, !tbaa !18
  %149 = fadd double %146, 1.000000e+00
  %150 = fmul double %148, %149
  %151 = call double @rec_HubbleConstant(ptr noundef nonnull %0, double noundef %146)
  %152 = load double, ptr %127, align 8, !tbaa !19
  %153 = load double, ptr %128, align 8, !tbaa !17
  %154 = call double @cube(double noundef %149) #6
  %155 = fmul double %153, %154
  %156 = call double @energy_injection_rate(ptr noundef nonnull %0, double noundef %146)
  %157 = call double @rec_Tmss(double noundef %147, double noundef %150, double noundef %151, double noundef %152, double noundef %155, double noundef %156)
  %158 = getelementptr inbounds double, ptr %4, i64 %.2273288362
  store double %157, ptr %158, align 8, !tbaa !15
  %159 = load double, ptr %137, align 8, !tbaa !15
  %160 = load double, ptr %0, align 8, !tbaa !18
  %161 = fmul double %149, %160
  %162 = fmul double %161, 0x3F16970041294F34
  %163 = load double, ptr %128, align 8, !tbaa !17
  %164 = call double @cube(double noundef %149) #6
  %165 = fmul double %163, %164
  %166 = fmul double %165, 0x3EB0C6F7A0B5ED8D
  %167 = load double, ptr %29, align 8, !tbaa !14
  %168 = load double, ptr %92, align 8, !tbaa !16
  %169 = trunc i64 %.2273288362 to i32
  call void @update_fminus_Saha(ptr noundef %17, ptr noundef nonnull %6, double noundef %159, double noundef %162, double noundef %166, ptr noundef %2, double noundef %167, double noundef %168, i32 noundef %169, double noundef %146, i32 noundef 0) #6
  %170 = load double, ptr %137, align 8, !tbaa !15
  %171 = load double, ptr %128, align 8, !tbaa !17
  %172 = load double, ptr %0, align 8, !tbaa !18
  %173 = call double @rec_saha_xe_H(double noundef %171, double noundef %172, double noundef %146) #6
  %174 = fsub double %170, %173
  %175 = call double @llvm.fabs.f64(double %174)
  store double %175, ptr %13, align 8, !tbaa !15
  %176 = add nuw nsw i64 %.2273288362, 1
  %177 = load i64, ptr %14, align 8, !tbaa !5
  %178 = icmp slt i64 %176, %177
  br i1 %178, label %129, label %.critedge, !llvm.loop !23

.critedge:                                        ; preds = %129, %133, %._crit_edge285
  %179 = phi i64 [ %125, %._crit_edge285 ], [ %177, %133 ], [ %177, %129 ]
  %.2273.lcssa = phi i64 [ %.1272.lcssa, %._crit_edge285 ], [ %176, %133 ], [ %176, %129 ]
  %.2.lcssa = phi double [ %.1.lcssa, %._crit_edge285 ], [ %146, %133 ], [ %146, %129 ]
  store double 0.000000e+00, ptr %13, align 8, !tbaa !15
  %180 = icmp slt i64 %.2273.lcssa, %179
  br i1 %180, label %.lr.ph298, label %._crit_edge299

.lr.ph298:                                        ; preds = %.critedge
  %181 = getelementptr inbounds i8, ptr %0, i64 72
  %182 = getelementptr inbounds i8, ptr %0, i64 80
  br label %183

183:                                              ; preds = %.lr.ph298, %183
  %.3274297 = phi i64 [ %.2273.lcssa, %.lr.ph298 ], [ %216, %183 ]
  %184 = load double, ptr %29, align 8, !tbaa !14
  %185 = fadd double %184, 1.000000e+00
  %186 = load double, ptr %92, align 8, !tbaa !16
  %187 = fneg double %186
  %188 = uitofp nneg i64 %.3274297 to double
  %189 = fmul double %188, %187
  %190 = call double @exp(double noundef %189) #6
  %191 = call double @llvm.fmuladd.f64(double %185, double %190, double -1.000000e+00)
  %192 = call double @rec_HubbleConstant(ptr noundef nonnull %0, double noundef %191)
  %193 = load double, ptr %181, align 8, !tbaa !17
  %194 = fadd double %191, 1.000000e+00
  %195 = call double @cube(double noundef %194) #6
  %196 = fmul double %193, %195
  %197 = fmul double %196, 0x3EB0C6F7A0B5ED8D
  %198 = load double, ptr %0, align 8, !tbaa !18
  %199 = fmul double %198, 0x3F16970041294F34
  %200 = fmul double %194, %199
  %201 = load double, ptr %29, align 8, !tbaa !14
  %202 = load double, ptr %92, align 8, !tbaa !16
  %203 = trunc i64 %.3274297 to i32
  %204 = call double @energy_injection_rate(ptr noundef nonnull %0, double noundef %191)
  %205 = call double @xe_PostSahaH(double noundef %197, double noundef %192, double noundef %200, ptr noundef %1, ptr noundef %2, double noundef %201, double noundef %202, ptr noundef %17, ptr noundef nonnull %6, i32 noundef %203, double noundef %191, ptr noundef nonnull %13, i32 noundef 1, double noundef %204) #6
  %206 = getelementptr inbounds double, ptr %3, i64 %.3274297
  store double %205, ptr %206, align 8, !tbaa !15
  %207 = load double, ptr %0, align 8, !tbaa !18
  %208 = fmul double %194, %207
  %209 = load double, ptr %182, align 8, !tbaa !19
  %210 = load double, ptr %181, align 8, !tbaa !17
  %211 = call double @cube(double noundef %194) #6
  %212 = fmul double %210, %211
  %213 = call double @energy_injection_rate(ptr noundef nonnull %0, double noundef %191)
  %214 = call double @rec_Tmss(double noundef %205, double noundef %208, double noundef %192, double noundef %209, double noundef %212, double noundef %213)
  %215 = getelementptr inbounds double, ptr %4, i64 %.3274297
  store double %214, ptr %215, align 8, !tbaa !15
  %216 = add nuw nsw i64 %.3274297, 1
  %217 = load i64, ptr %14, align 8, !tbaa !5
  %218 = icmp slt i64 %216, %217
  %219 = load double, ptr %13, align 8
  %220 = fcmp olt double %219, 5.000000e-05
  %221 = select i1 %218, i1 %220, i1 false
  br i1 %221, label %183, label %._crit_edge299, !llvm.loop !24

._crit_edge299:                                   ; preds = %183, %.critedge
  %.3274.lcssa = phi i64 [ %.2273.lcssa, %.critedge ], [ %216, %183 ]
  %.3.lcssa = phi double [ %.2.lcssa, %.critedge ], [ %191, %183 ]
  %222 = load double, ptr %29, align 8, !tbaa !14
  %223 = fadd double %222, 1.000000e+00
  %224 = load double, ptr %92, align 8, !tbaa !16
  %225 = fneg double %224
  %226 = add nsw i64 %.3274.lcssa, -3
  %227 = sitofp i64 %226 to double
  %228 = fmul double %227, %225
  %229 = call double @exp(double noundef %228) #6
  %230 = call double @llvm.fmuladd.f64(double %223, double %229, double -1.000000e+00)
  store double %230, ptr %9, align 8, !tbaa !15
  %231 = add nsw i64 %.3274.lcssa, -2
  %232 = getelementptr inbounds double, ptr %3, i64 %231
  %233 = load double, ptr %232, align 8, !tbaa !15
  %234 = getelementptr double, ptr %3, i64 %.3274.lcssa
  %235 = getelementptr i8, ptr %234, i64 -32
  %236 = load double, ptr %235, align 8, !tbaa !15
  %237 = fsub double %233, %236
  %238 = fmul double %237, 5.000000e-01
  %239 = load double, ptr %92, align 8, !tbaa !16
  %240 = fdiv double %238, %239
  store double %240, ptr %10, align 8, !tbaa !15
  %241 = load double, ptr %29, align 8, !tbaa !14
  %242 = fadd double %241, 1.000000e+00
  %243 = fneg double %239
  %244 = sitofp i64 %231 to double
  %245 = fmul double %244, %243
  %246 = call double @exp(double noundef %245) #6
  %247 = call double @llvm.fmuladd.f64(double %242, double %246, double -1.000000e+00)
  store double %247, ptr %7, align 8, !tbaa !15
  %248 = getelementptr i8, ptr %234, i64 -8
  %249 = load double, ptr %248, align 8, !tbaa !15
  %250 = getelementptr inbounds double, ptr %3, i64 %226
  %251 = load double, ptr %250, align 8, !tbaa !15
  %252 = fsub double %249, %251
  %253 = fmul double %252, 5.000000e-01
  %254 = load double, ptr %92, align 8, !tbaa !16
  %255 = fdiv double %253, %254
  store double %255, ptr %8, align 8, !tbaa !15
  %256 = load i64, ptr %14, align 8, !tbaa !5
  %257 = icmp slt i64 %.3274.lcssa, %256
  br i1 %257, label %.lr.ph305, label %.critedge2

.lr.ph305:                                        ; preds = %._crit_edge299
  %258 = getelementptr inbounds i8, ptr %0, i64 80
  %259 = getelementptr inbounds i8, ptr %0, i64 72
  br label %260

260:                                              ; preds = %.lr.ph305, %271
  %.4303 = phi double [ %.3.lcssa, %.lr.ph305 ], [ %283, %271 ]
  %.4275302 = phi i64 [ %.3274.lcssa, %.lr.ph305 ], [ %297, %271 ]
  %261 = add nsw i64 %.4275302, -1
  %262 = getelementptr inbounds double, ptr %4, i64 %261
  %263 = load double, ptr %262, align 8, !tbaa !15
  %264 = load double, ptr %0, align 8, !tbaa !18
  %265 = fdiv double %263, %264
  %266 = fadd double %.4303, 1.000000e+00
  %267 = fdiv double %265, %266
  %268 = fsub double 1.000000e+00, %267
  %269 = fcmp olt double %268, 5.000000e-04
  %270 = fcmp ogt double %.4303, 7.000000e+02
  %or.cond = and i1 %269, %270
  br i1 %or.cond, label %271, label %.critedge2.loopexit

271:                                              ; preds = %260
  %272 = getelementptr inbounds double, ptr %3, i64 %261
  %273 = load double, ptr %272, align 8, !tbaa !15
  %274 = getelementptr inbounds double, ptr %3, i64 %.4275302
  %275 = trunc i64 %261 to i32
  call void @rec_get_xe_next1(ptr noundef nonnull %0, double noundef %.4303, double noundef %273, ptr noundef %274, ptr noundef %1, i32 noundef 2, i32 noundef %275, ptr noundef %2, ptr noundef %17, ptr noundef nonnull %6, ptr noundef nonnull %7, ptr noundef nonnull %8, ptr noundef nonnull %9, ptr noundef nonnull %10)
  %276 = load double, ptr %29, align 8, !tbaa !14
  %277 = fadd double %276, 1.000000e+00
  %278 = load double, ptr %92, align 8, !tbaa !16
  %279 = fneg double %278
  %280 = uitofp nneg i64 %.4275302 to double
  %281 = fmul double %280, %279
  %282 = call double @exp(double noundef %281) #6
  %283 = call double @llvm.fmuladd.f64(double %277, double %282, double -1.000000e+00)
  %284 = load double, ptr %274, align 8, !tbaa !15
  %285 = load double, ptr %0, align 8, !tbaa !18
  %286 = fadd double %283, 1.000000e+00
  %287 = fmul double %285, %286
  %288 = call double @rec_HubbleConstant(ptr noundef nonnull %0, double noundef %283)
  %289 = load double, ptr %258, align 8, !tbaa !19
  %290 = load double, ptr %259, align 8, !tbaa !17
  %291 = call double @cube(double noundef %286) #6
  %292 = fmul double %290, %291
  %293 = fmul double %292, 0x3EB0C6F7A0B5ED8D
  %294 = call double @energy_injection_rate(ptr noundef nonnull %0, double noundef %283)
  %295 = call double @rec_Tmss(double noundef %284, double noundef %287, double noundef %288, double noundef %289, double noundef %293, double noundef %294)
  %296 = getelementptr inbounds double, ptr %4, i64 %.4275302
  store double %295, ptr %296, align 8, !tbaa !15
  %297 = add nuw nsw i64 %.4275302, 1
  %298 = load i64, ptr %14, align 8, !tbaa !5
  %299 = icmp slt i64 %297, %298
  br i1 %299, label %260, label %.critedge2.loopexit, !llvm.loop !25

.critedge2.loopexit:                              ; preds = %260, %271
  %.4275.lcssa.ph = phi i64 [ %297, %271 ], [ %.4275302, %260 ]
  %.4.lcssa.ph = phi double [ %283, %271 ], [ %.4303, %260 ]
  %.pre = add nsw i64 %.4275.lcssa.ph, -3
  %.phi.trans.insert = getelementptr inbounds double, ptr %3, i64 %.pre
  %.pre334 = load double, ptr %.phi.trans.insert, align 8, !tbaa !15
  %.pre335 = add nsw i64 %.4275.lcssa.ph, -2
  br label %.critedge2

.critedge2:                                       ; preds = %.critedge2.loopexit, %._crit_edge299
  %.pre-phi336 = phi i64 [ %.pre335, %.critedge2.loopexit ], [ %231, %._crit_edge299 ]
  %300 = phi double [ %.pre334, %.critedge2.loopexit ], [ %251, %._crit_edge299 ]
  %.pre-phi = phi i64 [ %.pre, %.critedge2.loopexit ], [ %226, %._crit_edge299 ]
  %.4275.lcssa = phi i64 [ %.4275.lcssa.ph, %.critedge2.loopexit ], [ %.3274.lcssa, %._crit_edge299 ]
  %.4.lcssa = phi double [ %.4.lcssa.ph, %.critedge2.loopexit ], [ %.3.lcssa, %._crit_edge299 ]
  %301 = getelementptr inbounds double, ptr %4, i64 %.pre-phi
  %302 = load double, ptr %301, align 8, !tbaa !15
  %303 = load double, ptr %0, align 8, !tbaa !18
  %304 = load double, ptr %9, align 8, !tbaa !15
  %305 = fadd double %304, 1.000000e+00
  %306 = fmul double %303, %305
  %307 = call double @rec_HubbleConstant(ptr noundef nonnull %0, double noundef %304)
  %308 = getelementptr inbounds i8, ptr %0, i64 80
  %309 = load double, ptr %308, align 8, !tbaa !19
  %310 = getelementptr inbounds i8, ptr %0, i64 72
  %311 = load double, ptr %310, align 8, !tbaa !17
  %312 = load double, ptr %9, align 8, !tbaa !15
  %313 = fadd double %312, 1.000000e+00
  %314 = call double @cube(double noundef %313) #6
  %315 = fmul double %311, %314
  %316 = load double, ptr %9, align 8, !tbaa !15
  %317 = call double @energy_injection_rate(ptr noundef nonnull %0, double noundef %316)
  %318 = call double @rec_dTmdlna(double noundef %300, double noundef %302, double noundef %306, double noundef %307, double noundef %309, double noundef %315, double noundef %317)
  store double %318, ptr %12, align 8, !tbaa !15
  %319 = getelementptr inbounds double, ptr %3, i64 %.pre-phi336
  %320 = load double, ptr %319, align 8, !tbaa !15
  %321 = getelementptr inbounds double, ptr %4, i64 %.pre-phi336
  %322 = load double, ptr %321, align 8, !tbaa !15
  %323 = load double, ptr %0, align 8, !tbaa !18
  %324 = load double, ptr %7, align 8, !tbaa !15
  %325 = fadd double %324, 1.000000e+00
  %326 = fmul double %323, %325
  %327 = call double @rec_HubbleConstant(ptr noundef nonnull %0, double noundef %324)
  %328 = load double, ptr %308, align 8, !tbaa !19
  %329 = load double, ptr %310, align 8, !tbaa !17
  %330 = load double, ptr %7, align 8, !tbaa !15
  %331 = fadd double %330, 1.000000e+00
  %332 = call double @cube(double noundef %331) #6
  %333 = fmul double %329, %332
  %334 = load double, ptr %7, align 8, !tbaa !15
  %335 = call double @energy_injection_rate(ptr noundef nonnull %0, double noundef %334)
  %336 = call double @rec_dTmdlna(double noundef %320, double noundef %322, double noundef %326, double noundef %327, double noundef %328, double noundef %333, double noundef %335)
  store double %336, ptr %11, align 8, !tbaa !15
  %337 = load i64, ptr %14, align 8, !tbaa !5
  %338 = icmp slt i64 %.4275.lcssa, %337
  %339 = fcmp ogt double %.4.lcssa, 7.000000e+02
  %340 = select i1 %338, i1 %339, i1 false
  br i1 %340, label %.lr.ph313, label %.preheader279

.preheader279:                                    ; preds = %.lr.ph313, %.critedge2
  %341 = phi i64 [ %337, %.critedge2 ], [ %362, %.lr.ph313 ]
  %.5276.lcssa = phi i64 [ %.4275.lcssa, %.critedge2 ], [ %361, %.lr.ph313 ]
  %.5.lcssa = phi double [ %.4.lcssa, %.critedge2 ], [ %360, %.lr.ph313 ]
  %342 = icmp slt i64 %.5276.lcssa, %341
  %343 = fcmp ogt double %.5.lcssa, 2.000000e+01
  %344 = select i1 %342, i1 %343, i1 false
  br i1 %344, label %.lr.ph318, label %.preheader

.lr.ph313:                                        ; preds = %.critedge2, %.lr.ph313
  %.5312 = phi double [ %360, %.lr.ph313 ], [ %.4.lcssa, %.critedge2 ]
  %.5276311 = phi i64 [ %361, %.lr.ph313 ], [ %.4275.lcssa, %.critedge2 ]
  %345 = add nsw i64 %.5276311, -1
  %346 = getelementptr inbounds double, ptr %3, i64 %345
  %347 = load double, ptr %346, align 8, !tbaa !15
  %348 = getelementptr inbounds double, ptr %4, i64 %345
  %349 = load double, ptr %348, align 8, !tbaa !15
  %350 = getelementptr inbounds double, ptr %3, i64 %.5276311
  %351 = getelementptr inbounds double, ptr %4, i64 %.5276311
  %352 = trunc i64 %345 to i32
  call void @rec_get_xe_next2(ptr noundef nonnull %0, double noundef %.5312, double noundef %347, double noundef %349, ptr noundef %350, ptr noundef %351, ptr noundef %1, i32 noundef 2, i32 noundef %352, ptr noundef %2, ptr noundef %17, ptr noundef nonnull %6, ptr noundef nonnull %7, ptr noundef nonnull %8, ptr noundef nonnull %11, ptr noundef nonnull %9, ptr noundef nonnull %10, ptr noundef nonnull %12)
  %353 = load double, ptr %29, align 8, !tbaa !14
  %354 = fadd double %353, 1.000000e+00
  %355 = load double, ptr %92, align 8, !tbaa !16
  %356 = fneg double %355
  %357 = uitofp nneg i64 %.5276311 to double
  %358 = fmul double %357, %356
  %359 = call double @exp(double noundef %358) #6
  %360 = call double @llvm.fmuladd.f64(double %354, double %359, double -1.000000e+00)
  %361 = add nuw nsw i64 %.5276311, 1
  %362 = load i64, ptr %14, align 8, !tbaa !5
  %363 = icmp slt i64 %361, %362
  %364 = fcmp ogt double %360, 7.000000e+02
  %365 = select i1 %363, i1 %364, i1 false
  br i1 %365, label %.lr.ph313, label %.preheader279, !llvm.loop !26

.preheader:                                       ; preds = %.lr.ph318, %.preheader279
  %366 = phi i64 [ %341, %.preheader279 ], [ %385, %.lr.ph318 ]
  %.6277.lcssa = phi i64 [ %.5276.lcssa, %.preheader279 ], [ %384, %.lr.ph318 ]
  %.6.lcssa = phi double [ %.5.lcssa, %.preheader279 ], [ %383, %.lr.ph318 ]
  %367 = icmp slt i64 %.6277.lcssa, %366
  br i1 %367, label %.lr.ph323, label %._crit_edge324

.lr.ph318:                                        ; preds = %.preheader279, %.lr.ph318
  %.6317 = phi double [ %383, %.lr.ph318 ], [ %.5.lcssa, %.preheader279 ]
  %.6277316 = phi i64 [ %384, %.lr.ph318 ], [ %.5276.lcssa, %.preheader279 ]
  %368 = add nsw i64 %.6277316, -1
  %369 = getelementptr inbounds double, ptr %3, i64 %368
  %370 = load double, ptr %369, align 8, !tbaa !15
  %371 = getelementptr inbounds double, ptr %4, i64 %368
  %372 = load double, ptr %371, align 8, !tbaa !15
  %373 = getelementptr inbounds double, ptr %3, i64 %.6277316
  %374 = getelementptr inbounds double, ptr %4, i64 %.6277316
  %375 = trunc i64 %368 to i32
  call void @rec_get_xe_next2(ptr noundef nonnull %0, double noundef %.6317, double noundef %370, double noundef %372, ptr noundef %373, ptr noundef %374, ptr noundef %1, i32 noundef 3, i32 noundef %375, ptr noundef %2, ptr noundef %17, ptr noundef nonnull %6, ptr noundef nonnull %7, ptr noundef nonnull %8, ptr noundef nonnull %11, ptr noundef nonnull %9, ptr noundef nonnull %10, ptr noundef nonnull %12)
  %376 = load double, ptr %29, align 8, !tbaa !14
  %377 = fadd double %376, 1.000000e+00
  %378 = load double, ptr %92, align 8, !tbaa !16
  %379 = fneg double %378
  %380 = sitofp i64 %.6277316 to double
  %381 = fmul double %380, %379
  %382 = call double @exp(double noundef %381) #6
  %383 = call double @llvm.fmuladd.f64(double %377, double %382, double -1.000000e+00)
  %384 = add nsw i64 %.6277316, 1
  %385 = load i64, ptr %14, align 8, !tbaa !5
  %386 = icmp slt i64 %384, %385
  %387 = fcmp ogt double %383, 2.000000e+01
  %388 = select i1 %386, i1 %387, i1 false
  br i1 %388, label %.lr.ph318, label %.preheader, !llvm.loop !27

.lr.ph323:                                        ; preds = %.preheader, %.lr.ph323
  %.7322 = phi double [ %404, %.lr.ph323 ], [ %.6.lcssa, %.preheader ]
  %.7278321 = phi i64 [ %405, %.lr.ph323 ], [ %.6277.lcssa, %.preheader ]
  %389 = add nsw i64 %.7278321, -1
  %390 = getelementptr inbounds double, ptr %3, i64 %389
  %391 = load double, ptr %390, align 8, !tbaa !15
  %392 = getelementptr inbounds double, ptr %4, i64 %389
  %393 = load double, ptr %392, align 8, !tbaa !15
  %394 = getelementptr inbounds double, ptr %3, i64 %.7278321
  %395 = getelementptr inbounds double, ptr %4, i64 %.7278321
  %396 = trunc i64 %389 to i32
  call void @rec_get_xe_next2(ptr noundef nonnull %0, double noundef %.7322, double noundef %391, double noundef %393, ptr noundef %394, ptr noundef %395, ptr noundef %1, i32 noundef 4, i32 noundef %396, ptr noundef %2, ptr noundef %17, ptr noundef nonnull %6, ptr noundef nonnull %7, ptr noundef nonnull %8, ptr noundef nonnull %11, ptr noundef nonnull %9, ptr noundef nonnull %10, ptr noundef nonnull %12)
  %397 = load double, ptr %29, align 8, !tbaa !14
  %398 = fadd double %397, 1.000000e+00
  %399 = load double, ptr %92, align 8, !tbaa !16
  %400 = fneg double %399
  %401 = sitofp i64 %.7278321 to double
  %402 = fmul double %401, %400
  %403 = call double @exp(double noundef %402) #6
  %404 = call double @llvm.fmuladd.f64(double %398, double %403, double -1.000000e+00)
  %405 = add nsw i64 %.7278321, 1
  %406 = load i64, ptr %14, align 8, !tbaa !5
  %407 = icmp slt i64 %405, %406
  br i1 %407, label %.lr.ph323, label %._crit_edge324, !llvm.loop !28

._crit_edge324:                                   ; preds = %.lr.ph323, %.preheader
  call void @free_2D_array(ptr noundef %17, i32 noundef 311) #6
  %408 = load ptr, ptr %6, align 16, !tbaa !12
  call void @free(ptr noundef %408) #6
  %409 = load ptr, ptr %24, align 8, !tbaa !12
  call void @free(ptr noundef %409) #6
  %410 = load ptr, ptr %28, align 16, !tbaa !12
  call void @free(ptr noundef %410) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %13) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %12) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %11) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %10) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %9) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %8) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %7) #6
  call void @llvm.lifetime.end.p0(i64 24, ptr nonnull %6) #6
  ret void
}

declare ptr @create_2D_array(i32 noundef, i32 noundef) local_unnamed_addr #4

declare ptr @create_1D_array(i32 noundef) local_unnamed_addr #4

declare double @rec_sahaHeII(double noundef, double noundef, double noundef, double noundef, ptr noundef) local_unnamed_addr #4

declare double @xe_PostSahaHe(double noundef, double noundef, double noundef, double noundef, double noundef, ptr noundef) local_unnamed_addr #4

declare double @cube(double noundef) local_unnamed_addr #4

declare void @update_fminus_Saha(ptr noundef, ptr noundef, double noundef, double noundef, double noundef, ptr noundef, double noundef, double noundef, i32 noundef, double noundef, i32 noundef) local_unnamed_addr #4

declare double @rec_saha_xe_H(double noundef, double noundef, double noundef) local_unnamed_addr #4

; Function Attrs: mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.fabs.f64(double) #3

declare double @xe_PostSahaH(double noundef, double noundef, double noundef, ptr noundef, ptr noundef, double noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, i32 noundef, double noundef) local_unnamed_addr #4

declare void @free_2D_array(ptr noundef, i32 noundef) local_unnamed_addr #4

; Function Attrs: mustprogress nounwind willreturn allockind("free") memory(argmem: readwrite, inaccessiblemem: readwrite)
declare void @free(ptr allocptr nocapture noundef) local_unnamed_addr #5

attributes #0 = { mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { mustprogress nofree nounwind willreturn memory(write) "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #4 = { "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #5 = { mustprogress nounwind willreturn allockind("free") memory(argmem: readwrite, inaccessiblemem: readwrite) "alloc-family"="malloc" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #6 = { nounwind }

!llvm.module.flags = !{!0, !1, !2, !3}
!llvm.ident = !{!4}

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !10, i64 112}
!6 = !{!"", !7, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !7, i64 40, !7, i64 48, !7, i64 56, !7, i64 64, !7, i64 72, !7, i64 80, !7, i64 88, !7, i64 96, !7, i64 104, !10, i64 112, !7, i64 120, !11, i64 128, !7, i64 136, !7, i64 144, !7, i64 152, !7, i64 160, !7, i64 168, !7, i64 176, !7, i64 184}
!7 = !{!"double", !8, i64 0}
!8 = !{!"omnipotent char", !9, i64 0}
!9 = !{!"Simple C/C++ TBAA"}
!10 = !{!"long", !8, i64 0}
!11 = !{!"short", !8, i64 0}
!12 = !{!13, !13, i64 0}
!13 = !{!"any pointer", !8, i64 0}
!14 = !{!6, !7, i64 88}
!15 = !{!7, !7, i64 0}
!16 = !{!6, !7, i64 104}
!17 = !{!6, !7, i64 72}
!18 = !{!6, !7, i64 0}
!19 = !{!6, !7, i64 80}
!20 = distinct !{!20, !21}
!21 = !{!"llvm.loop.mustprogress"}
!22 = distinct !{!22, !21}
!23 = distinct !{!23, !21}
!24 = distinct !{!24, !21}
!25 = distinct !{!25, !21}
!26 = distinct !{!26, !21}
!27 = distinct !{!27, !21}
!28 = distinct !{!28, !21}


; Effective O3 passes: BDCEPass,CorrelatedValuePropagationPass,EarlyCSEPass,GVNPass,GlobalOptPass,InferFunctionAttrsPass,InstCombinePass,InstSimplifyPass,JumpThreadingPass,LCSSAPass,LICMPass,LoopRotatePass,LoopSimplifyPass,LoopUnrollPass,ReassociatePass,SROAPass,SimplifyCFGPass,TailCallElimPass