; 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.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 }

@.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.137 = external hidden unnamed_addr constant [59 x i8], align 1
@__func__.thermodynamics_reionization_function = external hidden unnamed_addr constant [37 x i8], align 1
@.str.138 = external hidden unnamed_addr constant [38 x i8], align 1
@.str.139 = external hidden unnamed_addr constant [10 x i8], align 1
@.str.140 = external hidden unnamed_addr constant [1 x i8], align 1
@.str.141 = external hidden unnamed_addr constant [8 x i8], align 1
@.str.142 = external hidden unnamed_addr constant [8 x i8], align 1
@.str.143 = external hidden unnamed_addr constant [10 x i8], align 1
@.str.144 = external hidden unnamed_addr constant [7 x i8], align 1
@.str.145 = external hidden unnamed_addr constant [41 x i8], align 1
@.str.146 = external hidden unnamed_addr constant [7 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 uwtable
define dso_local i32 @thermodynamics_reionization_function(double noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3) #4 {
  %5 = alloca i32, align 4
  %6 = alloca double, align 8
  %7 = alloca ptr, align 8
  %8 = alloca ptr, align 8
  %9 = alloca ptr, align 8
  %10 = alloca double, align 8
  %11 = alloca i32, align 4
  %12 = alloca double, align 8
  %13 = alloca i32, align 4
  %14 = alloca double, align 8
  %15 = alloca double, align 8
  %16 = alloca double, align 8
  %17 = alloca double, align 8
  %18 = alloca double, align 8
  %19 = alloca i32, align 4
  %20 = alloca [2048 x i8], align 16
  %21 = alloca [2048 x i8], align 16
  %22 = alloca [2048 x i8], align 16
  %23 = alloca double, align 8
  %24 = alloca double, align 8
  %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 double, align 8
  %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
  store double %0, ptr %6, align 8, !tbaa !5
  store ptr %1, ptr %7, align 8, !tbaa !9
  store ptr %2, ptr %8, align 8, !tbaa !9
  store ptr %3, ptr %9, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %10) #5
  call void @llvm.lifetime.start.p0(i64 4, ptr %11) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #5
  call void @llvm.lifetime.start.p0(i64 4, ptr %13) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %17) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %18) #5
  %34 = load ptr, ptr %7, align 8, !tbaa !9
  %35 = getelementptr inbounds %struct.thermo, ptr %34, i32 0, i32 2
  %36 = load i32, ptr %35, align 4, !tbaa !11
  %37 = icmp eq i32 %36, 1
  br i1 %37, label %43, label %38

38:                                               ; preds = %4
  %39 = load ptr, ptr %7, align 8, !tbaa !9
  %40 = getelementptr inbounds %struct.thermo, ptr %39, i32 0, i32 2
  %41 = load i32, ptr %40, align 4, !tbaa !11
  %42 = icmp eq i32 %41, 3
  br i1 %42, label %43, label %261

43:                                               ; preds = %38, %4
  %44 = load double, ptr %6, align 8, !tbaa !5
  %45 = load ptr, ptr %8, align 8, !tbaa !9
  %46 = getelementptr inbounds %struct.reionization, ptr %45, i32 0, i32 25
  %47 = load ptr, ptr %46, align 8, !tbaa !15
  %48 = load ptr, ptr %8, align 8, !tbaa !9
  %49 = getelementptr inbounds %struct.reionization, ptr %48, i32 0, i32 24
  %50 = load i32, ptr %49, align 8, !tbaa !17
  %51 = sext i32 %50 to i64
  %52 = getelementptr inbounds double, ptr %47, i64 %51
  %53 = load double, ptr %52, align 8, !tbaa !5
  %54 = fcmp ogt double %44, %53
  br i1 %54, label %55, label %66

55:                                               ; preds = %43
  %56 = load ptr, ptr %8, align 8, !tbaa !9
  %57 = getelementptr inbounds %struct.reionization, ptr %56, i32 0, i32 25
  %58 = load ptr, ptr %57, align 8, !tbaa !15
  %59 = load ptr, ptr %8, align 8, !tbaa !9
  %60 = getelementptr inbounds %struct.reionization, ptr %59, i32 0, i32 15
  %61 = load i32, ptr %60, align 4, !tbaa !18
  %62 = sext i32 %61 to i64
  %63 = getelementptr inbounds double, ptr %58, i64 %62
  %64 = load double, ptr %63, align 8, !tbaa !5
  %65 = load ptr, ptr %9, align 8, !tbaa !9
  store double %64, ptr %65, align 8, !tbaa !5
  br label %260

66:                                               ; preds = %43
  %67 = load ptr, ptr %8, align 8, !tbaa !9
  %68 = getelementptr inbounds %struct.reionization, ptr %67, i32 0, i32 25
  %69 = load ptr, ptr %68, align 8, !tbaa !15
  %70 = load ptr, ptr %8, align 8, !tbaa !9
  %71 = getelementptr inbounds %struct.reionization, ptr %70, i32 0, i32 12
  %72 = load i32, ptr %71, align 8, !tbaa !19
  %73 = sext i32 %72 to i64
  %74 = getelementptr inbounds double, ptr %69, i64 %73
  %75 = load double, ptr %74, align 8, !tbaa !5
  %76 = fadd double 1.000000e+00, %75
  %77 = load ptr, ptr %8, align 8, !tbaa !9
  %78 = getelementptr inbounds %struct.reionization, ptr %77, i32 0, i32 25
  %79 = load ptr, ptr %78, align 8, !tbaa !15
  %80 = load ptr, ptr %8, align 8, !tbaa !9
  %81 = getelementptr inbounds %struct.reionization, ptr %80, i32 0, i32 13
  %82 = load i32, ptr %81, align 4, !tbaa !20
  %83 = sext i32 %82 to i64
  %84 = getelementptr inbounds double, ptr %79, i64 %83
  %85 = load double, ptr %84, align 8, !tbaa !5
  %86 = call double @pow(double noundef %76, double noundef %85) #5
  %87 = load double, ptr %6, align 8, !tbaa !5
  %88 = fadd double 1.000000e+00, %87
  %89 = load ptr, ptr %8, align 8, !tbaa !9
  %90 = getelementptr inbounds %struct.reionization, ptr %89, i32 0, i32 25
  %91 = load ptr, ptr %90, align 8, !tbaa !15
  %92 = load ptr, ptr %8, align 8, !tbaa !9
  %93 = getelementptr inbounds %struct.reionization, ptr %92, i32 0, i32 13
  %94 = load i32, ptr %93, align 4, !tbaa !20
  %95 = sext i32 %94 to i64
  %96 = getelementptr inbounds double, ptr %91, i64 %95
  %97 = load double, ptr %96, align 8, !tbaa !5
  %98 = call double @pow(double noundef %88, double noundef %97) #5
  %99 = fsub double %86, %98
  %100 = load ptr, ptr %8, align 8, !tbaa !9
  %101 = getelementptr inbounds %struct.reionization, ptr %100, i32 0, i32 25
  %102 = load ptr, ptr %101, align 8, !tbaa !15
  %103 = load ptr, ptr %8, align 8, !tbaa !9
  %104 = getelementptr inbounds %struct.reionization, ptr %103, i32 0, i32 13
  %105 = load i32, ptr %104, align 4, !tbaa !20
  %106 = sext i32 %105 to i64
  %107 = getelementptr inbounds double, ptr %102, i64 %106
  %108 = load double, ptr %107, align 8, !tbaa !5
  %109 = load ptr, ptr %8, align 8, !tbaa !9
  %110 = getelementptr inbounds %struct.reionization, ptr %109, i32 0, i32 25
  %111 = load ptr, ptr %110, align 8, !tbaa !15
  %112 = load ptr, ptr %8, align 8, !tbaa !9
  %113 = getelementptr inbounds %struct.reionization, ptr %112, i32 0, i32 12
  %114 = load i32, ptr %113, align 8, !tbaa !19
  %115 = sext i32 %114 to i64
  %116 = getelementptr inbounds double, ptr %111, i64 %115
  %117 = load double, ptr %116, align 8, !tbaa !5
  %118 = fadd double 1.000000e+00, %117
  %119 = load ptr, ptr %8, align 8, !tbaa !9
  %120 = getelementptr inbounds %struct.reionization, ptr %119, i32 0, i32 25
  %121 = load ptr, ptr %120, align 8, !tbaa !15
  %122 = load ptr, ptr %8, align 8, !tbaa !9
  %123 = getelementptr inbounds %struct.reionization, ptr %122, i32 0, i32 13
  %124 = load i32, ptr %123, align 4, !tbaa !20
  %125 = sext i32 %124 to i64
  %126 = getelementptr inbounds double, ptr %121, i64 %125
  %127 = load double, ptr %126, align 8, !tbaa !5
  %128 = fsub double %127, 1.000000e+00
  %129 = call double @pow(double noundef %118, double noundef %128) #5
  %130 = fmul double %108, %129
  %131 = fdiv double %99, %130
  %132 = load ptr, ptr %8, align 8, !tbaa !9
  %133 = getelementptr inbounds %struct.reionization, ptr %132, i32 0, i32 25
  %134 = load ptr, ptr %133, align 8, !tbaa !15
  %135 = load ptr, ptr %8, align 8, !tbaa !9
  %136 = getelementptr inbounds %struct.reionization, ptr %135, i32 0, i32 14
  %137 = load i32, ptr %136, align 8, !tbaa !21
  %138 = sext i32 %137 to i64
  %139 = getelementptr inbounds double, ptr %134, i64 %138
  %140 = load double, ptr %139, align 8, !tbaa !5
  %141 = fdiv double %131, %140
  store double %141, ptr %10, align 8, !tbaa !5
  %142 = load ptr, ptr %7, align 8, !tbaa !9
  %143 = getelementptr inbounds %struct.thermo, ptr %142, i32 0, i32 2
  %144 = load i32, ptr %143, align 4, !tbaa !11
  %145 = icmp eq i32 %144, 1
  br i1 %145, label %146, label %182

146:                                              ; preds = %66
  %147 = load ptr, ptr %8, align 8, !tbaa !9
  %148 = getelementptr inbounds %struct.reionization, ptr %147, i32 0, i32 25
  %149 = load ptr, ptr %148, align 8, !tbaa !15
  %150 = load ptr, ptr %8, align 8, !tbaa !9
  %151 = getelementptr inbounds %struct.reionization, ptr %150, i32 0, i32 16
  %152 = load i32, ptr %151, align 8, !tbaa !22
  %153 = sext i32 %152 to i64
  %154 = getelementptr inbounds double, ptr %149, i64 %153
  %155 = load double, ptr %154, align 8, !tbaa !5
  %156 = load ptr, ptr %8, align 8, !tbaa !9
  %157 = getelementptr inbounds %struct.reionization, ptr %156, i32 0, i32 25
  %158 = load ptr, ptr %157, align 8, !tbaa !15
  %159 = load ptr, ptr %8, align 8, !tbaa !9
  %160 = getelementptr inbounds %struct.reionization, ptr %159, i32 0, i32 15
  %161 = load i32, ptr %160, align 4, !tbaa !18
  %162 = sext i32 %161 to i64
  %163 = getelementptr inbounds double, ptr %158, i64 %162
  %164 = load double, ptr %163, align 8, !tbaa !5
  %165 = fsub double %155, %164
  %166 = load double, ptr %10, align 8, !tbaa !5
  %167 = call double @tanh(double noundef %166) #5
  %168 = fadd double %167, 1.000000e+00
  %169 = fmul double %165, %168
  %170 = fdiv double %169, 2.000000e+00
  %171 = load ptr, ptr %8, align 8, !tbaa !9
  %172 = getelementptr inbounds %struct.reionization, ptr %171, i32 0, i32 25
  %173 = load ptr, ptr %172, align 8, !tbaa !15
  %174 = load ptr, ptr %8, align 8, !tbaa !9
  %175 = getelementptr inbounds %struct.reionization, ptr %174, i32 0, i32 15
  %176 = load i32, ptr %175, align 4, !tbaa !18
  %177 = sext i32 %176 to i64
  %178 = getelementptr inbounds double, ptr %173, i64 %177
  %179 = load double, ptr %178, align 8, !tbaa !5
  %180 = fadd double %170, %179
  %181 = load ptr, ptr %9, align 8, !tbaa !9
  store double %180, ptr %181, align 8, !tbaa !5
  br label %215

