; LLVM IR from CLASS-SZ/class_sz :: class-sz/source (thermodynamics)

; 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.growTable = type { ptr, i64, i64, i32, [2048 x i8] }
%struct.thermo = type { double, i32, i32, i32, double, double, i16, i16, double, double, double, double, i32, ptr, ptr, double, i32, ptr, ptr, double, i32, ptr, ptr, double, i16, double, double, double, double, double, double, double, double, double, double, double, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, ptr, ptr, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i16, i16, i16, [2048 x i8] }
%struct.reionization = type { i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, double, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, i32, i32 }
%struct.recombination = type { i32, i32, i32, i32, i32, i32, i32, i32, ptr, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i16, double, double, double, double, double, double, double }
%struct.precision = type { double, double, double, double, double, double, double, double, double, double, double, double, double, [256 x i8], double, i32, double, i32, i32, double, i32, double, i32, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i32, [256 x i8], [256 x i8], [256 x i8], double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i32, double, double, i32, i32, i32, i32, i32, i32, i32, i32, double, double, double, double, double, double, double, double, i32, double, double, i32, double, double, i32, double, i32, double, double, i32, double, double, double, i32, double, double, double, double, double, i32, double, double, double, double, double, double, double, i32, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i32, double, double, double, double, i32, double, double, i32, i32, i32, double, double, [2048 x i8] }
%struct.background = type { double, double, double, double, double, double, double, i32, double, double, double, double, i16, double, double, double, double, double, double, double, double, double, i16, double, double, ptr, i32, i32, double, i32, ptr, ptr, double, ptr, double, ptr, double, ptr, double, ptr, ptr, ptr, double, double, double, double, i32, ptr, double, double, double, double, double, double, double, double, double, double, double, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, ptr, ptr, ptr, ptr, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, i16, i16, i16, i16, i16, i16, [2048 x i8], i16, [2048 x i8] }

@.str = external hidden unnamed_addr constant [18 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [13 x i8], align 1
@.str.10 = external hidden unnamed_addr constant [27 x i8], align 1
@.str.34 = external hidden unnamed_addr constant [31 x i8], align 1
@.str.37 = external hidden unnamed_addr constant [35 x i8], align 1
@.str.185 = external hidden unnamed_addr constant [12 x i8], align 1
@__func__.thermodynamics_reionization_sample = external hidden unnamed_addr constant [35 x i8], align 1
@.str.186 = external hidden unnamed_addr constant [17 x i8], align 1
@.str.187 = external hidden unnamed_addr constant [73 x i8], align 1
@.str.188 = external hidden unnamed_addr constant [22 x i8], align 1
@.str.189 = external hidden unnamed_addr constant [54 x i8], align 1
@.str.190 = external hidden unnamed_addr constant [34 x i8], align 1
@.str.191 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.192 = external hidden unnamed_addr constant [32 x i8], align 1
@.str.193 = external hidden unnamed_addr constant [72 x i8], align 1
@.str.194 = external hidden unnamed_addr constant [113 x i8], align 1
@.str.195 = external hidden unnamed_addr constant [37 x i8], align 1
@.str.196 = external hidden unnamed_addr constant [64 x i8], align 1
@.str.197 = external hidden unnamed_addr constant [39 x i8], align 1
@.str.198 = external hidden unnamed_addr constant [39 x i8], align 1
@.str.199 = external hidden unnamed_addr constant [59 x i8], align 1
@.str.200 = external hidden unnamed_addr constant [26 x i8], align 1
@.str.201 = external hidden unnamed_addr constant [34 x i8], align 1
@.str.202 = external hidden unnamed_addr constant [54 x i8], align 1
@.str.203 = external hidden unnamed_addr constant [61 x i8], align 1
@.str.204 = external hidden unnamed_addr constant [17 x i8], align 1
@.str.205 = external hidden unnamed_addr constant [93 x i8], align 1
@.str.206 = external hidden unnamed_addr constant [81 x i8], align 1
@.str.207 = external hidden unnamed_addr constant [168 x i8], align 1
@.str.208 = external hidden unnamed_addr constant [217 x i8], align 1

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

; Function Attrs: nounwind
declare double @pow(double noundef, double noundef) #1

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

declare void @class_protect_sprintf(ptr noundef, ptr noundef, ...) #3

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

; Function Attrs: nounwind allocsize(0)
declare noalias ptr @malloc(i64 noundef) #4

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

declare i32 @background_tau_of_z(ptr noundef, double noundef, ptr noundef) #3

declare i32 @background_at_tau(ptr noundef, double noundef, i16 noundef signext, i16 noundef signext, ptr noundef, ptr noundef) #3

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

; Function Attrs: nounwind uwtable
declare dso_local i32 @thermodynamics_energy_injection(ptr noundef, ptr noundef, ptr noundef, double noundef, ptr noundef, ptr noundef) #5

; Function Attrs: nounwind uwtable
declare dso_local i32 @thermodynamics_reionization_function(double noundef, ptr noundef, ptr noundef, ptr noundef) #5

; Function Attrs: nounwind uwtable
define dso_local i32 @thermodynamics_reionization_sample(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, ptr noundef %4, ptr noundef %5) #5 {
  %7 = alloca i32, align 4
  %8 = alloca ptr, align 8
  %9 = alloca ptr, align 8
  %10 = alloca ptr, align 8
  %11 = alloca ptr, align 8
  %12 = alloca ptr, align 8
  %13 = alloca ptr, align 8
  %14 = alloca %struct.growTable, align 8
  %15 = alloca ptr, align 8
  %16 = alloca ptr, align 8
  %17 = alloca ptr, align 8
  %18 = alloca i32, align 4
  %19 = alloca i32, align 4
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca double, align 8
  %23 = alloca double, align 8
  %24 = alloca double, align 8
  %25 = alloca double, align 8
  %26 = alloca double, align 8
  %27 = alloca double, align 8
  %28 = alloca double, align 8
  %29 = alloca double, align 8
  %30 = alloca double, align 8
  %31 = alloca double, align 8
  %32 = alloca double, align 8
  %33 = alloca double, align 8
  %34 = alloca double, align 8
  %35 = alloca double, align 8
  %36 = alloca double, align 8
  %37 = alloca double, align 8
  %38 = alloca i32, align 4
  %39 = alloca double, align 8
  %40 = alloca i32, align 4
  %41 = alloca [2048 x i8], align 16
  %42 = alloca i32, align 4
  %43 = alloca [2048 x i8], align 16
  %44 = alloca [2048 x i8], align 16
  %45 = alloca [2048 x i8], align 16
  %46 = alloca [2048 x i8], align 16
  %47 = alloca [2048 x i8], align 16
  %48 = alloca [2048 x i8], align 16
  %49 = alloca [2048 x i8], align 16
  %50 = alloca [2048 x i8], align 16
  %51 = alloca [2048 x i8], align 16
  %52 = alloca [2048 x i8], align 16
  %53 = alloca [2048 x i8], align 16
  %54 = alloca [2048 x i8], align 16
  %55 = alloca [2048 x i8], align 16
  %56 = alloca [2048 x i8], align 16
  %57 = alloca [2048 x i8], align 16
  %58 = alloca [2048 x i8], align 16
  %59 = alloca [2048 x i8], align 16
  %60 = alloca [2048 x i8], align 16
  %61 = alloca [2048 x i8], align 16
  %62 = alloca [2048 x i8], align 16
  %63 = alloca [2048 x i8], align 16
  %64 = alloca i32, align 4
  %65 = alloca [2048 x i8], align 16
  %66 = alloca [2048 x i8], align 16
  %67 = alloca [2048 x i8], align 16
  %68 = alloca [2048 x i8], align 16
  %69 = alloca [2048 x i8], align 16
  %70 = alloca [2048 x i8], align 16
  %71 = alloca [2048 x i8], align 16
  %72 = alloca [2048 x i8], align 16
  %73 = alloca [2048 x i8], align 16
  %74 = alloca [2048 x i8], align 16
  store ptr %0, ptr %8, align 8, !tbaa !5
  store ptr %1, ptr %9, align 8, !tbaa !5
  store ptr %2, ptr %10, align 8, !tbaa !5
  store ptr %3, ptr %11, align 8, !tbaa !5
  store ptr %4, ptr %12, align 8, !tbaa !5
  store ptr %5, ptr %13, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 2080, ptr %14) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %17) #7
  call void @llvm.lifetime.start.p0(i64 4, ptr %18) #7
  call void @llvm.lifetime.start.p0(i64 4, ptr %19) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %20) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %22) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %23) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %24) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %25) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %26) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %27) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %28) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %29) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %30) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %31) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %32) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %33) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %34) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %35) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %36) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %37) #7
  call void @llvm.lifetime.start.p0(i64 4, ptr %38) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr %39) #7
  %75 = load ptr, ptr %10, align 8, !tbaa !5
  %76 = getelementptr inbounds %struct.thermo, ptr %75, i32 0, i32 0
  %77 = load double, ptr %76, align 8, !tbaa !9
  store double %77, ptr %29, align 8, !tbaa !14
  %78 = load ptr, ptr %12, align 8, !tbaa !5
  %79 = getelementptr inbounds %struct.reionization, ptr %78, i32 0, i32 8
  %80 = load i32, ptr %79, align 8, !tbaa !15
  %81 = sext i32 %80 to i64
  %82 = mul i64 %81, 8
  %83 = call noalias ptr @malloc(i64 noundef %82) #8
  store ptr %83, ptr %17, align 8, !tbaa !5
  %84 = load ptr, ptr %17, align 8, !tbaa !5
  %85 = icmp eq ptr %84, null
  br i1 %85, label %86, label %99

86:                                               ; preds = %6
  call void @llvm.lifetime.start.p0(i64 4, ptr %40) #7
  %87 = load ptr, ptr %12, align 8, !tbaa !5
  %88 = getelementptr inbounds %struct.reionization, ptr %87, i32 0, i32 8
  %89 = load i32, ptr %88, align 8, !tbaa !15
  %90 = sext i32 %89 to i64
  %91 = mul i64 %90, 8
  %92 = trunc i64 %91 to i32
  store i32 %92, ptr %40, align 4, !tbaa !17
  call void @llvm.lifetime.start.p0(i64 2048, ptr %41) #7
  %93 = getelementptr inbounds [2048 x i8], ptr %41, i64 0, i64 0
  %94 = load i32, ptr %40, align 4, !tbaa !17
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %93, ptr noundef @.str.37, ptr noundef @.str.185, i32 noundef %94)
  %95 = load ptr, ptr %10, align 8, !tbaa !5
  %96 = getelementptr inbounds %struct.thermo, ptr %95, i32 0, i32 93
  %97 = getelementptr inbounds [2048 x i8], ptr %96, i64 0, i64 0
  %98 = getelementptr inbounds [2048 x i8], ptr %41, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %97, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2784, ptr noundef %98)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %41) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %40) #7
  br label %1323

99:                                               ; preds = %6
  %100 = call i32 @gt_init(ptr noundef %14)
  %101 = icmp eq i32 %100, 1
  br i1 %101, label %102, label %110

102:                                              ; preds = %99
  call void @llvm.lifetime.start.p0(i64 2048, ptr %43) #7
  %103 = getelementptr inbounds [2048 x i8], ptr %43, i64 0, i64 0
  %104 = getelementptr inbounds %struct.growTable, ptr %14, i32 0, i32 4
  %105 = getelementptr inbounds [2048 x i8], ptr %104, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %103, ptr noundef @.str, ptr noundef @.str.186, ptr noundef %105)
  %106 = load ptr, ptr %10, align 8, !tbaa !5
  %107 = getelementptr inbounds %struct.thermo, ptr %106, i32 0, i32 93
  %108 = getelementptr inbounds [2048 x i8], ptr %107, i64 0, i64 0
  %109 = getelementptr inbounds [2048 x i8], ptr %43, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %108, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2789, ptr noundef %109)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %43) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

110:                                              ; preds = %99
  store i32 0, ptr %18, align 4, !tbaa !17
  br label %111

111:                                              ; preds = %156, %110
  %112 = load ptr, ptr %11, align 8, !tbaa !5
  %113 = getelementptr inbounds %struct.recombination, ptr %112, i32 0, i32 8
  %114 = load ptr, ptr %113, align 8, !tbaa !18
  %115 = load i32, ptr %18, align 4, !tbaa !17
  %116 = load ptr, ptr %11, align 8, !tbaa !5
  %117 = getelementptr inbounds %struct.recombination, ptr %116, i32 0, i32 6
  %118 = load i32, ptr %117, align 8, !tbaa !20
  %119 = mul nsw i32 %115, %118
  %120 = load ptr, ptr %11, align 8, !tbaa !5
  %121 = getelementptr inbounds %struct.recombination, ptr %120, i32 0, i32 0
  %122 = load i32, ptr %121, align 8, !tbaa !21
  %123 = add nsw i32 %119, %122
  %124 = sext i32 %123 to i64
  %125 = getelementptr inbounds double, ptr %114, i64 %124
  %126 = load double, ptr %125, align 8, !tbaa !14
  %127 = load ptr, ptr %12, align 8, !tbaa !5
  %128 = getelementptr inbounds %struct.reionization, ptr %127, i32 0, i32 25
  %129 = load ptr, ptr %128, align 8, !tbaa !22
  %130 = load ptr, ptr %12, align 8, !tbaa !5
  %131 = getelementptr inbounds %struct.reionization, ptr %130, i32 0, i32 24
  %132 = load i32, ptr %131, align 8, !tbaa !23
  %133 = sext i32 %132 to i64
  %134 = getelementptr inbounds double, ptr %129, i64 %133
  %135 = load double, ptr %134, align 8, !tbaa !14
  %136 = fcmp olt double %126, %135
  br i1 %136, label %137, label %157

137:                                              ; preds = %111
  %138 = load i32, ptr %18, align 4, !tbaa !17
  %139 = add nsw i32 %138, 1
  store i32 %139, ptr %18, align 4, !tbaa !17
  %140 = load i32, ptr %18, align 4, !tbaa !17
  %141 = load ptr, ptr %8, align 8, !tbaa !5
  %142 = getelementptr inbounds %struct.precision, ptr %141, i32 0, i32 15
  %143 = load i32, ptr %142, align 8, !tbaa !24
  %144 = icmp eq i32 %140, %143
  br i1 %144, label %145, label %156

145:                                              ; preds = %137
  call void @llvm.lifetime.start.p0(i64 2048, ptr %44) #7
  %146 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  %147 = load ptr, ptr %8, align 8, !tbaa !5
  %148 = getelementptr inbounds %struct.precision, ptr %147, i32 0, i32 44
  %149 = load double, ptr %148, align 8, !tbaa !26
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %146, ptr noundef @.str.187, double noundef %149)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %45) #7
  %150 = getelementptr inbounds [2048 x i8], ptr %45, i64 0, i64 0
  %151 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %150, ptr noundef @.str.10, ptr noundef @.str.188, ptr noundef %151)
  %152 = load ptr, ptr %10, align 8, !tbaa !5
  %153 = getelementptr inbounds %struct.thermo, ptr %152, i32 0, i32 93
  %154 = getelementptr inbounds [2048 x i8], ptr %153, i64 0, i64 0
  %155 = getelementptr inbounds [2048 x i8], ptr %45, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %154, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2799, ptr noundef %155)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %45) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr %44) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

156:                                              ; preds = %137
  br label %111, !llvm.loop !27

157:                                              ; preds = %111
  %158 = load ptr, ptr %11, align 8, !tbaa !5
  %159 = getelementptr inbounds %struct.recombination, ptr %158, i32 0, i32 8
  %160 = load ptr, ptr %159, align 8, !tbaa !18
  %161 = load i32, ptr %18, align 4, !tbaa !17
  %162 = load ptr, ptr %11, align 8, !tbaa !5
  %163 = getelementptr inbounds %struct.recombination, ptr %162, i32 0, i32 6
  %164 = load i32, ptr %163, align 8, !tbaa !20
  %165 = mul nsw i32 %161, %164
  %166 = load ptr, ptr %11, align 8, !tbaa !5
  %167 = getelementptr inbounds %struct.recombination, ptr %166, i32 0, i32 0
  %168 = load i32, ptr %167, align 8, !tbaa !21
  %169 = add nsw i32 %165, %168
  %170 = sext i32 %169 to i64
  %171 = getelementptr inbounds double, ptr %160, i64 %170
  %172 = load double, ptr %171, align 8, !tbaa !14
  store double %172, ptr %22, align 8, !tbaa !14
  %173 = load double, ptr %22, align 8, !tbaa !14
  %174 = load ptr, ptr %17, align 8, !tbaa !5
  %175 = load ptr, ptr %12, align 8, !tbaa !5
  %176 = getelementptr inbounds %struct.reionization, ptr %175, i32 0, i32 0
  %177 = load i32, ptr %176, align 8, !tbaa !29
  %178 = sext i32 %177 to i64
  %179 = getelementptr inbounds double, ptr %174, i64 %178
  store double %173, ptr %179, align 8, !tbaa !14
  %180 = load i32, ptr %18, align 4, !tbaa !17
  %181 = load ptr, ptr %12, align 8, !tbaa !5
  %182 = getelementptr inbounds %struct.reionization, ptr %181, i32 0, i32 27
  store i32 %180, ptr %182, align 4, !tbaa !30
  %183 = load double, ptr %22, align 8, !tbaa !14
  %184 = load ptr, ptr %10, align 8, !tbaa !5
  %185 = load ptr, ptr %12, align 8, !tbaa !5
  %186 = call i32 @thermodynamics_reionization_function(double noundef %183, ptr noundef %184, ptr noundef %185, ptr noundef %24)
  %187 = icmp eq i32 %186, 1
  br i1 %187, label %188, label %197

188:                                              ; preds = %157
  call void @llvm.lifetime.start.p0(i64 2048, ptr %46) #7
  %189 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  %190 = load ptr, ptr %10, align 8, !tbaa !5
  %191 = getelementptr inbounds %struct.thermo, ptr %190, i32 0, i32 93
  %192 = getelementptr inbounds [2048 x i8], ptr %191, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %189, ptr noundef @.str, ptr noundef @.str.189, ptr noundef %192)
  %193 = load ptr, ptr %10, align 8, !tbaa !5
  %194 = getelementptr inbounds %struct.thermo, ptr %193, i32 0, i32 93
  %195 = getelementptr inbounds [2048 x i8], ptr %194, i64 0, i64 0
  %196 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %195, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2810, ptr noundef %196)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %46) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

197:                                              ; preds = %157
  %198 = load double, ptr %24, align 8, !tbaa !14
  %199 = load ptr, ptr %17, align 8, !tbaa !5
  %200 = load ptr, ptr %12, align 8, !tbaa !5
  %201 = getelementptr inbounds %struct.reionization, ptr %200, i32 0, i32 1
  %202 = load i32, ptr %201, align 4, !tbaa !31
  %203 = sext i32 %202 to i64
  %204 = getelementptr inbounds double, ptr %199, i64 %203
  store double %198, ptr %204, align 8, !tbaa !14
  %205 = load ptr, ptr %9, align 8, !tbaa !5
  %206 = load double, ptr %22, align 8, !tbaa !14
  %207 = call i32 @background_tau_of_z(ptr noundef %205, double noundef %206, ptr noundef %36)
  %208 = icmp eq i32 %207, 1
  br i1 %208, label %209, label %218

209:                                              ; preds = %197
  call void @llvm.lifetime.start.p0(i64 2048, ptr %47) #7
  %210 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  %211 = load ptr, ptr %9, align 8, !tbaa !5
  %212 = getelementptr inbounds %struct.background, ptr %211, i32 0, i32 149
  %213 = getelementptr inbounds [2048 x i8], ptr %212, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %210, ptr noundef @.str, ptr noundef @.str.190, ptr noundef %213)
  %214 = load ptr, ptr %10, align 8, !tbaa !5
  %215 = getelementptr inbounds %struct.thermo, ptr %214, i32 0, i32 93
  %216 = getelementptr inbounds [2048 x i8], ptr %215, i64 0, i64 0
  %217 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %216, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2820, ptr noundef %217)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %47) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

218:                                              ; preds = %197
  %219 = load ptr, ptr %9, align 8, !tbaa !5
  %220 = load double, ptr %36, align 8, !tbaa !14
  %221 = load ptr, ptr %9, align 8, !tbaa !5
  %222 = getelementptr inbounds %struct.background, ptr %221, i32 0, i32 141
  %223 = load i16, ptr %222, align 8, !tbaa !32
  %224 = load ptr, ptr %9, align 8, !tbaa !5
  %225 = getelementptr inbounds %struct.background, ptr %224, i32 0, i32 144
  %226 = load i16, ptr %225, align 2, !tbaa !34
  %227 = load ptr, ptr %13, align 8, !tbaa !5
  %228 = call i32 @background_at_tau(ptr noundef %219, double noundef %220, i16 noundef signext %223, i16 noundef signext %226, ptr noundef %38, ptr noundef %227)
  %229 = icmp eq i32 %228, 1
  br i1 %229, label %230, label %239

230:                                              ; preds = %218
  call void @llvm.lifetime.start.p0(i64 2048, ptr %48) #7
  %231 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  %232 = load ptr, ptr %9, align 8, !tbaa !5
  %233 = getelementptr inbounds %struct.background, ptr %232, i32 0, i32 149
  %234 = getelementptr inbounds [2048 x i8], ptr %233, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %231, ptr noundef @.str, ptr noundef @.str.191, ptr noundef %234)
  %235 = load ptr, ptr %10, align 8, !tbaa !5
  %236 = getelementptr inbounds %struct.thermo, ptr %235, i32 0, i32 93
  %237 = getelementptr inbounds [2048 x i8], ptr %236, i64 0, i64 0
  %238 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %237, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2829, ptr noundef %238)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %48) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