182:                                              ; preds = %66
  %183 = load ptr, ptr %8, align 8, !tbaa !9
  %184 = getelementptr inbounds %struct.reionization, ptr %183, i32 0, i32 25
  %185 = load ptr, ptr %184, align 8, !tbaa !15
  %186 = load ptr, ptr %8, align 8, !tbaa !9
  %187 = getelementptr inbounds %struct.reionization, ptr %186, i32 0, i32 16
  %188 = load i32, ptr %187, align 8, !tbaa !22
  %189 = sext i32 %188 to i64
  %190 = getelementptr inbounds double, ptr %185, i64 %189
  %191 = load double, ptr %190, align 8, !tbaa !5
  %192 = load ptr, ptr %8, align 8, !tbaa !9
  %193 = getelementptr inbounds %struct.reionization, ptr %192, i32 0, i32 25
  %194 = load ptr, ptr %193, align 8, !tbaa !15
  %195 = load ptr, ptr %8, align 8, !tbaa !9
  %196 = getelementptr inbounds %struct.reionization, ptr %195, i32 0, i32 15
  %197 = load i32, ptr %196, align 4, !tbaa !18
  %198 = sext i32 %197 to i64
  %199 = getelementptr inbounds double, ptr %194, i64 %198
  %200 = load double, ptr %199, align 8, !tbaa !5
  %201 = fsub double %191, %200
  %202 = load double, ptr %10, align 8, !tbaa !5
  %203 = call double @tanh(double noundef %202) #5
  %204 = load ptr, ptr %8, align 8, !tbaa !9
  %205 = getelementptr inbounds %struct.reionization, ptr %204, i32 0, i32 25
  %206 = load ptr, ptr %205, align 8, !tbaa !15
  %207 = load ptr, ptr %8, align 8, !tbaa !9
  %208 = getelementptr inbounds %struct.reionization, ptr %207, i32 0, i32 15
  %209 = load i32, ptr %208, align 4, !tbaa !18
  %210 = sext i32 %209 to i64
  %211 = getelementptr inbounds double, ptr %206, i64 %210
  %212 = load double, ptr %211, align 8, !tbaa !5
  %213 = call double @llvm.fmuladd.f64(double %201, double %203, double %212)
  %214 = load ptr, ptr %9, align 8, !tbaa !9
  store double %213, ptr %214, align 8, !tbaa !5
  br label %215

215:                                              ; preds = %182, %146
  %216 = load ptr, ptr %7, align 8, !tbaa !9
  %217 = getelementptr inbounds %struct.thermo, ptr %216, i32 0, i32 2
  %218 = load i32, ptr %217, align 4, !tbaa !11
  %219 = icmp eq i32 %218, 1
  br i1 %219, label %220, label %259

220:                                              ; preds = %215
  %221 = load ptr, ptr %8, align 8, !tbaa !9
  %222 = getelementptr inbounds %struct.reionization, ptr %221, i32 0, i32 25
  %223 = load ptr, ptr %222, align 8, !tbaa !15
  %224 = load ptr, ptr %8, align 8, !tbaa !9
  %225 = getelementptr inbounds %struct.reionization, ptr %224, i32 0, i32 18
  %226 = load i32, ptr %225, align 8, !tbaa !23
  %227 = sext i32 %226 to i64
  %228 = getelementptr inbounds double, ptr %223, i64 %227
  %229 = load double, ptr %228, align 8, !tbaa !5
  %230 = load double, ptr %6, align 8, !tbaa !5
  %231 = fsub double %229, %230
  %232 = load ptr, ptr %8, align 8, !tbaa !9
  %233 = getelementptr inbounds %struct.reionization, ptr %232, i32 0, i32 25
  %234 = load ptr, ptr %233, align 8, !tbaa !15
  %235 = load ptr, ptr %8, align 8, !tbaa !9
  %236 = getelementptr inbounds %struct.reionization, ptr %235, i32 0, i32 19
  %237 = load i32, ptr %236, align 4, !tbaa !24
  %238 = sext i32 %237 to i64
  %239 = getelementptr inbounds double, ptr %234, i64 %238
  %240 = load double, ptr %239, align 8, !tbaa !5
  %241 = fdiv double %231, %240
  store double %241, ptr %10, align 8, !tbaa !5
  %242 = load ptr, ptr %8, align 8, !tbaa !9
  %243 = getelementptr inbounds %struct.reionization, ptr %242, i32 0, i32 25
  %244 = load ptr, ptr %243, align 8, !tbaa !15
  %245 = load ptr, ptr %8, align 8, !tbaa !9
  %246 = getelementptr inbounds %struct.reionization, ptr %245, i32 0, i32 17
  %247 = load i32, ptr %246, align 4, !tbaa !25
  %248 = sext i32 %247 to i64
  %249 = getelementptr inbounds double, ptr %244, i64 %248
  %250 = load double, ptr %249, align 8, !tbaa !5
  %251 = load double, ptr %10, align 8, !tbaa !5
  %252 = call double @tanh(double noundef %251) #5
  %253 = fadd double %252, 1.000000e+00
  %254 = fmul double %250, %253
  %255 = fdiv double %254, 2.000000e+00
  %256 = load ptr, ptr %9, align 8, !tbaa !9
  %257 = load double, ptr %256, align 8, !tbaa !5
  %258 = fadd double %257, %255
  store double %258, ptr %256, align 8, !tbaa !5
  br label %259

259:                                              ; preds = %220, %215
  br label %260

260:                                              ; preds = %259, %55
  store i32 0, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %876

261:                                              ; preds = %38
  %262 = load ptr, ptr %7, align 8, !tbaa !9
  %263 = getelementptr inbounds %struct.thermo, ptr %262, i32 0, i32 2
  %264 = load i32, ptr %263, align 4, !tbaa !11
  %265 = icmp eq i32 %264, 2
  br i1 %265, label %266, label %467

266:                                              ; preds = %261
  %267 = load double, ptr %6, align 8, !tbaa !5
  %268 = load ptr, ptr %8, align 8, !tbaa !9
  %269 = getelementptr inbounds %struct.reionization, ptr %268, i32 0, i32 25
  %270 = load ptr, ptr %269, align 8, !tbaa !15
  %271 = load ptr, ptr %8, align 8, !tbaa !9
  %272 = getelementptr inbounds %struct.reionization, ptr %271, i32 0, i32 21
  %273 = load i32, ptr %272, align 4, !tbaa !26
  %274 = load ptr, ptr %8, align 8, !tbaa !9
  %275 = getelementptr inbounds %struct.reionization, ptr %274, i32 0, i32 20
  %276 = load i32, ptr %275, align 8, !tbaa !27
  %277 = add nsw i32 %273, %276
  %278 = sub nsw i32 %277, 1
  %279 = sext i32 %278 to i64
  %280 = getelementptr inbounds double, ptr %270, i64 %279
  %281 = load double, ptr %280, align 8, !tbaa !5
  %282 = fcmp ogt double %267, %281
  br i1 %282, label %283, label %299

283:                                              ; preds = %266
  %284 = load ptr, ptr %8, align 8, !tbaa !9
  %285 = getelementptr inbounds %struct.reionization, ptr %284, i32 0, i32 25
  %286 = load ptr, ptr %285, align 8, !tbaa !15
  %287 = load ptr, ptr %8, align 8, !tbaa !9
  %288 = getelementptr inbounds %struct.reionization, ptr %287, i32 0, i32 22
  %289 = load i32, ptr %288, align 8, !tbaa !28
  %290 = load ptr, ptr %8, align 8, !tbaa !9
  %291 = getelementptr inbounds %struct.reionization, ptr %290, i32 0, i32 20
  %292 = load i32, ptr %291, align 8, !tbaa !27
  %293 = add nsw i32 %289, %292
  %294 = sub nsw i32 %293, 1
  %295 = sext i32 %294 to i64
  %296 = getelementptr inbounds double, ptr %286, i64 %295
  %297 = load double, ptr %296, align 8, !tbaa !5
  %298 = load ptr, ptr %9, align 8, !tbaa !9
  store double %297, ptr %298, align 8, !tbaa !5
  br label %466

299:                                              ; preds = %266
  %300 = load double, ptr %6, align 8, !tbaa !5
  %301 = load ptr, ptr %8, align 8, !tbaa !9
  %302 = getelementptr inbounds %struct.reionization, ptr %301, i32 0, i32 25
  %303 = load ptr, ptr %302, align 8, !tbaa !15
  %304 = load ptr, ptr %8, align 8, !tbaa !9
  %305 = getelementptr inbounds %struct.reionization, ptr %304, i32 0, i32 21
  %306 = load i32, ptr %305, align 4, !tbaa !26
  %307 = sext i32 %306 to i64
  %308 = getelementptr inbounds double, ptr %303, i64 %307
  %309 = load double, ptr %308, align 8, !tbaa !5
  %310 = fcmp olt double %300, %309
  br i1 %310, label %311, label %322

311:                                              ; preds = %299
  %312 = load ptr, ptr %8, align 8, !tbaa !9
  %313 = getelementptr inbounds %struct.reionization, ptr %312, i32 0, i32 25
  %314 = load ptr, ptr %313, align 8, !tbaa !15
  %315 = load ptr, ptr %8, align 8, !tbaa !9
  %316 = getelementptr inbounds %struct.reionization, ptr %315, i32 0, i32 22
  %317 = load i32, ptr %316, align 8, !tbaa !28
  %318 = sext i32 %317 to i64
  %319 = getelementptr inbounds double, ptr %314, i64 %318
  %320 = load double, ptr %319, align 8, !tbaa !5
  %321 = load ptr, ptr %9, align 8, !tbaa !9
  store double %320, ptr %321, align 8, !tbaa !5
  br label %465

322:                                              ; preds = %299
  store i32 0, ptr %11, align 4, !tbaa !29
  br label %323

323:                                              ; preds = %338, %322
  %324 = load ptr, ptr %8, align 8, !tbaa !9
  %325 = getelementptr inbounds %struct.reionization, ptr %324, i32 0, i32 25
  %326 = load ptr, ptr %325, align 8, !tbaa !15
  %327 = load ptr, ptr %8, align 8, !tbaa !9
  %328 = getelementptr inbounds %struct.reionization, ptr %327, i32 0, i32 21
  %329 = load i32, ptr %328, align 4, !tbaa !26
  %330 = load i32, ptr %11, align 4, !tbaa !29
  %331 = add nsw i32 %329, %330
  %332 = add nsw i32 %331, 1
  %333 = sext i32 %332 to i64
  %334 = getelementptr inbounds double, ptr %326, i64 %333
  %335 = load double, ptr %334, align 8, !tbaa !5
  %336 = load double, ptr %6, align 8, !tbaa !5
  %337 = fcmp olt double %335, %336
  br i1 %337, label %338, label %341

338:                                              ; preds = %323
  %339 = load i32, ptr %11, align 4, !tbaa !29
  %340 = add nsw i32 %339, 1
  store i32 %340, ptr %11, align 4, !tbaa !29
  br label %323, !llvm.loop !30