239:                                              ; preds = %218
  %240 = load double, ptr %22, align 8, !tbaa !14
  %241 = fadd double 1.000000e+00, %240
  %242 = load double, ptr %22, align 8, !tbaa !14
  %243 = fadd double 1.000000e+00, %242
  %244 = fmul double %241, %243
  %245 = load ptr, ptr %10, align 8, !tbaa !5
  %246 = getelementptr inbounds %struct.thermo, ptr %245, i32 0, i32 89
  %247 = load double, ptr %246, align 8, !tbaa !35
  %248 = fmul double %244, %247
  %249 = load double, ptr %24, align 8, !tbaa !14
  %250 = fmul double %248, %249
  %251 = fmul double %250, 0x3A15151E37FBF236
  %252 = fmul double %251, 0x449A22FF260E3826
  %253 = load ptr, ptr %17, align 8, !tbaa !5
  %254 = load ptr, ptr %12, align 8, !tbaa !5
  %255 = getelementptr inbounds %struct.reionization, ptr %254, i32 0, i32 5
  %256 = load i32, ptr %255, align 4, !tbaa !36
  %257 = sext i32 %256 to i64
  %258 = getelementptr inbounds double, ptr %253, i64 %257
  store double %252, ptr %258, align 8, !tbaa !14
  %259 = load ptr, ptr %13, align 8, !tbaa !5
  %260 = load ptr, ptr %9, align 8, !tbaa !5
  %261 = getelementptr inbounds %struct.background, ptr %260, i32 0, i32 60
  %262 = load i32, ptr %261, align 4, !tbaa !37
  %263 = sext i32 %262 to i64
  %264 = getelementptr inbounds double, ptr %259, i64 %263
  %265 = load double, ptr %264, align 8, !tbaa !14
  %266 = fcmp oeq double %265, 0.000000e+00
  br i1 %266, label %267, label %275

267:                                              ; preds = %239
  call void @llvm.lifetime.start.p0(i64 2048, ptr %49) #7
  %268 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %268, ptr noundef @.str.34)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %50) #7
  %269 = getelementptr inbounds [2048 x i8], ptr %50, i64 0, i64 0
  %270 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %269, ptr noundef @.str.10, ptr noundef @.str.192, ptr noundef %270)
  %271 = load ptr, ptr %10, align 8, !tbaa !5
  %272 = getelementptr inbounds %struct.thermo, ptr %271, i32 0, i32 93
  %273 = getelementptr inbounds [2048 x i8], ptr %272, i64 0, i64 0
  %274 = getelementptr inbounds [2048 x i8], ptr %50, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %273, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2835, ptr noundef %274)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %50) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr %49) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

275:                                              ; preds = %239
  %276 = load ptr, ptr %17, align 8, !tbaa !5
  %277 = load ptr, ptr %12, align 8, !tbaa !5
  %278 = getelementptr inbounds %struct.reionization, ptr %277, i32 0, i32 5
  %279 = load i32, ptr %278, align 4, !tbaa !36
  %280 = sext i32 %279 to i64
  %281 = getelementptr inbounds double, ptr %276, i64 %280
  %282 = load double, ptr %281, align 8, !tbaa !14
  %283 = load ptr, ptr %13, align 8, !tbaa !5
  %284 = load ptr, ptr %9, align 8, !tbaa !5
  %285 = getelementptr inbounds %struct.background, ptr %284, i32 0, i32 60
  %286 = load i32, ptr %285, align 4, !tbaa !37
  %287 = sext i32 %286 to i64
  %288 = getelementptr inbounds double, ptr %283, i64 %287
  %289 = load double, ptr %288, align 8, !tbaa !14
  %290 = fdiv double %282, %289
  %291 = load ptr, ptr %17, align 8, !tbaa !5
  %292 = load ptr, ptr %12, align 8, !tbaa !5
  %293 = getelementptr inbounds %struct.reionization, ptr %292, i32 0, i32 6
  %294 = load i32, ptr %293, align 8, !tbaa !38
  %295 = sext i32 %294 to i64
  %296 = getelementptr inbounds double, ptr %291, i64 %295
  store double %290, ptr %296, align 8, !tbaa !14
  %297 = load ptr, ptr %17, align 8, !tbaa !5
  %298 = load ptr, ptr %12, align 8, !tbaa !5
  %299 = getelementptr inbounds %struct.reionization, ptr %298, i32 0, i32 6
  %300 = load i32, ptr %299, align 8, !tbaa !38
  %301 = sext i32 %300 to i64
  %302 = getelementptr inbounds double, ptr %297, i64 %301
  %303 = load double, ptr %302, align 8, !tbaa !14
  store double %303, ptr %26, align 8, !tbaa !14
  %304 = load ptr, ptr %17, align 8, !tbaa !5
  %305 = load ptr, ptr %12, align 8, !tbaa !5
  %306 = getelementptr inbounds %struct.reionization, ptr %305, i32 0, i32 5
  %307 = load i32, ptr %306, align 4, !tbaa !36
  %308 = sext i32 %307 to i64
  %309 = getelementptr inbounds double, ptr %304, i64 %308
  %310 = load double, ptr %309, align 8, !tbaa !14
  store double %310, ptr %33, align 8, !tbaa !14
  %311 = load ptr, ptr %11, align 8, !tbaa !5
  %312 = getelementptr inbounds %struct.recombination, ptr %311, i32 0, i32 8
  %313 = load ptr, ptr %312, align 8, !tbaa !18
  %314 = load i32, ptr %18, align 4, !tbaa !17
  %315 = load ptr, ptr %11, align 8, !tbaa !5
  %316 = getelementptr inbounds %struct.recombination, ptr %315, i32 0, i32 6
  %317 = load i32, ptr %316, align 8, !tbaa !20
  %318 = mul nsw i32 %314, %317
  %319 = load ptr, ptr %11, align 8, !tbaa !5
  %320 = getelementptr inbounds %struct.recombination, ptr %319, i32 0, i32 2
  %321 = load i32, ptr %320, align 8, !tbaa !39
  %322 = add nsw i32 %318, %321
  %323 = sext i32 %322 to i64
  %324 = getelementptr inbounds double, ptr %313, i64 %323
  %325 = load double, ptr %324, align 8, !tbaa !14
  store double %325, ptr %28, align 8, !tbaa !14
  %326 = load double, ptr %28, align 8, !tbaa !14
  %327 = load ptr, ptr %17, align 8, !tbaa !5
  %328 = load ptr, ptr %12, align 8, !tbaa !5
  %329 = getelementptr inbounds %struct.reionization, ptr %328, i32 0, i32 2
  %330 = load i32, ptr %329, align 8, !tbaa !40
  %331 = sext i32 %330 to i64
  %332 = getelementptr inbounds double, ptr %327, i64 %331
  store double %326, ptr %332, align 8, !tbaa !14
  %333 = load double, ptr %29, align 8, !tbaa !14
  %334 = call double @llvm.fmuladd.f64(double 0xBFE7F14DA42541E8, double %333, double 1.000000e+00)
  %335 = load double, ptr %24, align 8, !tbaa !14
  %336 = load double, ptr %29, align 8, !tbaa !14
  %337 = fsub double 1.000000e+00, %336
  %338 = call double @llvm.fmuladd.f64(double %335, double %337, double %334)
  %339 = fmul double 0x3D39D63187BB40C4, %338
  %340 = load double, ptr %28, align 8, !tbaa !14
  %341 = fmul double %339, %340
  %342 = load ptr, ptr %17, align 8, !tbaa !5
  %343 = load ptr, ptr %12, align 8, !tbaa !5
  %344 = getelementptr inbounds %struct.reionization, ptr %343, i32 0, i32 3
  %345 = load i32, ptr %344, align 4, !tbaa !41
  %346 = sext i32 %345 to i64
  %347 = getelementptr inbounds double, ptr %342, i64 %346
  store double %341, ptr %347, align 8, !tbaa !14
  %348 = load ptr, ptr %17, align 8, !tbaa !5
  %349 = load ptr, ptr %12, align 8, !tbaa !5
  %350 = getelementptr inbounds %struct.reionization, ptr %349, i32 0, i32 3
  %351 = load i32, ptr %350, align 4, !tbaa !41
  %352 = sext i32 %351 to i64
  %353 = getelementptr inbounds double, ptr %348, i64 %352
  %354 = load double, ptr %353, align 8, !tbaa !14
  %355 = fmul double 0x3FFAAAAAAAAAAAAB, %354
  %356 = load ptr, ptr %17, align 8, !tbaa !5
  %357 = load ptr, ptr %12, align 8, !tbaa !5
  %358 = getelementptr inbounds %struct.reionization, ptr %357, i32 0, i32 4
  %359 = load i32, ptr %358, align 8, !tbaa !42
  %360 = sext i32 %359 to i64
  %361 = getelementptr inbounds double, ptr %356, i64 %360
  store double %355, ptr %361, align 8, !tbaa !14
  %362 = load ptr, ptr %17, align 8, !tbaa !5
  %363 = load ptr, ptr %12, align 8, !tbaa !5
  %364 = getelementptr inbounds %struct.reionization, ptr %363, i32 0, i32 8
  %365 = load i32, ptr %364, align 8, !tbaa !15
  %366 = sext i32 %365 to i64
  %367 = mul i64 8, %366
  %368 = call i32 @gt_add(ptr noundef %14, i64 noundef -1, ptr noundef %362, i64 noundef %367)
  %369 = icmp eq i32 %368, 1
  br i1 %369, label %370, label %378

370:                                              ; preds = %275
  call void @llvm.lifetime.start.p0(i64 2048, ptr %51) #7
  %371 = getelementptr inbounds [2048 x i8], ptr %51, i64 0, i64 0
  %372 = getelementptr inbounds %struct.growTable, ptr %14, i32 0, i32 4
  %373 = getelementptr inbounds [2048 x i8], ptr %372, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %371, ptr noundef @.str, ptr noundef @.str.193, ptr noundef %373)
  %374 = load ptr, ptr %10, align 8, !tbaa !5
  %375 = getelementptr inbounds %struct.thermo, ptr %374, i32 0, i32 93
  %376 = getelementptr inbounds [2048 x i8], ptr %375, i64 0, i64 0
  %377 = getelementptr inbounds [2048 x i8], ptr %51, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %376, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2858, ptr noundef %377)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %51) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

378:                                              ; preds = %275
  store i32 1, ptr %19, align 4, !tbaa !17
  %379 = load ptr, ptr %11, align 8, !tbaa !5
  %380 = getelementptr inbounds %struct.recombination, ptr %379, i32 0, i32 8
  %381 = load ptr, ptr %380, align 8, !tbaa !18
  %382 = load i32, ptr %18, align 4, !tbaa !17
  %383 = load ptr, ptr %11, align 8, !tbaa !5
  %384 = getelementptr inbounds %struct.recombination, ptr %383, i32 0, i32 6
  %385 = load i32, ptr %384, align 8, !tbaa !20
  %386 = mul nsw i32 %382, %385
  %387 = load ptr, ptr %11, align 8, !tbaa !5
  %388 = getelementptr inbounds %struct.recombination, ptr %387, i32 0, i32 0
  %389 = load i32, ptr %388, align 8, !tbaa !21
  %390 = add nsw i32 %386, %389
  %391 = sext i32 %390 to i64
  %392 = getelementptr inbounds double, ptr %381, i64 %391
  %393 = load double, ptr %392, align 8, !tbaa !14
  %394 = load ptr, ptr %11, align 8, !tbaa !5
  %395 = getelementptr inbounds %struct.recombination, ptr %394, i32 0, i32 8
  %396 = load ptr, ptr %395, align 8, !tbaa !18
  %397 = load i32, ptr %18, align 4, !tbaa !17
  %398 = sub nsw i32 %397, 1
  %399 = load ptr, ptr %11, align 8, !tbaa !5
  %400 = getelementptr inbounds %struct.recombination, ptr %399, i32 0, i32 6
  %401 = load i32, ptr %400, align 8, !tbaa !20
  %402 = mul nsw i32 %398, %401
  %403 = load ptr, ptr %11, align 8, !tbaa !5
  %404 = getelementptr inbounds %struct.recombination, ptr %403, i32 0, i32 0
  %405 = load i32, ptr %404, align 8, !tbaa !21
  %406 = add nsw i32 %402, %405
  %407 = sext i32 %406 to i64
  %408 = getelementptr inbounds double, ptr %396, i64 %407
  %409 = load double, ptr %408, align 8, !tbaa !14
  %410 = fsub double %393, %409
  store double %410, ptr %21, align 8, !tbaa !14
  %411 = load double, ptr %21, align 8, !tbaa !14
  store double %411, ptr %20, align 8, !tbaa !14
  br label %412

412:                                              ; preds = %683, %378
  %413 = load double, ptr %22, align 8, !tbaa !14
  %414 = fcmp ogt double %413, 0.000000e+00
  br i1 %414, label %415, label %684

415:                                              ; preds = %412
  %416 = load double, ptr %20, align 8, !tbaa !14
  %417 = load ptr, ptr %8, align 8, !tbaa !5
  %418 = getelementptr inbounds %struct.precision, ptr %417, i32 0, i32 175
  %419 = load double, ptr %418, align 8, !tbaa !43
  %420 = fcmp olt double %416, %419
  br i1 %420, label %421, label %429

421:                                              ; preds = %415
  call void @llvm.lifetime.start.p0(i64 2048, ptr %52) #7
  %422 = getelementptr inbounds [2048 x i8], ptr %52, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %422, ptr noundef @.str.194)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %53) #7
  %423 = getelementptr inbounds [2048 x i8], ptr %53, i64 0, i64 0
  %424 = getelementptr inbounds [2048 x i8], ptr %52, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %423, ptr noundef @.str.10, ptr noundef @.str.195, ptr noundef %424)
  %425 = load ptr, ptr %10, align 8, !tbaa !5
  %426 = getelementptr inbounds %struct.thermo, ptr %425, i32 0, i32 93
  %427 = getelementptr inbounds [2048 x i8], ptr %426, i64 0, i64 0
  %428 = getelementptr inbounds [2048 x i8], ptr %53, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %427, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2875, ptr noundef %428)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %53) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr %52) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

429:                                              ; preds = %415
  %430 = load double, ptr %22, align 8, !tbaa !14
  %431 = load double, ptr %20, align 8, !tbaa !14
  %432 = fsub double %430, %431
  store double %432, ptr %23, align 8, !tbaa !14
  %433 = load double, ptr %23, align 8, !tbaa !14
  %434 = fcmp olt double %433, 0.000000e+00
  br i1 %434, label %435, label %436

435:                                              ; preds = %429
  store double 0.000000e+00, ptr %23, align 8, !tbaa !14
  br label %436

436:                                              ; preds = %435, %429
  %437 = load double, ptr %23, align 8, !tbaa !14
  %438 = load ptr, ptr %10, align 8, !tbaa !5
  %439 = load ptr, ptr %12, align 8, !tbaa !5
  %440 = call i32 @thermodynamics_reionization_function(double noundef %437, ptr noundef %438, ptr noundef %439, ptr noundef %25)
  %441 = icmp eq i32 %440, 1
  br i1 %441, label %442, label %451

442:                                              ; preds = %436
  call void @llvm.lifetime.start.p0(i64 2048, ptr %54) #7
  %443 = getelementptr inbounds [2048 x i8], ptr %54, i64 0, i64 0
  %444 = load ptr, ptr %10, align 8, !tbaa !5
  %445 = getelementptr inbounds %struct.thermo, ptr %444, i32 0, i32 93
  %446 = getelementptr inbounds [2048 x i8], ptr %445, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %443, ptr noundef @.str, ptr noundef @.str.196, ptr noundef %446)
  %447 = load ptr, ptr %10, align 8, !tbaa !5
  %448 = getelementptr inbounds %struct.thermo, ptr %447, i32 0, i32 93
  %449 = getelementptr inbounds [2048 x i8], ptr %448, i64 0, i64 0
  %450 = getelementptr inbounds [2048 x i8], ptr %54, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %449, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2883, ptr noundef %450)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %54) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

451:                                              ; preds = %436
  %452 = load ptr, ptr %9, align 8, !tbaa !5
  %453 = load double, ptr %23, align 8, !tbaa !14
  %454 = call i32 @background_tau_of_z(ptr noundef %452, double noundef %453, ptr noundef %36)
  %455 = icmp eq i32 %454, 1
  br i1 %455, label %456, label %465

456:                                              ; preds = %451
  call void @llvm.lifetime.start.p0(i64 2048, ptr %55) #7
  %457 = getelementptr inbounds [2048 x i8], ptr %55, i64 0, i64 0
  %458 = load ptr, ptr %9, align 8, !tbaa !5
  %459 = getelementptr inbounds %struct.background, ptr %458, i32 0, i32 149
  %460 = getelementptr inbounds [2048 x i8], ptr %459, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %457, ptr noundef @.str, ptr noundef @.str.197, ptr noundef %460)
  %461 = load ptr, ptr %10, align 8, !tbaa !5
  %462 = getelementptr inbounds %struct.thermo, ptr %461, i32 0, i32 93
  %463 = getelementptr inbounds [2048 x i8], ptr %462, i64 0, i64 0
  %464 = getelementptr inbounds [2048 x i8], ptr %55, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %463, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2889, ptr noundef %464)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %55) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

465:                                              ; preds = %451
  %466 = load ptr, ptr %9, align 8, !tbaa !5
  %467 = load double, ptr %36, align 8, !tbaa !14
  %468 = load ptr, ptr %9, align 8, !tbaa !5
  %469 = getelementptr inbounds %struct.background, ptr %468, i32 0, i32 141
  %470 = load i16, ptr %469, align 8, !tbaa !32
  %471 = load ptr, ptr %9, align 8, !tbaa !5
  %472 = getelementptr inbounds %struct.background, ptr %471, i32 0, i32 144
  %473 = load i16, ptr %472, align 2, !tbaa !34
  %474 = load ptr, ptr %13, align 8, !tbaa !5
  %475 = call i32 @background_at_tau(ptr noundef %466, double noundef %467, i16 noundef signext %470, i16 noundef signext %473, ptr noundef %38, ptr noundef %474)
  %476 = icmp eq i32 %475, 1
  br i1 %476, label %477, label %486

477:                                              ; preds = %465
  call void @llvm.lifetime.start.p0(i64 2048, ptr %56) #7
  %478 = getelementptr inbounds [2048 x i8], ptr %56, i64 0, i64 0
  %479 = load ptr, ptr %9, align 8, !tbaa !5
  %480 = getelementptr inbounds %struct.background, ptr %479, i32 0, i32 149
  %481 = getelementptr inbounds [2048 x i8], ptr %480, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %478, ptr noundef @.str, ptr noundef @.str.191, ptr noundef %481)
  %482 = load ptr, ptr %10, align 8, !tbaa !5
  %483 = getelementptr inbounds %struct.thermo, ptr %482, i32 0, i32 93
  %484 = getelementptr inbounds [2048 x i8], ptr %483, i64 0, i64 0
  %485 = getelementptr inbounds [2048 x i8], ptr %56, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %484, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2898, ptr noundef %485)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %56) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

486:                                              ; preds = %465
  %487 = load ptr, ptr %13, align 8, !tbaa !5
  %488 = load ptr, ptr %9, align 8, !tbaa !5
  %489 = getelementptr inbounds %struct.background, ptr %488, i32 0, i32 60
  %490 = load i32, ptr %489, align 4, !tbaa !37
  %491 = sext i32 %490 to i64
  %492 = getelementptr inbounds double, ptr %487, i64 %491
  %493 = load double, ptr %492, align 8, !tbaa !14
  %494 = fcmp oeq double %493, 0.000000e+00
  br i1 %494, label %495, label %503

495:                                              ; preds = %486
  call void @llvm.lifetime.start.p0(i64 2048, ptr %57) #7
  %496 = getelementptr inbounds [2048 x i8], ptr %57, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %496, ptr noundef @.str.34)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %58) #7
  %497 = getelementptr inbounds [2048 x i8], ptr %58, i64 0, i64 0
  %498 = getelementptr inbounds [2048 x i8], ptr %57, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %497, ptr noundef @.str.10, ptr noundef @.str.192, ptr noundef %498)
  %499 = load ptr, ptr %10, align 8, !tbaa !5
  %500 = getelementptr inbounds %struct.thermo, ptr %499, i32 0, i32 93
  %501 = getelementptr inbounds [2048 x i8], ptr %500, i64 0, i64 0
  %502 = getelementptr inbounds [2048 x i8], ptr %58, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %501, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2902, ptr noundef %502)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %58) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr %57) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

503:                                              ; preds = %486
  %504 = load double, ptr %23, align 8, !tbaa !14
  %505 = fadd double 1.000000e+00, %504
  %506 = load double, ptr %23, align 8, !tbaa !14
  %507 = fadd double 1.000000e+00, %506
  %508 = fmul double %505, %507
  %509 = load ptr, ptr %10, align 8, !tbaa !5
  %510 = getelementptr inbounds %struct.thermo, ptr %509, i32 0, i32 89
  %511 = load double, ptr %510, align 8, !tbaa !35
  %512 = fmul double %508, %511
  %513 = load double, ptr %25, align 8, !tbaa !14
  %514 = fmul double %512, %513
  %515 = fmul double %514, 0x3A15151E37FBF236
  %516 = fmul double %515, 0x449A22FF260E3826
  %517 = load ptr, ptr %13, align 8, !tbaa !5
  %518 = load ptr, ptr %9, align 8, !tbaa !5
  %519 = getelementptr inbounds %struct.background, ptr %518, i32 0, i32 60
  %520 = load i32, ptr %519, align 4, !tbaa !37
  %521 = sext i32 %520 to i64
  %522 = getelementptr inbounds double, ptr %517, i64 %521
  %523 = load double, ptr %522, align 8, !tbaa !14
  %524 = fdiv double %516, %523
  store double %524, ptr %27, align 8, !tbaa !14
  %525 = load double, ptr %23, align 8, !tbaa !14
  %526 = fadd double 1.000000e+00, %525
  %527 = load double, ptr %23, align 8, !tbaa !14
  %528 = fadd double 1.000000e+00, %527
  %529 = fmul double %526, %528
  %530 = load ptr, ptr %10, align 8, !tbaa !5
  %531 = getelementptr inbounds %struct.thermo, ptr %530, i32 0, i32 89
  %532 = load double, ptr %531, align 8, !tbaa !35
  %533 = fmul double %529, %532
  %534 = load double, ptr %25, align 8, !tbaa !14
  %535 = fmul double %533, %534
  %536 = fmul double %535, 0x3A15151E37FBF236
  %537 = fmul double %536, 0x449A22FF260E3826
  store double %537, ptr %34, align 8, !tbaa !14
  %538 = load double, ptr %26, align 8, !tbaa !14
  %539 = fcmp oeq double %538, 0.000000e+00
  br i1 %539, label %543, label %540

540:                                              ; preds = %503
  %541 = load double, ptr %33, align 8, !tbaa !14
  %542 = fcmp oeq double %541, 0.000000e+00
  br i1 %542, label %543, label %551

543:                                              ; preds = %540, %503
  call void @llvm.lifetime.start.p0(i64 2048, ptr %59) #7
  %544 = getelementptr inbounds [2048 x i8], ptr %59, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %544, ptr noundef @.str.34)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %60) #7
  %545 = getelementptr inbounds [2048 x i8], ptr %60, i64 0, i64 0
  %546 = getelementptr inbounds [2048 x i8], ptr %59, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %545, ptr noundef @.str.10, ptr noundef @.str.198, ptr noundef %546)
  %547 = load ptr, ptr %10, align 8, !tbaa !5
  %548 = getelementptr inbounds %struct.thermo, ptr %547, i32 0, i32 93
  %549 = getelementptr inbounds [2048 x i8], ptr %548, i64 0, i64 0
  %550 = getelementptr inbounds [2048 x i8], ptr %60, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %549, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2910, ptr noundef %550)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %60) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr %59) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

551:                                              ; preds = %540
  %552 = load double, ptr %27, align 8, !tbaa !14
  %553 = load double, ptr %26, align 8, !tbaa !14
  %554 = fsub double %552, %553
  %555 = load double, ptr %26, align 8, !tbaa !14
  %556 = fdiv double %554, %555
  %557 = call double @llvm.fabs.f64(double %556)
  %558 = load double, ptr %34, align 8, !tbaa !14
  %559 = load double, ptr %33, align 8, !tbaa !14
  %560 = fsub double %558, %559
  %561 = load double, ptr %33, align 8, !tbaa !14
  %562 = fdiv double %560, %561
  %563 = call double @llvm.fabs.f64(double %562)
  %564 = fadd double %557, %563
  store double %564, ptr %39, align 8, !tbaa !14
  %565 = load double, ptr %39, align 8, !tbaa !14
  %566 = load ptr, ptr %8, align 8, !tbaa !5
  %567 = getelementptr inbounds %struct.precision, ptr %566, i32 0, i32 45
  %568 = load double, ptr %567, align 8, !tbaa !44
  %569 = fcmp olt double %565, %568
  br i1 %569, label %570, label %674

570:                                              ; preds = %551
  %571 = load double, ptr %23, align 8, !tbaa !14
  store double %571, ptr %22, align 8, !tbaa !14
  %572 = load double, ptr %25, align 8, !tbaa !14
  store double %572, ptr %24, align 8, !tbaa !14
  %573 = load double, ptr %27, align 8, !tbaa !14
  store double %573, ptr %26, align 8, !tbaa !14
  %574 = load double, ptr %34, align 8, !tbaa !14
  store double %574, ptr %33, align 8, !tbaa !14
  %575 = load double, ptr %26, align 8, !tbaa !14
  %576 = fcmp oeq double %575, 0.000000e+00
  br i1 %576, label %580, label %577

577:                                              ; preds = %570
  %578 = load double, ptr %33, align 8, !tbaa !14
  %579 = fcmp oeq double %578, 0.000000e+00
  br i1 %579, label %580, label %590

580:                                              ; preds = %577, %570
  call void @llvm.lifetime.start.p0(i64 2048, ptr %61) #7
  %581 = getelementptr inbounds [2048 x i8], ptr %61, i64 0, i64 0
  %582 = load double, ptr %26, align 8, !tbaa !14
  %583 = load double, ptr %33, align 8, !tbaa !14
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %581, ptr noundef @.str.199, double noundef %582, double noundef %583)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %62) #7
  %584 = getelementptr inbounds [2048 x i8], ptr %62, i64 0, i64 0
  %585 = getelementptr inbounds [2048 x i8], ptr %61, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %584, ptr noundef @.str.10, ptr noundef @.str.198, ptr noundef %585)
  %586 = load ptr, ptr %10, align 8, !tbaa !5
  %587 = getelementptr inbounds %struct.thermo, ptr %586, i32 0, i32 93
  %588 = getelementptr inbounds [2048 x i8], ptr %587, i64 0, i64 0
  %589 = getelementptr inbounds [2048 x i8], ptr %62, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %588, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2925, ptr noundef %589)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %62) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr %61) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

590:                                              ; preds = %577
  %591 = load double, ptr %22, align 8, !tbaa !14
  %592 = load ptr, ptr %17, align 8, !tbaa !5
  %593 = load ptr, ptr %12, align 8, !tbaa !5
  %594 = getelementptr inbounds %struct.reionization, ptr %593, i32 0, i32 0
  %595 = load i32, ptr %594, align 8, !tbaa !29
  %596 = sext i32 %595 to i64
  %597 = getelementptr inbounds double, ptr %592, i64 %596
  store double %591, ptr %597, align 8, !tbaa !14
  %598 = load double, ptr %24, align 8, !tbaa !14
  %599 = load ptr, ptr %17, align 8, !tbaa !5
  %600 = load ptr, ptr %12, align 8, !tbaa !5
  %601 = getelementptr inbounds %struct.reionization, ptr %600, i32 0, i32 1
  %602 = load i32, ptr %601, align 4, !tbaa !31
  %603 = sext i32 %602 to i64
  %604 = getelementptr inbounds double, ptr %599, i64 %603
  store double %598, ptr %604, align 8, !tbaa !14
  %605 = load double, ptr %26, align 8, !tbaa !14
  %606 = load ptr, ptr %17, align 8, !tbaa !5
  %607 = load ptr, ptr %12, align 8, !tbaa !5
  %608 = getelementptr inbounds %struct.reionization, ptr %607, i32 0, i32 6
  %609 = load i32, ptr %608, align 8, !tbaa !38
  %610 = sext i32 %609 to i64
  %611 = getelementptr inbounds double, ptr %606, i64 %610
  store double %605, ptr %611, align 8, !tbaa !14
  %612 = load double, ptr %26, align 8, !tbaa !14
  %613 = load ptr, ptr %13, align 8, !tbaa !5
  %614 = load ptr, ptr %9, align 8, !tbaa !5
  %615 = getelementptr inbounds %struct.background, ptr %614, i32 0, i32 60
  %616 = load i32, ptr %615, align 4, !tbaa !37
  %617 = sext i32 %616 to i64
  %618 = getelementptr inbounds double, ptr %613, i64 %617
  %619 = load double, ptr %618, align 8, !tbaa !14
  %620 = fmul double %612, %619
  %621 = load ptr, ptr %17, align 8, !tbaa !5
  %622 = load ptr, ptr %12, align 8, !tbaa !5
  %623 = getelementptr inbounds %struct.reionization, ptr %622, i32 0, i32 5
  %624 = load i32, ptr %623, align 4, !tbaa !36
  %625 = sext i32 %624 to i64
  %626 = getelementptr inbounds double, ptr %621, i64 %625
  store double %620, ptr %626, align 8, !tbaa !14
  %627 = load ptr, ptr %17, align 8, !tbaa !5
  %628 = load ptr, ptr %12, align 8, !tbaa !5
  %629 = getelementptr inbounds %struct.reionization, ptr %628, i32 0, i32 8
  %630 = load i32, ptr %629, align 8, !tbaa !15
  %631 = sext i32 %630 to i64
  %632 = mul i64 8, %631
  %633 = call i32 @gt_add(ptr noundef %14, i64 noundef -1, ptr noundef %627, i64 noundef %632)
  %634 = icmp eq i32 %633, 1
  br i1 %634, label %635, label %643

635:                                              ; preds = %590
  call void @llvm.lifetime.start.p0(i64 2048, ptr %63) #7
  %636 = getelementptr inbounds [2048 x i8], ptr %63, i64 0, i64 0
  %637 = getelementptr inbounds %struct.growTable, ptr %14, i32 0, i32 4
  %638 = getelementptr inbounds [2048 x i8], ptr %637, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %636, ptr noundef @.str, ptr noundef @.str.193, ptr noundef %638)
  %639 = load ptr, ptr %10, align 8, !tbaa !5
  %640 = getelementptr inbounds %struct.thermo, ptr %639, i32 0, i32 93
  %641 = getelementptr inbounds [2048 x i8], ptr %640, i64 0, i64 0
  %642 = getelementptr inbounds [2048 x i8], ptr %63, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %641, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2934, ptr noundef %642)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %63) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

643:                                              ; preds = %590
  %644 = load i32, ptr %19, align 4, !tbaa !17
  %645 = add nsw i32 %644, 1
  store i32 %645, ptr %19, align 4, !tbaa !17
  %646 = load ptr, ptr %8, align 8, !tbaa !5
  %647 = getelementptr inbounds %struct.precision, ptr %646, i32 0, i32 45
  %648 = load double, ptr %647, align 8, !tbaa !44
  %649 = load double, ptr %39, align 8, !tbaa !14
  %650 = fdiv double %648, %649
  %651 = fmul double 9.000000e-01, %650
  %652 = fcmp olt double %651, 5.000000e+00
  br i1 %652, label %653, label %660

653:                                              ; preds = %643
  %654 = load ptr, ptr %8, align 8, !tbaa !5
  %655 = getelementptr inbounds %struct.precision, ptr %654, i32 0, i32 45
  %656 = load double, ptr %655, align 8, !tbaa !44
  %657 = load double, ptr %39, align 8, !tbaa !14
  %658 = fdiv double %656, %657
  %659 = fmul double 9.000000e-01, %658
  br label %661

660:                                              ; preds = %643
  br label %661

661:                                              ; preds = %660, %653
  %662 = phi double [ %659, %653 ], [ 5.000000e+00, %660 ]
  %663 = load double, ptr %20, align 8, !tbaa !14
  %664 = fmul double %662, %663
  store double %664, ptr %20, align 8, !tbaa !14
  %665 = load double, ptr %20, align 8, !tbaa !14
  %666 = load double, ptr %21, align 8, !tbaa !14
  %667 = fcmp olt double %665, %666
  br i1 %667, label %668, label %670

668:                                              ; preds = %661
  %669 = load double, ptr %20, align 8, !tbaa !14
  br label %672

670:                                              ; preds = %661
  %671 = load double, ptr %21, align 8, !tbaa !14
  br label %672

672:                                              ; preds = %670, %668
  %673 = phi double [ %669, %668 ], [ %671, %670 ]
  store double %673, ptr %20, align 8, !tbaa !14
  br label %683

674:                                              ; preds = %551
  %675 = load ptr, ptr %8, align 8, !tbaa !5
  %676 = getelementptr inbounds %struct.precision, ptr %675, i32 0, i32 45
  %677 = load double, ptr %676, align 8, !tbaa !44
  %678 = load double, ptr %39, align 8, !tbaa !14
  %679 = fdiv double %677, %678
  %680 = fmul double 9.000000e-01, %679
  %681 = load double, ptr %20, align 8, !tbaa !14
  %682 = fmul double %680, %681
  store double %682, ptr %20, align 8, !tbaa !14
  br label %683

683:                                              ; preds = %674, %672
  br label %412, !llvm.loop !45

684:                                              ; preds = %412
  %685 = load ptr, ptr %12, align 8, !tbaa !5
  %686 = getelementptr inbounds %struct.reionization, ptr %685, i32 0, i32 8
  %687 = load i32, ptr %686, align 8, !tbaa !15
  %688 = load i32, ptr %19, align 4, !tbaa !17
  %689 = mul nsw i32 %687, %688
  %690 = sext i32 %689 to i64
  %691 = mul i64 %690, 8
  %692 = call noalias ptr @malloc(i64 noundef %691) #8
  %693 = load ptr, ptr %12, align 8, !tbaa !5
  %694 = getelementptr inbounds %struct.reionization, ptr %693, i32 0, i32 10
  store ptr %692, ptr %694, align 8, !tbaa !46
  %695 = load ptr, ptr %12, align 8, !tbaa !5
  %696 = getelementptr inbounds %struct.reionization, ptr %695, i32 0, i32 10
  %697 = load ptr, ptr %696, align 8, !tbaa !46
  %698 = icmp eq ptr %697, null
  br i1 %698, label %699, label %714

699:                                              ; preds = %684
  call void @llvm.lifetime.start.p0(i64 4, ptr %64) #7
  %700 = load ptr, ptr %12, align 8, !tbaa !5
  %701 = getelementptr inbounds %struct.reionization, ptr %700, i32 0, i32 8
  %702 = load i32, ptr %701, align 8, !tbaa !15
  %703 = load i32, ptr %19, align 4, !tbaa !17
  %704 = mul nsw i32 %702, %703
  %705 = sext i32 %704 to i64
  %706 = mul i64 %705, 8
  %707 = trunc i64 %706 to i32
  store i32 %707, ptr %64, align 4, !tbaa !17
  call void @llvm.lifetime.start.p0(i64 2048, ptr %65) #7
  %708 = getelementptr inbounds [2048 x i8], ptr %65, i64 0, i64 0
  %709 = load i32, ptr %64, align 4, !tbaa !17
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %708, ptr noundef @.str.37, ptr noundef @.str.200, i32 noundef %709)
  %710 = load ptr, ptr %10, align 8, !tbaa !5
  %711 = getelementptr inbounds %struct.thermo, ptr %710, i32 0, i32 93
  %712 = getelementptr inbounds [2048 x i8], ptr %711, i64 0, i64 0
  %713 = getelementptr inbounds [2048 x i8], ptr %65, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %712, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2948, ptr noundef %713)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %65) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %64) #7
  br label %1323

714:                                              ; preds = %684
  %715 = load i32, ptr %19, align 4, !tbaa !17
  %716 = load ptr, ptr %12, align 8, !tbaa !5
  %717 = getelementptr inbounds %struct.reionization, ptr %716, i32 0, i32 9
  store i32 %715, ptr %717, align 4, !tbaa !47
  %718 = call i32 @gt_getPtr(ptr noundef %14, ptr noundef %15)
  %719 = icmp eq i32 %718, 1
  br i1 %719, label %720, label %728

720:                                              ; preds = %714
  call void @llvm.lifetime.start.p0(i64 2048, ptr %66) #7
  %721 = getelementptr inbounds [2048 x i8], ptr %66, i64 0, i64 0
  %722 = getelementptr inbounds %struct.growTable, ptr %14, i32 0, i32 4
  %723 = getelementptr inbounds [2048 x i8], ptr %722, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %721, ptr noundef @.str, ptr noundef @.str.201, ptr noundef %723)
  %724 = load ptr, ptr %10, align 8, !tbaa !5
  %725 = getelementptr inbounds %struct.thermo, ptr %724, i32 0, i32 93
  %726 = getelementptr inbounds [2048 x i8], ptr %725, i64 0, i64 0
  %727 = getelementptr inbounds [2048 x i8], ptr %66, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %726, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2955, ptr noundef %727)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %66) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

728:                                              ; preds = %714
  store i32 0, ptr %18, align 4, !tbaa !17
  br label %729

729:                                              ; preds = %785, %728
  %730 = load i32, ptr %18, align 4, !tbaa !17
  %731 = load ptr, ptr %12, align 8, !tbaa !5
  %732 = getelementptr inbounds %struct.reionization, ptr %731, i32 0, i32 9
  %733 = load i32, ptr %732, align 4, !tbaa !47
  %734 = icmp slt i32 %730, %733
  br i1 %734, label %735, label %788

735:                                              ; preds = %729
  %736 = load ptr, ptr %12, align 8, !tbaa !5
  %737 = getelementptr inbounds %struct.reionization, ptr %736, i32 0, i32 10
  %738 = load ptr, ptr %737, align 8, !tbaa !46
  %739 = load i32, ptr %18, align 4, !tbaa !17
  %740 = load ptr, ptr %12, align 8, !tbaa !5
  %741 = getelementptr inbounds %struct.reionization, ptr %740, i32 0, i32 8
  %742 = load i32, ptr %741, align 8, !tbaa !15
  %743 = mul nsw i32 %739, %742
  %744 = sext i32 %743 to i64
  %745 = getelementptr inbounds double, ptr %738, i64 %744
  %746 = load ptr, ptr %15, align 8, !tbaa !5
  %747 = load ptr, ptr %12, align 8, !tbaa !5
  %748 = getelementptr inbounds %struct.reionization, ptr %747, i32 0, i32 9
  %749 = load i32, ptr %748, align 4, !tbaa !47
  %750 = load i32, ptr %18, align 4, !tbaa !17
  %751 = sub nsw i32 %749, %750
  %752 = sub nsw i32 %751, 1
  %753 = load ptr, ptr %12, align 8, !tbaa !5
  %754 = getelementptr inbounds %struct.reionization, ptr %753, i32 0, i32 8
  %755 = load i32, ptr %754, align 8, !tbaa !15
  %756 = mul nsw i32 %752, %755
  %757 = sext i32 %756 to i64
  %758 = getelementptr inbounds double, ptr %746, i64 %757
  %759 = load ptr, ptr %12, align 8, !tbaa !5
  %760 = getelementptr inbounds %struct.reionization, ptr %759, i32 0, i32 8
  %761 = load i32, ptr %760, align 8, !tbaa !15
  %762 = sext i32 %761 to i64
  %763 = mul i64 %762, 8
  call void @llvm.memcpy.p0.p0.i64(ptr align 8 %745, ptr align 8 %758, i64 %763, i1 false)
  store ptr %745, ptr %16, align 8, !tbaa !5
  %764 = load ptr, ptr %16, align 8, !tbaa !5
  %765 = load ptr, ptr %12, align 8, !tbaa !5
  %766 = getelementptr inbounds %struct.reionization, ptr %765, i32 0, i32 10
  %767 = load ptr, ptr %766, align 8, !tbaa !46
  %768 = load i32, ptr %18, align 4, !tbaa !17
  %769 = load ptr, ptr %12, align 8, !tbaa !5
  %770 = getelementptr inbounds %struct.reionization, ptr %769, i32 0, i32 8
  %771 = load i32, ptr %770, align 8, !tbaa !15
  %772 = mul nsw i32 %768, %771
  %773 = sext i32 %772 to i64
  %774 = getelementptr inbounds double, ptr %767, i64 %773
  %775 = icmp ne ptr %764, %774
  br i1 %775, label %776, label %784

776:                                              ; preds = %735
  call void @llvm.lifetime.start.p0(i64 2048, ptr %67) #7
  %777 = getelementptr inbounds [2048 x i8], ptr %67, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %777, ptr noundef @.str.202)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %68) #7
  %778 = getelementptr inbounds [2048 x i8], ptr %68, i64 0, i64 0
  %779 = getelementptr inbounds [2048 x i8], ptr %67, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %778, ptr noundef @.str.10, ptr noundef @.str.203, ptr noundef %779)
  %780 = load ptr, ptr %10, align 8, !tbaa !5
  %781 = getelementptr inbounds %struct.thermo, ptr %780, i32 0, i32 93
  %782 = getelementptr inbounds [2048 x i8], ptr %781, i64 0, i64 0
  %783 = getelementptr inbounds [2048 x i8], ptr %68, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %782, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2962, ptr noundef %783)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %68) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr %67) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

784:                                              ; preds = %735
  br label %785

785:                                              ; preds = %784
  %786 = load i32, ptr %18, align 4, !tbaa !17
  %787 = add nsw i32 %786, 1
  store i32 %787, ptr %18, align 4, !tbaa !17
  br label %729, !llvm.loop !48

788:                                              ; preds = %729
  %789 = call i32 @gt_free(ptr noundef %14)
  %790 = icmp eq i32 %789, 1
  br i1 %790, label %791, label %799

791:                                              ; preds = %788
  call void @llvm.lifetime.start.p0(i64 2048, ptr %69) #7
  %792 = getelementptr inbounds [2048 x i8], ptr %69, i64 0, i64 0
  %793 = getelementptr inbounds %struct.growTable, ptr %14, i32 0, i32 4
  %794 = getelementptr inbounds [2048 x i8], ptr %793, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %792, ptr noundef @.str, ptr noundef @.str.204, ptr noundef %794)
  %795 = load ptr, ptr %10, align 8, !tbaa !5
  %796 = getelementptr inbounds %struct.thermo, ptr %795, i32 0, i32 93
  %797 = getelementptr inbounds [2048 x i8], ptr %796, i64 0, i64 0
  %798 = getelementptr inbounds [2048 x i8], ptr %69, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %797, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2969, ptr noundef %798)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %69) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

799:                                              ; preds = %788
  %800 = load ptr, ptr %17, align 8, !tbaa !5
  call void @free(ptr noundef %800) #7
  %801 = load ptr, ptr %12, align 8, !tbaa !5
  %802 = getelementptr inbounds %struct.reionization, ptr %801, i32 0, i32 9
  %803 = load i32, ptr %802, align 4, !tbaa !47
  %804 = sub nsw i32 %803, 1
  store i32 %804, ptr %18, align 4, !tbaa !17
  br label %805

805:                                              ; preds = %1251, %799
  %806 = load i32, ptr %18, align 4, !tbaa !17
  %807 = icmp sgt i32 %806, 0
  br i1 %807, label %808, label %1254

808:                                              ; preds = %805
  %809 = load ptr, ptr %12, align 8, !tbaa !5
  %810 = getelementptr inbounds %struct.reionization, ptr %809, i32 0, i32 10
  %811 = load ptr, ptr %810, align 8, !tbaa !46
  %812 = load i32, ptr %18, align 4, !tbaa !17
  %813 = load ptr, ptr %12, align 8, !tbaa !5
  %814 = getelementptr inbounds %struct.reionization, ptr %813, i32 0, i32 8
  %815 = load i32, ptr %814, align 8, !tbaa !15
  %816 = mul nsw i32 %812, %815
  %817 = load ptr, ptr %12, align 8, !tbaa !5
  %818 = getelementptr inbounds %struct.reionization, ptr %817, i32 0, i32 0
  %819 = load i32, ptr %818, align 8, !tbaa !29
  %820 = add nsw i32 %816, %819
  %821 = sext i32 %820 to i64
  %822 = getelementptr inbounds double, ptr %811, i64 %821
  %823 = load double, ptr %822, align 8, !tbaa !14
  store double %823, ptr %22, align 8, !tbaa !14
  %824 = load ptr, ptr %9, align 8, !tbaa !5
  %825 = load double, ptr %22, align 8, !tbaa !14
  %826 = call i32 @background_tau_of_z(ptr noundef %824, double noundef %825, ptr noundef %36)
  %827 = icmp eq i32 %826, 1
  br i1 %827, label %828, label %837