341:                                              ; preds = %323
  %342 = load i32, ptr %11, align 4, !tbaa !29
  %343 = load ptr, ptr %8, align 8, !tbaa !9
  %344 = getelementptr inbounds %struct.reionization, ptr %343, i32 0, i32 20
  %345 = load i32, ptr %344, align 8, !tbaa !27
  %346 = sub nsw i32 %345, 2
  %347 = icmp eq i32 %342, %346
  br i1 %347, label %348, label %385

348:                                              ; preds = %341
  %349 = load ptr, ptr %8, align 8, !tbaa !9
  %350 = getelementptr inbounds %struct.reionization, ptr %349, i32 0, i32 25
  %351 = load ptr, ptr %350, align 8, !tbaa !15
  %352 = load ptr, ptr %8, align 8, !tbaa !9
  %353 = getelementptr inbounds %struct.reionization, ptr %352, i32 0, i32 21
  %354 = load i32, ptr %353, align 4, !tbaa !26
  %355 = load i32, ptr %11, align 4, !tbaa !29
  %356 = add nsw i32 %354, %355
  %357 = sext i32 %356 to i64
  %358 = getelementptr inbounds double, ptr %351, i64 %357
  %359 = load double, ptr %358, align 8, !tbaa !5
  %360 = load ptr, ptr %8, align 8, !tbaa !9
  %361 = getelementptr inbounds %struct.reionization, ptr %360, i32 0, i32 25
  %362 = load ptr, ptr %361, align 8, !tbaa !15
  %363 = load ptr, ptr %8, align 8, !tbaa !9
  %364 = getelementptr inbounds %struct.reionization, ptr %363, i32 0, i32 21
  %365 = load i32, ptr %364, align 4, !tbaa !26
  %366 = load i32, ptr %11, align 4, !tbaa !29
  %367 = add nsw i32 %365, %366
  %368 = sext i32 %367 to i64
  %369 = getelementptr inbounds double, ptr %362, i64 %368
  %370 = load double, ptr %369, align 8, !tbaa !5
  %371 = load ptr, ptr %8, align 8, !tbaa !9
  %372 = getelementptr inbounds %struct.reionization, ptr %371, i32 0, i32 25
  %373 = load ptr, ptr %372, align 8, !tbaa !15
  %374 = load ptr, ptr %8, align 8, !tbaa !9
  %375 = getelementptr inbounds %struct.reionization, ptr %374, i32 0, i32 21
  %376 = load i32, ptr %375, align 4, !tbaa !26
  %377 = load i32, ptr %11, align 4, !tbaa !29
  %378 = add nsw i32 %376, %377
  %379 = sub nsw i32 %378, 1
  %380 = sext i32 %379 to i64
  %381 = getelementptr inbounds double, ptr %373, i64 %380
  %382 = load double, ptr %381, align 8, !tbaa !5
  %383 = fsub double %370, %382
  %384 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %383, double %359)
  store double %384, ptr %12, align 8, !tbaa !5
  br label %411

385:                                              ; preds = %341
  %386 = load ptr, ptr %8, align 8, !tbaa !9
  %387 = getelementptr inbounds %struct.reionization, ptr %386, i32 0, i32 25
  %388 = load ptr, ptr %387, align 8, !tbaa !15
  %389 = load ptr, ptr %8, align 8, !tbaa !9
  %390 = getelementptr inbounds %struct.reionization, ptr %389, i32 0, i32 21
  %391 = load i32, ptr %390, align 4, !tbaa !26
  %392 = load i32, ptr %11, align 4, !tbaa !29
  %393 = add nsw i32 %391, %392
  %394 = add nsw i32 %393, 1
  %395 = sext i32 %394 to i64
  %396 = getelementptr inbounds double, ptr %388, i64 %395
  %397 = load double, ptr %396, align 8, !tbaa !5
  %398 = load ptr, ptr %8, align 8, !tbaa !9
  %399 = getelementptr inbounds %struct.reionization, ptr %398, i32 0, i32 25
  %400 = load ptr, ptr %399, align 8, !tbaa !15
  %401 = load ptr, ptr %8, align 8, !tbaa !9
  %402 = getelementptr inbounds %struct.reionization, ptr %401, i32 0, i32 21
  %403 = load i32, ptr %402, align 4, !tbaa !26
  %404 = load i32, ptr %11, align 4, !tbaa !29
  %405 = add nsw i32 %403, %404
  %406 = sext i32 %405 to i64
  %407 = getelementptr inbounds double, ptr %400, i64 %406
  %408 = load double, ptr %407, align 8, !tbaa !5
  %409 = fadd double %397, %408
  %410 = fmul double 5.000000e-01, %409
  store double %410, ptr %12, align 8, !tbaa !5
  br label %411

411:                                              ; preds = %385, %348
  %412 = load ptr, ptr %8, align 8, !tbaa !9
  %413 = getelementptr inbounds %struct.reionization, ptr %412, i32 0, i32 25
  %414 = load ptr, ptr %413, align 8, !tbaa !15
  %415 = load ptr, ptr %8, align 8, !tbaa !9
  %416 = getelementptr inbounds %struct.reionization, ptr %415, i32 0, i32 22
  %417 = load i32, ptr %416, align 8, !tbaa !28
  %418 = load i32, ptr %11, align 4, !tbaa !29
  %419 = add nsw i32 %417, %418
  %420 = sext i32 %419 to i64
  %421 = getelementptr inbounds double, ptr %414, i64 %420
  %422 = load double, ptr %421, align 8, !tbaa !5
  %423 = load double, ptr %6, align 8, !tbaa !5
  %424 = load double, ptr %12, align 8, !tbaa !5
  %425 = fsub double %423, %424
  %426 = load ptr, ptr %8, align 8, !tbaa !9
  %427 = getelementptr inbounds %struct.reionization, ptr %426, i32 0, i32 25
  %428 = load ptr, ptr %427, align 8, !tbaa !15
  %429 = load ptr, ptr %8, align 8, !tbaa !9
  %430 = getelementptr inbounds %struct.reionization, ptr %429, i32 0, i32 23
  %431 = load i32, ptr %430, align 4, !tbaa !32
  %432 = sext i32 %431 to i64
  %433 = getelementptr inbounds double, ptr %428, i64 %432
  %434 = load double, ptr %433, align 8, !tbaa !5
  %435 = fdiv double %425, %434
  %436 = call double @tanh(double noundef %435) #5
  %437 = fadd double %436, 1.000000e+00
  %438 = fmul double 5.000000e-01, %437
  %439 = load ptr, ptr %8, align 8, !tbaa !9
  %440 = getelementptr inbounds %struct.reionization, ptr %439, i32 0, i32 25
  %441 = load ptr, ptr %440, align 8, !tbaa !15
  %442 = load ptr, ptr %8, align 8, !tbaa !9
  %443 = getelementptr inbounds %struct.reionization, ptr %442, i32 0, i32 22
  %444 = load i32, ptr %443, align 8, !tbaa !28
  %445 = load i32, ptr %11, align 4, !tbaa !29
  %446 = add nsw i32 %444, %445
  %447 = add nsw i32 %446, 1
  %448 = sext i32 %447 to i64
  %449 = getelementptr inbounds double, ptr %441, i64 %448
  %450 = load double, ptr %449, align 8, !tbaa !5
  %451 = load ptr, ptr %8, align 8, !tbaa !9
  %452 = getelementptr inbounds %struct.reionization, ptr %451, i32 0, i32 25
  %453 = load ptr, ptr %452, align 8, !tbaa !15
  %454 = load ptr, ptr %8, align 8, !tbaa !9
  %455 = getelementptr inbounds %struct.reionization, ptr %454, i32 0, i32 22
  %456 = load i32, ptr %455, align 8, !tbaa !28
  %457 = load i32, ptr %11, align 4, !tbaa !29
  %458 = add nsw i32 %456, %457
  %459 = sext i32 %458 to i64
  %460 = getelementptr inbounds double, ptr %453, i64 %459
  %461 = load double, ptr %460, align 8, !tbaa !5
  %462 = fsub double %450, %461
  %463 = call double @llvm.fmuladd.f64(double %438, double %462, double %422)
  %464 = load ptr, ptr %9, align 8, !tbaa !9
  store double %463, ptr %464, align 8, !tbaa !5
  br label %465

465:                                              ; preds = %411, %311
  br label %466

466:                                              ; preds = %465, %283
  store i32 0, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %876

467:                                              ; preds = %261
  %468 = load ptr, ptr %7, align 8, !tbaa !9
  %469 = getelementptr inbounds %struct.thermo, ptr %468, i32 0, i32 2
  %470 = load i32, ptr %469, align 4, !tbaa !11
  %471 = icmp eq i32 %470, 4
  br i1 %471, label %472, label %636

472:                                              ; preds = %467
  %473 = load double, ptr %6, align 8, !tbaa !5
  %474 = load ptr, ptr %8, align 8, !tbaa !9
  %475 = getelementptr inbounds %struct.reionization, ptr %474, i32 0, i32 25
  %476 = load ptr, ptr %475, align 8, !tbaa !15
  %477 = load ptr, ptr %8, align 8, !tbaa !9
  %478 = getelementptr inbounds %struct.reionization, ptr %477, i32 0, i32 21
  %479 = load i32, ptr %478, align 4, !tbaa !26
  %480 = load ptr, ptr %8, align 8, !tbaa !9
  %481 = getelementptr inbounds %struct.reionization, ptr %480, i32 0, i32 20
  %482 = load i32, ptr %481, align 8, !tbaa !27
  %483 = add nsw i32 %479, %482
  %484 = sub nsw i32 %483, 1
  %485 = sext i32 %484 to i64
  %486 = getelementptr inbounds double, ptr %476, i64 %485
  %487 = load double, ptr %486, align 8, !tbaa !5
  %488 = fcmp ogt double %473, %487
  br i1 %488, label %489, label %505

489:                                              ; preds = %472
  %490 = load ptr, ptr %8, align 8, !tbaa !9
  %491 = getelementptr inbounds %struct.reionization, ptr %490, i32 0, i32 25
  %492 = load ptr, ptr %491, align 8, !tbaa !15
  %493 = load ptr, ptr %8, align 8, !tbaa !9
  %494 = getelementptr inbounds %struct.reionization, ptr %493, i32 0, i32 22
  %495 = load i32, ptr %494, align 8, !tbaa !28
  %496 = load ptr, ptr %8, align 8, !tbaa !9
  %497 = getelementptr inbounds %struct.reionization, ptr %496, i32 0, i32 20
  %498 = load i32, ptr %497, align 8, !tbaa !27
  %499 = add nsw i32 %495, %498
  %500 = sub nsw i32 %499, 1
  %501 = sext i32 %500 to i64
  %502 = getelementptr inbounds double, ptr %492, i64 %501
  %503 = load double, ptr %502, align 8, !tbaa !5
  %504 = load ptr, ptr %9, align 8, !tbaa !9
  store double %503, ptr %504, align 8, !tbaa !5
  br label %635

505:                                              ; preds = %472
  %506 = load double, ptr %6, align 8, !tbaa !5
  %507 = load ptr, ptr %8, align 8, !tbaa !9
  %508 = getelementptr inbounds %struct.reionization, ptr %507, i32 0, i32 25
  %509 = load ptr, ptr %508, align 8, !tbaa !15
  %510 = load ptr, ptr %8, align 8, !tbaa !9
  %511 = getelementptr inbounds %struct.reionization, ptr %510, i32 0, i32 21
  %512 = load i32, ptr %511, align 4, !tbaa !26
  %513 = sext i32 %512 to i64
  %514 = getelementptr inbounds double, ptr %509, i64 %513
  %515 = load double, ptr %514, align 8, !tbaa !5
  %516 = fcmp ogt double %506, %515
  br i1 %516, label %517, label %623