828:                                              ; preds = %808
  call void @llvm.lifetime.start.p0(i64 2048, ptr %70) #7
  %829 = getelementptr inbounds [2048 x i8], ptr %70, i64 0, i64 0
  %830 = load ptr, ptr %9, align 8, !tbaa !5
  %831 = getelementptr inbounds %struct.background, ptr %830, i32 0, i32 149
  %832 = getelementptr inbounds [2048 x i8], ptr %831, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %829, ptr noundef @.str, ptr noundef @.str.190, ptr noundef %832)
  %833 = load ptr, ptr %10, align 8, !tbaa !5
  %834 = getelementptr inbounds %struct.thermo, ptr %833, i32 0, i32 93
  %835 = getelementptr inbounds [2048 x i8], ptr %834, i64 0, i64 0
  %836 = getelementptr inbounds [2048 x i8], ptr %70, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %835, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2982, ptr noundef %836)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %70) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

837:                                              ; preds = %808
  %838 = load ptr, ptr %9, align 8, !tbaa !5
  %839 = load double, ptr %36, align 8, !tbaa !14
  %840 = load ptr, ptr %9, align 8, !tbaa !5
  %841 = getelementptr inbounds %struct.background, ptr %840, i32 0, i32 142
  %842 = load i16, ptr %841, align 2, !tbaa !49
  %843 = load ptr, ptr %9, align 8, !tbaa !5
  %844 = getelementptr inbounds %struct.background, ptr %843, i32 0, i32 144
  %845 = load i16, ptr %844, align 2, !tbaa !34
  %846 = load ptr, ptr %13, align 8, !tbaa !5
  %847 = call i32 @background_at_tau(ptr noundef %838, double noundef %839, i16 noundef signext %842, i16 noundef signext %845, ptr noundef %38, ptr noundef %846)
  %848 = icmp eq i32 %847, 1
  br i1 %848, label %849, label %858

849:                                              ; preds = %837
  call void @llvm.lifetime.start.p0(i64 2048, ptr %71) #7
  %850 = getelementptr inbounds [2048 x i8], ptr %71, i64 0, i64 0
  %851 = load ptr, ptr %9, align 8, !tbaa !5
  %852 = getelementptr inbounds %struct.background, ptr %851, i32 0, i32 149
  %853 = getelementptr inbounds [2048 x i8], ptr %852, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %850, ptr noundef @.str, ptr noundef @.str.205, ptr noundef %853)
  %854 = load ptr, ptr %10, align 8, !tbaa !5
  %855 = getelementptr inbounds %struct.thermo, ptr %854, i32 0, i32 93
  %856 = getelementptr inbounds [2048 x i8], ptr %855, i64 0, i64 0
  %857 = getelementptr inbounds [2048 x i8], ptr %71, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %856, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 2991, ptr noundef %857)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %71) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

858:                                              ; preds = %837
  %859 = load ptr, ptr %12, align 8, !tbaa !5
  %860 = getelementptr inbounds %struct.reionization, ptr %859, i32 0, i32 10
  %861 = load ptr, ptr %860, align 8, !tbaa !46
  %862 = load i32, ptr %18, align 4, !tbaa !17
  %863 = load ptr, ptr %12, align 8, !tbaa !5
  %864 = getelementptr inbounds %struct.reionization, ptr %863, i32 0, i32 8
  %865 = load i32, ptr %864, align 8, !tbaa !15
  %866 = mul nsw i32 %862, %865
  %867 = load ptr, ptr %12, align 8, !tbaa !5
  %868 = getelementptr inbounds %struct.reionization, ptr %867, i32 0, i32 0
  %869 = load i32, ptr %868, align 8, !tbaa !29
  %870 = add nsw i32 %866, %869
  %871 = sext i32 %870 to i64
  %872 = getelementptr inbounds double, ptr %861, i64 %871
  %873 = load double, ptr %872, align 8, !tbaa !14
  %874 = load ptr, ptr %12, align 8, !tbaa !5
  %875 = getelementptr inbounds %struct.reionization, ptr %874, i32 0, i32 10
  %876 = load ptr, ptr %875, align 8, !tbaa !46
  %877 = load i32, ptr %18, align 4, !tbaa !17
  %878 = sub nsw i32 %877, 1
  %879 = load ptr, ptr %12, align 8, !tbaa !5
  %880 = getelementptr inbounds %struct.reionization, ptr %879, i32 0, i32 8
  %881 = load i32, ptr %880, align 8, !tbaa !15
  %882 = mul nsw i32 %878, %881
  %883 = load ptr, ptr %12, align 8, !tbaa !5
  %884 = getelementptr inbounds %struct.reionization, ptr %883, i32 0, i32 0
  %885 = load i32, ptr %884, align 8, !tbaa !29
  %886 = add nsw i32 %882, %885
  %887 = sext i32 %886 to i64
  %888 = getelementptr inbounds double, ptr %876, i64 %887
  %889 = load double, ptr %888, align 8, !tbaa !14
  %890 = fsub double %873, %889
  store double %890, ptr %20, align 8, !tbaa !14
  %891 = load double, ptr %22, align 8, !tbaa !14
  %892 = fadd double 1.000000e+00, %891
  %893 = load double, ptr %22, align 8, !tbaa !14
  %894 = fadd double 1.000000e+00, %893
  %895 = fmul double %892, %894
  %896 = load ptr, ptr %10, align 8, !tbaa !5
  %897 = getelementptr inbounds %struct.thermo, ptr %896, i32 0, i32 89
  %898 = load double, ptr %897, align 8, !tbaa !35
  %899 = fmul double %895, %898
  %900 = load ptr, ptr %12, align 8, !tbaa !5
  %901 = getelementptr inbounds %struct.reionization, ptr %900, i32 0, i32 10
  %902 = load ptr, ptr %901, align 8, !tbaa !46
  %903 = load i32, ptr %18, align 4, !tbaa !17
  %904 = load ptr, ptr %12, align 8, !tbaa !5
  %905 = getelementptr inbounds %struct.reionization, ptr %904, i32 0, i32 8
  %906 = load i32, ptr %905, align 8, !tbaa !15
  %907 = mul nsw i32 %903, %906
  %908 = load ptr, ptr %12, align 8, !tbaa !5
  %909 = getelementptr inbounds %struct.reionization, ptr %908, i32 0, i32 1
  %910 = load i32, ptr %909, align 4, !tbaa !31
  %911 = add nsw i32 %907, %910
  %912 = sext i32 %911 to i64
  %913 = getelementptr inbounds double, ptr %902, i64 %912
  %914 = load double, ptr %913, align 8, !tbaa !14
  %915 = fmul double %899, %914
  %916 = fmul double %915, 0x3A15151E37FBF236
  %917 = fmul double %916, 0x449A22FF260E3826
  store double %917, ptr %31, align 8, !tbaa !14
  %918 = load ptr, ptr %10, align 8, !tbaa !5
  %919 = getelementptr inbounds %struct.thermo, ptr %918, i32 0, i32 0
  %920 = load double, ptr %919, align 8, !tbaa !9
  %921 = call double @llvm.fmuladd.f64(double 0xBFE7F14DA42541E8, double %920, double 1.000000e+00)
  %922 = load ptr, ptr %12, align 8, !tbaa !5
  %923 = getelementptr inbounds %struct.reionization, ptr %922, i32 0, i32 10
  %924 = load ptr, ptr %923, align 8, !tbaa !46
  %925 = load i32, ptr %18, align 4, !tbaa !17
  %926 = load ptr, ptr %12, align 8, !tbaa !5
  %927 = getelementptr inbounds %struct.reionization, ptr %926, i32 0, i32 8
  %928 = load i32, ptr %927, align 8, !tbaa !15
  %929 = mul nsw i32 %925, %928
  %930 = load ptr, ptr %12, align 8, !tbaa !5
  %931 = getelementptr inbounds %struct.reionization, ptr %930, i32 0, i32 1
  %932 = load i32, ptr %931, align 4, !tbaa !31
  %933 = add nsw i32 %929, %932
  %934 = sext i32 %933 to i64
  %935 = getelementptr inbounds double, ptr %924, i64 %934
  %936 = load double, ptr %935, align 8, !tbaa !14
  %937 = load ptr, ptr %10, align 8, !tbaa !5
  %938 = getelementptr inbounds %struct.thermo, ptr %937, i32 0, i32 0
  %939 = load double, ptr %938, align 8, !tbaa !9
  %940 = fsub double 1.000000e+00, %939
  %941 = call double @llvm.fmuladd.f64(double %936, double %940, double %921)
  %942 = fdiv double 0x3A60930446E176B5, %941
  store double %942, ptr %32, align 8, !tbaa !14
  %943 = load double, ptr %22, align 8, !tbaa !14
  %944 = fadd double 1.000000e+00, %943
  %945 = fdiv double 2.000000e+00, %944
  %946 = load ptr, ptr %12, align 8, !tbaa !5
  %947 = getelementptr inbounds %struct.reionization, ptr %946, i32 0, i32 10
  %948 = load ptr, ptr %947, align 8, !tbaa !46
  %949 = load i32, ptr %18, align 4, !tbaa !17
  %950 = load ptr, ptr %12, align 8, !tbaa !5
  %951 = getelementptr inbounds %struct.reionization, ptr %950, i32 0, i32 8
  %952 = load i32, ptr %951, align 8, !tbaa !15
  %953 = mul nsw i32 %949, %952
  %954 = load ptr, ptr %12, align 8, !tbaa !5
  %955 = getelementptr inbounds %struct.reionization, ptr %954, i32 0, i32 2
  %956 = load i32, ptr %955, align 8, !tbaa !40
  %957 = add nsw i32 %953, %956
  %958 = sext i32 %957 to i64
  %959 = getelementptr inbounds double, ptr %948, i64 %958
  %960 = load double, ptr %959, align 8, !tbaa !14
  %961 = load double, ptr %32, align 8, !tbaa !14
  %962 = fmul double 2.000000e+00, %961
  %963 = fdiv double %962, 0x39B279DC93EB3F4E
  %964 = fmul double %963, 4.000000e+00
  %965 = load ptr, ptr %13, align 8, !tbaa !5
  %966 = load ptr, ptr %9, align 8, !tbaa !5
  %967 = getelementptr inbounds %struct.background, ptr %966, i32 0, i32 62
  %968 = load i32, ptr %967, align 4, !tbaa !50
  %969 = sext i32 %968 to i64
  %970 = getelementptr inbounds double, ptr %965, i64 %969
  %971 = load double, ptr %970, align 8, !tbaa !14
  %972 = fmul double %964, %971
  %973 = fdiv double %972, 3.000000e+00
  %974 = load ptr, ptr %13, align 8, !tbaa !5
  %975 = load ptr, ptr %9, align 8, !tbaa !5
  %976 = getelementptr inbounds %struct.background, ptr %975, i32 0, i32 63
  %977 = load i32, ptr %976, align 8, !tbaa !51
  %978 = sext i32 %977 to i64
  %979 = getelementptr inbounds double, ptr %974, i64 %978
  %980 = load double, ptr %979, align 8, !tbaa !14
  %981 = fdiv double %973, %980
  %982 = load double, ptr %31, align 8, !tbaa !14
  %983 = fmul double %981, %982
  %984 = load ptr, ptr %9, align 8, !tbaa !5
  %985 = getelementptr inbounds %struct.background, ptr %984, i32 0, i32 2
  %986 = load double, ptr %985, align 8, !tbaa !52
  %987 = load double, ptr %22, align 8, !tbaa !14
  %988 = fadd double 1.000000e+00, %987
  %989 = load ptr, ptr %12, align 8, !tbaa !5
  %990 = getelementptr inbounds %struct.reionization, ptr %989, i32 0, i32 10
  %991 = load ptr, ptr %990, align 8, !tbaa !46
  %992 = load i32, ptr %18, align 4, !tbaa !17
  %993 = load ptr, ptr %12, align 8, !tbaa !5
  %994 = getelementptr inbounds %struct.reionization, ptr %993, i32 0, i32 8
  %995 = load i32, ptr %994, align 8, !tbaa !15
  %996 = mul nsw i32 %992, %995
  %997 = load ptr, ptr %12, align 8, !tbaa !5
  %998 = getelementptr inbounds %struct.reionization, ptr %997, i32 0, i32 2
  %999 = load i32, ptr %998, align 8, !tbaa !40
  %1000 = add nsw i32 %996, %999
  %1001 = sext i32 %1000 to i64
  %1002 = getelementptr inbounds double, ptr %991, i64 %1001
  %1003 = load double, ptr %1002, align 8, !tbaa !14
  %1004 = fneg double %1003
  %1005 = call double @llvm.fmuladd.f64(double %986, double %988, double %1004)
  %1006 = fmul double %983, %1005
  %1007 = load ptr, ptr %13, align 8, !tbaa !5
  %1008 = load ptr, ptr %9, align 8, !tbaa !5
  %1009 = getelementptr inbounds %struct.background, ptr %1008, i32 0, i32 60
  %1010 = load i32, ptr %1009, align 4, !tbaa !37
  %1011 = sext i32 %1010 to i64
  %1012 = getelementptr inbounds double, ptr %1007, i64 %1011
  %1013 = load double, ptr %1012, align 8, !tbaa !14
  %1014 = fdiv double %1006, %1013
  %1015 = fneg double %1014
  %1016 = call double @llvm.fmuladd.f64(double %945, double %960, double %1015)
  store double %1016, ptr %30, align 8, !tbaa !14
  %1017 = load ptr, ptr %11, align 8, !tbaa !5
  %1018 = getelementptr inbounds %struct.recombination, ptr %1017, i32 0, i32 28
  %1019 = load double, ptr %1018, align 8, !tbaa !53
  %1020 = fcmp ogt double %1019, 0.000000e+00
  br i1 %1020, label %1021, label %1126

1021:                                             ; preds = %858
  %1022 = load ptr, ptr %8, align 8, !tbaa !5
  %1023 = load ptr, ptr %9, align 8, !tbaa !5
  %1024 = load ptr, ptr %11, align 8, !tbaa !5
  %1025 = load double, ptr %22, align 8, !tbaa !14
  %1026 = load ptr, ptr %10, align 8, !tbaa !5
  %1027 = getelementptr inbounds %struct.thermo, ptr %1026, i32 0, i32 93
  %1028 = getelementptr inbounds [2048 x i8], ptr %1027, i64 0, i64 0
  %1029 = call i32 @thermodynamics_energy_injection(ptr noundef %1022, ptr noundef %1023, ptr noundef %1024, double noundef %1025, ptr noundef %35, ptr noundef %1028)
  %1030 = icmp eq i32 %1029, 1
  br i1 %1030, label %1031, label %1040

1031:                                             ; preds = %1021
  call void @llvm.lifetime.start.p0(i64 2048, ptr %72) #7
  %1032 = getelementptr inbounds [2048 x i8], ptr %72, i64 0, i64 0
  %1033 = load ptr, ptr %10, align 8, !tbaa !5
  %1034 = getelementptr inbounds %struct.thermo, ptr %1033, i32 0, i32 93
  %1035 = getelementptr inbounds [2048 x i8], ptr %1034, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1032, ptr noundef @.str, ptr noundef @.str.206, ptr noundef %1035)
  %1036 = load ptr, ptr %10, align 8, !tbaa !5
  %1037 = getelementptr inbounds %struct.thermo, ptr %1036, i32 0, i32 93
  %1038 = getelementptr inbounds [2048 x i8], ptr %1037, i64 0, i64 0
  %1039 = getelementptr inbounds [2048 x i8], ptr %72, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1038, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 3010, ptr noundef %1039)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %72) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

1040:                                             ; preds = %1021
  %1041 = load ptr, ptr %12, align 8, !tbaa !5
  %1042 = getelementptr inbounds %struct.reionization, ptr %1041, i32 0, i32 10
  %1043 = load ptr, ptr %1042, align 8, !tbaa !46
  %1044 = load i32, ptr %18, align 4, !tbaa !17
  %1045 = load ptr, ptr %12, align 8, !tbaa !5
  %1046 = getelementptr inbounds %struct.reionization, ptr %1045, i32 0, i32 8
  %1047 = load i32, ptr %1046, align 8, !tbaa !15
  %1048 = mul nsw i32 %1044, %1047
  %1049 = load ptr, ptr %12, align 8, !tbaa !5
  %1050 = getelementptr inbounds %struct.reionization, ptr %1049, i32 0, i32 1
  %1051 = load i32, ptr %1050, align 4, !tbaa !31
  %1052 = add nsw i32 %1048, %1051
  %1053 = sext i32 %1052 to i64
  %1054 = getelementptr inbounds double, ptr %1043, i64 %1053
  %1055 = load double, ptr %1054, align 8, !tbaa !14
  store double %1055, ptr %24, align 8, !tbaa !14
  %1056 = load double, ptr %24, align 8, !tbaa !14
  %1057 = fcmp olt double %1056, 1.000000e+00
  br i1 %1057, label %1058, label %1076

1058:                                             ; preds = %1040
  %1059 = load double, ptr %24, align 8, !tbaa !14
  %1060 = call double @pow(double noundef %1059, double noundef 3.001340e-01) #7
  %1061 = fsub double 1.000000e+00, %1060
  %1062 = call double @pow(double noundef %1061, double noundef 1.510350e+00) #7
  %1063 = fsub double 1.000000e+00, %1062
  %1064 = fmul double 0x3FEFE640A6B93CCD, %1063
  %1065 = fcmp olt double %1064, 1.000000e+00
  br i1 %1065, label %1066, label %1073

1066:                                             ; preds = %1058
  %1067 = load double, ptr %24, align 8, !tbaa !14
  %1068 = call double @pow(double noundef %1067, double noundef 3.001340e-01) #7
  %1069 = fsub double 1.000000e+00, %1068
  %1070 = call double @pow(double noundef %1069, double noundef 1.510350e+00) #7
  %1071 = fsub double 1.000000e+00, %1070
  %1072 = fmul double 0x3FEFE640A6B93CCD, %1071
  br label %1074

1073:                                             ; preds = %1058
  br label %1074

1074:                                             ; preds = %1073, %1066
  %1075 = phi double [ %1072, %1066 ], [ 1.000000e+00, %1073 ]
  store double %1075, ptr %37, align 8, !tbaa !14
  br label %1077

1076:                                             ; preds = %1040
  store double 1.000000e+00, ptr %37, align 8, !tbaa !14
  br label %1077

1077:                                             ; preds = %1076, %1074
  %1078 = load double, ptr %35, align 8, !tbaa !14
  %1079 = fmul double 0xC4A4733969CB9B70, %1078
  %1080 = load double, ptr %37, align 8, !tbaa !14
  %1081 = fmul double %1079, %1080
  %1082 = load ptr, ptr %11, align 8, !tbaa !5
  %1083 = getelementptr inbounds %struct.recombination, ptr %1082, i32 0, i32 21
  %1084 = load double, ptr %1083, align 8, !tbaa !54
  %1085 = load double, ptr %22, align 8, !tbaa !14
  %1086 = fadd double 1.000000e+00, %1085
  %1087 = call double @pow(double noundef %1086, double noundef 3.000000e+00) #7
  %1088 = fmul double %1084, %1087
  %1089 = fdiv double %1081, %1088
  %1090 = load ptr, ptr %11, align 8, !tbaa !5
  %1091 = getelementptr inbounds %struct.recombination, ptr %1090, i32 0, i32 14
  %1092 = load double, ptr %1091, align 8, !tbaa !55
  %1093 = fadd double 1.000000e+00, %1092
  %1094 = load ptr, ptr %12, align 8, !tbaa !5
  %1095 = getelementptr inbounds %struct.reionization, ptr %1094, i32 0, i32 10
  %1096 = load ptr, ptr %1095, align 8, !tbaa !46
  %1097 = load i32, ptr %18, align 4, !tbaa !17
  %1098 = load ptr, ptr %12, align 8, !tbaa !5
  %1099 = getelementptr inbounds %struct.reionization, ptr %1098, i32 0, i32 8
  %1100 = load i32, ptr %1099, align 8, !tbaa !15
  %1101 = mul nsw i32 %1097, %1100
  %1102 = load ptr, ptr %12, align 8, !tbaa !5
  %1103 = getelementptr inbounds %struct.reionization, ptr %1102, i32 0, i32 1
  %1104 = load i32, ptr %1103, align 4, !tbaa !31
  %1105 = add nsw i32 %1101, %1104
  %1106 = sext i32 %1105 to i64
  %1107 = getelementptr inbounds double, ptr %1096, i64 %1106
  %1108 = load double, ptr %1107, align 8, !tbaa !14
  %1109 = fadd double %1093, %1108
  %1110 = fdiv double %1089, %1109
  %1111 = load ptr, ptr %13, align 8, !tbaa !5
  %1112 = load ptr, ptr %9, align 8, !tbaa !5
  %1113 = getelementptr inbounds %struct.background, ptr %1112, i32 0, i32 60
  %1114 = load i32, ptr %1113, align 4, !tbaa !37
  %1115 = sext i32 %1114 to i64
  %1116 = getelementptr inbounds double, ptr %1111, i64 %1115
  %1117 = load double, ptr %1116, align 8, !tbaa !14
  %1118 = fmul double %1117, 0x41B1DE784A000000
  %1119 = fdiv double %1118, 0x449A22FF260E3826
  %1120 = load double, ptr %22, align 8, !tbaa !14
  %1121 = fadd double 1.000000e+00, %1120
  %1122 = fmul double %1119, %1121
  %1123 = fdiv double %1110, %1122
  %1124 = load double, ptr %30, align 8, !tbaa !14
  %1125 = fadd double %1124, %1123
  store double %1125, ptr %30, align 8, !tbaa !14
  br label %1126

1126:                                             ; preds = %1077, %858
  %1127 = load ptr, ptr %12, align 8, !tbaa !5
  %1128 = getelementptr inbounds %struct.reionization, ptr %1127, i32 0, i32 10
  %1129 = load ptr, ptr %1128, align 8, !tbaa !46
  %1130 = load i32, ptr %18, align 4, !tbaa !17
  %1131 = load ptr, ptr %12, align 8, !tbaa !5
  %1132 = getelementptr inbounds %struct.reionization, ptr %1131, i32 0, i32 8
  %1133 = load i32, ptr %1132, align 8, !tbaa !15
  %1134 = mul nsw i32 %1130, %1133
  %1135 = load ptr, ptr %12, align 8, !tbaa !5
  %1136 = getelementptr inbounds %struct.reionization, ptr %1135, i32 0, i32 2
  %1137 = load i32, ptr %1136, align 8, !tbaa !40
  %1138 = add nsw i32 %1134, %1137
  %1139 = sext i32 %1138 to i64
  %1140 = getelementptr inbounds double, ptr %1129, i64 %1139
  %1141 = load double, ptr %1140, align 8, !tbaa !14
  %1142 = load double, ptr %30, align 8, !tbaa !14
  %1143 = load double, ptr %20, align 8, !tbaa !14
  %1144 = fneg double %1142
  %1145 = call double @llvm.fmuladd.f64(double %1144, double %1143, double %1141)
  %1146 = load ptr, ptr %12, align 8, !tbaa !5
  %1147 = getelementptr inbounds %struct.reionization, ptr %1146, i32 0, i32 10
  %1148 = load ptr, ptr %1147, align 8, !tbaa !46
  %1149 = load i32, ptr %18, align 4, !tbaa !17
  %1150 = sub nsw i32 %1149, 1
  %1151 = load ptr, ptr %12, align 8, !tbaa !5
  %1152 = getelementptr inbounds %struct.reionization, ptr %1151, i32 0, i32 8
  %1153 = load i32, ptr %1152, align 8, !tbaa !15
  %1154 = mul nsw i32 %1150, %1153
  %1155 = load ptr, ptr %12, align 8, !tbaa !5
  %1156 = getelementptr inbounds %struct.reionization, ptr %1155, i32 0, i32 2
  %1157 = load i32, ptr %1156, align 8, !tbaa !40
  %1158 = add nsw i32 %1154, %1157
  %1159 = sext i32 %1158 to i64
  %1160 = getelementptr inbounds double, ptr %1148, i64 %1159
  store double %1145, ptr %1160, align 8, !tbaa !14
  %1161 = load double, ptr %32, align 8, !tbaa !14
  %1162 = fmul double 0x4373F4D4EACDD756, %1161
  %1163 = fdiv double 0x3B30B0E7F1528C6D, %1162
  %1164 = load ptr, ptr %12, align 8, !tbaa !5
  %1165 = getelementptr inbounds %struct.reionization, ptr %1164, i32 0, i32 10
  %1166 = load ptr, ptr %1165, align 8, !tbaa !46
  %1167 = load i32, ptr %18, align 4, !tbaa !17
  %1168 = sub nsw i32 %1167, 1
  %1169 = load ptr, ptr %12, align 8, !tbaa !5
  %1170 = getelementptr inbounds %struct.reionization, ptr %1169, i32 0, i32 8
  %1171 = load i32, ptr %1170, align 8, !tbaa !15
  %1172 = mul nsw i32 %1168, %1171
  %1173 = load ptr, ptr %12, align 8, !tbaa !5
  %1174 = getelementptr inbounds %struct.reionization, ptr %1173, i32 0, i32 2
  %1175 = load i32, ptr %1174, align 8, !tbaa !40
  %1176 = add nsw i32 %1172, %1175
  %1177 = sext i32 %1176 to i64
  %1178 = getelementptr inbounds double, ptr %1166, i64 %1177
  %1179 = load double, ptr %1178, align 8, !tbaa !14
  %1180 = fmul double %1163, %1179
  %1181 = load ptr, ptr %12, align 8, !tbaa !5
  %1182 = getelementptr inbounds %struct.reionization, ptr %1181, i32 0, i32 10
  %1183 = load ptr, ptr %1182, align 8, !tbaa !46
  %1184 = load i32, ptr %18, align 4, !tbaa !17
  %1185 = sub nsw i32 %1184, 1
  %1186 = load ptr, ptr %12, align 8, !tbaa !5
  %1187 = getelementptr inbounds %struct.reionization, ptr %1186, i32 0, i32 8
  %1188 = load i32, ptr %1187, align 8, !tbaa !15
  %1189 = mul nsw i32 %1185, %1188
  %1190 = load ptr, ptr %12, align 8, !tbaa !5
  %1191 = getelementptr inbounds %struct.reionization, ptr %1190, i32 0, i32 3
  %1192 = load i32, ptr %1191, align 4, !tbaa !41
  %1193 = add nsw i32 %1189, %1192
  %1194 = sext i32 %1193 to i64
  %1195 = getelementptr inbounds double, ptr %1183, i64 %1194
  store double %1180, ptr %1195, align 8, !tbaa !14
  %1196 = load ptr, ptr %12, align 8, !tbaa !5
  %1197 = getelementptr inbounds %struct.reionization, ptr %1196, i32 0, i32 10
  %1198 = load ptr, ptr %1197, align 8, !tbaa !46
  %1199 = load i32, ptr %18, align 4, !tbaa !17
  %1200 = sub nsw i32 %1199, 1
  %1201 = load ptr, ptr %12, align 8, !tbaa !5
  %1202 = getelementptr inbounds %struct.reionization, ptr %1201, i32 0, i32 8
  %1203 = load i32, ptr %1202, align 8, !tbaa !15
  %1204 = mul nsw i32 %1200, %1203
  %1205 = load ptr, ptr %12, align 8, !tbaa !5
  %1206 = getelementptr inbounds %struct.reionization, ptr %1205, i32 0, i32 3
  %1207 = load i32, ptr %1206, align 4, !tbaa !41
  %1208 = add nsw i32 %1204, %1207
  %1209 = sext i32 %1208 to i64
  %1210 = getelementptr inbounds double, ptr %1198, i64 %1209
  %1211 = load double, ptr %1210, align 8, !tbaa !14
  %1212 = load double, ptr %22, align 8, !tbaa !14
  %1213 = fadd double 1.000000e+00, %1212
  %1214 = fdiv double %1213, 3.000000e+00
  %1215 = load double, ptr %30, align 8, !tbaa !14
  %1216 = fmul double %1214, %1215
  %1217 = load ptr, ptr %12, align 8, !tbaa !5
  %1218 = getelementptr inbounds %struct.reionization, ptr %1217, i32 0, i32 10
  %1219 = load ptr, ptr %1218, align 8, !tbaa !46
  %1220 = load i32, ptr %18, align 4, !tbaa !17
  %1221 = sub nsw i32 %1220, 1
  %1222 = load ptr, ptr %12, align 8, !tbaa !5
  %1223 = getelementptr inbounds %struct.reionization, ptr %1222, i32 0, i32 8
  %1224 = load i32, ptr %1223, align 8, !tbaa !15
  %1225 = mul nsw i32 %1221, %1224
  %1226 = load ptr, ptr %12, align 8, !tbaa !5
  %1227 = getelementptr inbounds %struct.reionization, ptr %1226, i32 0, i32 2
  %1228 = load i32, ptr %1227, align 8, !tbaa !40
  %1229 = add nsw i32 %1225, %1228
  %1230 = sext i32 %1229 to i64
  %1231 = getelementptr inbounds double, ptr %1219, i64 %1230
  %1232 = load double, ptr %1231, align 8, !tbaa !14
  %1233 = fdiv double %1216, %1232
  %1234 = fadd double 1.000000e+00, %1233
  %1235 = fmul double %1211, %1234
  %1236 = load ptr, ptr %12, align 8, !tbaa !5
  %1237 = getelementptr inbounds %struct.reionization, ptr %1236, i32 0, i32 10
  %1238 = load ptr, ptr %1237, align 8, !tbaa !46
  %1239 = load i32, ptr %18, align 4, !tbaa !17
  %1240 = sub nsw i32 %1239, 1
  %1241 = load ptr, ptr %12, align 8, !tbaa !5
  %1242 = getelementptr inbounds %struct.reionization, ptr %1241, i32 0, i32 8
  %1243 = load i32, ptr %1242, align 8, !tbaa !15
  %1244 = mul nsw i32 %1240, %1243
  %1245 = load ptr, ptr %12, align 8, !tbaa !5
  %1246 = getelementptr inbounds %struct.reionization, ptr %1245, i32 0, i32 4
  %1247 = load i32, ptr %1246, align 8, !tbaa !42
  %1248 = add nsw i32 %1244, %1247
  %1249 = sext i32 %1248 to i64
  %1250 = getelementptr inbounds double, ptr %1238, i64 %1249
  store double %1235, ptr %1250, align 8, !tbaa !14
  br label %1251

1251:                                             ; preds = %1126
  %1252 = load i32, ptr %18, align 4, !tbaa !17
  %1253 = add nsw i32 %1252, -1
  store i32 %1253, ptr %18, align 4, !tbaa !17
  br label %805, !llvm.loop !56

1254:                                             ; preds = %805
  %1255 = load ptr, ptr %12, align 8, !tbaa !5
  %1256 = getelementptr inbounds %struct.reionization, ptr %1255, i32 0, i32 10
  %1257 = load ptr, ptr %1256, align 8, !tbaa !46
  %1258 = load ptr, ptr %12, align 8, !tbaa !5
  %1259 = getelementptr inbounds %struct.reionization, ptr %1258, i32 0, i32 8
  %1260 = load i32, ptr %1259, align 8, !tbaa !15
  %1261 = load ptr, ptr %12, align 8, !tbaa !5
  %1262 = getelementptr inbounds %struct.reionization, ptr %1261, i32 0, i32 9
  %1263 = load i32, ptr %1262, align 4, !tbaa !47
  %1264 = load ptr, ptr %12, align 8, !tbaa !5
  %1265 = getelementptr inbounds %struct.reionization, ptr %1264, i32 0, i32 0
  %1266 = load i32, ptr %1265, align 8, !tbaa !29
  %1267 = load ptr, ptr %12, align 8, !tbaa !5
  %1268 = getelementptr inbounds %struct.reionization, ptr %1267, i32 0, i32 6
  %1269 = load i32, ptr %1268, align 8, !tbaa !38
  %1270 = load ptr, ptr %12, align 8, !tbaa !5
  %1271 = getelementptr inbounds %struct.reionization, ptr %1270, i32 0, i32 7
  %1272 = load i32, ptr %1271, align 4, !tbaa !57
  %1273 = load ptr, ptr %10, align 8, !tbaa !5
  %1274 = getelementptr inbounds %struct.thermo, ptr %1273, i32 0, i32 93
  %1275 = getelementptr inbounds [2048 x i8], ptr %1274, i64 0, i64 0
  %1276 = call i32 @array_spline(ptr noundef %1257, i32 noundef %1260, i32 noundef %1263, i32 noundef %1266, i32 noundef %1269, i32 noundef %1272, i16 noundef signext 1, ptr noundef %1275)
  %1277 = icmp eq i32 %1276, 1
  br i1 %1277, label %1278, label %1287

1278:                                             ; preds = %1254
  call void @llvm.lifetime.start.p0(i64 2048, ptr %73) #7
  %1279 = getelementptr inbounds [2048 x i8], ptr %73, i64 0, i64 0
  %1280 = load ptr, ptr %10, align 8, !tbaa !5
  %1281 = getelementptr inbounds %struct.thermo, ptr %1280, i32 0, i32 93
  %1282 = getelementptr inbounds [2048 x i8], ptr %1281, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1279, ptr noundef @.str, ptr noundef @.str.207, ptr noundef %1282)
  %1283 = load ptr, ptr %10, align 8, !tbaa !5
  %1284 = getelementptr inbounds %struct.thermo, ptr %1283, i32 0, i32 93
  %1285 = getelementptr inbounds [2048 x i8], ptr %1284, i64 0, i64 0
  %1286 = getelementptr inbounds [2048 x i8], ptr %73, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1285, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 3058, ptr noundef %1286)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %73) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

1287:                                             ; preds = %1254
  %1288 = load ptr, ptr %12, align 8, !tbaa !5
  %1289 = getelementptr inbounds %struct.reionization, ptr %1288, i32 0, i32 10
  %1290 = load ptr, ptr %1289, align 8, !tbaa !46
  %1291 = load ptr, ptr %12, align 8, !tbaa !5
  %1292 = getelementptr inbounds %struct.reionization, ptr %1291, i32 0, i32 8
  %1293 = load i32, ptr %1292, align 8, !tbaa !15
  %1294 = load ptr, ptr %12, align 8, !tbaa !5
  %1295 = getelementptr inbounds %struct.reionization, ptr %1294, i32 0, i32 9
  %1296 = load i32, ptr %1295, align 4, !tbaa !47
  %1297 = load ptr, ptr %12, align 8, !tbaa !5
  %1298 = getelementptr inbounds %struct.reionization, ptr %1297, i32 0, i32 0
  %1299 = load i32, ptr %1298, align 8, !tbaa !29
  %1300 = load ptr, ptr %12, align 8, !tbaa !5
  %1301 = getelementptr inbounds %struct.reionization, ptr %1300, i32 0, i32 6
  %1302 = load i32, ptr %1301, align 8, !tbaa !38
  %1303 = load ptr, ptr %12, align 8, !tbaa !5
  %1304 = getelementptr inbounds %struct.reionization, ptr %1303, i32 0, i32 7
  %1305 = load i32, ptr %1304, align 4, !tbaa !57
  %1306 = load ptr, ptr %12, align 8, !tbaa !5
  %1307 = getelementptr inbounds %struct.reionization, ptr %1306, i32 0, i32 11
  %1308 = load ptr, ptr %10, align 8, !tbaa !5
  %1309 = getelementptr inbounds %struct.thermo, ptr %1308, i32 0, i32 93
  %1310 = getelementptr inbounds [2048 x i8], ptr %1309, i64 0, i64 0
  %1311 = call i32 @array_integrate_all_spline(ptr noundef %1290, i32 noundef %1293, i32 noundef %1296, i32 noundef %1299, i32 noundef %1302, i32 noundef %1305, ptr noundef %1307, ptr noundef %1310)
  %1312 = icmp eq i32 %1311, 1
  br i1 %1312, label %1313, label %1322

1313:                                             ; preds = %1287
  call void @llvm.lifetime.start.p0(i64 2048, ptr %74) #7
  %1314 = getelementptr inbounds [2048 x i8], ptr %74, i64 0, i64 0
  %1315 = load ptr, ptr %10, align 8, !tbaa !5
  %1316 = getelementptr inbounds %struct.thermo, ptr %1315, i32 0, i32 93
  %1317 = getelementptr inbounds [2048 x i8], ptr %1316, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1314, ptr noundef @.str, ptr noundef @.str.208, ptr noundef %1317)
  %1318 = load ptr, ptr %10, align 8, !tbaa !5
  %1319 = getelementptr inbounds %struct.thermo, ptr %1318, i32 0, i32 93
  %1320 = getelementptr inbounds [2048 x i8], ptr %1319, i64 0, i64 0
  %1321 = getelementptr inbounds [2048 x i8], ptr %74, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1320, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_sample, i32 noundef 3070, ptr noundef %1321)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %74) #7
  store i32 1, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

1322:                                             ; preds = %1287
  store i32 0, ptr %7, align 4
  store i32 1, ptr %42, align 4
  br label %1323

1323:                                             ; preds = %1322, %1313, %1278, %1031, %849, %828, %791, %776, %720, %699, %635, %580, %543, %495, %477, %456, %442, %421, %370, %267, %230, %209, %188, %145, %102, %86
  call void @llvm.lifetime.end.p0(i64 8, ptr %39) #7
  call void @llvm.lifetime.end.p0(i64 4, ptr %38) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %37) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %36) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %35) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %34) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %33) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %32) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %31) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %30) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %29) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %28) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %27) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %26) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %25) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %22) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %20) #7
  call void @llvm.lifetime.end.p0(i64 4, ptr %19) #7
  call void @llvm.lifetime.end.p0(i64 4, ptr %18) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %17) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #7
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #7
  call void @llvm.lifetime.end.p0(i64 2080, ptr %14) #7
  %1324 = load i32, ptr %7, align 4
  ret i32 %1324
}

declare i32 @gt_init(ptr noundef) #3

declare i32 @gt_add(ptr noundef, i64 noundef, ptr noundef, i64 noundef) #3

declare i32 @gt_getPtr(ptr noundef, ptr noundef) #3

; Function Attrs: nocallback nofree nounwind willreturn memory(argmem: readwrite)
declare void @llvm.memcpy.p0.p0.i64(ptr noalias nocapture writeonly, ptr noalias nocapture readonly, i64, i1 immarg) #6

declare i32 @gt_free(ptr noundef) #3

declare i32 @array_spline(ptr noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, i16 noundef signext, ptr noundef) #3