517:                                              ; preds = %505
  %518 = load ptr, ptr %8, align 8, !tbaa !9
  %519 = getelementptr inbounds %struct.reionization, ptr %518, i32 0, i32 25
  %520 = load ptr, ptr %519, align 8, !tbaa !15
  %521 = load ptr, ptr %8, align 8, !tbaa !9
  %522 = getelementptr inbounds %struct.reionization, ptr %521, i32 0, i32 22
  %523 = load i32, ptr %522, align 8, !tbaa !28
  %524 = load ptr, ptr %8, align 8, !tbaa !9
  %525 = getelementptr inbounds %struct.reionization, ptr %524, i32 0, i32 20
  %526 = load i32, ptr %525, align 8, !tbaa !27
  %527 = add nsw i32 %523, %526
  %528 = sub nsw i32 %527, 1
  %529 = sext i32 %528 to i64
  %530 = getelementptr inbounds double, ptr %520, i64 %529
  %531 = load double, ptr %530, align 8, !tbaa !5
  %532 = load ptr, ptr %9, align 8, !tbaa !9
  store double %531, ptr %532, align 8, !tbaa !5
  store i32 1, ptr %13, align 4, !tbaa !29
  br label %533

533:                                              ; preds = %619, %517
  %534 = load i32, ptr %13, align 4, !tbaa !29
  %535 = load ptr, ptr %8, align 8, !tbaa !9
  %536 = getelementptr inbounds %struct.reionization, ptr %535, i32 0, i32 20
  %537 = load i32, ptr %536, align 8, !tbaa !27
  %538 = sub nsw i32 %537, 1
  %539 = icmp slt i32 %534, %538
  br i1 %539, label %540, label %622

540:                                              ; preds = %533
  %541 = load ptr, ptr %8, align 8, !tbaa !9
  %542 = getelementptr inbounds %struct.reionization, ptr %541, i32 0, i32 25
  %543 = load ptr, ptr %542, align 8, !tbaa !15
  %544 = load ptr, ptr %8, align 8, !tbaa !9
  %545 = getelementptr inbounds %struct.reionization, ptr %544, i32 0, i32 21
  %546 = load i32, ptr %545, align 4, !tbaa !26
  %547 = load ptr, ptr %8, align 8, !tbaa !9
  %548 = getelementptr inbounds %struct.reionization, ptr %547, i32 0, i32 20
  %549 = load i32, ptr %548, align 8, !tbaa !27
  %550 = add nsw i32 %546, %549
  %551 = sub nsw i32 %550, 1
  %552 = load i32, ptr %13, align 4, !tbaa !29
  %553 = sub nsw i32 %551, %552
  %554 = sext i32 %553 to i64
  %555 = getelementptr inbounds double, ptr %543, i64 %554
  %556 = load double, ptr %555, align 8, !tbaa !5
  store double %556, ptr %14, align 8, !tbaa !5
  %557 = load ptr, ptr %8, align 8, !tbaa !9
  %558 = getelementptr inbounds %struct.reionization, ptr %557, i32 0, i32 25
  %559 = load ptr, ptr %558, align 8, !tbaa !15
  %560 = load ptr, ptr %8, align 8, !tbaa !9
  %561 = getelementptr inbounds %struct.reionization, ptr %560, i32 0, i32 22
  %562 = load i32, ptr %561, align 8, !tbaa !28
  %563 = load ptr, ptr %8, align 8, !tbaa !9
  %564 = getelementptr inbounds %struct.reionization, ptr %563, i32 0, i32 20
  %565 = load i32, ptr %564, align 8, !tbaa !27
  %566 = add nsw i32 %562, %565
  %567 = sub nsw i32 %566, 1
  %568 = load i32, ptr %13, align 4, !tbaa !29
  %569 = sub nsw i32 %567, %568
  %570 = sext i32 %569 to i64
  %571 = getelementptr inbounds double, ptr %559, i64 %570
  %572 = load double, ptr %571, align 8, !tbaa !5
  %573 = load ptr, ptr %8, align 8, !tbaa !9
  %574 = getelementptr inbounds %struct.reionization, ptr %573, i32 0, i32 25
  %575 = load ptr, ptr %574, align 8, !tbaa !15
  %576 = load ptr, ptr %8, align 8, !tbaa !9
  %577 = getelementptr inbounds %struct.reionization, ptr %576, i32 0, i32 22
  %578 = load i32, ptr %577, align 8, !tbaa !28
  %579 = load ptr, ptr %8, align 8, !tbaa !9
  %580 = getelementptr inbounds %struct.reionization, ptr %579, i32 0, i32 20
  %581 = load i32, ptr %580, align 8, !tbaa !27
  %582 = add nsw i32 %578, %581
  %583 = load i32, ptr %13, align 4, !tbaa !29
  %584 = sub nsw i32 %582, %583
  %585 = sext i32 %584 to i64
  %586 = getelementptr inbounds double, ptr %575, i64 %585
  %587 = load double, ptr %586, align 8, !tbaa !5
  %588 = fsub double %572, %587
  store double %588, ptr %15, align 8, !tbaa !5
  store double 0.000000e+00, ptr %16, align 8, !tbaa !5
  %589 = load ptr, ptr %8, align 8, !tbaa !9
  %590 = getelementptr inbounds %struct.reionization, ptr %589, i32 0, i32 25
  %591 = load ptr, ptr %590, align 8, !tbaa !15
  %592 = load ptr, ptr %8, align 8, !tbaa !9
  %593 = getelementptr inbounds %struct.reionization, ptr %592, i32 0, i32 23
  %594 = load i32, ptr %593, align 4, !tbaa !32
  %595 = sext i32 %594 to i64
  %596 = getelementptr inbounds double, ptr %591, i64 %595
  %597 = load double, ptr %596, align 8, !tbaa !5
  store double %597, ptr %17, align 8, !tbaa !5
  %598 = load double, ptr %6, align 8, !tbaa !5
  %599 = load double, ptr %14, align 8, !tbaa !5
  %600 = load double, ptr %15, align 8, !tbaa !5
  %601 = load double, ptr %16, align 8, !tbaa !5
  %602 = load double, ptr %17, align 8, !tbaa !5
  %603 = call i32 @thermodynamics_tanh(double noundef %598, double noundef %599, double noundef %600, double noundef %601, double noundef %602, ptr noundef %18)
  %604 = icmp eq i32 %603, 1
  br i1 %604, label %605, label %614

605:                                              ; preds = %540
  call void @llvm.lifetime.start.p0(i64 2048, ptr %20) #5
  %606 = getelementptr inbounds [2048 x i8], ptr %20, i64 0, i64 0
  %607 = load ptr, ptr %7, align 8, !tbaa !9
  %608 = getelementptr inbounds %struct.thermo, ptr %607, i32 0, i32 93
  %609 = getelementptr inbounds [2048 x i8], ptr %608, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %606, ptr noundef @.str, ptr noundef @.str.137, ptr noundef %609)
  %610 = load ptr, ptr %7, align 8, !tbaa !9
  %611 = getelementptr inbounds %struct.thermo, ptr %610, i32 0, i32 93
  %612 = getelementptr inbounds [2048 x i8], ptr %611, i64 0, i64 0
  %613 = getelementptr inbounds [2048 x i8], ptr %20, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %612, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_function, i32 noundef 2089, ptr noundef %613)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %20) #5
  store i32 1, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %876

614:                                              ; preds = %540
  %615 = load double, ptr %18, align 8, !tbaa !5
  %616 = load ptr, ptr %9, align 8, !tbaa !9
  %617 = load double, ptr %616, align 8, !tbaa !5
  %618 = fadd double %617, %615
  store double %618, ptr %616, align 8, !tbaa !5
  br label %619

619:                                              ; preds = %614
  %620 = load i32, ptr %13, align 4, !tbaa !29
  %621 = add nsw i32 %620, 1
  store i32 %621, ptr %13, align 4, !tbaa !29
  br label %533, !llvm.loop !33

622:                                              ; preds = %533
  br label %634

623:                                              ; preds = %505
  %624 = load ptr, ptr %8, align 8, !tbaa !9
  %625 = getelementptr inbounds %struct.reionization, ptr %624, i32 0, i32 25
  %626 = load ptr, ptr %625, align 8, !tbaa !15
  %627 = load ptr, ptr %8, align 8, !tbaa !9
  %628 = getelementptr inbounds %struct.reionization, ptr %627, i32 0, i32 22
  %629 = load i32, ptr %628, align 8, !tbaa !28
  %630 = sext i32 %629 to i64
  %631 = getelementptr inbounds double, ptr %626, i64 %630
  %632 = load double, ptr %631, align 8, !tbaa !5
  %633 = load ptr, ptr %9, align 8, !tbaa !9
  store double %632, ptr %633, align 8, !tbaa !5
  br label %634

634:                                              ; preds = %623, %622
  br label %635

635:                                              ; preds = %634, %489
  store i32 0, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %876

636:                                              ; preds = %467
  %637 = load ptr, ptr %7, align 8, !tbaa !9
  %638 = getelementptr inbounds %struct.thermo, ptr %637, i32 0, i32 2
  %639 = load i32, ptr %638, align 4, !tbaa !11
  %640 = icmp eq i32 %639, 5
  br i1 %640, label %641, label %865

641:                                              ; preds = %636
  %642 = load double, ptr %6, align 8, !tbaa !5
  %643 = load ptr, ptr %8, align 8, !tbaa !9
  %644 = getelementptr inbounds %struct.reionization, ptr %643, i32 0, i32 25
  %645 = load ptr, ptr %644, align 8, !tbaa !15
  %646 = load ptr, ptr %8, align 8, !tbaa !9
  %647 = getelementptr inbounds %struct.reionization, ptr %646, i32 0, i32 21
  %648 = load i32, ptr %647, align 4, !tbaa !26
  %649 = load ptr, ptr %8, align 8, !tbaa !9
  %650 = getelementptr inbounds %struct.reionization, ptr %649, i32 0, i32 20
  %651 = load i32, ptr %650, align 8, !tbaa !27
  %652 = add nsw i32 %648, %651
  %653 = sub nsw i32 %652, 1
  %654 = sext i32 %653 to i64
  %655 = getelementptr inbounds double, ptr %645, i64 %654
  %656 = load double, ptr %655, align 8, !tbaa !5
  %657 = fcmp ogt double %642, %656
  br i1 %657, label %658, label %681