declare i32 @array_integrate_all_spline(ptr noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, ptr noundef, ptr noundef) #3

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { 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 #2 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #3 = { "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 #4 = { nounwind allocsize(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 #5 = { 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 #6 = { nocallback nofree nounwind willreturn memory(argmem: readwrite) }
attributes #7 = { nounwind }
attributes #8 = { nounwind allocsize(0) }

!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, !11, i64 0}
!10 = !{!"thermo", !11, i64 0, !12, i64 8, !12, i64 12, !12, i64 16, !11, i64 24, !11, i64 32, !13, i64 40, !13, i64 42, !11, i64 48, !11, i64 56, !11, i64 64, !11, i64 72, !12, i64 80, !6, i64 88, !6, i64 96, !11, i64 104, !12, i64 112, !6, i64 120, !6, i64 128, !11, i64 136, !12, i64 144, !6, i64 152, !6, i64 160, !11, i64 168, !13, i64 176, !11, i64 184, !11, i64 192, !11, i64 200, !11, i64 208, !11, i64 216, !11, i64 224, !11, i64 232, !11, i64 240, !11, i64 248, !11, i64 256, !11, i64 264, !12, i64 272, !12, i64 276, !12, i64 280, !12, i64 284, !12, i64 288, !12, i64 292, !12, i64 296, !12, i64 300, !12, i64 304, !12, i64 308, !12, i64 312, !12, i64 316, !12, i64 320, !12, i64 324, !12, i64 328, !12, i64 332, !12, i64 336, !12, i64 340, !12, i64 344, !12, i64 348, !12, i64 352, !12, i64 356, !12, i64 360, !12, i64 364, !12, i64 368, !12, i64 372, !12, i64 376, !6, i64 384, !6, i64 392, !6, i64 400, !11, i64 408, !11, i64 416, !11, i64 424, !11, i64 432, !11, i64 440, !11, i64 448, !11, i64 456, !11, i64 464, !11, i64 472, !11, i64 480, !11, i64 488, !11, i64 496, !11, i64 504, !11, i64 512, !11, i64 520, !11, i64 528, !11, i64 536, !11, i64 544, !11, i64 552, !11, i64 560, !11, i64 568, !11, i64 576, !11, i64 584, !11, i64 592, !13, i64 600, !13, i64 602, !13, i64 604, !7, i64 606}
!11 = !{!"double", !7, i64 0}
!12 = !{!"int", !7, i64 0}
!13 = !{!"short", !7, i64 0}
!14 = !{!11, !11, i64 0}
!15 = !{!16, !12, i64 32}
!16 = !{!"reionization", !12, i64 0, !12, i64 4, !12, i64 8, !12, i64 12, !12, i64 16, !12, i64 20, !12, i64 24, !12, i64 28, !12, i64 32, !12, i64 36, !6, i64 40, !11, i64 48, !12, i64 56, !12, i64 60, !12, i64 64, !12, i64 68, !12, i64 72, !12, i64 76, !12, i64 80, !12, i64 84, !12, i64 88, !12, i64 92, !12, i64 96, !12, i64 100, !12, i64 104, !6, i64 112, !12, i64 120, !12, i64 124}
!17 = !{!12, !12, i64 0}
!18 = !{!19, !6, i64 32}
!19 = !{!"recombination", !12, i64 0, !12, i64 4, !12, i64 8, !12, i64 12, !12, i64 16, !12, i64 20, !12, i64 24, !12, i64 28, !6, 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, !11, i64 112, !11, i64 120, !11, i64 128, !11, i64 136, !11, i64 144, !11, i64 152, !11, i64 160, !11, i64 168, !11, i64 176, !11, i64 184, !11, i64 192, !13, i64 200, !11, i64 208, !11, i64 216, !11, i64 224, !11, i64 232, !11, i64 240, !11, i64 248, !11, i64 256}
!20 = !{!19, !12, i64 24}
!21 = !{!19, !12, i64 0}
!22 = !{!16, !6, i64 112}
!23 = !{!16, !12, i64 104}
!24 = !{!25, !12, i64 368}
!25 = !{!"precision", !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, !7, i64 104, !11, i64 360, !12, i64 368, !11, i64 376, !12, i64 384, !12, i64 388, !11, i64 392, !12, i64 400, !11, i64 408, !12, i64 416, !11, i64 424, !11, i64 432, !11, i64 440, !11, i64 448, !11, i64 456, !11, i64 464, !11, i64 472, !11, i64 480, !11, i64 488, !11, i64 496, !11, i64 504, !11, i64 512, !11, i64 520, !11, i64 528, !11, i64 536, !11, i64 544, !11, i64 552, !11, i64 560, !11, i64 568, !11, i64 576, !11, i64 584, !11, i64 592, !11, i64 600, !11, i64 608, !11, i64 616, !12, i64 624, !7, i64 628, !7, i64 884, !7, i64 1140, !11, i64 1400, !11, i64 1408, !11, i64 1416, !11, i64 1424, !11, i64 1432, !11, i64 1440, !11, i64 1448, !11, i64 1456, !11, i64 1464, !11, i64 1472, !11, i64 1480, !11, i64 1488, !11, i64 1496, !11, i64 1504, !11, i64 1512, !11, i64 1520, !12, i64 1528, !11, i64 1536, !11, i64 1544, !12, i64 1552, !12, i64 1556, !12, i64 1560, !12, i64 1564, !12, i64 1568, !12, i64 1572, !12, i64 1576, !12, i64 1580, !11, i64 1584, !11, i64 1592, !11, i64 1600, !11, i64 1608, !11, i64 1616, !11, i64 1624, !11, i64 1632, !11, i64 1640, !12, i64 1648, !11, i64 1656, !11, i64 1664, !12, i64 1672, !11, i64 1680, !11, i64 1688, !12, i64 1696, !11, i64 1704, !12, i64 1712, !11, i64 1720, !11, i64 1728, !12, i64 1736, !11, i64 1744, !11, i64 1752, !11, i64 1760, !12, i64 1768, !11, i64 1776, !11, i64 1784, !11, i64 1792, !11, i64 1800, !11, i64 1808, !12, i64 1816, !11, i64 1824, !11, i64 1832, !11, i64 1840, !11, i64 1848, !11, i64 1856, !11, i64 1864, !11, i64 1872, !12, i64 1880, !11, i64 1888, !11, i64 1896, !11, i64 1904, !11, i64 1912, !11, i64 1920, !11, i64 1928, !11, i64 1936, !11, i64 1944, !11, i64 1952, !11, i64 1960, !11, i64 1968, !11, i64 1976, !11, i64 1984, !11, i64 1992, !11, i64 2000, !11, i64 2008, !11, i64 2016, !11, i64 2024, !11, i64 2032, !11, i64 2040, !11, i64 2048, !11, i64 2056, !11, i64 2064, !11, i64 2072, !11, i64 2080, !11, i64 2088, !11, i64 2096, !11, i64 2104, !11, i64 2112, !11, i64 2120, !11, i64 2128, !11, i64 2136, !11, i64 2144, !11, i64 2152, !11, i64 2160, !11, i64 2168, !11, i64 2176, !11, i64 2184, !11, i64 2192, !11, i64 2200, !11, i64 2208, !11, i64 2216, !11, i64 2224, !11, i64 2232, !11, i64 2240, !11, i64 2248, !12, i64 2256, !11, i64 2264, !11, i64 2272, !11, i64 2280, !11, i64 2288, !12, i64 2296, !11, i64 2304, !11, i64 2312, !12, i64 2320, !12, i64 2324, !12, i64 2328, !11, i64 2336, !11, i64 2344, !7, i64 2352}
!26 = !{!25, !11, i64 592}
!27 = distinct !{!27, !28}
!28 = !{!"llvm.loop.mustprogress"}
!29 = !{!16, !12, i64 0}
!30 = !{!16, !12, i64 124}
!31 = !{!16, !12, i64 4}
!32 = !{!33, !13, i64 840}
!33 = !{!"background", !11, i64 0, !11, i64 8, !11, i64 16, !11, i64 24, !11, i64 32, !11, i64 40, !11, i64 48, !12, i64 56, !11, i64 64, !11, i64 72, !11, i64 80, !11, i64 88, !13, i64 96, !11, i64 104, !11, i64 112, !11, i64 120, !11, i64 128, !11, i64 136, !11, i64 144, !11, i64 152, !11, i64 160, !11, i64 168, !13, i64 176, !11, i64 184, !11, i64 192, !6, i64 200, !12, i64 208, !12, i64 212, !11, i64 216, !12, i64 224, !6, i64 232, !6, i64 240, !11, i64 248, !6, i64 256, !11, i64 264, !6, i64 272, !11, i64 280, !6, i64 288, !11, i64 296, !6, i64 304, !6, i64 312, !6, i64 320, !11, i64 328, !11, i64 336, !11, i64 344, !11, i64 352, !12, i64 360, !6, i64 368, !11, i64 376, !11, i64 384, !11, i64 392, !11, i64 400, !11, i64 408, !11, i64 416, !11, i64 424, !11, i64 432, !11, i64 440, !11, i64 448, !11, i64 456, !12, i64 464, !12, i64 468, !12, i64 472, !12, i64 476, !12, i64 480, !12, i64 484, !12, i64 488, !12, i64 492, !12, i64 496, !12, i64 500, !12, i64 504, !12, i64 508, !12, i64 512, !12, i64 516, !12, i64 520, !12, i64 524, !12, i64 528, !12, i64 532, !12, i64 536, !12, i64 540, !12, i64 544, !12, i64 548, !12, i64 552, !12, i64 556, !12, i64 560, !12, i64 564, !12, i64 568, !12, i64 572, !12, i64 576, !12, i64 580, !12, i64 584, !12, i64 588, !12, i64 592, !12, i64 596, !12, i64 600, !12, i64 604, !12, i64 608, !12, i64 612, !12, i64 616, !12, i64 620, !12, i64 624, !12, i64 628, !6, i64 632, !6, i64 640, !6, i64 648, !6, i64 656, !6, i64 664, !12, i64 672, !12, i64 676, !12, i64 680, !12, i64 684, !12, i64 688, !12, i64 692, !12, i64 696, !12, i64 700, !12, i64 704, !12, i64 708, !12, i64 712, !12, i64 716, !12, i64 720, !13, i64 724, !13, i64 726, !13, i64 728, !13, i64 730, !13, i64 732, !13, i64 734, !13, i64 736, !13, i64 738, !13, i64 740, !13, i64 742, !13, i64 744, !6, i64 752, !6, i64 760, !6, i64 768, !6, i64 776, !6, i64 784, !6, i64 792, !6, i64 800, !6, i64 808, !6, i64 816, !6, i64 824, !6, i64 832, !13, i64 840, !13, i64 842, !13, i64 844, !13, i64 846, !13, i64 848, !13, i64 850, !7, i64 852, !13, i64 2900, !7, i64 2902}
!34 = !{!33, !13, i64 846}
!35 = !{!10, !11, i64 592}
!36 = !{!16, !12, i64 20}
!37 = !{!33, !12, i64 468}
!38 = !{!16, !12, i64 24}
!39 = !{!19, !12, i64 8}
!40 = !{!16, !12, i64 8}
!41 = !{!16, !12, i64 12}
!42 = !{!16, !12, i64 16}
!43 = !{!25, !11, i64 2344}
!44 = !{!25, !11, i64 600}
!45 = distinct !{!45, !28}
!46 = !{!16, !6, i64 40}
!47 = !{!16, !12, i64 36}
!48 = distinct !{!48, !28}
!49 = !{!33, !13, i64 842}
!50 = !{!33, !12, i64 476}
!51 = !{!33, !12, i64 480}
!52 = !{!33, !11, i64 16}
!53 = !{!19, !11, i64 192}
!54 = !{!19, !11, i64 136}
!55 = !{!19, !11, i64 80}
!56 = distinct !{!56, !28}
!57 = !{!16, !12, i64 28}


; 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"

%struct.growTable = type { ptr, i64, i64, i32, [2048 x i8] }

@.str = external hidden unnamed_addr constant [18 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [13 x i8], align 1
@.str.10 = external hidden unnamed_addr constant [27 x i8], align 1
@.str.34 = external hidden unnamed_addr constant [31 x i8], align 1
@.str.37 = external hidden unnamed_addr constant [35 x i8], align 1
@.str.185 = external hidden unnamed_addr constant [12 x i8], align 1
@__func__.thermodynamics_reionization_sample = external hidden unnamed_addr constant [35 x i8], align 1
@.str.186 = external hidden unnamed_addr constant [17 x i8], align 1
@.str.187 = external hidden unnamed_addr constant [73 x i8], align 1
@.str.188 = external hidden unnamed_addr constant [22 x i8], align 1
@.str.189 = external hidden unnamed_addr constant [54 x i8], align 1
@.str.190 = external hidden unnamed_addr constant [34 x i8], align 1
@.str.191 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.192 = external hidden unnamed_addr constant [32 x i8], align 1
@.str.193 = external hidden unnamed_addr constant [72 x i8], align 1
@.str.194 = external hidden unnamed_addr constant [113 x i8], align 1
@.str.195 = external hidden unnamed_addr constant [37 x i8], align 1
@.str.196 = external hidden unnamed_addr constant [64 x i8], align 1
@.str.197 = external hidden unnamed_addr constant [39 x i8], align 1
@.str.198 = external hidden unnamed_addr constant [39 x i8], align 1
@.str.199 = external hidden unnamed_addr constant [59 x i8], align 1
@.str.200 = external hidden unnamed_addr constant [26 x i8], align 1
@.str.201 = external hidden unnamed_addr constant [34 x i8], align 1
@.str.202 = external hidden unnamed_addr constant [54 x i8], align 1
@.str.203 = external hidden unnamed_addr constant [61 x i8], align 1
@.str.204 = external hidden unnamed_addr constant [17 x i8], align 1
@.str.205 = external hidden unnamed_addr constant [93 x i8], align 1
@.str.206 = external hidden unnamed_addr constant [81 x i8], align 1
@.str.207 = external hidden unnamed_addr constant [168 x i8], align 1
@.str.208 = external hidden unnamed_addr constant [217 x i8], align 1

; 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 nofree nounwind willreturn memory(write)
declare double @pow(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) #2

declare void @class_protect_sprintf(ptr noundef, ptr noundef, ...) local_unnamed_addr #3

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

; Function Attrs: mustprogress nofree nounwind willreturn allockind("alloc,uninitialized") allocsize(0) memory(inaccessiblemem: readwrite)
declare noalias noundef ptr @malloc(i64 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

declare i32 @background_tau_of_z(ptr noundef, double noundef, ptr noundef) local_unnamed_addr #3

declare i32 @background_at_tau(ptr noundef, double noundef, i16 noundef signext, i16 noundef signext, ptr noundef, ptr noundef) local_unnamed_addr #3

; Function Attrs: nounwind uwtable
declare dso_local i32 @thermodynamics_energy_injection(ptr noundef, ptr noundef, ptr noundef, double noundef, ptr noundef, ptr noundef) local_unnamed_addr #6

; Function Attrs: nounwind uwtable
declare dso_local i32 @thermodynamics_reionization_function(double noundef, ptr noundef, ptr noundef, ptr noundef) local_unnamed_addr #6

; Function Attrs: nounwind uwtable
define dso_local range(i32 0, 2) i32 @thermodynamics_reionization_sample(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, ptr noundef %4, ptr noundef %5) local_unnamed_addr #6 {
  %7 = alloca %struct.growTable, align 8
  %8 = alloca ptr, align 8
  %9 = alloca double, align 8
  %10 = alloca double, align 8
  %11 = alloca double, align 8
  %12 = alloca double, align 8
  %13 = alloca i32, align 4
  %14 = alloca [2048 x i8], align 16
  %15 = alloca [2048 x i8], align 16
  %16 = alloca [2048 x i8], align 16
  %17 = alloca [2048 x i8], align 16
  %18 = alloca [2048 x i8], align 16
  %19 = alloca [2048 x i8], align 16
  %20 = alloca [2048 x i8], align 16
  %21 = alloca [2048 x i8], align 16
  %22 = alloca [2048 x i8], align 16
  %23 = alloca [2048 x i8], align 16
  %24 = alloca [2048 x i8], align 16
  %25 = alloca [2048 x i8], align 16
  %26 = alloca [2048 x i8], align 16
  %27 = alloca [2048 x i8], align 16
  %28 = alloca [2048 x i8], align 16
  %29 = alloca [2048 x i8], align 16
  %30 = alloca [2048 x i8], align 16
  %31 = alloca [2048 x i8], align 16
  %32 = alloca [2048 x i8], align 16
  %33 = alloca [2048 x i8], align 16
  %34 = alloca [2048 x i8], align 16
  %35 = alloca [2048 x i8], align 16
  %36 = alloca [2048 x i8], align 16
  %37 = alloca [2048 x i8], align 16
  %38 = alloca [2048 x i8], align 16
  %39 = alloca [2048 x i8], align 16
  %40 = alloca [2048 x i8], align 16
  %41 = alloca [2048 x i8], align 16
  %42 = alloca [2048 x i8], align 16
  %43 = alloca [2048 x i8], align 16
  %44 = alloca [2048 x i8], align 16
  %45 = alloca [2048 x i8], align 16
  call void @llvm.lifetime.start.p0(i64 2080, ptr nonnull %7) #9
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %8) #9
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %9) #9
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %10) #9
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %11) #9
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %12) #9
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %13) #9
  %46 = load double, ptr %2, align 8, !tbaa !5
  %47 = getelementptr inbounds i8, ptr %4, i64 32
  %48 = load i32, ptr %47, align 8, !tbaa !13
  %49 = sext i32 %48 to i64
  %50 = shl nsw i64 %49, 3
  %51 = tail call noalias ptr @malloc(i64 noundef %50) #10
  %52 = icmp eq ptr %51, null
  br i1 %52, label %53, label %56

53:                                               ; preds = %6
  %54 = shl i32 %48, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %14) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %14, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.185, i32 noundef %54) #9
  %55 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %55, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2784, ptr noundef nonnull %14) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %14) #9
  br label %567

56:                                               ; preds = %6
  %57 = call i32 @gt_init(ptr noundef nonnull %7) #9
  %58 = icmp eq i32 %57, 1
  br i1 %58, label %74, label %.preheader363

.preheader363:                                    ; preds = %56
  %59 = getelementptr inbounds i8, ptr %3, i64 32
  %60 = load ptr, ptr %59, align 8, !tbaa !15
  %61 = getelementptr inbounds i8, ptr %3, i64 24
  %62 = load i32, ptr %61, align 8, !tbaa !17
  %63 = load i32, ptr %3, align 8, !tbaa !18
  %64 = getelementptr inbounds i8, ptr %4, i64 112
  %65 = load ptr, ptr %64, align 8, !tbaa !19
  %66 = getelementptr inbounds i8, ptr %4, i64 104
  %67 = load i32, ptr %66, align 8, !tbaa !20
  %68 = sext i32 %67 to i64
  %69 = getelementptr inbounds double, ptr %65, i64 %68
  %70 = load double, ptr %69, align 8, !tbaa !21
  %71 = getelementptr inbounds i8, ptr %0, i64 368
  %72 = sext i32 %62 to i64
  %73 = sext i32 %63 to i64
  %invariant.gep = getelementptr double, ptr %60, i64 %73
  br label %77

74:                                               ; preds = %56
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %15) #9
  %75 = getelementptr inbounds i8, ptr %7, i64 28
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %15, ptr noundef nonnull @.str, ptr noundef nonnull @.str.186, ptr noundef nonnull %75) #9
  %76 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %76, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2789, ptr noundef nonnull %15) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %15) #9
  br label %567

77:                                               ; preds = %.preheader363, %81
  %indvars.iv = phi i64 [ 0, %.preheader363 ], [ %indvars.iv.next, %81 ]
  %78 = mul nsw i64 %indvars.iv, %72
  %gep = getelementptr double, ptr %invariant.gep, i64 %78
  %79 = load double, ptr %gep, align 8, !tbaa !21
  %80 = fcmp olt double %79, %70
  br i1 %80, label %81, label %89

81:                                               ; preds = %77
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %82 = load i32, ptr %71, align 8, !tbaa !22
  %83 = trunc nsw i64 %indvars.iv.next to i32
  %84 = icmp eq i32 %82, %83
  br i1 %84, label %85, label %77, !llvm.loop !24

85:                                               ; preds = %81
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %16) #9
  %86 = getelementptr inbounds i8, ptr %0, i64 592
  %87 = load double, ptr %86, align 8, !tbaa !26
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %16, ptr noundef nonnull @.str.187, double noundef %87) #9
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %17) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %17, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.188, ptr noundef nonnull %16) #9
  %88 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %88, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2799, ptr noundef nonnull %17) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %17) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %16) #9
  br label %567

89:                                               ; preds = %77
  %90 = trunc nuw nsw i64 %indvars.iv to i32
  %91 = load i32, ptr %4, align 8, !tbaa !27
  %92 = sext i32 %91 to i64
  %93 = getelementptr inbounds double, ptr %51, i64 %92
  store double %79, ptr %93, align 8, !tbaa !21
  %94 = getelementptr inbounds i8, ptr %4, i64 124
  store i32 %90, ptr %94, align 4, !tbaa !28
  %95 = call i32 @thermodynamics_reionization_function(double noundef %79, ptr noundef nonnull %2, ptr noundef nonnull %4, ptr noundef nonnull %9)
  %96 = icmp eq i32 %95, 1
  br i1 %96, label %97, label %99

97:                                               ; preds = %89
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %18) #9
  %98 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %18, ptr noundef nonnull @.str, ptr noundef nonnull @.str.189, ptr noundef nonnull %98) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %98, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2810, ptr noundef nonnull %18) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %18) #9
  br label %567

99:                                               ; preds = %89
  %100 = load double, ptr %9, align 8, !tbaa !21
  %101 = getelementptr inbounds i8, ptr %4, i64 4
  %102 = load i32, ptr %101, align 4, !tbaa !29
  %103 = sext i32 %102 to i64
  %104 = getelementptr inbounds double, ptr %51, i64 %103
  store double %100, ptr %104, align 8, !tbaa !21
  %105 = call i32 @background_tau_of_z(ptr noundef %1, double noundef %79, ptr noundef nonnull %12) #9
  %106 = icmp eq i32 %105, 1
  br i1 %106, label %107, label %110

107:                                              ; preds = %99
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %19) #9
  %108 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %19, ptr noundef nonnull @.str, ptr noundef nonnull @.str.190, ptr noundef nonnull %108) #9
  %109 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %109, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2820, ptr noundef nonnull %19) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %19) #9
  br label %567

110:                                              ; preds = %99
  %111 = load double, ptr %12, align 8, !tbaa !21
  %112 = getelementptr inbounds i8, ptr %1, i64 840
  %113 = load i16, ptr %112, align 8, !tbaa !30
  %114 = getelementptr inbounds i8, ptr %1, i64 846
  %115 = load i16, ptr %114, align 2, !tbaa !32
  %116 = call i32 @background_at_tau(ptr noundef %1, double noundef %111, i16 noundef signext %113, i16 noundef signext %115, ptr noundef nonnull %13, ptr noundef %5) #9
  %117 = icmp eq i32 %116, 1
  br i1 %117, label %118, label %121

118:                                              ; preds = %110
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %20) #9
  %119 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %20, ptr noundef nonnull @.str, ptr noundef nonnull @.str.191, ptr noundef nonnull %119) #9
  %120 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %120, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2829, ptr noundef nonnull %20) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %20) #9
  br label %567

121:                                              ; preds = %110
  %122 = fadd double %79, 1.000000e+00
  %123 = fmul double %122, %122
  %124 = getelementptr inbounds i8, ptr %2, i64 592
  %125 = load double, ptr %124, align 8, !tbaa !33
  %126 = fmul double %123, %125
  %127 = load double, ptr %9, align 8, !tbaa !21
  %128 = fmul double %126, %127
  %129 = fmul double %128, 0x3A15151E37FBF236
  %130 = fmul double %129, 0x449A22FF260E3826
  %131 = getelementptr inbounds i8, ptr %4, i64 20
  %132 = load i32, ptr %131, align 4, !tbaa !34
  %133 = sext i32 %132 to i64
  %134 = getelementptr inbounds double, ptr %51, i64 %133
  store double %130, ptr %134, align 8, !tbaa !21
  %135 = getelementptr inbounds i8, ptr %1, i64 468
  %136 = load i32, ptr %135, align 4, !tbaa !35
  %137 = sext i32 %136 to i64
  %138 = getelementptr inbounds double, ptr %5, i64 %137
  %139 = load double, ptr %138, align 8, !tbaa !21
  %140 = fcmp oeq double %139, 0.000000e+00
  br i1 %140, label %141, label %143

141:                                              ; preds = %121
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %21) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %21, ptr noundef nonnull @.str.34) #9
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %22) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %22, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.192, ptr noundef nonnull %21) #9
  %142 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %142, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2835, ptr noundef nonnull %22) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %22) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %21) #9
  br label %567

143:                                              ; preds = %121
  %144 = fdiv double %130, %139
  %145 = getelementptr inbounds i8, ptr %4, i64 24
  %146 = load i32, ptr %145, align 8, !tbaa !36
  %147 = sext i32 %146 to i64
  %148 = getelementptr inbounds double, ptr %51, i64 %147
  store double %144, ptr %148, align 8, !tbaa !21
  %149 = load double, ptr %134, align 8, !tbaa !21
  %150 = load ptr, ptr %59, align 8, !tbaa !15
  %151 = load i32, ptr %61, align 8, !tbaa !17
  %152 = mul nsw i32 %151, %90
  %153 = getelementptr inbounds i8, ptr %3, i64 8
  %154 = load i32, ptr %153, align 8, !tbaa !37
  %155 = add nsw i32 %152, %154
  %156 = sext i32 %155 to i64
  %157 = getelementptr inbounds double, ptr %150, i64 %156
  %158 = load double, ptr %157, align 8, !tbaa !21
  %159 = getelementptr inbounds i8, ptr %4, i64 8
  %160 = load i32, ptr %159, align 8, !tbaa !38
  %161 = sext i32 %160 to i64
  %162 = getelementptr inbounds double, ptr %51, i64 %161
  store double %158, ptr %162, align 8, !tbaa !21
  %163 = call double @llvm.fmuladd.f64(double %46, double 0xBFE7F14DA42541E8, double 1.000000e+00)
  %164 = fsub double 1.000000e+00, %46
  %165 = call double @llvm.fmuladd.f64(double %127, double %164, double %163)
  %166 = fmul double %165, 0x3D39D63187BB40C4
  %167 = fmul double %166, %158
  %168 = getelementptr inbounds i8, ptr %4, i64 12
  %169 = load i32, ptr %168, align 4, !tbaa !39
  %170 = sext i32 %169 to i64
  %171 = getelementptr inbounds double, ptr %51, i64 %170
  store double %167, ptr %171, align 8, !tbaa !21
  %172 = fmul double %167, 0x3FFAAAAAAAAAAAAB
  %173 = getelementptr inbounds i8, ptr %4, i64 16
  %174 = load i32, ptr %173, align 8, !tbaa !40
  %175 = sext i32 %174 to i64
  %176 = getelementptr inbounds double, ptr %51, i64 %175
  store double %172, ptr %176, align 8, !tbaa !21
  %177 = load i32, ptr %47, align 8, !tbaa !13
  %178 = sext i32 %177 to i64
  %179 = shl nsw i64 %178, 3
  %180 = call i32 @gt_add(ptr noundef nonnull %7, i64 noundef -1, ptr noundef nonnull %51, i64 noundef %179) #9
  %181 = icmp eq i32 %180, 1
  br i1 %181, label %182, label %185

182:                                              ; preds = %143
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %23) #9
  %183 = getelementptr inbounds i8, ptr %7, i64 28
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %23, ptr noundef nonnull @.str, ptr noundef nonnull @.str.193, ptr noundef nonnull %183) #9
  %184 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %184, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2858, ptr noundef nonnull %23) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %23) #9
  br label %567

185:                                              ; preds = %143
  %186 = load ptr, ptr %59, align 8, !tbaa !15
  %187 = load i32, ptr %61, align 8, !tbaa !17
  %188 = mul nsw i32 %187, %90
  %189 = load i32, ptr %3, align 8, !tbaa !18
  %190 = add nsw i32 %188, %189
  %191 = sext i32 %190 to i64
  %192 = getelementptr inbounds double, ptr %186, i64 %191
  %193 = load double, ptr %192, align 8, !tbaa !21
  %194 = add nsw i32 %90, -1
  %195 = mul nsw i32 %187, %194
  %196 = add nsw i32 %195, %189
  %197 = sext i32 %196 to i64
  %198 = getelementptr inbounds double, ptr %186, i64 %197
  %199 = load double, ptr %198, align 8, !tbaa !21
  %200 = fsub double %193, %199
  %201 = fcmp ogt double %79, 0.000000e+00
  br i1 %201, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %185
  %202 = getelementptr inbounds i8, ptr %0, i64 2344
  %203 = getelementptr inbounds i8, ptr %0, i64 600
  %204 = insertelement <2 x double> poison, double %149, i64 0
  %205 = insertelement <2 x double> %204, double %144, i64 1
  br label %206

206:                                              ; preds = %.lr.ph, %309
  %.0336389 = phi i32 [ 1, %.lr.ph ], [ %.1337, %309 ]
  %.0339388 = phi double [ %200, %.lr.ph ], [ %.1340, %309 ]
  %.0343386 = phi double [ %79, %.lr.ph ], [ %.1344, %309 ]
  %207 = phi <2 x double> [ %205, %.lr.ph ], [ %310, %309 ]
  %208 = load double, ptr %202, align 8, !tbaa !41
  %209 = fcmp olt double %.0339388, %208
  br i1 %209, label %210, label %212

210:                                              ; preds = %206
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %24) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %24, ptr noundef nonnull @.str.194) #9
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %25) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %25, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.195, ptr noundef nonnull %24) #9
  %211 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %211, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2875, ptr noundef nonnull %25) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %25) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %24) #9
  br label %567

212:                                              ; preds = %206
  %213 = fsub double %.0343386, %.0339388
  %214 = fcmp olt double %213, 0.000000e+00
  %.0345 = select i1 %214, double 0.000000e+00, double %213
  %215 = call i32 @thermodynamics_reionization_function(double noundef %.0345, ptr noundef nonnull %2, ptr noundef nonnull %4, ptr noundef nonnull %10)
  %216 = icmp eq i32 %215, 1
  br i1 %216, label %217, label %219

217:                                              ; preds = %212
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %26) #9
  %218 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %26, ptr noundef nonnull @.str, ptr noundef nonnull @.str.196, ptr noundef nonnull %218) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %218, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2883, ptr noundef nonnull %26) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %26) #9
  br label %567

219:                                              ; preds = %212
  %220 = call i32 @background_tau_of_z(ptr noundef nonnull %1, double noundef %.0345, ptr noundef nonnull %12) #9
  %221 = icmp eq i32 %220, 1
  br i1 %221, label %222, label %225

222:                                              ; preds = %219
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %27) #9
  %223 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %27, ptr noundef nonnull @.str, ptr noundef nonnull @.str.197, ptr noundef nonnull %223) #9
  %224 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %224, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2889, ptr noundef nonnull %27) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %27) #9
  br label %567

225:                                              ; preds = %219
  %226 = load double, ptr %12, align 8, !tbaa !21
  %227 = load i16, ptr %112, align 8, !tbaa !30
  %228 = load i16, ptr %114, align 2, !tbaa !32
  %229 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %226, i16 noundef signext %227, i16 noundef signext %228, ptr noundef nonnull %13, ptr noundef nonnull %5) #9
  %230 = icmp eq i32 %229, 1
  br i1 %230, label %231, label %234

231:                                              ; preds = %225
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %28) #9
  %232 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %28, ptr noundef nonnull @.str, ptr noundef nonnull @.str.191, ptr noundef nonnull %232) #9
  %233 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %233, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2898, ptr noundef nonnull %28) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %28) #9
  br label %567

234:                                              ; preds = %225
  %235 = load i32, ptr %135, align 4, !tbaa !35
  %236 = sext i32 %235 to i64
  %237 = getelementptr inbounds double, ptr %5, i64 %236
  %238 = load double, ptr %237, align 8, !tbaa !21
  %239 = fcmp oeq double %238, 0.000000e+00
  br i1 %239, label %240, label %242

240:                                              ; preds = %234
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %29) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %29, ptr noundef nonnull @.str.34) #9
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %30) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %30, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.192, ptr noundef nonnull %29) #9
  %241 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %241, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2902, ptr noundef nonnull %30) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %30) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %29) #9
  br label %567

242:                                              ; preds = %234
  %243 = fadd double %.0345, 1.000000e+00
  %244 = fmul double %243, %243
  %245 = load double, ptr %124, align 8, !tbaa !33
  %246 = fmul double %244, %245
  %247 = load double, ptr %10, align 8, !tbaa !21
  %248 = fmul double %246, %247
  %249 = fmul double %248, 0x3A15151E37FBF236
  %250 = fmul double %249, 0x449A22FF260E3826
  %251 = fdiv double %250, %238
  %252 = extractelement <2 x double> %207, i64 1
  %253 = fcmp oeq double %252, 0.000000e+00
  %254 = extractelement <2 x double> %207, i64 0
  %255 = fcmp oeq double %254, 0.000000e+00
  %or.cond = select i1 %253, i1 true, i1 %255
  br i1 %or.cond, label %256, label %258

256:                                              ; preds = %242
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %31) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %31, ptr noundef nonnull @.str.34) #9
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %32) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %32, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.198, ptr noundef nonnull %31) #9
  %257 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %257, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2910, ptr noundef nonnull %32) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %32) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %31) #9
  br label %567

258:                                              ; preds = %242
  %259 = insertelement <2 x double> poison, double %250, i64 0
  %260 = insertelement <2 x double> %259, double %251, i64 1
  %261 = fsub <2 x double> %260, %207
  %262 = fdiv <2 x double> %261, %207
  %263 = call <2 x double> @llvm.fabs.v2f64(<2 x double> %262)
  %shift = shufflevector <2 x double> %263, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %264 = fadd <2 x double> %263, %shift
  %265 = extractelement <2 x double> %264, i64 0
  %266 = load double, ptr %203, align 8, !tbaa !42
  %267 = fcmp olt double %265, %266
  br i1 %267, label %268, label %305

268:                                              ; preds = %258
  store double %247, ptr %9, align 8, !tbaa !21
  %269 = fcmp oeq double %251, 0.000000e+00
  %270 = fcmp oeq double %250, 0.000000e+00
  %or.cond3 = or i1 %270, %269
  br i1 %or.cond3, label %271, label %273

271:                                              ; preds = %268
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %33) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %33, ptr noundef nonnull @.str.199, double noundef %251, double noundef %250) #9
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %34) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %34, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.198, ptr noundef nonnull %33) #9
  %272 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %272, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2925, ptr noundef nonnull %34) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %34) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %33) #9
  br label %567

273:                                              ; preds = %268
  %274 = load i32, ptr %4, align 8, !tbaa !27
  %275 = sext i32 %274 to i64
  %276 = getelementptr inbounds double, ptr %51, i64 %275
  store double %.0345, ptr %276, align 8, !tbaa !21
  %277 = load i32, ptr %101, align 4, !tbaa !29
  %278 = sext i32 %277 to i64
  %279 = getelementptr inbounds double, ptr %51, i64 %278
  store double %247, ptr %279, align 8, !tbaa !21
  %280 = load i32, ptr %145, align 8, !tbaa !36
  %281 = sext i32 %280 to i64
  %282 = getelementptr inbounds double, ptr %51, i64 %281
  store double %251, ptr %282, align 8, !tbaa !21
  %283 = fmul double %238, %251
  %284 = load i32, ptr %131, align 4, !tbaa !34
  %285 = sext i32 %284 to i64
  %286 = getelementptr inbounds double, ptr %51, i64 %285
  store double %283, ptr %286, align 8, !tbaa !21
  %287 = load i32, ptr %47, align 8, !tbaa !13
  %288 = sext i32 %287 to i64
  %289 = shl nsw i64 %288, 3
  %290 = call i32 @gt_add(ptr noundef nonnull %7, i64 noundef -1, ptr noundef nonnull %51, i64 noundef %289) #9
  %291 = icmp eq i32 %290, 1
  br i1 %291, label %292, label %295

292:                                              ; preds = %273
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %35) #9
  %293 = getelementptr inbounds i8, ptr %7, i64 28
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %35, ptr noundef nonnull @.str, ptr noundef nonnull @.str.193, ptr noundef nonnull %293) #9
  %294 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %294, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2934, ptr noundef nonnull %35) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %35) #9
  br label %567

295:                                              ; preds = %273
  %296 = add nsw i32 %.0336389, 1
  %297 = load double, ptr %203, align 8, !tbaa !42
  %298 = fdiv double %297, %265
  %299 = fmul double %298, 9.000000e-01
  %300 = fcmp olt double %299, 5.000000e+00
  %301 = select i1 %300, double %299, double 5.000000e+00
  %302 = fmul double %.0339388, %301
  %303 = fcmp olt double %302, %200
  %304 = select i1 %303, double %302, double %200
  br label %309

305:                                              ; preds = %258
  %306 = fdiv double %266, %265
  %307 = fmul double %306, 9.000000e-01
  %308 = fmul double %.0339388, %307
  br label %309

309:                                              ; preds = %305, %295
  %.1344 = phi double [ %.0345, %295 ], [ %.0343386, %305 ]
  %.1340 = phi double [ %304, %295 ], [ %308, %305 ]
  %.1337 = phi i32 [ %296, %295 ], [ %.0336389, %305 ]
  %310 = phi <2 x double> [ %260, %295 ], [ %207, %305 ]
  %311 = fcmp ogt double %.1344, 0.000000e+00
  br i1 %311, label %206, label %._crit_edge, !llvm.loop !43

._crit_edge:                                      ; preds = %309, %185
  %.0336.lcssa = phi i32 [ 1, %185 ], [ %.1337, %309 ]
  %312 = load i32, ptr %47, align 8, !tbaa !13
  %313 = mul nsw i32 %312, %.0336.lcssa
  %314 = sext i32 %313 to i64
  %315 = shl nsw i64 %314, 3
  %316 = call noalias ptr @malloc(i64 noundef %315) #10
  %317 = getelementptr inbounds i8, ptr %4, i64 40
  store ptr %316, ptr %317, align 8, !tbaa !44
  %318 = icmp eq ptr %316, null
  br i1 %318, label %319, label %322

319:                                              ; preds = %._crit_edge
  %320 = shl i32 %313, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %36) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %36, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.200, i32 noundef %320) #9
  %321 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %321, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2948, ptr noundef nonnull %36) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %36) #9
  br label %567

322:                                              ; preds = %._crit_edge
  %323 = getelementptr inbounds i8, ptr %4, i64 36
  store i32 %.0336.lcssa, ptr %323, align 4, !tbaa !45
  %324 = call i32 @gt_getPtr(ptr noundef nonnull %7, ptr noundef nonnull %8) #9
  %325 = icmp eq i32 %324, 1
  br i1 %325, label %328, label %.preheader

.preheader:                                       ; preds = %322
  %326 = load i32, ptr %323, align 4, !tbaa !45
  %327 = icmp sgt i32 %326, 0
  br i1 %327, label %.lr.ph391.preheader, label %._crit_edge392

.lr.ph391.preheader:                              ; preds = %.preheader
  %.pre = load ptr, ptr %317, align 8, !tbaa !44
  %.pre409 = load i32, ptr %47, align 8, !tbaa !13
  br label %.lr.ph391

328:                                              ; preds = %322
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %37) #9
  %329 = getelementptr inbounds i8, ptr %7, i64 28
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %37, ptr noundef nonnull @.str, ptr noundef nonnull @.str.201, ptr noundef nonnull %329) #9
  %330 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %330, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2955, ptr noundef nonnull %37) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %37) #9
  br label %567

331:                                              ; preds = %.lr.ph391
  %332 = add nuw nsw i32 %.1390, 1
  %333 = load i32, ptr %323, align 4, !tbaa !45
  %334 = icmp slt i32 %332, %333
  br i1 %334, label %.lr.ph391, label %._crit_edge392, !llvm.loop !46

.lr.ph391:                                        ; preds = %.lr.ph391.preheader, %331
  %335 = phi i32 [ %350, %331 ], [ %.pre409, %.lr.ph391.preheader ]
  %336 = phi ptr [ %349, %331 ], [ %.pre, %.lr.ph391.preheader ]
  %337 = phi i32 [ %333, %331 ], [ %326, %.lr.ph391.preheader ]
  %.1390 = phi i32 [ %332, %331 ], [ 0, %.lr.ph391.preheader ]
  %338 = mul nsw i32 %335, %.1390
  %339 = sext i32 %338 to i64
  %340 = getelementptr inbounds double, ptr %336, i64 %339
  %341 = load ptr, ptr %8, align 8, !tbaa !47
  %342 = xor i32 %.1390, -1
  %343 = add i32 %337, %342
  %344 = mul nsw i32 %335, %343
  %345 = sext i32 %344 to i64
  %346 = getelementptr inbounds double, ptr %341, i64 %345
  %347 = sext i32 %335 to i64
  %348 = shl nsw i64 %347, 3
  call void @llvm.memcpy.p0.p0.i64(ptr align 8 %340, ptr align 8 %346, i64 %348, i1 false)
  %349 = load ptr, ptr %317, align 8, !tbaa !44
  %350 = load i32, ptr %47, align 8, !tbaa !13
  %351 = mul nsw i32 %350, %.1390
  %352 = sext i32 %351 to i64
  %353 = getelementptr inbounds double, ptr %349, i64 %352
  %.not = icmp eq ptr %340, %353
  br i1 %.not, label %331, label %354

354:                                              ; preds = %.lr.ph391
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %38) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %38, ptr noundef nonnull @.str.202) #9
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %39) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %39, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.203, ptr noundef nonnull %38) #9
  %355 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %355, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2962, ptr noundef nonnull %39) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %39) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %38) #9
  br label %567

._crit_edge392:                                   ; preds = %331, %.preheader
  %356 = call i32 @gt_free(ptr noundef nonnull %7) #9
  %357 = icmp eq i32 %356, 1
  br i1 %357, label %358, label %361

358:                                              ; preds = %._crit_edge392
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %40) #9
  %359 = getelementptr inbounds i8, ptr %7, i64 28
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %40, ptr noundef nonnull @.str, ptr noundef nonnull @.str.204, ptr noundef nonnull %359) #9
  %360 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %360, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2969, ptr noundef nonnull %40) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %40) #9
  br label %567

361:                                              ; preds = %._crit_edge392
  call void @free(ptr noundef %51) #9
  %362 = load i32, ptr %323, align 4, !tbaa !45
  %363 = icmp sgt i32 %362, 1
  %.pre413 = load ptr, ptr %317, align 8, !tbaa !44
  %.pre414 = load i32, ptr %47, align 8, !tbaa !13
  br i1 %363, label %.lr.ph397, label %._crit_edge398

.lr.ph397:                                        ; preds = %361
  %364 = getelementptr inbounds i8, ptr %1, i64 842
  %365 = getelementptr inbounds i8, ptr %1, i64 476
  %366 = getelementptr inbounds i8, ptr %1, i64 480
  %367 = getelementptr inbounds i8, ptr %1, i64 16
  %368 = getelementptr inbounds i8, ptr %3, i64 192
  %369 = getelementptr inbounds i8, ptr %2, i64 606
  %370 = getelementptr inbounds i8, ptr %3, i64 136
  %371 = getelementptr inbounds i8, ptr %3, i64 80
  br label %372

372:                                              ; preds = %.lr.ph397, %516
  %373 = phi i32 [ %.pre414, %.lr.ph397 ], [ %518, %516 ]
  %374 = phi ptr [ %.pre413, %.lr.ph397 ], [ %519, %516 ]
  %.2.in394 = phi i32 [ %362, %.lr.ph397 ], [ %.2395, %516 ]
  %.2395 = add nsw i32 %.2.in394, -1
  %375 = mul nsw i32 %373, %.2395
  %376 = load i32, ptr %4, align 8, !tbaa !27
  %377 = add nsw i32 %375, %376
  %378 = sext i32 %377 to i64
  %379 = getelementptr inbounds double, ptr %374, i64 %378
  %380 = load double, ptr %379, align 8, !tbaa !21
  %381 = call i32 @background_tau_of_z(ptr noundef %1, double noundef %380, ptr noundef nonnull %12) #9
  %382 = icmp eq i32 %381, 1
  br i1 %382, label %383, label %385

383:                                              ; preds = %372
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %41) #9
  %384 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %41, ptr noundef nonnull @.str, ptr noundef nonnull @.str.190, ptr noundef nonnull %384) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %369, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2982, ptr noundef nonnull %41) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %41) #9
  br label %567

385:                                              ; preds = %372
  %386 = load double, ptr %12, align 8, !tbaa !21
  %387 = load i16, ptr %364, align 2, !tbaa !48
  %388 = load i16, ptr %114, align 2, !tbaa !32
  %389 = call i32 @background_at_tau(ptr noundef %1, double noundef %386, i16 noundef signext %387, i16 noundef signext %388, ptr noundef nonnull %13, ptr noundef %5) #9
  %390 = icmp eq i32 %389, 1
  br i1 %390, label %391, label %393

391:                                              ; preds = %385
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %42) #9
  %392 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %42, ptr noundef nonnull @.str, ptr noundef nonnull @.str.205, ptr noundef nonnull %392) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %369, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 2991, ptr noundef nonnull %42) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %42) #9
  br label %567

393:                                              ; preds = %385
  %394 = load ptr, ptr %317, align 8, !tbaa !44
  %395 = load i32, ptr %47, align 8, !tbaa !13
  %396 = mul nsw i32 %395, %.2395
  %397 = load i32, ptr %4, align 8, !tbaa !27
  %398 = add nsw i32 %396, %397
  %399 = sext i32 %398 to i64
  %400 = getelementptr inbounds double, ptr %394, i64 %399
  %401 = load double, ptr %400, align 8, !tbaa !21
  %402 = add nsw i32 %.2.in394, -2
  %403 = mul nsw i32 %395, %402
  %404 = add nsw i32 %403, %397
  %405 = sext i32 %404 to i64
  %406 = getelementptr inbounds double, ptr %394, i64 %405
  %407 = load double, ptr %406, align 8, !tbaa !21
  %408 = fsub double %401, %407
  %409 = fadd double %380, 1.000000e+00
  %410 = fmul double %409, %409
  %411 = load double, ptr %124, align 8, !tbaa !33
  %412 = fmul double %410, %411
  %413 = load i32, ptr %101, align 4, !tbaa !29
  %414 = add nsw i32 %413, %396
  %415 = sext i32 %414 to i64
  %416 = getelementptr inbounds double, ptr %394, i64 %415
  %417 = load double, ptr %416, align 8, !tbaa !21
  %418 = fmul double %412, %417
  %419 = fmul double %418, 0x3A15151E37FBF236
  %420 = fmul double %419, 0x449A22FF260E3826
  %421 = load double, ptr %2, align 8, !tbaa !5
  %422 = call double @llvm.fmuladd.f64(double %421, double 0xBFE7F14DA42541E8, double 1.000000e+00)
  %423 = fsub double 1.000000e+00, %421
  %424 = call double @llvm.fmuladd.f64(double %417, double %423, double %422)
  %425 = fdiv double 0x3A60930446E176B5, %424
  %426 = fdiv double 2.000000e+00, %409
  %427 = load i32, ptr %159, align 8, !tbaa !38
  %428 = add nsw i32 %427, %396
  %429 = sext i32 %428 to i64
  %430 = getelementptr inbounds double, ptr %394, i64 %429
  %431 = load double, ptr %430, align 8, !tbaa !21
  %432 = fmul double %425, 2.000000e+00
  %433 = fdiv double %432, 0x39B279DC93EB3F4E
  %434 = fmul double %433, 4.000000e+00
  %435 = load i32, ptr %365, align 4, !tbaa !49
  %436 = sext i32 %435 to i64
  %437 = getelementptr inbounds double, ptr %5, i64 %436
  %438 = load double, ptr %437, align 8, !tbaa !21
  %439 = fmul double %438, %434
  %440 = load i32, ptr %366, align 8, !tbaa !50
  %441 = sext i32 %440 to i64
  %442 = getelementptr inbounds double, ptr %5, i64 %441
  %443 = load double, ptr %442, align 8, !tbaa !21
  %444 = load double, ptr %367, align 8, !tbaa !51
  %445 = fneg double %431
  %446 = call double @llvm.fmuladd.f64(double %444, double %409, double %445)
  %447 = load i32, ptr %135, align 4, !tbaa !35
  %448 = sext i32 %447 to i64
  %449 = getelementptr inbounds double, ptr %5, i64 %448
  %450 = load double, ptr %449, align 8, !tbaa !21
  %451 = fdiv double %439, -3.000000e+00
  %452 = fdiv double %451, %443
  %453 = fmul double %420, %452
  %454 = fmul double %446, %453
  %455 = fdiv double %454, %450
  %456 = call double @llvm.fmuladd.f64(double %426, double %431, double %455)
  %457 = load double, ptr %368, align 8, !tbaa !52
  %458 = fcmp ogt double %457, 0.000000e+00
  br i1 %458, label %459, label %516

459:                                              ; preds = %393
  %460 = call i32 @thermodynamics_energy_injection(ptr noundef %0, ptr noundef nonnull %1, ptr noundef nonnull %3, double noundef %380, ptr noundef nonnull %11, ptr noundef nonnull %369)
  %461 = icmp eq i32 %460, 1
  br i1 %461, label %462, label %463