658:                                              ; preds = %641
  %659 = load ptr, ptr %8, align 8, !tbaa !9
  %660 = getelementptr inbounds %struct.reionization, ptr %659, i32 0, i32 25
  %661 = load ptr, ptr %660, align 8, !tbaa !15
  %662 = load ptr, ptr %8, align 8, !tbaa !9
  %663 = getelementptr inbounds %struct.reionization, ptr %662, i32 0, i32 22
  %664 = load i32, ptr %663, align 8, !tbaa !28
  %665 = load ptr, ptr %8, align 8, !tbaa !9
  %666 = getelementptr inbounds %struct.reionization, ptr %665, i32 0, i32 20
  %667 = load i32, ptr %666, align 8, !tbaa !27
  %668 = add nsw i32 %664, %667
  %669 = sub nsw i32 %668, 1
  %670 = sext i32 %669 to i64
  %671 = getelementptr inbounds double, ptr %661, i64 %670
  %672 = load double, ptr %671, align 8, !tbaa !5
  %673 = load ptr, ptr %9, align 8, !tbaa !9
  store double %672, ptr %673, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 2048, ptr %21) #5
  %674 = getelementptr inbounds [2048 x i8], ptr %21, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %674, ptr noundef @.str.138)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %22) #5
  %675 = getelementptr inbounds [2048 x i8], ptr %22, i64 0, i64 0
  %676 = getelementptr inbounds [2048 x i8], ptr %21, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %675, ptr noundef @.str.139, ptr noundef %676)
  %677 = load ptr, ptr %7, align 8, !tbaa !9
  %678 = getelementptr inbounds %struct.thermo, ptr %677, i32 0, i32 93
  %679 = getelementptr inbounds [2048 x i8], ptr %678, i64 0, i64 0
  %680 = getelementptr inbounds [2048 x i8], ptr %22, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %679, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_function, i32 noundef 2113, ptr noundef %680)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %22) #5
  store i32 1, ptr %5, align 4
  store i32 1, ptr %19, align 4
  call void @llvm.lifetime.end.p0(i64 2048, ptr %21) #5
  br label %876

681:                                              ; preds = %641
  store i32 0, ptr %11, align 4, !tbaa !29
  br label %682

682:                                              ; preds = %697, %681
  %683 = load ptr, ptr %8, align 8, !tbaa !9
  %684 = getelementptr inbounds %struct.reionization, ptr %683, i32 0, i32 25
  %685 = load ptr, ptr %684, align 8, !tbaa !15
  %686 = load ptr, ptr %8, align 8, !tbaa !9
  %687 = getelementptr inbounds %struct.reionization, ptr %686, i32 0, i32 21
  %688 = load i32, ptr %687, align 4, !tbaa !26
  %689 = load i32, ptr %11, align 4, !tbaa !29
  %690 = add nsw i32 %688, %689
  %691 = add nsw i32 %690, 1
  %692 = sext i32 %691 to i64
  %693 = getelementptr inbounds double, ptr %685, i64 %692
  %694 = load double, ptr %693, align 8, !tbaa !5
  %695 = load double, ptr %6, align 8, !tbaa !5
  %696 = fcmp olt double %694, %695
  br i1 %696, label %697, label %700

697:                                              ; preds = %682
  %698 = load i32, ptr %11, align 4, !tbaa !29
  %699 = add nsw i32 %698, 1
  store i32 %699, ptr %11, align 4, !tbaa !29
  br label %682, !llvm.loop !34

700:                                              ; preds = %682
  call void @llvm.lifetime.start.p0(i64 8, ptr %23) #5
  %701 = load ptr, ptr %8, align 8, !tbaa !9
  %702 = getelementptr inbounds %struct.reionization, ptr %701, i32 0, i32 25
  %703 = load ptr, ptr %702, align 8, !tbaa !15
  %704 = load ptr, ptr %8, align 8, !tbaa !9
  %705 = getelementptr inbounds %struct.reionization, ptr %704, i32 0, i32 21
  %706 = load i32, ptr %705, align 4, !tbaa !26
  %707 = load i32, ptr %11, align 4, !tbaa !29
  %708 = add nsw i32 %706, %707
  %709 = sext i32 %708 to i64
  %710 = getelementptr inbounds double, ptr %703, i64 %709
  %711 = load double, ptr %710, align 8, !tbaa !5
  store double %711, ptr %23, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %24) #5
  %712 = load ptr, ptr %8, align 8, !tbaa !9
  %713 = getelementptr inbounds %struct.reionization, ptr %712, i32 0, i32 25
  %714 = load ptr, ptr %713, align 8, !tbaa !15
  %715 = load ptr, ptr %8, align 8, !tbaa !9
  %716 = getelementptr inbounds %struct.reionization, ptr %715, i32 0, i32 21
  %717 = load i32, ptr %716, align 4, !tbaa !26
  %718 = load i32, ptr %11, align 4, !tbaa !29
  %719 = add nsw i32 %717, %718
  %720 = add nsw i32 %719, 1
  %721 = sext i32 %720 to i64
  %722 = getelementptr inbounds double, ptr %714, i64 %721
  %723 = load double, ptr %722, align 8, !tbaa !5
  store double %723, ptr %24, align 8, !tbaa !5
  %724 = load double, ptr %6, align 8, !tbaa !5
  %725 = load double, ptr %23, align 8, !tbaa !5
  %726 = fcmp olt double %724, %725
  br i1 %726, label %727, label %735

727:                                              ; preds = %700
  call void @llvm.lifetime.start.p0(i64 2048, ptr %25) #5
  %728 = getelementptr inbounds [2048 x i8], ptr %25, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %728, ptr noundef @.str.140)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %26) #5
  %729 = getelementptr inbounds [2048 x i8], ptr %26, i64 0, i64 0
  %730 = getelementptr inbounds [2048 x i8], ptr %25, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %729, ptr noundef @.str.10, ptr noundef @.str.141, ptr noundef %730)
  %731 = load ptr, ptr %7, align 8, !tbaa !9
  %732 = getelementptr inbounds %struct.thermo, ptr %731, i32 0, i32 93
  %733 = getelementptr inbounds [2048 x i8], ptr %732, i64 0, i64 0
  %734 = getelementptr inbounds [2048 x i8], ptr %26, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %733, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_function, i32 noundef 2126, ptr noundef %734)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %26) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr %25) #5
  store i32 1, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %861

735:                                              ; preds = %700
  %736 = load double, ptr %6, align 8, !tbaa !5
  %737 = load double, ptr %24, align 8, !tbaa !5
  %738 = fcmp ogt double %736, %737
  br i1 %738, label %739, label %747

739:                                              ; preds = %735
  call void @llvm.lifetime.start.p0(i64 2048, ptr %27) #5
  %740 = getelementptr inbounds [2048 x i8], ptr %27, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %740, ptr noundef @.str.140)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %28) #5
  %741 = getelementptr inbounds [2048 x i8], ptr %28, i64 0, i64 0
  %742 = getelementptr inbounds [2048 x i8], ptr %27, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %741, ptr noundef @.str.10, ptr noundef @.str.142, ptr noundef %742)
  %743 = load ptr, ptr %7, align 8, !tbaa !9
  %744 = getelementptr inbounds %struct.thermo, ptr %743, i32 0, i32 93
  %745 = getelementptr inbounds [2048 x i8], ptr %744, i64 0, i64 0
  %746 = getelementptr inbounds [2048 x i8], ptr %28, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %745, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_function, i32 noundef 2130, ptr noundef %746)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %28) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr %27) #5
  store i32 1, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %861

747:                                              ; preds = %735
  call void @llvm.lifetime.start.p0(i64 8, ptr %29) #5
  %748 = load double, ptr %6, align 8, !tbaa !5
  %749 = load ptr, ptr %8, align 8, !tbaa !9
  %750 = getelementptr inbounds %struct.reionization, ptr %749, i32 0, i32 25
  %751 = load ptr, ptr %750, align 8, !tbaa !15
  %752 = load ptr, ptr %8, align 8, !tbaa !9
  %753 = getelementptr inbounds %struct.reionization, ptr %752, i32 0, i32 21
  %754 = load i32, ptr %753, align 4, !tbaa !26
  %755 = load i32, ptr %11, align 4, !tbaa !29
  %756 = add nsw i32 %754, %755
  %757 = sext i32 %756 to i64
  %758 = getelementptr inbounds double, ptr %751, i64 %757
  %759 = load double, ptr %758, align 8, !tbaa !5
  %760 = fsub double %748, %759
  %761 = load ptr, ptr %8, align 8, !tbaa !9
  %762 = getelementptr inbounds %struct.reionization, ptr %761, i32 0, i32 25
  %763 = load ptr, ptr %762, align 8, !tbaa !15
  %764 = load ptr, ptr %8, align 8, !tbaa !9
  %765 = getelementptr inbounds %struct.reionization, ptr %764, i32 0, i32 21
  %766 = load i32, ptr %765, align 4, !tbaa !26
  %767 = load i32, ptr %11, align 4, !tbaa !29
  %768 = add nsw i32 %766, %767
  %769 = add nsw i32 %768, 1
  %770 = sext i32 %769 to i64
  %771 = getelementptr inbounds double, ptr %763, i64 %770
  %772 = load double, ptr %771, align 8, !tbaa !5
  %773 = load ptr, ptr %8, align 8, !tbaa !9
  %774 = getelementptr inbounds %struct.reionization, ptr %773, i32 0, i32 25
  %775 = load ptr, ptr %774, align 8, !tbaa !15
  %776 = load ptr, ptr %8, align 8, !tbaa !9
  %777 = getelementptr inbounds %struct.reionization, ptr %776, i32 0, i32 21
  %778 = load i32, ptr %777, align 4, !tbaa !26
  %779 = load i32, ptr %11, align 4, !tbaa !29
  %780 = add nsw i32 %778, %779
  %781 = sext i32 %780 to i64
  %782 = getelementptr inbounds double, ptr %775, i64 %781
  %783 = load double, ptr %782, align 8, !tbaa !5
  %784 = fsub double %772, %783
  %785 = fdiv double %760, %784
  store double %785, ptr %29, align 8, !tbaa !5
  %786 = load ptr, ptr %8, align 8, !tbaa !9
  %787 = getelementptr inbounds %struct.reionization, ptr %786, i32 0, i32 25
  %788 = load ptr, ptr %787, align 8, !tbaa !15
  %789 = load ptr, ptr %8, align 8, !tbaa !9
  %790 = getelementptr inbounds %struct.reionization, ptr %789, i32 0, i32 22
  %791 = load i32, ptr %790, align 8, !tbaa !28
  %792 = load i32, ptr %11, align 4, !tbaa !29
  %793 = add nsw i32 %791, %792
  %794 = sext i32 %793 to i64
  %795 = getelementptr inbounds double, ptr %788, i64 %794
  %796 = load double, ptr %795, align 8, !tbaa !5
  %797 = load double, ptr %29, align 8, !tbaa !5
  %798 = load ptr, ptr %8, align 8, !tbaa !9
  %799 = getelementptr inbounds %struct.reionization, ptr %798, i32 0, i32 25
  %800 = load ptr, ptr %799, align 8, !tbaa !15
  %801 = load ptr, ptr %8, align 8, !tbaa !9
  %802 = getelementptr inbounds %struct.reionization, ptr %801, i32 0, i32 22
  %803 = load i32, ptr %802, align 8, !tbaa !28
  %804 = load i32, ptr %11, align 4, !tbaa !29
  %805 = add nsw i32 %803, %804
  %806 = add nsw i32 %805, 1
  %807 = sext i32 %806 to i64
  %808 = getelementptr inbounds double, ptr %800, i64 %807
  %809 = load double, ptr %808, align 8, !tbaa !5
  %810 = load ptr, ptr %8, align 8, !tbaa !9
  %811 = getelementptr inbounds %struct.reionization, ptr %810, i32 0, i32 25
  %812 = load ptr, ptr %811, align 8, !tbaa !15
  %813 = load ptr, ptr %8, align 8, !tbaa !9
  %814 = getelementptr inbounds %struct.reionization, ptr %813, i32 0, i32 22
  %815 = load i32, ptr %814, align 8, !tbaa !28
  %816 = load i32, ptr %11, align 4, !tbaa !29
  %817 = add nsw i32 %815, %816
  %818 = sext i32 %817 to i64
  %819 = getelementptr inbounds double, ptr %812, i64 %818
  %820 = load double, ptr %819, align 8, !tbaa !5
  %821 = fsub double %809, %820
  %822 = call double @llvm.fmuladd.f64(double %797, double %821, double %796)
  %823 = load ptr, ptr %9, align 8, !tbaa !9
  store double %822, ptr %823, align 8, !tbaa !5
  %824 = load ptr, ptr %9, align 8, !tbaa !9
  %825 = load double, ptr %824, align 8, !tbaa !5
  %826 = fcmp olt double %825, 0.000000e+00
  br i1 %826, label %827, label %859

827:                                              ; preds = %747
  call void @llvm.lifetime.start.p0(i64 2048, ptr %30) #5
  %828 = getelementptr inbounds [2048 x i8], ptr %30, i64 0, i64 0
  %829 = load double, ptr %29, align 8, !tbaa !5
  %830 = load ptr, ptr %8, align 8, !tbaa !9
  %831 = getelementptr inbounds %struct.reionization, ptr %830, i32 0, i32 25
  %832 = load ptr, ptr %831, align 8, !tbaa !15
  %833 = load ptr, ptr %8, align 8, !tbaa !9
  %834 = getelementptr inbounds %struct.reionization, ptr %833, i32 0, i32 22
  %835 = load i32, ptr %834, align 8, !tbaa !28
  %836 = load i32, ptr %11, align 4, !tbaa !29
  %837 = add nsw i32 %835, %836
  %838 = sext i32 %837 to i64
  %839 = getelementptr inbounds double, ptr %832, i64 %838
  %840 = load double, ptr %839, align 8, !tbaa !5
  %841 = load ptr, ptr %8, align 8, !tbaa !9
  %842 = getelementptr inbounds %struct.reionization, ptr %841, i32 0, i32 25
  %843 = load ptr, ptr %842, align 8, !tbaa !15
  %844 = load ptr, ptr %8, align 8, !tbaa !9
  %845 = getelementptr inbounds %struct.reionization, ptr %844, i32 0, i32 22
  %846 = load i32, ptr %845, align 8, !tbaa !28
  %847 = load i32, ptr %11, align 4, !tbaa !29
  %848 = add nsw i32 %846, %847
  %849 = add nsw i32 %848, 1
  %850 = sext i32 %849 to i64
  %851 = getelementptr inbounds double, ptr %843, i64 %850
  %852 = load double, ptr %851, align 8, !tbaa !5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %828, ptr noundef @.str.143, double noundef %829, double noundef %840, double noundef %852)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %31) #5
  %853 = getelementptr inbounds [2048 x i8], ptr %31, i64 0, i64 0
  %854 = getelementptr inbounds [2048 x i8], ptr %30, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %853, ptr noundef @.str.10, ptr noundef @.str.144, ptr noundef %854)
  %855 = load ptr, ptr %7, align 8, !tbaa !9
  %856 = getelementptr inbounds %struct.thermo, ptr %855, i32 0, i32 93
  %857 = getelementptr inbounds [2048 x i8], ptr %856, i64 0, i64 0
  %858 = getelementptr inbounds [2048 x i8], ptr %31, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %857, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_function, i32 noundef 2145, ptr noundef %858)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %31) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr %30) #5
  store i32 1, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %860

859:                                              ; preds = %747
  store i32 0, ptr %19, align 4
  br label %860

860:                                              ; preds = %859, %827
  call void @llvm.lifetime.end.p0(i64 8, ptr %29) #5
  br label %861

861:                                              ; preds = %860, %739, %727
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #5
  %862 = load i32, ptr %19, align 4
  switch i32 %862, label %876 [
    i32 0, label %863
  ]

863:                                              ; preds = %861
  br label %864

864:                                              ; preds = %863
  store i32 0, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %876

865:                                              ; preds = %636
  call void @llvm.lifetime.start.p0(i64 2048, ptr %32) #5
  %866 = getelementptr inbounds [2048 x i8], ptr %32, i64 0, i64 0
  %867 = load ptr, ptr %7, align 8, !tbaa !9
  %868 = getelementptr inbounds %struct.thermo, ptr %867, i32 0, i32 2
  %869 = load i32, ptr %868, align 4, !tbaa !11
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %866, ptr noundef @.str.145, i32 noundef %869)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %33) #5
  %870 = getelementptr inbounds [2048 x i8], ptr %33, i64 0, i64 0
  %871 = getelementptr inbounds [2048 x i8], ptr %32, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %870, ptr noundef @.str.10, ptr noundef @.str.146, ptr noundef %871)
  %872 = load ptr, ptr %7, align 8, !tbaa !9
  %873 = getelementptr inbounds %struct.thermo, ptr %872, i32 0, i32 93
  %874 = getelementptr inbounds [2048 x i8], ptr %873, i64 0, i64 0
  %875 = getelementptr inbounds [2048 x i8], ptr %33, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %874, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_reionization_function, i32 noundef 2155, ptr noundef %875)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %33) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr %32) #5
  store i32 1, ptr %5, align 4
  store i32 1, ptr %19, align 4
  br label %876

876:                                              ; preds = %865, %864, %861, %658, %635, %605, %466, %260
  call void @llvm.lifetime.end.p0(i64 8, ptr %18) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %17) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %14) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %13) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %11) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %10) #5
  %877 = load i32, ptr %5, align 4
  ret i32 %877
}

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

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

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 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 #5 = { nounwind }

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"double", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"any pointer", !7, i64 0}
!11 = !{!12, !13, i64 12}
!12 = !{!"thermo", !6, i64 0, !13, i64 8, !13, i64 12, !13, i64 16, !6, i64 24, !6, i64 32, !14, i64 40, !14, i64 42, !6, i64 48, !6, i64 56, !6, i64 64, !6, i64 72, !13, i64 80, !10, i64 88, !10, i64 96, !6, i64 104, !13, i64 112, !10, i64 120, !10, i64 128, !6, i64 136, !13, i64 144, !10, i64 152, !10, i64 160, !6, i64 168, !14, i64 176, !6, i64 184, !6, i64 192, !6, i64 200, !6, i64 208, !6, i64 216, !6, i64 224, !6, i64 232, !6, i64 240, !6, i64 248, !6, i64 256, !6, i64 264, !13, i64 272, !13, i64 276, !13, i64 280, !13, i64 284, !13, i64 288, !13, i64 292, !13, i64 296, !13, i64 300, !13, i64 304, !13, i64 308, !13, i64 312, !13, i64 316, !13, i64 320, !13, i64 324, !13, i64 328, !13, i64 332, !13, i64 336, !13, i64 340, !13, i64 344, !13, i64 348, !13, i64 352, !13, i64 356, !13, i64 360, !13, i64 364, !13, i64 368, !13, i64 372, !13, i64 376, !10, i64 384, !10, i64 392, !10, i64 400, !6, i64 408, !6, i64 416, !6, i64 424, !6, i64 432, !6, i64 440, !6, i64 448, !6, i64 456, !6, i64 464, !6, i64 472, !6, i64 480, !6, i64 488, !6, i64 496, !6, i64 504, !6, i64 512, !6, i64 520, !6, i64 528, !6, i64 536, !6, i64 544, !6, i64 552, !6, i64 560, !6, i64 568, !6, i64 576, !6, i64 584, !6, i64 592, !14, i64 600, !14, i64 602, !14, i64 604, !7, i64 606}
!13 = !{!"int", !7, i64 0}
!14 = !{!"short", !7, i64 0}
!15 = !{!16, !10, i64 112}
!16 = !{!"reionization", !13, i64 0, !13, i64 4, !13, i64 8, !13, i64 12, !13, i64 16, !13, i64 20, !13, i64 24, !13, i64 28, !13, i64 32, !13, i64 36, !10, i64 40, !6, i64 48, !13, i64 56, !13, i64 60, !13, i64 64, !13, i64 68, !13, i64 72, !13, i64 76, !13, i64 80, !13, i64 84, !13, i64 88, !13, i64 92, !13, i64 96, !13, i64 100, !13, i64 104, !10, i64 112, !13, i64 120, !13, i64 124}
!17 = !{!16, !13, i64 104}
!18 = !{!16, !13, i64 68}
!19 = !{!16, !13, i64 56}
!20 = !{!16, !13, i64 60}
!21 = !{!16, !13, i64 64}
!22 = !{!16, !13, i64 72}
!23 = !{!16, !13, i64 80}
!24 = !{!16, !13, i64 84}
!25 = !{!16, !13, i64 76}
!26 = !{!16, !13, i64 92}
!27 = !{!16, !13, i64 88}
!28 = !{!16, !13, i64 96}
!29 = !{!13, !13, i64 0}
!30 = distinct !{!30, !31}
!31 = !{!"llvm.loop.mustprogress"}
!32 = !{!16, !13, i64 100}
!33 = distinct !{!33, !31}
!34 = distinct !{!34, !31}


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

@.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.137 = external hidden unnamed_addr constant [59 x i8], align 1
@__func__.thermodynamics_reionization_function = external hidden unnamed_addr constant [37 x i8], align 1
@.str.138 = external hidden unnamed_addr constant [38 x i8], align 1
@.str.139 = external hidden unnamed_addr constant [10 x i8], align 1
@.str.140 = external hidden unnamed_addr constant [1 x i8], align 1
@.str.141 = external hidden unnamed_addr constant [8 x i8], align 1
@.str.143 = external hidden unnamed_addr constant [10 x i8], align 1
@.str.144 = external hidden unnamed_addr constant [7 x i8], align 1
@.str.145 = external hidden unnamed_addr constant [41 x i8], align 1
@.str.146 = external hidden unnamed_addr constant [7 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: nounwind uwtable
define dso_local range(i32 0, 2) i32 @thermodynamics_reionization_function(double noundef %0, ptr noundef %1, ptr nocapture noundef readonly %2, ptr nocapture noundef %3) local_unnamed_addr #4 {
  %5 = alloca double, align 8
  %6 = alloca [2048 x i8], align 16
  %7 = alloca [2048 x i8], align 16
  %8 = alloca [2048 x i8], align 16
  %9 = alloca [2048 x i8], align 16
  %10 = alloca [2048 x i8], align 16
  %11 = alloca [2048 x i8], align 16
  %12 = alloca [2048 x i8], align 16
  %13 = alloca [2048 x i8], align 16
  %14 = alloca [2048 x i8], align 16
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %5) #5
  %15 = getelementptr inbounds i8, ptr %1, i64 12
  %16 = load i32, ptr %15, align 4, !tbaa !5
  switch i32 %16, label %336 [
    i32 1, label %17
    i32 3, label %17
    i32 2, label %133
    i32 4, label %208
    i32 5, label %285
  ]

17:                                               ; preds = %4, %4
  %18 = getelementptr inbounds i8, ptr %2, i64 112
  %19 = load ptr, ptr %18, align 8, !tbaa !13
  %20 = getelementptr inbounds i8, ptr %2, i64 104
  %21 = load i32, ptr %20, align 8, !tbaa !15
  %22 = sext i32 %21 to i64
  %23 = getelementptr inbounds double, ptr %19, i64 %22
  %24 = load double, ptr %23, align 8, !tbaa !16
  %25 = fcmp olt double %24, %0
  br i1 %25, label %26, label %32

26:                                               ; preds = %17
  %27 = getelementptr inbounds i8, ptr %2, i64 68
  %28 = load i32, ptr %27, align 4, !tbaa !17
  %29 = sext i32 %28 to i64
  %30 = getelementptr inbounds double, ptr %19, i64 %29
  %31 = load double, ptr %30, align 8, !tbaa !16
  store double %31, ptr %3, align 8, !tbaa !16
  br label %.thread