462:                                              ; preds = %459
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %43) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %43, ptr noundef nonnull @.str, ptr noundef nonnull @.str.206, ptr noundef nonnull %369) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %369, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 3010, ptr noundef nonnull %43) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %43) #9
  br label %567

463:                                              ; preds = %459
  %464 = load ptr, ptr %317, align 8, !tbaa !44
  %465 = load i32, ptr %47, align 8, !tbaa !13
  %466 = mul nsw i32 %465, %.2395
  %467 = load i32, ptr %101, align 4, !tbaa !29
  %468 = add nsw i32 %466, %467
  %469 = sext i32 %468 to i64
  %470 = getelementptr inbounds double, ptr %464, i64 %469
  %471 = load double, ptr %470, align 8, !tbaa !21
  store double %471, ptr %9, align 8, !tbaa !21
  %472 = fcmp olt double %471, 1.000000e+00
  br i1 %472, label %473, label %487

473:                                              ; preds = %463
  %474 = call double @pow(double noundef %471, double noundef 3.001340e-01) #9
  %475 = fsub double 1.000000e+00, %474
  %476 = call double @pow(double noundef %475, double noundef 1.510350e+00) #9
  %477 = fsub double 1.000000e+00, %476
  %478 = fmul double %477, 0x3FEFE640A6B93CCD
  %479 = fcmp olt double %478, 1.000000e+00
  br i1 %479, label %480, label %487

480:                                              ; preds = %473
  %481 = load double, ptr %9, align 8, !tbaa !21
  %482 = call double @pow(double noundef %481, double noundef 3.001340e-01) #9
  %483 = fsub double 1.000000e+00, %482
  %484 = call double @pow(double noundef %483, double noundef 1.510350e+00) #9
  %485 = fsub double 1.000000e+00, %484
  %486 = fmul double %485, 0x3FEFE640A6B93CCD
  br label %487

487:                                              ; preds = %463, %480, %473
  %.0338 = phi double [ %486, %480 ], [ 1.000000e+00, %473 ], [ 1.000000e+00, %463 ]
  %488 = load double, ptr %11, align 8, !tbaa !21
  %489 = fmul double %488, 0x44A4733969CB9B70
  %490 = fmul double %.0338, %489
  %491 = load double, ptr %370, align 8, !tbaa !53
  %492 = call double @pow(double noundef %409, double noundef 3.000000e+00) #9
  %493 = fmul double %491, %492
  %494 = fdiv double %490, %493
  %495 = load double, ptr %371, align 8, !tbaa !54
  %496 = fadd double %495, 1.000000e+00
  %497 = load ptr, ptr %317, align 8, !tbaa !44
  %498 = load i32, ptr %47, align 8, !tbaa !13
  %499 = mul nsw i32 %498, %.2395
  %500 = load i32, ptr %101, align 4, !tbaa !29
  %501 = add nsw i32 %499, %500
  %502 = sext i32 %501 to i64
  %503 = getelementptr inbounds double, ptr %497, i64 %502
  %504 = load double, ptr %503, align 8, !tbaa !21
  %505 = fadd double %496, %504
  %506 = fdiv double %494, %505
  %507 = load i32, ptr %135, align 4, !tbaa !35
  %508 = sext i32 %507 to i64
  %509 = getelementptr inbounds double, ptr %5, i64 %508
  %510 = load double, ptr %509, align 8, !tbaa !21
  %511 = fmul double %510, 0x41B1DE784A000000
  %512 = fdiv double %511, 0x449A22FF260E3826
  %513 = fmul double %409, %512
  %514 = fdiv double %506, %513
  %515 = fsub double %456, %514
  %.pre412 = load i32, ptr %159, align 8, !tbaa !38
  %.pre416 = add nsw i32 %499, %.pre412
  %.pre418 = sext i32 %.pre416 to i64
  %.pre420 = mul nsw i32 %498, %402
  br label %516

516:                                              ; preds = %487, %393
  %.pre-phi421 = phi i32 [ %.pre420, %487 ], [ %403, %393 ]
  %.pre-phi419 = phi i64 [ %.pre418, %487 ], [ %429, %393 ]
  %517 = phi i32 [ %.pre412, %487 ], [ %427, %393 ]
  %518 = phi i32 [ %498, %487 ], [ %395, %393 ]
  %519 = phi ptr [ %497, %487 ], [ %394, %393 ]
  %.0346 = phi double [ %515, %487 ], [ %456, %393 ]
  %520 = getelementptr inbounds double, ptr %519, i64 %.pre-phi419
  %521 = load double, ptr %520, align 8, !tbaa !21
  %522 = fneg double %.0346
  %523 = call double @llvm.fmuladd.f64(double %522, double %408, double %521)
  %524 = add nsw i32 %.pre-phi421, %517
  %525 = sext i32 %524 to i64
  %526 = getelementptr inbounds double, ptr %519, i64 %525
  store double %523, ptr %526, align 8, !tbaa !21
  %527 = fmul double %425, 0x4373F4D4EACDD756
  %528 = fdiv double 0x3B30B0E7F1528C6D, %527
  %529 = fmul double %528, %523
  %530 = load i32, ptr %168, align 4, !tbaa !39
  %531 = add nsw i32 %530, %.pre-phi421
  %532 = sext i32 %531 to i64
  %533 = getelementptr inbounds double, ptr %519, i64 %532
  store double %529, ptr %533, align 8, !tbaa !21
  %534 = fdiv double %409, 3.000000e+00
  %535 = fmul double %534, %.0346
  %536 = load double, ptr %526, align 8, !tbaa !21
  %537 = fdiv double %535, %536
  %538 = fadd double %537, 1.000000e+00
  %539 = fmul double %529, %538
  %540 = load i32, ptr %173, align 8, !tbaa !40
  %541 = add nsw i32 %540, %.pre-phi421
  %542 = sext i32 %541 to i64
  %543 = getelementptr inbounds double, ptr %519, i64 %542
  store double %539, ptr %543, align 8, !tbaa !21
  %544 = icmp ugt i32 %.2.in394, 2
  br i1 %544, label %372, label %._crit_edge398.loopexit, !llvm.loop !55

._crit_edge398.loopexit:                          ; preds = %516
  %.pre415 = load i32, ptr %323, align 4, !tbaa !45
  br label %._crit_edge398

._crit_edge398:                                   ; preds = %._crit_edge398.loopexit, %361
  %545 = phi i32 [ %.pre415, %._crit_edge398.loopexit ], [ %362, %361 ]
  %546 = phi i32 [ %518, %._crit_edge398.loopexit ], [ %.pre414, %361 ]
  %547 = phi ptr [ %519, %._crit_edge398.loopexit ], [ %.pre413, %361 ]
  %548 = load i32, ptr %4, align 8, !tbaa !27
  %549 = load i32, ptr %145, align 8, !tbaa !36
  %550 = getelementptr inbounds i8, ptr %4, i64 28
  %551 = load i32, ptr %550, align 4, !tbaa !56
  %552 = getelementptr inbounds i8, ptr %2, i64 606
  %553 = call i32 @array_spline(ptr noundef %547, i32 noundef %546, i32 noundef %545, i32 noundef %548, i32 noundef %549, i32 noundef %551, i16 noundef signext 1, ptr noundef nonnull %552) #9
  %554 = icmp eq i32 %553, 1
  br i1 %554, label %555, label %556

555:                                              ; preds = %._crit_edge398
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %44) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %44, ptr noundef nonnull @.str, ptr noundef nonnull @.str.207, ptr noundef nonnull %552) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %552, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 3058, ptr noundef nonnull %44) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %44) #9
  br label %567

556:                                              ; preds = %._crit_edge398
  %557 = load ptr, ptr %317, align 8, !tbaa !44
  %558 = load i32, ptr %47, align 8, !tbaa !13
  %559 = load i32, ptr %323, align 4, !tbaa !45
  %560 = load i32, ptr %4, align 8, !tbaa !27
  %561 = load i32, ptr %145, align 8, !tbaa !36
  %562 = load i32, ptr %550, align 4, !tbaa !56
  %563 = getelementptr inbounds i8, ptr %4, i64 48
  %564 = call i32 @array_integrate_all_spline(ptr noundef %557, i32 noundef %558, i32 noundef %559, i32 noundef %560, i32 noundef %561, i32 noundef %562, ptr noundef nonnull %563, ptr noundef nonnull %552) #9
  %565 = icmp eq i32 %564, 1
  br i1 %565, label %566, label %567

566:                                              ; preds = %556
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %45) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %45, ptr noundef nonnull @.str, ptr noundef nonnull @.str.208, ptr noundef nonnull %552) #9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %552, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_sample, i32 noundef 3070, ptr noundef nonnull %45) #9
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %45) #9
  br label %567

567:                                              ; preds = %556, %566, %555, %462, %391, %383, %358, %354, %328, %319, %292, %271, %256, %240, %231, %222, %217, %210, %182, %141, %118, %107, %97, %85, %74, %53
  %.0 = phi i32 [ 1, %53 ], [ 1, %74 ], [ 1, %85 ], [ 1, %97 ], [ 1, %107 ], [ 1, %118 ], [ 1, %141 ], [ 1, %182 ], [ 1, %210 ], [ 1, %217 ], [ 1, %222 ], [ 1, %231 ], [ 1, %240 ], [ 1, %256 ], [ 1, %271 ], [ 1, %292 ], [ 1, %319 ], [ 1, %328 ], [ 1, %354 ], [ 1, %358 ], [ 1, %383 ], [ 1, %391 ], [ 1, %462 ], [ 1, %555 ], [ 1, %566 ], [ 0, %556 ]
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %13) #9
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %12) #9
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %11) #9
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %10) #9
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %9) #9
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %8) #9
  call void @llvm.lifetime.end.p0(i64 2080, ptr nonnull %7) #9
  ret i32 %.0
}

declare i32 @gt_init(ptr noundef) local_unnamed_addr #3

declare i32 @gt_add(ptr noundef, i64 noundef, ptr noundef, i64 noundef) local_unnamed_addr #3

declare i32 @gt_getPtr(ptr noundef, ptr noundef) local_unnamed_addr #3

; Function Attrs: mustprogress nocallback nofree nounwind willreturn memory(argmem: readwrite)
declare void @llvm.memcpy.p0.p0.i64(ptr noalias nocapture writeonly, ptr noalias nocapture readonly, i64, i1 immarg) #7

declare i32 @gt_free(ptr noundef) local_unnamed_addr #3

declare i32 @array_spline(ptr noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, i16 noundef signext, ptr noundef) local_unnamed_addr #3

declare i32 @array_integrate_all_spline(ptr noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, i32 noundef, ptr noundef, ptr noundef) local_unnamed_addr #3

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare <2 x double> @llvm.fabs.v2f64(<2 x double>) #8

attributes #0 = { mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { 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 #2 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #3 = { "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 #4 = { mustprogress nofree nounwind willreturn allockind("alloc,uninitialized") allocsize(0) memory(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 #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 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 #7 = { mustprogress nocallback nofree nounwind willreturn memory(argmem: readwrite) }
attributes #8 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #9 = { nounwind }
attributes #10 = { nounwind allocsize(0) }

!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, !7, i64 0}
!6 = !{!"thermo", !7, i64 0, !10, i64 8, !10, i64 12, !10, i64 16, !7, i64 24, !7, i64 32, !11, i64 40, !11, i64 42, !7, i64 48, !7, i64 56, !7, i64 64, !7, i64 72, !10, i64 80, !12, i64 88, !12, i64 96, !7, i64 104, !10, i64 112, !12, i64 120, !12, i64 128, !7, i64 136, !10, i64 144, !12, i64 152, !12, i64 160, !7, i64 168, !11, i64 176, !7, i64 184, !7, i64 192, !7, i64 200, !7, i64 208, !7, i64 216, !7, i64 224, !7, i64 232, !7, i64 240, !7, i64 248, !7, i64 256, !7, i64 264, !10, i64 272, !10, i64 276, !10, i64 280, !10, i64 284, !10, i64 288, !10, i64 292, !10, i64 296, !10, i64 300, !10, i64 304, !10, i64 308, !10, i64 312, !10, i64 316, !10, i64 320, !10, i64 324, !10, i64 328, !10, i64 332, !10, i64 336, !10, i64 340, !10, i64 344, !10, i64 348, !10, i64 352, !10, i64 356, !10, i64 360, !10, i64 364, !10, i64 368, !10, i64 372, !10, i64 376, !12, i64 384, !12, i64 392, !12, i64 400, !7, i64 408, !7, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !7, i64 464, !7, i64 472, !7, i64 480, !7, i64 488, !7, i64 496, !7, i64 504, !7, i64 512, !7, i64 520, !7, i64 528, !7, i64 536, !7, i64 544, !7, i64 552, !7, i64 560, !7, i64 568, !7, i64 576, !7, i64 584, !7, i64 592, !11, i64 600, !11, i64 602, !11, i64 604, !8, i64 606}
!7 = !{!"double", !8, i64 0}
!8 = !{!"omnipotent char", !9, i64 0}
!9 = !{!"Simple C/C++ TBAA"}
!10 = !{!"int", !8, i64 0}
!11 = !{!"short", !8, i64 0}
!12 = !{!"any pointer", !8, i64 0}
!13 = !{!14, !10, i64 32}
!14 = !{!"reionization", !10, i64 0, !10, i64 4, !10, i64 8, !10, i64 12, !10, i64 16, !10, i64 20, !10, i64 24, !10, i64 28, !10, i64 32, !10, i64 36, !12, i64 40, !7, i64 48, !10, i64 56, !10, i64 60, !10, i64 64, !10, i64 68, !10, i64 72, !10, i64 76, !10, i64 80, !10, i64 84, !10, i64 88, !10, i64 92, !10, i64 96, !10, i64 100, !10, i64 104, !12, i64 112, !10, i64 120, !10, i64 124}
!15 = !{!16, !12, i64 32}
!16 = !{!"recombination", !10, i64 0, !10, i64 4, !10, i64 8, !10, i64 12, !10, i64 16, !10, i64 20, !10, i64 24, !10, i64 28, !12, 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, !7, i64 112, !7, i64 120, !7, i64 128, !7, i64 136, !7, i64 144, !7, i64 152, !7, i64 160, !7, i64 168, !7, i64 176, !7, i64 184, !7, i64 192, !11, i64 200, !7, i64 208, !7, i64 216, !7, i64 224, !7, i64 232, !7, i64 240, !7, i64 248, !7, i64 256}
!17 = !{!16, !10, i64 24}
!18 = !{!16, !10, i64 0}
!19 = !{!14, !12, i64 112}
!20 = !{!14, !10, i64 104}
!21 = !{!7, !7, i64 0}
!22 = !{!23, !10, i64 368}
!23 = !{!"precision", !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, !8, i64 104, !7, i64 360, !10, i64 368, !7, i64 376, !10, i64 384, !10, i64 388, !7, i64 392, !10, i64 400, !7, i64 408, !10, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !7, i64 464, !7, i64 472, !7, i64 480, !7, i64 488, !7, i64 496, !7, i64 504, !7, i64 512, !7, i64 520, !7, i64 528, !7, i64 536, !7, i64 544, !7, i64 552, !7, i64 560, !7, i64 568, !7, i64 576, !7, i64 584, !7, i64 592, !7, i64 600, !7, i64 608, !7, i64 616, !10, i64 624, !8, i64 628, !8, i64 884, !8, i64 1140, !7, i64 1400, !7, i64 1408, !7, i64 1416, !7, i64 1424, !7, i64 1432, !7, i64 1440, !7, i64 1448, !7, i64 1456, !7, i64 1464, !7, i64 1472, !7, i64 1480, !7, i64 1488, !7, i64 1496, !7, i64 1504, !7, i64 1512, !7, i64 1520, !10, i64 1528, !7, i64 1536, !7, i64 1544, !10, i64 1552, !10, i64 1556, !10, i64 1560, !10, i64 1564, !10, i64 1568, !10, i64 1572, !10, i64 1576, !10, i64 1580, !7, i64 1584, !7, i64 1592, !7, i64 1600, !7, i64 1608, !7, i64 1616, !7, i64 1624, !7, i64 1632, !7, i64 1640, !10, i64 1648, !7, i64 1656, !7, i64 1664, !10, i64 1672, !7, i64 1680, !7, i64 1688, !10, i64 1696, !7, i64 1704, !10, i64 1712, !7, i64 1720, !7, i64 1728, !10, i64 1736, !7, i64 1744, !7, i64 1752, !7, i64 1760, !10, i64 1768, !7, i64 1776, !7, i64 1784, !7, i64 1792, !7, i64 1800, !7, i64 1808, !10, i64 1816, !7, i64 1824, !7, i64 1832, !7, i64 1840, !7, i64 1848, !7, i64 1856, !7, i64 1864, !7, i64 1872, !10, i64 1880, !7, i64 1888, !7, i64 1896, !7, i64 1904, !7, i64 1912, !7, i64 1920, !7, i64 1928, !7, i64 1936, !7, i64 1944, !7, i64 1952, !7, i64 1960, !7, i64 1968, !7, i64 1976, !7, i64 1984, !7, i64 1992, !7, i64 2000, !7, i64 2008, !7, i64 2016, !7, i64 2024, !7, i64 2032, !7, i64 2040, !7, i64 2048, !7, i64 2056, !7, i64 2064, !7, i64 2072, !7, i64 2080, !7, i64 2088, !7, i64 2096, !7, i64 2104, !7, i64 2112, !7, i64 2120, !7, i64 2128, !7, i64 2136, !7, i64 2144, !7, i64 2152, !7, i64 2160, !7, i64 2168, !7, i64 2176, !7, i64 2184, !7, i64 2192, !7, i64 2200, !7, i64 2208, !7, i64 2216, !7, i64 2224, !7, i64 2232, !7, i64 2240, !7, i64 2248, !10, i64 2256, !7, i64 2264, !7, i64 2272, !7, i64 2280, !7, i64 2288, !10, i64 2296, !7, i64 2304, !7, i64 2312, !10, i64 2320, !10, i64 2324, !10, i64 2328, !7, i64 2336, !7, i64 2344, !8, i64 2352}
!24 = distinct !{!24, !25}
!25 = !{!"llvm.loop.mustprogress"}
!26 = !{!23, !7, i64 592}
!27 = !{!14, !10, i64 0}
!28 = !{!14, !10, i64 124}
!29 = !{!14, !10, i64 4}
!30 = !{!31, !11, i64 840}
!31 = !{!"background", !7, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !7, i64 40, !7, i64 48, !10, i64 56, !7, i64 64, !7, i64 72, !7, i64 80, !7, i64 88, !11, i64 96, !7, i64 104, !7, i64 112, !7, i64 120, !7, i64 128, !7, i64 136, !7, i64 144, !7, i64 152, !7, i64 160, !7, i64 168, !11, i64 176, !7, i64 184, !7, i64 192, !12, i64 200, !10, i64 208, !10, i64 212, !7, i64 216, !10, i64 224, !12, i64 232, !12, i64 240, !7, i64 248, !12, i64 256, !7, i64 264, !12, i64 272, !7, i64 280, !12, i64 288, !7, i64 296, !12, i64 304, !12, i64 312, !12, i64 320, !7, i64 328, !7, i64 336, !7, i64 344, !7, i64 352, !10, i64 360, !12, i64 368, !7, i64 376, !7, i64 384, !7, i64 392, !7, i64 400, !7, i64 408, !7, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !10, i64 464, !10, i64 468, !10, i64 472, !10, i64 476, !10, i64 480, !10, i64 484, !10, i64 488, !10, i64 492, !10, i64 496, !10, i64 500, !10, i64 504, !10, i64 508, !10, i64 512, !10, i64 516, !10, i64 520, !10, i64 524, !10, i64 528, !10, i64 532, !10, i64 536, !10, i64 540, !10, i64 544, !10, i64 548, !10, i64 552, !10, i64 556, !10, i64 560, !10, i64 564, !10, i64 568, !10, i64 572, !10, i64 576, !10, i64 580, !10, i64 584, !10, i64 588, !10, i64 592, !10, i64 596, !10, i64 600, !10, i64 604, !10, i64 608, !10, i64 612, !10, i64 616, !10, i64 620, !10, i64 624, !10, i64 628, !12, i64 632, !12, i64 640, !12, i64 648, !12, i64 656, !12, i64 664, !10, i64 672, !10, i64 676, !10, i64 680, !10, i64 684, !10, i64 688, !10, i64 692, !10, i64 696, !10, i64 700, !10, i64 704, !10, i64 708, !10, i64 712, !10, i64 716, !10, i64 720, !11, i64 724, !11, i64 726, !11, i64 728, !11, i64 730, !11, i64 732, !11, i64 734, !11, i64 736, !11, i64 738, !11, i64 740, !11, i64 742, !11, i64 744, !12, i64 752, !12, i64 760, !12, i64 768, !12, i64 776, !12, i64 784, !12, i64 792, !12, i64 800, !12, i64 808, !12, i64 816, !12, i64 824, !12, i64 832, !11, i64 840, !11, i64 842, !11, i64 844, !11, i64 846, !11, i64 848, !11, i64 850, !8, i64 852, !11, i64 2900, !8, i64 2902}
!32 = !{!31, !11, i64 846}
!33 = !{!6, !7, i64 592}
!34 = !{!14, !10, i64 20}
!35 = !{!31, !10, i64 468}
!36 = !{!14, !10, i64 24}
!37 = !{!16, !10, i64 8}
!38 = !{!14, !10, i64 8}
!39 = !{!14, !10, i64 12}
!40 = !{!14, !10, i64 16}
!41 = !{!23, !7, i64 2344}
!42 = !{!23, !7, i64 600}
!43 = distinct !{!43, !25}
!44 = !{!14, !12, i64 40}
!45 = !{!14, !10, i64 36}
!46 = distinct !{!46, !25}
!47 = !{!12, !12, i64 0}
!48 = !{!31, !11, i64 842}
!49 = !{!31, !10, i64 476}
!50 = !{!31, !10, i64 480}
!51 = !{!31, !7, i64 16}
!52 = !{!16, !7, i64 192}
!53 = !{!16, !7, i64 136}
!54 = !{!16, !7, i64 80}
!55 = distinct !{!55, !25}
!56 = !{!14, !10, i64 28}


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