32:                                               ; preds = %17
  %33 = getelementptr inbounds i8, ptr %2, i64 56
  %34 = load i32, ptr %33, align 8, !tbaa !18
  %35 = sext i32 %34 to i64
  %36 = getelementptr inbounds double, ptr %19, i64 %35
  %37 = load double, ptr %36, align 8, !tbaa !16
  %38 = fadd double %37, 1.000000e+00
  %39 = getelementptr inbounds i8, ptr %2, i64 60
  %40 = load i32, ptr %39, align 4, !tbaa !19
  %41 = sext i32 %40 to i64
  %42 = getelementptr inbounds double, ptr %19, i64 %41
  %43 = load double, ptr %42, align 8, !tbaa !16
  %44 = tail call double @pow(double noundef %38, double noundef %43) #5
  %45 = fadd double %0, 1.000000e+00
  %46 = load ptr, ptr %18, align 8, !tbaa !13
  %47 = load i32, ptr %39, align 4, !tbaa !19
  %48 = sext i32 %47 to i64
  %49 = getelementptr inbounds double, ptr %46, i64 %48
  %50 = load double, ptr %49, align 8, !tbaa !16
  %51 = tail call double @pow(double noundef %45, double noundef %50) #5
  %52 = fsub double %44, %51
  %53 = load ptr, ptr %18, align 8, !tbaa !13
  %54 = load i32, ptr %39, align 4, !tbaa !19
  %55 = sext i32 %54 to i64
  %56 = getelementptr inbounds double, ptr %53, i64 %55
  %57 = load double, ptr %56, align 8, !tbaa !16
  %58 = load i32, ptr %33, align 8, !tbaa !18
  %59 = sext i32 %58 to i64
  %60 = getelementptr inbounds double, ptr %53, i64 %59
  %61 = load double, ptr %60, align 8, !tbaa !16
  %62 = fadd double %61, 1.000000e+00
  %63 = fadd double %57, -1.000000e+00
  %64 = tail call double @pow(double noundef %62, double noundef %63) #5
  %65 = fmul double %57, %64
  %66 = fdiv double %52, %65
  %67 = load ptr, ptr %18, align 8, !tbaa !13
  %68 = getelementptr inbounds i8, ptr %2, i64 64
  %69 = load i32, ptr %68, align 8, !tbaa !20
  %70 = sext i32 %69 to i64
  %71 = getelementptr inbounds double, ptr %67, i64 %70
  %72 = load double, ptr %71, align 8, !tbaa !16
  %73 = fdiv double %66, %72
  %74 = load i32, ptr %15, align 4, !tbaa !5
  %75 = icmp eq i32 %74, 1
  %76 = getelementptr inbounds i8, ptr %2, i64 72
  %77 = load i32, ptr %76, align 8, !tbaa !21
  %78 = sext i32 %77 to i64
  %79 = getelementptr inbounds double, ptr %67, i64 %78
  %80 = load double, ptr %79, align 8, !tbaa !16
  %81 = getelementptr inbounds i8, ptr %2, i64 68
  %82 = load i32, ptr %81, align 4, !tbaa !17
  %83 = sext i32 %82 to i64
  %84 = getelementptr inbounds double, ptr %67, i64 %83
  %85 = load double, ptr %84, align 8, !tbaa !16
  %86 = fsub double %80, %85
  %87 = tail call double @tanh(double noundef %73) #5
  br i1 %75, label %88, label %98

88:                                               ; preds = %32
  %89 = fadd double %87, 1.000000e+00
  %90 = fmul double %86, %89
  %91 = fmul double %90, 5.000000e-01
  %92 = load ptr, ptr %18, align 8, !tbaa !13
  %93 = load i32, ptr %81, align 4, !tbaa !17
  %94 = sext i32 %93 to i64
  %95 = getelementptr inbounds double, ptr %92, i64 %94
  %96 = load double, ptr %95, align 8, !tbaa !16
  %97 = fadd double %91, %96
  br label %105

98:                                               ; preds = %32
  %99 = load ptr, ptr %18, align 8, !tbaa !13
  %100 = load i32, ptr %81, align 4, !tbaa !17
  %101 = sext i32 %100 to i64
  %102 = getelementptr inbounds double, ptr %99, i64 %101
  %103 = load double, ptr %102, align 8, !tbaa !16
  %104 = tail call double @llvm.fmuladd.f64(double %86, double %87, double %103)
  br label %105

105:                                              ; preds = %98, %88
  %106 = phi ptr [ %99, %98 ], [ %92, %88 ]
  %storemerge = phi double [ %104, %98 ], [ %97, %88 ]
  store double %storemerge, ptr %3, align 8, !tbaa !16
  %107 = load i32, ptr %15, align 4, !tbaa !5
  %108 = icmp eq i32 %107, 1
  br i1 %108, label %109, label %.thread

109:                                              ; preds = %105
  %110 = getelementptr inbounds i8, ptr %2, i64 80
  %111 = load i32, ptr %110, align 8, !tbaa !22
  %112 = sext i32 %111 to i64
  %113 = getelementptr inbounds double, ptr %106, i64 %112
  %114 = load double, ptr %113, align 8, !tbaa !16
  %115 = fsub double %114, %0
  %116 = getelementptr inbounds i8, ptr %2, i64 84
  %117 = load i32, ptr %116, align 4, !tbaa !23
  %118 = sext i32 %117 to i64
  %119 = getelementptr inbounds double, ptr %106, i64 %118
  %120 = load double, ptr %119, align 8, !tbaa !16
  %121 = fdiv double %115, %120
  %122 = getelementptr inbounds i8, ptr %2, i64 76
  %123 = load i32, ptr %122, align 4, !tbaa !24
  %124 = sext i32 %123 to i64
  %125 = getelementptr inbounds double, ptr %106, i64 %124
  %126 = load double, ptr %125, align 8, !tbaa !16
  %127 = tail call double @tanh(double noundef %121) #5
  %128 = fadd double %127, 1.000000e+00
  %129 = fmul double %126, %128
  %130 = fmul double %129, 5.000000e-01
  %131 = load double, ptr %3, align 8, !tbaa !16
  %132 = fadd double %131, %130
  store double %132, ptr %3, align 8, !tbaa !16
  br label %.thread

133:                                              ; preds = %4
  %134 = getelementptr inbounds i8, ptr %2, i64 112
  %135 = load ptr, ptr %134, align 8, !tbaa !13
  %136 = getelementptr inbounds i8, ptr %2, i64 92
  %137 = load i32, ptr %136, align 4, !tbaa !25
  %138 = getelementptr inbounds i8, ptr %2, i64 88
  %139 = load i32, ptr %138, align 8, !tbaa !26
  %140 = add nsw i32 %139, %137
  %141 = sext i32 %140 to i64
  %142 = getelementptr double, ptr %135, i64 %141
  %143 = getelementptr i8, ptr %142, i64 -8
  %144 = load double, ptr %143, align 8, !tbaa !16
  %145 = fcmp olt double %144, %0
  br i1 %145, label %146, label %154

146:                                              ; preds = %133
  %147 = getelementptr inbounds i8, ptr %2, i64 96
  %148 = load i32, ptr %147, align 8, !tbaa !27
  %149 = add nsw i32 %148, %139
  %150 = sext i32 %149 to i64
  %151 = getelementptr double, ptr %135, i64 %150
  %152 = getelementptr i8, ptr %151, i64 -8
  %153 = load double, ptr %152, align 8, !tbaa !16
  store double %153, ptr %3, align 8, !tbaa !16
  br label %.thread

154:                                              ; preds = %133
  %155 = sext i32 %137 to i64
  %156 = getelementptr inbounds double, ptr %135, i64 %155
  %157 = load double, ptr %156, align 8, !tbaa !16
  %158 = fcmp ogt double %157, %0
  br i1 %158, label %159, label %.preheader.preheader

.preheader.preheader:                             ; preds = %154
  %invariant.gep231 = getelementptr double, ptr %135, i64 %155
  br label %.preheader

159:                                              ; preds = %154
  %160 = getelementptr inbounds i8, ptr %2, i64 96
  %161 = load i32, ptr %160, align 8, !tbaa !27
  %162 = sext i32 %161 to i64
  %163 = getelementptr inbounds double, ptr %135, i64 %162
  %164 = load double, ptr %163, align 8, !tbaa !16
  store double %164, ptr %3, align 8, !tbaa !16
  br label %.thread

.preheader:                                       ; preds = %.preheader.preheader, %.preheader
  %indvars.iv225 = phi i64 [ 0, %.preheader.preheader ], [ %indvars.iv.next226, %.preheader ]
  %gep232 = getelementptr double, ptr %invariant.gep231, i64 %indvars.iv225
  %165 = getelementptr i8, ptr %gep232, i64 8
  %166 = load double, ptr %165, align 8, !tbaa !16
  %167 = fcmp olt double %166, %0
  %indvars.iv.next226 = add nuw nsw i64 %indvars.iv225, 1
  br i1 %167, label %.preheader, label %168, !llvm.loop !28

168:                                              ; preds = %.preheader
  %169 = trunc nuw nsw i64 %indvars.iv225 to i32
  %170 = add nsw i32 %139, -2
  %171 = icmp eq i32 %170, %169
  %172 = load double, ptr %gep232, align 8, !tbaa !16
  br i1 %171, label %173, label %178

173:                                              ; preds = %168
  %174 = getelementptr i8, ptr %gep232, i64 -8
  %175 = load double, ptr %174, align 8, !tbaa !16
  %176 = fsub double %172, %175
  %177 = tail call double @llvm.fmuladd.f64(double %176, double 5.000000e-01, double %172)
  br label %181

178:                                              ; preds = %168
  %179 = fadd double %166, %172
  %180 = fmul double %179, 5.000000e-01
  br label %181

181:                                              ; preds = %178, %173
  %.0209 = phi double [ %177, %173 ], [ %180, %178 ]
  %182 = getelementptr inbounds i8, ptr %2, i64 96
  %183 = load i32, ptr %182, align 8, !tbaa !27
  %184 = add nsw i32 %183, %169
  %185 = sext i32 %184 to i64
  %186 = getelementptr inbounds double, ptr %135, i64 %185
  %187 = load double, ptr %186, align 8, !tbaa !16
  %188 = fsub double %0, %.0209
  %189 = getelementptr inbounds i8, ptr %2, i64 100
  %190 = load i32, ptr %189, align 4, !tbaa !30
  %191 = sext i32 %190 to i64
  %192 = getelementptr inbounds double, ptr %135, i64 %191
  %193 = load double, ptr %192, align 8, !tbaa !16
  %194 = fdiv double %188, %193
  %195 = tail call double @tanh(double noundef %194) #5
  %196 = fadd double %195, 1.000000e+00
  %197 = fmul double %196, 5.000000e-01
  %198 = load ptr, ptr %134, align 8, !tbaa !13
  %199 = load i32, ptr %182, align 8, !tbaa !27
  %200 = add nsw i32 %199, %169
  %201 = sext i32 %200 to i64
  %202 = getelementptr double, ptr %198, i64 %201
  %203 = getelementptr i8, ptr %202, i64 8
  %204 = load double, ptr %203, align 8, !tbaa !16
  %205 = load double, ptr %202, align 8, !tbaa !16
  %206 = fsub double %204, %205
  %207 = tail call double @llvm.fmuladd.f64(double %197, double %206, double %187)
  store double %207, ptr %3, align 8, !tbaa !16
  br label %.thread

208:                                              ; preds = %4
  %209 = getelementptr inbounds i8, ptr %2, i64 112
  %210 = load ptr, ptr %209, align 8, !tbaa !13
  %211 = getelementptr inbounds i8, ptr %2, i64 92
  %212 = load i32, ptr %211, align 4, !tbaa !25
  %213 = getelementptr inbounds i8, ptr %2, i64 88
  %214 = load i32, ptr %213, align 8, !tbaa !26
  %215 = add nsw i32 %214, %212
  %216 = sext i32 %215 to i64
  %217 = getelementptr double, ptr %210, i64 %216
  %218 = getelementptr i8, ptr %217, i64 -8
  %219 = load double, ptr %218, align 8, !tbaa !16
  %220 = fcmp olt double %219, %0
  br i1 %220, label %221, label %229

221:                                              ; preds = %208
  %222 = getelementptr inbounds i8, ptr %2, i64 96
  %223 = load i32, ptr %222, align 8, !tbaa !27
  %224 = add nsw i32 %223, %214
  %225 = sext i32 %224 to i64
  %226 = getelementptr double, ptr %210, i64 %225
  %227 = getelementptr i8, ptr %226, i64 -8
  %228 = load double, ptr %227, align 8, !tbaa !16
  store double %228, ptr %3, align 8, !tbaa !16
  br label %.thread

229:                                              ; preds = %208
  %230 = sext i32 %212 to i64
  %231 = getelementptr inbounds double, ptr %210, i64 %230
  %232 = load double, ptr %231, align 8, !tbaa !16
  %233 = fcmp olt double %232, %0
  %234 = getelementptr inbounds i8, ptr %2, i64 96
  %235 = load i32, ptr %234, align 8, !tbaa !27
  br i1 %233, label %236, label %281

236:                                              ; preds = %229
  %237 = add nsw i32 %235, %214
  %238 = sext i32 %237 to i64
  %239 = getelementptr double, ptr %210, i64 %238
  %240 = getelementptr i8, ptr %239, i64 -8
  %241 = load double, ptr %240, align 8, !tbaa !16
  store double %241, ptr %3, align 8, !tbaa !16
  %242 = icmp sgt i32 %214, 2
  br i1 %242, label %.lr.ph, label %.thread

.lr.ph:                                           ; preds = %236
  %243 = getelementptr inbounds i8, ptr %2, i64 100
  br label %244

244:                                              ; preds = %.lr.ph, %273
  %245 = phi i32 [ %214, %.lr.ph ], [ %278, %273 ]
  %.0210220 = phi i32 [ 1, %.lr.ph ], [ %277, %273 ]
  %246 = load ptr, ptr %209, align 8, !tbaa !13
  %247 = load i32, ptr %211, align 4, !tbaa !25
  %248 = xor i32 %.0210220, -1
  %249 = add i32 %245, %248
  %250 = add i32 %249, %247
  %251 = sext i32 %250 to i64
  %252 = getelementptr inbounds double, ptr %246, i64 %251
  %253 = load double, ptr %252, align 8, !tbaa !16
  %254 = load i32, ptr %234, align 8, !tbaa !27
  %255 = add nsw i32 %254, %245
  %256 = add i32 %255, %248
  %257 = sext i32 %256 to i64
  %258 = getelementptr inbounds double, ptr %246, i64 %257
  %259 = load double, ptr %258, align 8, !tbaa !16
  %260 = sub nsw i32 %255, %.0210220
  %261 = sext i32 %260 to i64
  %262 = getelementptr inbounds double, ptr %246, i64 %261
  %263 = load double, ptr %262, align 8, !tbaa !16
  %264 = fsub double %259, %263
  %265 = load i32, ptr %243, align 4, !tbaa !30
  %266 = sext i32 %265 to i64
  %267 = getelementptr inbounds double, ptr %246, i64 %266
  %268 = load double, ptr %267, align 8, !tbaa !16
  %269 = call i32 @thermodynamics_tanh(double noundef %0, double noundef %253, double noundef %264, double noundef 0.000000e+00, double noundef %268, ptr noundef nonnull %5)
  %270 = icmp eq i32 %269, 1
  br i1 %270, label %271, label %273

271:                                              ; preds = %244
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %6) #5
  %272 = getelementptr inbounds i8, ptr %1, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %6, ptr noundef nonnull @.str, ptr noundef nonnull @.str.137, ptr noundef nonnull %272) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %272, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_function, i32 noundef 2089, ptr noundef nonnull %6) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %6) #5
  br label %.thread

273:                                              ; preds = %244
  %274 = load double, ptr %5, align 8, !tbaa !16
  %275 = load double, ptr %3, align 8, !tbaa !16
  %276 = fadd double %274, %275
  store double %276, ptr %3, align 8, !tbaa !16
  %277 = add nuw nsw i32 %.0210220, 1
  %278 = load i32, ptr %213, align 8, !tbaa !26
  %279 = add nsw i32 %278, -1
  %280 = icmp slt i32 %277, %279
  br i1 %280, label %244, label %.thread, !llvm.loop !31

281:                                              ; preds = %229
  %282 = sext i32 %235 to i64
  %283 = getelementptr inbounds double, ptr %210, i64 %282
  %284 = load double, ptr %283, align 8, !tbaa !16
  store double %284, ptr %3, align 8, !tbaa !16
  br label %.thread

285:                                              ; preds = %4
  %286 = getelementptr inbounds i8, ptr %2, i64 112
  %287 = load ptr, ptr %286, align 8, !tbaa !13
  %288 = getelementptr inbounds i8, ptr %2, i64 92
  %289 = load i32, ptr %288, align 4, !tbaa !25
  %290 = getelementptr inbounds i8, ptr %2, i64 88
  %291 = load i32, ptr %290, align 8, !tbaa !26
  %292 = add nsw i32 %291, %289
  %293 = sext i32 %292 to i64
  %294 = getelementptr double, ptr %287, i64 %293
  %295 = getelementptr i8, ptr %294, i64 -8
  %296 = load double, ptr %295, align 8, !tbaa !16
  %297 = fcmp olt double %296, %0
  br i1 %297, label %299, label %.preheader216.preheader

.preheader216.preheader:                          ; preds = %285
  %298 = sext i32 %289 to i64
  %invariant.gep = getelementptr double, ptr %287, i64 %298
  br label %.preheader216

299:                                              ; preds = %285
  %300 = getelementptr inbounds i8, ptr %2, i64 96
  %301 = load i32, ptr %300, align 8, !tbaa !27
  %302 = add nsw i32 %301, %291
  %303 = sext i32 %302 to i64
  %304 = getelementptr double, ptr %287, i64 %303
  %305 = getelementptr i8, ptr %304, i64 -8
  %306 = load double, ptr %305, align 8, !tbaa !16
  store double %306, ptr %3, align 8, !tbaa !16
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %7) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %7, ptr noundef nonnull @.str.138) #5
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %8) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %8, ptr noundef nonnull @.str.139, ptr noundef nonnull %7) #5
  %307 = getelementptr inbounds i8, ptr %1, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %307, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_function, i32 noundef 2113, ptr noundef nonnull %8) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %8) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %7) #5
  br label %.thread

.preheader216:                                    ; preds = %.preheader216.preheader, %.preheader216
  %indvars.iv = phi i64 [ 0, %.preheader216.preheader ], [ %indvars.iv.next, %.preheader216 ]
  %gep = getelementptr double, ptr %invariant.gep, i64 %indvars.iv
  %308 = getelementptr i8, ptr %gep, i64 8
  %309 = load double, ptr %308, align 8, !tbaa !16
  %310 = fcmp olt double %309, %0
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  br i1 %310, label %.preheader216, label %311, !llvm.loop !32

311:                                              ; preds = %.preheader216
  %312 = load double, ptr %gep, align 8, !tbaa !16
  %313 = fcmp ogt double %312, %0
  br i1 %313, label %314, label %316

314:                                              ; preds = %311
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %9) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %9, ptr noundef nonnull @.str.140) #5
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %10) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %10, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.141, ptr noundef nonnull %9) #5
  %315 = getelementptr inbounds i8, ptr %1, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %315, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_function, i32 noundef 2126, ptr noundef nonnull %10) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %10) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %9) #5
  br label %.thread

316:                                              ; preds = %311
  %317 = trunc nuw nsw i64 %indvars.iv to i32
  %318 = fsub double %0, %312
  %319 = fsub double %309, %312
  %320 = fdiv double %318, %319
  %321 = getelementptr inbounds i8, ptr %2, i64 96
  %322 = load i32, ptr %321, align 8, !tbaa !27
  %323 = add nsw i32 %322, %317
  %324 = sext i32 %323 to i64
  %325 = getelementptr inbounds double, ptr %287, i64 %324
  %326 = load double, ptr %325, align 8, !tbaa !16
  %327 = getelementptr i8, ptr %325, i64 8
  %328 = load double, ptr %327, align 8, !tbaa !16
  %329 = fsub double %328, %326
  %330 = tail call double @llvm.fmuladd.f64(double %320, double %329, double %326)
  store double %330, ptr %3, align 8, !tbaa !16
  %331 = fcmp uge double %330, 0.000000e+00
  br i1 %331, label %.thread, label %332

332:                                              ; preds = %316
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %11) #5
  %333 = load double, ptr %325, align 8, !tbaa !16
  %334 = load double, ptr %327, align 8, !tbaa !16
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %11, ptr noundef nonnull @.str.143, double noundef %320, double noundef %333, double noundef %334) #5
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %12) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %12, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.144, ptr noundef nonnull %11) #5
  %335 = getelementptr inbounds i8, ptr %1, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %335, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_function, i32 noundef 2145, ptr noundef nonnull %12) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %12) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %11) #5
  br label %.thread

336:                                              ; preds = %4
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %13) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %13, ptr noundef nonnull @.str.145, i32 noundef %16) #5
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %14) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %14, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.146, ptr noundef nonnull %13) #5
  %337 = getelementptr inbounds i8, ptr %1, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %337, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_reionization_function, i32 noundef 2155, ptr noundef nonnull %14) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %14) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %13) #5
  br label %.thread

.thread:                                          ; preds = %273, %236, %332, %314, %316, %221, %281, %146, %181, %159, %26, %109, %105, %336, %299, %271
  %.2 = phi i32 [ 1, %271 ], [ 1, %299 ], [ 1, %336 ], [ 0, %105 ], [ 0, %109 ], [ 0, %26 ], [ 0, %159 ], [ 0, %181 ], [ 0, %146 ], [ 0, %281 ], [ 0, %221 ], [ 0, %316 ], [ 1, %314 ], [ 1, %332 ], [ 0, %236 ], [ 0, %273 ]
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %5) #5
  ret i32 %.2
}

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

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

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 = { 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 #5 = { nounwind }

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !10, i64 12}
!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, !12, i64 112}
!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 = !{!14, !10, i64 104}
!16 = !{!7, !7, i64 0}
!17 = !{!14, !10, i64 68}
!18 = !{!14, !10, i64 56}
!19 = !{!14, !10, i64 60}
!20 = !{!14, !10, i64 64}
!21 = !{!14, !10, i64 72}
!22 = !{!14, !10, i64 80}
!23 = !{!14, !10, i64 84}
!24 = !{!14, !10, i64 76}
!25 = !{!14, !10, i64 92}
!26 = !{!14, !10, i64 88}
!27 = !{!14, !10, i64 96}
!28 = distinct !{!28, !29}
!29 = !{!"llvm.loop.mustprogress"}
!30 = !{!14, !10, i64 100}
!31 = distinct !{!31, !29}
!32 = distinct !{!32, !29}


; Effective O3 passes: EarlyCSEPass,GVNPass,GlobalOptPass,IPSCCPPass,IndVarSimplifyPass,InferFunctionAttrsPass,InstCombinePass,InstSimplifyPass,JumpThreadingPass,LCSSAPass,LICMPass,LoopRotatePass,LoopSimplifyPass,LoopUnrollPass,PostOrderFunctionAttrsPass,PromotePass,ReassociatePass,SROAPass,SimplifyCFGPass,TailCallElimPass