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

@.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.37 = external hidden unnamed_addr constant [35 x i8], align 1
@.str.38 = external hidden unnamed_addr constant [9 x i8], align 1
@__func__.thermodynamics_helium_from_bbn = external hidden unnamed_addr constant [31 x i8], align 1
@.str.104 = external hidden unnamed_addr constant [42 x i8], align 1
@.str.105 = external hidden unnamed_addr constant [90 x i8], align 1
@.str.106 = external hidden unnamed_addr constant [2 x i8], align 1
@.str.107 = external hidden unnamed_addr constant [43 x i8], align 1
@.str.108 = external hidden unnamed_addr constant [3 x i8], align 1
@.str.109 = external hidden unnamed_addr constant [4 x i8], align 1
@.str.110 = external hidden unnamed_addr constant [6 x i8], align 1
@.str.111 = external hidden unnamed_addr constant [71 x i8], align 1
@.str.112 = external hidden unnamed_addr constant [50 x i8], align 1
@.str.113 = external hidden unnamed_addr constant [7 x i8], align 1
@.str.114 = external hidden unnamed_addr constant [7 x i8], align 1
@.str.115 = external hidden unnamed_addr constant [4 x i8], align 1
@.str.116 = external hidden unnamed_addr constant [6 x i8], align 1
@.str.117 = external hidden unnamed_addr constant [14 x i8], align 1
@.str.118 = external hidden unnamed_addr constant [16 x i8], align 1
@.str.119 = external hidden unnamed_addr constant [12 x i8], align 1
@.str.120 = external hidden unnamed_addr constant [67 x i8], align 1
@.str.121 = external hidden unnamed_addr constant [121 x i8], align 1
@.str.122 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.123 = external hidden unnamed_addr constant [248 x i8], align 1
@.str.124 = external hidden unnamed_addr constant [20 x i8], align 1
@.str.125 = external hidden unnamed_addr constant [247 x i8], align 1
@.str.126 = external hidden unnamed_addr constant [31 x i8], align 1
@.str.127 = external hidden unnamed_addr constant [255 x i8], align 1
@.str.128 = external hidden unnamed_addr constant [22 x i8], align 1
@.str.129 = external hidden unnamed_addr constant [254 x i8], align 1
@.str.130 = external hidden unnamed_addr constant [33 x i8], align 1
@.str.131 = external hidden unnamed_addr constant [140 x i8], align 1
@.str.132 = external hidden unnamed_addr constant [103 x i8], align 1
@.str.133 = external hidden unnamed_addr constant [138 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

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

; Function Attrs: nounwind uwtable
define dso_local i32 @thermodynamics_helium_from_bbn(ptr noundef %0, ptr noundef %1, ptr noundef %2) #4 {
  %4 = alloca i32, align 4
  %5 = alloca ptr, align 8
  %6 = alloca ptr, align 8
  %7 = alloca ptr, align 8
  %8 = alloca ptr, align 8
  %9 = alloca [8192 x i8], align 16
  %10 = alloca ptr, align 8
  %11 = alloca i32, align 4
  %12 = alloca i32, align 4
  %13 = alloca ptr, align 8
  %14 = alloca ptr, align 8
  %15 = alloca ptr, align 8
  %16 = alloca ptr, align 8
  %17 = alloca ptr, align 8
  %18 = alloca ptr, align 8
  %19 = alloca i32, align 4
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca i32, align 4
  %23 = alloca double, align 8
  %24 = alloca double, align 8
  %25 = alloca double, align 8
  %26 = alloca ptr, align 8
  %27 = alloca i32, align 4
  %28 = alloca [2048 x i8], align 16
  %29 = alloca i32, align 4
  %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 i32, align 4
  %36 = alloca [2048 x i8], align 16
  %37 = alloca i32, align 4
  %38 = alloca [2048 x i8], align 16
  %39 = alloca i32, align 4
  %40 = alloca [2048 x i8], align 16
  %41 = alloca i32, align 4
  %42 = alloca [2048 x i8], align 16
  %43 = alloca i32, align 4
  %44 = alloca [2048 x i8], align 16
  %45 = alloca i32, align 4
  %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
  store ptr %0, ptr %5, align 8, !tbaa !5
  store ptr %1, ptr %6, align 8, !tbaa !5
  store ptr %2, ptr %7, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %8) #6
  call void @llvm.lifetime.start.p0(i64 8192, ptr %9) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %10) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %11) #6
  store i32 0, ptr %11, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 4, ptr %12) #6
  store i32 0, ptr %12, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %13) #6
  store ptr null, ptr %13, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #6
  store ptr null, ptr %14, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #6
  store ptr null, ptr %15, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #6
  store ptr null, ptr %16, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %17) #6
  store ptr null, ptr %17, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %18) #6
  store ptr null, ptr %18, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr %19) #6
  store i32 0, ptr %19, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %20) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %22) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %23) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %24) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %25) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %26) #6
  %61 = load ptr, ptr %6, align 8, !tbaa !5
  %62 = getelementptr inbounds %struct.background, ptr %61, i32 0, i32 99
  %63 = load i32, ptr %62, align 8, !tbaa !11
  %64 = sext i32 %63 to i64
  %65 = mul i64 %64, 8
  %66 = call noalias ptr @malloc(i64 noundef %65) #7
  store ptr %66, ptr %26, align 8, !tbaa !5
  %67 = load ptr, ptr %26, align 8, !tbaa !5
  %68 = icmp eq ptr %67, null
  br i1 %68, label %69, label %82

69:                                               ; preds = %3
  call void @llvm.lifetime.start.p0(i64 4, ptr %27) #6
  %70 = load ptr, ptr %6, align 8, !tbaa !5
  %71 = getelementptr inbounds %struct.background, ptr %70, i32 0, i32 99
  %72 = load i32, ptr %71, align 8, !tbaa !11
  %73 = sext i32 %72 to i64
  %74 = mul i64 %73, 8
  %75 = trunc i64 %74 to i32
  store i32 %75, ptr %27, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 2048, ptr %28) #6
  %76 = getelementptr inbounds [2048 x i8], ptr %28, i64 0, i64 0
  %77 = load i32, ptr %27, align 4, !tbaa !9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %76, ptr noundef @.str.37, ptr noundef @.str.38, i32 noundef %77)
  %78 = load ptr, ptr %6, align 8, !tbaa !5
  %79 = getelementptr inbounds %struct.background, ptr %78, i32 0, i32 149
  %80 = getelementptr inbounds [2048 x i8], ptr %79, i64 0, i64 0
  %81 = getelementptr inbounds [2048 x i8], ptr %28, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %80, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1575, ptr noundef %81)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %28) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %27) #6
  br label %567

82:                                               ; preds = %3
  %83 = load ptr, ptr %6, align 8, !tbaa !5
  %84 = getelementptr inbounds %struct.background, ptr %83, i32 0, i32 2
  %85 = load double, ptr %84, align 8, !tbaa !15
  %86 = fmul double 0x3DD7AFE2C34414DC, %85
  %87 = fdiv double 1.000000e-01, %86
  %88 = fsub double %87, 1.000000e+00
  store double %88, ptr %24, align 8, !tbaa !16
  %89 = load ptr, ptr %6, align 8, !tbaa !5
  %90 = load double, ptr %24, align 8, !tbaa !16
  %91 = call i32 @background_tau_of_z(ptr noundef %89, double noundef %90, ptr noundef %25)
  %92 = icmp eq i32 %91, 1
  br i1 %92, label %93, label %102

93:                                               ; preds = %82
  call void @llvm.lifetime.start.p0(i64 2048, ptr %30) #6
  %94 = getelementptr inbounds [2048 x i8], ptr %30, i64 0, i64 0
  %95 = load ptr, ptr %6, align 8, !tbaa !5
  %96 = getelementptr inbounds %struct.background, ptr %95, i32 0, i32 149
  %97 = getelementptr inbounds [2048 x i8], ptr %96, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %94, ptr noundef @.str, ptr noundef @.str.104, ptr noundef %97)
  %98 = load ptr, ptr %7, align 8, !tbaa !5
  %99 = getelementptr inbounds %struct.thermo, ptr %98, i32 0, i32 93
  %100 = getelementptr inbounds [2048 x i8], ptr %99, i64 0, i64 0
  %101 = getelementptr inbounds [2048 x i8], ptr %30, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %100, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1584, ptr noundef %101)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %30) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

102:                                              ; preds = %82
  %103 = load ptr, ptr %6, align 8, !tbaa !5
  %104 = load double, ptr %25, align 8, !tbaa !16
  %105 = load ptr, ptr %6, align 8, !tbaa !5
  %106 = getelementptr inbounds %struct.background, ptr %105, i32 0, i32 143
  %107 = load i16, ptr %106, align 4, !tbaa !17
  %108 = load ptr, ptr %6, align 8, !tbaa !5
  %109 = getelementptr inbounds %struct.background, ptr %108, i32 0, i32 144
  %110 = load i16, ptr %109, align 2, !tbaa !18
  %111 = load ptr, ptr %26, align 8, !tbaa !5
  %112 = call i32 @background_at_tau(ptr noundef %103, double noundef %104, i16 noundef signext %107, i16 noundef signext %110, ptr noundef %22, ptr noundef %111)
  %113 = icmp eq i32 %112, 1
  br i1 %113, label %114, label %123

114:                                              ; preds = %102
  call void @llvm.lifetime.start.p0(i64 2048, ptr %31) #6
  %115 = getelementptr inbounds [2048 x i8], ptr %31, i64 0, i64 0
  %116 = load ptr, ptr %6, align 8, !tbaa !5
  %117 = getelementptr inbounds %struct.background, ptr %116, i32 0, i32 149
  %118 = getelementptr inbounds [2048 x i8], ptr %117, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %115, ptr noundef @.str, ptr noundef @.str.105, ptr noundef %118)
  %119 = load ptr, ptr %7, align 8, !tbaa !5
  %120 = getelementptr inbounds %struct.thermo, ptr %119, i32 0, i32 93
  %121 = getelementptr inbounds [2048 x i8], ptr %120, i64 0, i64 0
  %122 = getelementptr inbounds [2048 x i8], ptr %31, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %121, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1593, ptr noundef %122)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %31) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

123:                                              ; preds = %102
  %124 = load ptr, ptr %26, align 8, !tbaa !5
  %125 = load ptr, ptr %6, align 8, !tbaa !5
  %126 = getelementptr inbounds %struct.background, ptr %125, i32 0, i32 87
  %127 = load i32, ptr %126, align 8, !tbaa !19
  %128 = sext i32 %127 to i64
  %129 = getelementptr inbounds double, ptr %124, i64 %128
  %130 = load double, ptr %129, align 8, !tbaa !16
  %131 = load ptr, ptr %26, align 8, !tbaa !5
  %132 = load ptr, ptr %6, align 8, !tbaa !5
  %133 = getelementptr inbounds %struct.background, ptr %132, i32 0, i32 88
  %134 = load i32, ptr %133, align 4, !tbaa !20
  %135 = sext i32 %134 to i64
  %136 = getelementptr inbounds double, ptr %131, i64 %135
  %137 = load double, ptr %136, align 8, !tbaa !16
  %138 = load ptr, ptr %26, align 8, !tbaa !5
  %139 = load ptr, ptr %6, align 8, !tbaa !5
  %140 = getelementptr inbounds %struct.background, ptr %139, i32 0, i32 62
  %141 = load i32, ptr %140, align 4, !tbaa !21
  %142 = sext i32 %141 to i64
  %143 = getelementptr inbounds double, ptr %138, i64 %142
  %144 = load double, ptr %143, align 8, !tbaa !16
  %145 = fneg double %144
  %146 = call double @llvm.fmuladd.f64(double %130, double %137, double %145)
  %147 = call double @pow(double noundef 0x3FD745D1745D1746, double noundef 0x3FF5555555555555) #6
  %148 = fmul double 8.750000e-01, %147
  %149 = load ptr, ptr %26, align 8, !tbaa !5
  %150 = load ptr, ptr %6, align 8, !tbaa !5
  %151 = getelementptr inbounds %struct.background, ptr %150, i32 0, i32 62
  %152 = load i32, ptr %151, align 4, !tbaa !21
  %153 = sext i32 %152 to i64
  %154 = getelementptr inbounds double, ptr %149, i64 %153
  %155 = load double, ptr %154, align 8, !tbaa !16
  %156 = fmul double %148, %155
  %157 = fdiv double %146, %156
  store double %157, ptr %23, align 8, !tbaa !16
  %158 = load ptr, ptr %26, align 8, !tbaa !5
  call void @free(ptr noundef %158) #6
  %159 = load double, ptr %23, align 8, !tbaa !16
  %160 = fsub double %159, 3.046000e+00
  store double %160, ptr %20, align 8, !tbaa !16
  %161 = load ptr, ptr %5, align 8, !tbaa !5
  %162 = getelementptr inbounds %struct.precision, ptr %161, i32 0, i32 13
  %163 = getelementptr inbounds [256 x i8], ptr %162, i64 0, i64 0
  %164 = call noalias ptr @fopen(ptr noundef %163, ptr noundef @.str.106)
  store ptr %164, ptr %8, align 8, !tbaa !5
  %165 = load ptr, ptr %8, align 8, !tbaa !5
  %166 = icmp eq ptr %165, null
  br i1 %166, label %167, label %176

167:                                              ; preds = %123
  call void @llvm.lifetime.start.p0(i64 2048, ptr %32) #6
  %168 = getelementptr inbounds [2048 x i8], ptr %32, i64 0, i64 0
  %169 = load ptr, ptr %5, align 8, !tbaa !5
  %170 = getelementptr inbounds %struct.precision, ptr %169, i32 0, i32 13
  %171 = getelementptr inbounds [256 x i8], ptr %170, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %168, ptr noundef @.str.107, ptr noundef @.str.108, ptr noundef %171, ptr noundef @.str.109)
  %172 = load ptr, ptr %7, align 8, !tbaa !5
  %173 = getelementptr inbounds %struct.thermo, ptr %172, i32 0, i32 93
  %174 = getelementptr inbounds [2048 x i8], ptr %173, i64 0, i64 0
  %175 = getelementptr inbounds [2048 x i8], ptr %32, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %174, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1619, ptr noundef %175)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %32) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

176:                                              ; preds = %123
  br label %177

177:                                              ; preds = %372, %176
  %178 = getelementptr inbounds [8192 x i8], ptr %9, i64 0, i64 0
  %179 = load ptr, ptr %8, align 8, !tbaa !5
  %180 = call ptr @fgets(ptr noundef %178, i32 noundef 8191, ptr noundef %179)
  %181 = icmp ne ptr %180, null
  br i1 %181, label %182, label %373

182:                                              ; preds = %177
  %183 = getelementptr inbounds [8192 x i8], ptr %9, i64 0, i64 0
  store ptr %183, ptr %10, align 8, !tbaa !5
  br label %184

184:                                              ; preds = %190, %182
  %185 = load ptr, ptr %10, align 8, !tbaa !5
  %186 = getelementptr inbounds i8, ptr %185, i64 0
  %187 = load i8, ptr %186, align 1, !tbaa !22
  %188 = sext i8 %187 to i32
  %189 = icmp eq i32 %188, 32
  br i1 %189, label %190, label %193

190:                                              ; preds = %184
  %191 = load ptr, ptr %10, align 8, !tbaa !5
  %192 = getelementptr inbounds i8, ptr %191, i32 1
  store ptr %192, ptr %10, align 8, !tbaa !5
  br label %184, !llvm.loop !23

193:                                              ; preds = %184
  %194 = load ptr, ptr %10, align 8, !tbaa !5
  %195 = getelementptr inbounds i8, ptr %194, i64 0
  %196 = load i8, ptr %195, align 1, !tbaa !22
  %197 = sext i8 %196 to i32
  %198 = icmp sgt i32 %197, 39
  br i1 %198, label %199, label %372

199:                                              ; preds = %193
  %200 = load i32, ptr %11, align 4, !tbaa !9
  %201 = icmp eq i32 %200, 0
  br i1 %201, label %202, label %337

202:                                              ; preds = %199
  %203 = load i32, ptr %12, align 4, !tbaa !9
  %204 = icmp eq i32 %203, 0
  br i1 %204, label %205, label %337

205:                                              ; preds = %202
  %206 = getelementptr inbounds [8192 x i8], ptr %9, i64 0, i64 0
  %207 = call i32 (ptr, ptr, ...) @__isoc99_sscanf(ptr noundef %206, ptr noundef @.str.110, ptr noundef %11, ptr noundef %12) #6
  %208 = icmp ne i32 %207, 2
  br i1 %208, label %209, label %220

209:                                              ; preds = %205
  call void @llvm.lifetime.start.p0(i64 2048, ptr %33) #6
  %210 = getelementptr inbounds [2048 x i8], ptr %33, i64 0, i64 0
  %211 = load ptr, ptr %5, align 8, !tbaa !5
  %212 = getelementptr inbounds %struct.precision, ptr %211, i32 0, i32 13
  %213 = getelementptr inbounds [256 x i8], ptr %212, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %210, ptr noundef @.str.111, ptr noundef %213)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %34) #6
  %214 = getelementptr inbounds [2048 x i8], ptr %34, i64 0, i64 0
  %215 = getelementptr inbounds [2048 x i8], ptr %33, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %214, ptr noundef @.str.10, ptr noundef @.str.112, ptr noundef %215)
  %216 = load ptr, ptr %7, align 8, !tbaa !5
  %217 = getelementptr inbounds %struct.thermo, ptr %216, i32 0, i32 93
  %218 = getelementptr inbounds [2048 x i8], ptr %217, i64 0, i64 0
  %219 = getelementptr inbounds [2048 x i8], ptr %34, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %218, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1645, ptr noundef %219)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %34) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %33) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

220:                                              ; preds = %205
  %221 = load i32, ptr %11, align 4, !tbaa !9
  %222 = sext i32 %221 to i64
  %223 = mul i64 %222, 8
  %224 = call noalias ptr @malloc(i64 noundef %223) #7
  store ptr %224, ptr %13, align 8, !tbaa !5
  %225 = load ptr, ptr %13, align 8, !tbaa !5
  %226 = icmp eq ptr %225, null
  br i1 %226, label %227, label %238

227:                                              ; preds = %220
  call void @llvm.lifetime.start.p0(i64 4, ptr %35) #6
  %228 = load i32, ptr %11, align 4, !tbaa !9
  %229 = sext i32 %228 to i64
  %230 = mul i64 %229, 8
  %231 = trunc i64 %230 to i32
  store i32 %231, ptr %35, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 2048, ptr %36) #6
  %232 = getelementptr inbounds [2048 x i8], ptr %36, i64 0, i64 0
  %233 = load i32, ptr %35, align 4, !tbaa !9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %232, ptr noundef @.str.37, ptr noundef @.str.113, i32 noundef %233)
  %234 = load ptr, ptr %7, align 8, !tbaa !5
  %235 = getelementptr inbounds %struct.thermo, ptr %234, i32 0, i32 93
  %236 = getelementptr inbounds [2048 x i8], ptr %235, i64 0, i64 0
  %237 = getelementptr inbounds [2048 x i8], ptr %36, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %236, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1647, ptr noundef %237)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %36) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %35) #6
  br label %567

238:                                              ; preds = %220
  %239 = load i32, ptr %12, align 4, !tbaa !9
  %240 = sext i32 %239 to i64
  %241 = mul i64 %240, 8
  %242 = call noalias ptr @malloc(i64 noundef %241) #7
  store ptr %242, ptr %14, align 8, !tbaa !5
  %243 = load ptr, ptr %14, align 8, !tbaa !5
  %244 = icmp eq ptr %243, null
  br i1 %244, label %245, label %256

245:                                              ; preds = %238
  call void @llvm.lifetime.start.p0(i64 4, ptr %37) #6
  %246 = load i32, ptr %12, align 4, !tbaa !9
  %247 = sext i32 %246 to i64
  %248 = mul i64 %247, 8
  %249 = trunc i64 %248 to i32
  store i32 %249, ptr %37, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 2048, ptr %38) #6
  %250 = getelementptr inbounds [2048 x i8], ptr %38, i64 0, i64 0
  %251 = load i32, ptr %37, align 4, !tbaa !9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %250, ptr noundef @.str.37, ptr noundef @.str.114, i32 noundef %251)
  %252 = load ptr, ptr %7, align 8, !tbaa !5
  %253 = getelementptr inbounds %struct.thermo, ptr %252, i32 0, i32 93
  %254 = getelementptr inbounds [2048 x i8], ptr %253, i64 0, i64 0
  %255 = getelementptr inbounds [2048 x i8], ptr %38, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %254, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1648, ptr noundef %255)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %38) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %37) #6
  br label %567

256:                                              ; preds = %238
  %257 = load i32, ptr %11, align 4, !tbaa !9
  %258 = load i32, ptr %12, align 4, !tbaa !9
  %259 = mul nsw i32 %257, %258
  %260 = sext i32 %259 to i64
  %261 = mul i64 %260, 8
  %262 = call noalias ptr @malloc(i64 noundef %261) #7
  store ptr %262, ptr %15, align 8, !tbaa !5
  %263 = load ptr, ptr %15, align 8, !tbaa !5
  %264 = icmp eq ptr %263, null
  br i1 %264, label %265, label %278

265:                                              ; preds = %256
  call void @llvm.lifetime.start.p0(i64 4, ptr %39) #6
  %266 = load i32, ptr %11, align 4, !tbaa !9
  %267 = load i32, ptr %12, align 4, !tbaa !9
  %268 = mul nsw i32 %266, %267
  %269 = sext i32 %268 to i64
  %270 = mul i64 %269, 8
  %271 = trunc i64 %270 to i32
  store i32 %271, ptr %39, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 2048, ptr %40) #6
  %272 = getelementptr inbounds [2048 x i8], ptr %40, i64 0, i64 0
  %273 = load i32, ptr %39, align 4, !tbaa !9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %272, ptr noundef @.str.37, ptr noundef @.str.115, i32 noundef %273)
  %274 = load ptr, ptr %7, align 8, !tbaa !5
  %275 = getelementptr inbounds %struct.thermo, ptr %274, i32 0, i32 93
  %276 = getelementptr inbounds [2048 x i8], ptr %275, i64 0, i64 0
  %277 = getelementptr inbounds [2048 x i8], ptr %40, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %276, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1649, ptr noundef %277)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %40) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %39) #6
  br label %567

278:                                              ; preds = %256
  %279 = load i32, ptr %11, align 4, !tbaa !9
  %280 = load i32, ptr %12, align 4, !tbaa !9
  %281 = mul nsw i32 %279, %280
  %282 = sext i32 %281 to i64
  %283 = mul i64 %282, 8
  %284 = call noalias ptr @malloc(i64 noundef %283) #7
  store ptr %284, ptr %16, align 8, !tbaa !5
  %285 = load ptr, ptr %16, align 8, !tbaa !5
  %286 = icmp eq ptr %285, null
  br i1 %286, label %287, label %300

287:                                              ; preds = %278
  call void @llvm.lifetime.start.p0(i64 4, ptr %41) #6
  %288 = load i32, ptr %11, align 4, !tbaa !9
  %289 = load i32, ptr %12, align 4, !tbaa !9
  %290 = mul nsw i32 %288, %289
  %291 = sext i32 %290 to i64
  %292 = mul i64 %291, 8
  %293 = trunc i64 %292 to i32
  store i32 %293, ptr %41, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 2048, ptr %42) #6
  %294 = getelementptr inbounds [2048 x i8], ptr %42, i64 0, i64 0
  %295 = load i32, ptr %41, align 4, !tbaa !9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %294, ptr noundef @.str.37, ptr noundef @.str.116, i32 noundef %295)
  %296 = load ptr, ptr %7, align 8, !tbaa !5
  %297 = getelementptr inbounds %struct.thermo, ptr %296, i32 0, i32 93
  %298 = getelementptr inbounds [2048 x i8], ptr %297, i64 0, i64 0
  %299 = getelementptr inbounds [2048 x i8], ptr %42, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %298, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1650, ptr noundef %299)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %42) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %41) #6
  br label %567

300:                                              ; preds = %278
  %301 = load i32, ptr %11, align 4, !tbaa !9
  %302 = sext i32 %301 to i64
  %303 = mul i64 %302, 8
  %304 = call noalias ptr @malloc(i64 noundef %303) #7
  store ptr %304, ptr %17, align 8, !tbaa !5
  %305 = load ptr, ptr %17, align 8, !tbaa !5
  %306 = icmp eq ptr %305, null
  br i1 %306, label %307, label %318

307:                                              ; preds = %300
  call void @llvm.lifetime.start.p0(i64 4, ptr %43) #6
  %308 = load i32, ptr %11, align 4, !tbaa !9
  %309 = sext i32 %308 to i64
  %310 = mul i64 %309, 8
  %311 = trunc i64 %310 to i32
  store i32 %311, ptr %43, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 2048, ptr %44) #6
  %312 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  %313 = load i32, ptr %43, align 4, !tbaa !9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %312, ptr noundef @.str.37, ptr noundef @.str.117, i32 noundef %313)
  %314 = load ptr, ptr %7, align 8, !tbaa !5
  %315 = getelementptr inbounds %struct.thermo, ptr %314, i32 0, i32 93
  %316 = getelementptr inbounds [2048 x i8], ptr %315, i64 0, i64 0
  %317 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %316, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1651, ptr noundef %317)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %44) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %43) #6
  br label %567

318:                                              ; preds = %300
  %319 = load i32, ptr %11, align 4, !tbaa !9
  %320 = sext i32 %319 to i64
  %321 = mul i64 %320, 8
  %322 = call noalias ptr @malloc(i64 noundef %321) #7
  store ptr %322, ptr %18, align 8, !tbaa !5
  %323 = load ptr, ptr %18, align 8, !tbaa !5
  %324 = icmp eq ptr %323, null
  br i1 %324, label %325, label %336

325:                                              ; preds = %318
  call void @llvm.lifetime.start.p0(i64 4, ptr %45) #6
  %326 = load i32, ptr %11, align 4, !tbaa !9
  %327 = sext i32 %326 to i64
  %328 = mul i64 %327, 8
  %329 = trunc i64 %328 to i32
  store i32 %329, ptr %45, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 2048, ptr %46) #6
  %330 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  %331 = load i32, ptr %45, align 4, !tbaa !9
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %330, ptr noundef @.str.37, ptr noundef @.str.118, i32 noundef %331)
  %332 = load ptr, ptr %7, align 8, !tbaa !5
  %333 = getelementptr inbounds %struct.thermo, ptr %332, i32 0, i32 93
  %334 = getelementptr inbounds [2048 x i8], ptr %333, i64 0, i64 0
  %335 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %334, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1652, ptr noundef %335)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %46) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %45) #6
  br label %567

336:                                              ; preds = %318
  store i32 0, ptr %19, align 4, !tbaa !9
  br label %371

337:                                              ; preds = %202, %199
  %338 = getelementptr inbounds [8192 x i8], ptr %9, i64 0, i64 0
  %339 = load ptr, ptr %13, align 8, !tbaa !5
  %340 = load i32, ptr %19, align 4, !tbaa !9
  %341 = load i32, ptr %11, align 4, !tbaa !9
  %342 = srem i32 %340, %341
  %343 = sext i32 %342 to i64
  %344 = getelementptr inbounds double, ptr %339, i64 %343
  %345 = load ptr, ptr %14, align 8, !tbaa !5
  %346 = load i32, ptr %19, align 4, !tbaa !9
  %347 = load i32, ptr %11, align 4, !tbaa !9
  %348 = sdiv i32 %346, %347
  %349 = sext i32 %348 to i64
  %350 = getelementptr inbounds double, ptr %345, i64 %349
  %351 = load ptr, ptr %15, align 8, !tbaa !5
  %352 = load i32, ptr %19, align 4, !tbaa !9
  %353 = sext i32 %352 to i64
  %354 = getelementptr inbounds double, ptr %351, i64 %353
  %355 = call i32 (ptr, ptr, ...) @__isoc99_sscanf(ptr noundef %338, ptr noundef @.str.119, ptr noundef %344, ptr noundef %350, ptr noundef %354) #6
  %356 = icmp ne i32 %355, 3
  br i1 %356, label %357, label %368

357:                                              ; preds = %337
  call void @llvm.lifetime.start.p0(i64 2048, ptr %47) #6
  %358 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  %359 = load ptr, ptr %5, align 8, !tbaa !5
  %360 = getelementptr inbounds %struct.precision, ptr %359, i32 0, i32 13
  %361 = getelementptr inbounds [256 x i8], ptr %360, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %358, ptr noundef @.str.120, ptr noundef %361)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %48) #6
  %362 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  %363 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %362, ptr noundef @.str.10, ptr noundef @.str.121, ptr noundef %363)
  %364 = load ptr, ptr %7, align 8, !tbaa !5
  %365 = getelementptr inbounds %struct.thermo, ptr %364, i32 0, i32 93
  %366 = getelementptr inbounds [2048 x i8], ptr %365, i64 0, i64 0
  %367 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %366, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1665, ptr noundef %367)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %48) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %47) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

368:                                              ; preds = %337
  %369 = load i32, ptr %19, align 4, !tbaa !9
  %370 = add nsw i32 %369, 1
  store i32 %370, ptr %19, align 4, !tbaa !9
  br label %371

371:                                              ; preds = %368, %336
  br label %372

372:                                              ; preds = %371, %193
  br label %177, !llvm.loop !25

373:                                              ; preds = %177
  %374 = load ptr, ptr %8, align 8, !tbaa !5
  %375 = call i32 @fclose(ptr noundef %374)
  %376 = load ptr, ptr %14, align 8, !tbaa !5
  %377 = load i32, ptr %12, align 4, !tbaa !9
  %378 = load ptr, ptr %15, align 8, !tbaa !5
  %379 = load i32, ptr %11, align 4, !tbaa !9
  %380 = load ptr, ptr %16, align 8, !tbaa !5
  %381 = load ptr, ptr %7, align 8, !tbaa !5
  %382 = getelementptr inbounds %struct.thermo, ptr %381, i32 0, i32 93
  %383 = getelementptr inbounds [2048 x i8], ptr %382, i64 0, i64 0
  %384 = call i32 @array_spline_table_lines(ptr noundef %376, i32 noundef %377, ptr noundef %378, i32 noundef %379, ptr noundef %380, i16 noundef signext 0, ptr noundef %383)
  %385 = icmp eq i32 %384, 1
  br i1 %385, label %386, label %395

386:                                              ; preds = %373
  call void @llvm.lifetime.start.p0(i64 2048, ptr %49) #6
  %387 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  %388 = load ptr, ptr %7, align 8, !tbaa !5
  %389 = getelementptr inbounds %struct.thermo, ptr %388, i32 0, i32 93
  %390 = getelementptr inbounds [2048 x i8], ptr %389, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %387, ptr noundef @.str, ptr noundef @.str.122, ptr noundef %390)
  %391 = load ptr, ptr %7, align 8, !tbaa !5
  %392 = getelementptr inbounds %struct.thermo, ptr %391, i32 0, i32 93
  %393 = getelementptr inbounds [2048 x i8], ptr %392, i64 0, i64 0
  %394 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %393, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1682, ptr noundef %394)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %49) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

395:                                              ; preds = %373
  %396 = load ptr, ptr %6, align 8, !tbaa !5
  %397 = getelementptr inbounds %struct.background, ptr %396, i32 0, i32 3
  %398 = load double, ptr %397, align 8, !tbaa !26
  %399 = load ptr, ptr %6, align 8, !tbaa !5
  %400 = getelementptr inbounds %struct.background, ptr %399, i32 0, i32 42
  %401 = load double, ptr %400, align 8, !tbaa !27
  %402 = fmul double %398, %401
  %403 = load ptr, ptr %6, align 8, !tbaa !5
  %404 = getelementptr inbounds %struct.background, ptr %403, i32 0, i32 42
  %405 = load double, ptr %404, align 8, !tbaa !27
  %406 = fmul double %402, %405
  store double %406, ptr %21, align 8, !tbaa !16
  %407 = load double, ptr %21, align 8, !tbaa !16
  %408 = load ptr, ptr %13, align 8, !tbaa !5
  %409 = getelementptr inbounds double, ptr %408, i64 0
  %410 = load double, ptr %409, align 8, !tbaa !16
  %411 = fcmp olt double %407, %410
  br i1 %411, label %412, label %427

412:                                              ; preds = %395
  call void @llvm.lifetime.start.p0(i64 2048, ptr %50) #6
  %413 = getelementptr inbounds [2048 x i8], ptr %50, i64 0, i64 0
  %414 = load double, ptr %21, align 8, !tbaa !16
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %413, ptr noundef @.str.123, double noundef %414)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %51) #6
  %415 = getelementptr inbounds [2048 x i8], ptr %51, i64 0, i64 0
  %416 = getelementptr inbounds [2048 x i8], ptr %50, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %415, ptr noundef @.str.10, ptr noundef @.str.124, ptr noundef %416)
  %417 = load ptr, ptr %7, align 8, !tbaa !5
  %418 = getelementptr inbounds %struct.thermo, ptr %417, i32 0, i32 93
  %419 = getelementptr inbounds [2048 x i8], ptr %418, i64 0, i64 0
  %420 = getelementptr inbounds [2048 x i8], ptr %51, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %419, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1690, ptr noundef %420)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %51) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %50) #6
  %421 = load ptr, ptr %13, align 8, !tbaa !5
  call void @free(ptr noundef %421) #6
  %422 = load ptr, ptr %14, align 8, !tbaa !5
  call void @free(ptr noundef %422) #6
  %423 = load ptr, ptr %15, align 8, !tbaa !5
  call void @free(ptr noundef %423) #6
  %424 = load ptr, ptr %16, align 8, !tbaa !5
  call void @free(ptr noundef %424) #6
  %425 = load ptr, ptr %17, align 8, !tbaa !5
  call void @free(ptr noundef %425) #6
  %426 = load ptr, ptr %18, align 8, !tbaa !5
  call void @free(ptr noundef %426) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

427:                                              ; preds = %395
  %428 = load double, ptr %21, align 8, !tbaa !16
  %429 = load ptr, ptr %13, align 8, !tbaa !5
  %430 = load i32, ptr %11, align 4, !tbaa !9
  %431 = sub nsw i32 %430, 1
  %432 = sext i32 %431 to i64
  %433 = getelementptr inbounds double, ptr %429, i64 %432
  %434 = load double, ptr %433, align 8, !tbaa !16
  %435 = fcmp ogt double %428, %434
  br i1 %435, label %436, label %451

436:                                              ; preds = %427
  call void @llvm.lifetime.start.p0(i64 2048, ptr %52) #6
  %437 = getelementptr inbounds [2048 x i8], ptr %52, i64 0, i64 0
  %438 = load double, ptr %21, align 8, !tbaa !16
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %437, ptr noundef @.str.125, double noundef %438)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %53) #6
  %439 = getelementptr inbounds [2048 x i8], ptr %53, i64 0, i64 0
  %440 = getelementptr inbounds [2048 x i8], ptr %52, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %439, ptr noundef @.str.10, ptr noundef @.str.126, ptr noundef %440)
  %441 = load ptr, ptr %7, align 8, !tbaa !5
  %442 = getelementptr inbounds %struct.thermo, ptr %441, i32 0, i32 93
  %443 = getelementptr inbounds [2048 x i8], ptr %442, i64 0, i64 0
  %444 = getelementptr inbounds [2048 x i8], ptr %53, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %443, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1696, ptr noundef %444)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %53) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %52) #6
  %445 = load ptr, ptr %13, align 8, !tbaa !5
  call void @free(ptr noundef %445) #6
  %446 = load ptr, ptr %14, align 8, !tbaa !5
  call void @free(ptr noundef %446) #6
  %447 = load ptr, ptr %15, align 8, !tbaa !5
  call void @free(ptr noundef %447) #6
  %448 = load ptr, ptr %16, align 8, !tbaa !5
  call void @free(ptr noundef %448) #6
  %449 = load ptr, ptr %17, align 8, !tbaa !5
  call void @free(ptr noundef %449) #6
  %450 = load ptr, ptr %18, align 8, !tbaa !5
  call void @free(ptr noundef %450) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

451:                                              ; preds = %427
  %452 = load double, ptr %20, align 8, !tbaa !16
  %453 = load ptr, ptr %14, align 8, !tbaa !5
  %454 = getelementptr inbounds double, ptr %453, i64 0
  %455 = load double, ptr %454, align 8, !tbaa !16
  %456 = fcmp olt double %452, %455
  br i1 %456, label %457, label %472

457:                                              ; preds = %451
  call void @llvm.lifetime.start.p0(i64 2048, ptr %54) #6
  %458 = getelementptr inbounds [2048 x i8], ptr %54, i64 0, i64 0
  %459 = load double, ptr %20, align 8, !tbaa !16
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %458, ptr noundef @.str.127, double noundef %459)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %55) #6
  %460 = getelementptr inbounds [2048 x i8], ptr %55, i64 0, i64 0
  %461 = getelementptr inbounds [2048 x i8], ptr %54, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %460, ptr noundef @.str.10, ptr noundef @.str.128, ptr noundef %461)
  %462 = load ptr, ptr %7, align 8, !tbaa !5
  %463 = getelementptr inbounds %struct.thermo, ptr %462, i32 0, i32 93
  %464 = getelementptr inbounds [2048 x i8], ptr %463, i64 0, i64 0
  %465 = getelementptr inbounds [2048 x i8], ptr %55, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %464, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1702, ptr noundef %465)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %55) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %54) #6
  %466 = load ptr, ptr %13, align 8, !tbaa !5
  call void @free(ptr noundef %466) #6
  %467 = load ptr, ptr %14, align 8, !tbaa !5
  call void @free(ptr noundef %467) #6
  %468 = load ptr, ptr %15, align 8, !tbaa !5
  call void @free(ptr noundef %468) #6
  %469 = load ptr, ptr %16, align 8, !tbaa !5
  call void @free(ptr noundef %469) #6
  %470 = load ptr, ptr %17, align 8, !tbaa !5
  call void @free(ptr noundef %470) #6
  %471 = load ptr, ptr %18, align 8, !tbaa !5
  call void @free(ptr noundef %471) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

472:                                              ; preds = %451
  %473 = load double, ptr %20, align 8, !tbaa !16
  %474 = load ptr, ptr %14, align 8, !tbaa !5
  %475 = load i32, ptr %12, align 4, !tbaa !9
  %476 = sub nsw i32 %475, 1
  %477 = sext i32 %476 to i64
  %478 = getelementptr inbounds double, ptr %474, i64 %477
  %479 = load double, ptr %478, align 8, !tbaa !16
  %480 = fcmp ogt double %473, %479
  br i1 %480, label %481, label %496

481:                                              ; preds = %472
  call void @llvm.lifetime.start.p0(i64 2048, ptr %56) #6
  %482 = getelementptr inbounds [2048 x i8], ptr %56, i64 0, i64 0
  %483 = load double, ptr %20, align 8, !tbaa !16
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %482, ptr noundef @.str.129, double noundef %483)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %57) #6
  %484 = getelementptr inbounds [2048 x i8], ptr %57, 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.10, ptr noundef @.str.130, ptr noundef %485)
  %486 = load ptr, ptr %7, align 8, !tbaa !5
  %487 = getelementptr inbounds %struct.thermo, ptr %486, i32 0, i32 93
  %488 = getelementptr inbounds [2048 x i8], ptr %487, i64 0, i64 0
  %489 = getelementptr inbounds [2048 x i8], ptr %57, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %488, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1708, ptr noundef %489)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %57) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %56) #6
  %490 = load ptr, ptr %13, align 8, !tbaa !5
  call void @free(ptr noundef %490) #6
  %491 = load ptr, ptr %14, align 8, !tbaa !5
  call void @free(ptr noundef %491) #6
  %492 = load ptr, ptr %15, align 8, !tbaa !5
  call void @free(ptr noundef %492) #6
  %493 = load ptr, ptr %16, align 8, !tbaa !5
  call void @free(ptr noundef %493) #6
  %494 = load ptr, ptr %17, align 8, !tbaa !5
  call void @free(ptr noundef %494) #6
  %495 = load ptr, ptr %18, align 8, !tbaa !5
  call void @free(ptr noundef %495) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

496:                                              ; preds = %472
  %497 = load ptr, ptr %14, align 8, !tbaa !5
  %498 = load i32, ptr %12, align 4, !tbaa !9
  %499 = load ptr, ptr %15, align 8, !tbaa !5
  %500 = load ptr, ptr %16, align 8, !tbaa !5
  %501 = load i32, ptr %11, align 4, !tbaa !9
  %502 = load double, ptr %20, align 8, !tbaa !16
  %503 = load ptr, ptr %17, align 8, !tbaa !5
  %504 = load i32, ptr %11, align 4, !tbaa !9
  %505 = load ptr, ptr %7, align 8, !tbaa !5
  %506 = getelementptr inbounds %struct.thermo, ptr %505, i32 0, i32 93
  %507 = getelementptr inbounds [2048 x i8], ptr %506, i64 0, i64 0
  %508 = call i32 @array_interpolate_spline(ptr noundef %497, i32 noundef %498, ptr noundef %499, ptr noundef %500, i32 noundef %501, double noundef %502, ptr noundef %22, ptr noundef %503, i32 noundef %504, ptr noundef %507)
  %509 = icmp eq i32 %508, 1
  br i1 %509, label %510, label %519

510:                                              ; preds = %496
  call void @llvm.lifetime.start.p0(i64 2048, ptr %58) #6
  %511 = getelementptr inbounds [2048 x i8], ptr %58, i64 0, i64 0
  %512 = load ptr, ptr %7, align 8, !tbaa !5
  %513 = getelementptr inbounds %struct.thermo, ptr %512, i32 0, i32 93
  %514 = getelementptr inbounds [2048 x i8], ptr %513, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %511, ptr noundef @.str, ptr noundef @.str.131, ptr noundef %514)
  %515 = load ptr, ptr %7, align 8, !tbaa !5
  %516 = getelementptr inbounds %struct.thermo, ptr %515, i32 0, i32 93
  %517 = getelementptr inbounds [2048 x i8], ptr %516, i64 0, i64 0
  %518 = getelementptr inbounds [2048 x i8], ptr %58, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %517, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1722, ptr noundef %518)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %58) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

519:                                              ; preds = %496
  %520 = load ptr, ptr %13, align 8, !tbaa !5
  %521 = load i32, ptr %11, align 4, !tbaa !9
  %522 = load ptr, ptr %17, align 8, !tbaa !5
  %523 = load ptr, ptr %18, align 8, !tbaa !5
  %524 = load ptr, ptr %7, align 8, !tbaa !5
  %525 = getelementptr inbounds %struct.thermo, ptr %524, i32 0, i32 93
  %526 = getelementptr inbounds [2048 x i8], ptr %525, i64 0, i64 0
  %527 = call i32 @array_spline_table_lines(ptr noundef %520, i32 noundef %521, ptr noundef %522, i32 noundef 1, ptr noundef %523, i16 noundef signext 0, ptr noundef %526)
  %528 = icmp eq i32 %527, 1
  br i1 %528, label %529, label %538

529:                                              ; preds = %519
  call void @llvm.lifetime.start.p0(i64 2048, ptr %59) #6
  %530 = getelementptr inbounds [2048 x i8], ptr %59, i64 0, i64 0
  %531 = load ptr, ptr %7, align 8, !tbaa !5
  %532 = getelementptr inbounds %struct.thermo, ptr %531, i32 0, i32 93
  %533 = getelementptr inbounds [2048 x i8], ptr %532, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %530, ptr noundef @.str, ptr noundef @.str.132, ptr noundef %533)
  %534 = load ptr, ptr %7, align 8, !tbaa !5
  %535 = getelementptr inbounds %struct.thermo, ptr %534, i32 0, i32 93
  %536 = getelementptr inbounds [2048 x i8], ptr %535, i64 0, i64 0
  %537 = getelementptr inbounds [2048 x i8], ptr %59, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %536, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1733, ptr noundef %537)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %59) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

538:                                              ; preds = %519
  %539 = load ptr, ptr %13, align 8, !tbaa !5
  %540 = load i32, ptr %11, align 4, !tbaa !9
  %541 = load ptr, ptr %17, align 8, !tbaa !5
  %542 = load ptr, ptr %18, align 8, !tbaa !5
  %543 = load double, ptr %21, align 8, !tbaa !16
  %544 = load ptr, ptr %7, align 8, !tbaa !5
  %545 = getelementptr inbounds %struct.thermo, ptr %544, i32 0, i32 0
  %546 = load ptr, ptr %7, align 8, !tbaa !5
  %547 = getelementptr inbounds %struct.thermo, ptr %546, i32 0, i32 93
  %548 = getelementptr inbounds [2048 x i8], ptr %547, i64 0, i64 0
  %549 = call i32 @array_interpolate_spline(ptr noundef %539, i32 noundef %540, ptr noundef %541, ptr noundef %542, i32 noundef 1, double noundef %543, ptr noundef %22, ptr noundef %545, i32 noundef 1, ptr noundef %548)
  %550 = icmp eq i32 %549, 1
  br i1 %550, label %551, label %560

551:                                              ; preds = %538
  call void @llvm.lifetime.start.p0(i64 2048, ptr %60) #6
  %552 = getelementptr inbounds [2048 x i8], ptr %60, i64 0, i64 0
  %553 = load ptr, ptr %7, align 8, !tbaa !5
  %554 = getelementptr inbounds %struct.thermo, ptr %553, i32 0, i32 93
  %555 = getelementptr inbounds [2048 x i8], ptr %554, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %552, ptr noundef @.str, ptr noundef @.str.133, ptr noundef %555)
  %556 = load ptr, ptr %7, align 8, !tbaa !5
  %557 = getelementptr inbounds %struct.thermo, ptr %556, i32 0, i32 93
  %558 = getelementptr inbounds [2048 x i8], ptr %557, i64 0, i64 0
  %559 = getelementptr inbounds [2048 x i8], ptr %60, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %558, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_helium_from_bbn, i32 noundef 1747, ptr noundef %559)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %60) #6
  store i32 1, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

560:                                              ; preds = %538
  %561 = load ptr, ptr %13, align 8, !tbaa !5
  call void @free(ptr noundef %561) #6
  %562 = load ptr, ptr %14, align 8, !tbaa !5
  call void @free(ptr noundef %562) #6
  %563 = load ptr, ptr %15, align 8, !tbaa !5
  call void @free(ptr noundef %563) #6
  %564 = load ptr, ptr %16, align 8, !tbaa !5
  call void @free(ptr noundef %564) #6
  %565 = load ptr, ptr %17, align 8, !tbaa !5
  call void @free(ptr noundef %565) #6
  %566 = load ptr, ptr %18, align 8, !tbaa !5
  call void @free(ptr noundef %566) #6
  store i32 0, ptr %4, align 4
  store i32 1, ptr %29, align 4
  br label %567

567:                                              ; preds = %560, %551, %529, %510, %481, %457, %436, %412, %386, %357, %325, %307, %287, %265, %245, %227, %209, %167, %114, %93, %69
  call void @llvm.lifetime.end.p0(i64 8, ptr %26) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %25) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %22) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %20) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %19) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %18) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %17) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %14) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %13) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %12) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %11) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %10) #6
  call void @llvm.lifetime.end.p0(i64 8192, ptr %9) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %8) #6
  %568 = load i32, ptr %4, align 4
  ret i32 %568
}

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

; 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

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

declare noalias ptr @fopen(ptr noundef, ptr noundef) #3

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

; Function Attrs: nounwind
declare i32 @__isoc99_sscanf(ptr noundef, ptr noundef, ...) #1

declare i32 @fclose(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 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 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 #6 = { nounwind }
attributes #7 = { 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, !10, i64 0}
!10 = !{!"int", !7, i64 0}
!11 = !{!12, !10, i64 624}
!12 = !{!"background", !13, i64 0, !13, i64 8, !13, i64 16, !13, i64 24, !13, i64 32, !13, i64 40, !13, i64 48, !10, i64 56, !13, i64 64, !13, i64 72, !13, i64 80, !13, i64 88, !14, i64 96, !13, i64 104, !13, i64 112, !13, i64 120, !13, i64 128, !13, i64 136, !13, i64 144, !13, i64 152, !13, i64 160, !13, i64 168, !14, i64 176, !13, i64 184, !13, i64 192, !6, i64 200, !10, i64 208, !10, i64 212, !13, i64 216, !10, i64 224, !6, i64 232, !6, i64 240, !13, i64 248, !6, i64 256, !13, i64 264, !6, i64 272, !13, i64 280, !6, i64 288, !13, i64 296, !6, i64 304, !6, i64 312, !6, i64 320, !13, i64 328, !13, i64 336, !13, i64 344, !13, i64 352, !10, i64 360, !6, i64 368, !13, i64 376, !13, i64 384, !13, i64 392, !13, i64 400, !13, i64 408, !13, i64 416, !13, i64 424, !13, i64 432, !13, i64 440, !13, i64 448, !13, 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, !6, i64 632, !6, i64 640, !6, i64 648, !6, i64 656, !6, 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, !14, i64 724, !14, i64 726, !14, i64 728, !14, i64 730, !14, i64 732, !14, i64 734, !14, i64 736, !14, i64 738, !14, i64 740, !14, i64 742, !14, 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, !14, i64 840, !14, i64 842, !14, i64 844, !14, i64 846, !14, i64 848, !14, i64 850, !7, i64 852, !14, i64 2900, !7, i64 2902}
!13 = !{!"double", !7, i64 0}
!14 = !{!"short", !7, i64 0}
!15 = !{!12, !13, i64 16}
!16 = !{!13, !13, i64 0}
!17 = !{!12, !14, i64 844}
!18 = !{!12, !14, i64 846}
!19 = !{!12, !10, i64 576}
!20 = !{!12, !10, i64 580}
!21 = !{!12, !10, i64 476}
!22 = !{!7, !7, i64 0}
!23 = distinct !{!23, !24}
!24 = !{!"llvm.loop.mustprogress"}
!25 = distinct !{!25, !24}
!26 = !{!12, !13, i64 24}
!27 = !{!12, !13, i64 328}


; 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.37 = external hidden unnamed_addr constant [35 x i8], align 1
@.str.38 = external hidden unnamed_addr constant [9 x i8], align 1
@__func__.thermodynamics_helium_from_bbn = external hidden unnamed_addr constant [31 x i8], align 1
@.str.104 = external hidden unnamed_addr constant [42 x i8], align 1
@.str.105 = external hidden unnamed_addr constant [90 x i8], align 1
@.str.106 = external hidden unnamed_addr constant [2 x i8], align 1
@.str.107 = external hidden unnamed_addr constant [43 x i8], align 1
@.str.108 = external hidden unnamed_addr constant [3 x i8], align 1
@.str.109 = external hidden unnamed_addr constant [4 x i8], align 1
@.str.110 = external hidden unnamed_addr constant [6 x i8], align 1
@.str.111 = external hidden unnamed_addr constant [71 x i8], align 1
@.str.112 = external hidden unnamed_addr constant [50 x i8], align 1
@.str.113 = external hidden unnamed_addr constant [7 x i8], align 1
@.str.114 = external hidden unnamed_addr constant [7 x i8], align 1
@.str.115 = external hidden unnamed_addr constant [4 x i8], align 1
@.str.116 = external hidden unnamed_addr constant [6 x i8], align 1
@.str.117 = external hidden unnamed_addr constant [14 x i8], align 1
@.str.118 = external hidden unnamed_addr constant [16 x i8], align 1
@.str.119 = external hidden unnamed_addr constant [12 x i8], align 1
@.str.120 = external hidden unnamed_addr constant [67 x i8], align 1
@.str.121 = external hidden unnamed_addr constant [121 x i8], align 1
@.str.122 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.123 = external hidden unnamed_addr constant [248 x i8], align 1
@.str.124 = external hidden unnamed_addr constant [20 x i8], align 1
@.str.125 = external hidden unnamed_addr constant [247 x i8], align 1
@.str.126 = external hidden unnamed_addr constant [31 x i8], align 1
@.str.127 = external hidden unnamed_addr constant [255 x i8], align 1
@.str.128 = external hidden unnamed_addr constant [22 x i8], align 1
@.str.129 = external hidden unnamed_addr constant [254 x i8], align 1
@.str.130 = external hidden unnamed_addr constant [33 x i8], align 1
@.str.131 = external hidden unnamed_addr constant [140 x i8], align 1
@.str.132 = external hidden unnamed_addr constant [103 x i8], align 1
@.str.133 = external hidden unnamed_addr constant [138 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 nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.fmuladd.f64(double, double, double) #1

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

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

declare i32 @array_interpolate_spline(ptr noundef, i32 noundef, ptr noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, ptr noundef) local_unnamed_addr #2

; Function Attrs: nounwind uwtable
define dso_local range(i32 0, 2) i32 @thermodynamics_helium_from_bbn(ptr noundef %0, ptr noundef %1, ptr noundef %2) local_unnamed_addr #3 {
  %4 = alloca [8192 x i8], align 16
  %5 = alloca i32, align 4
  %6 = alloca i32, align 4
  %7 = alloca i32, align 4
  %8 = alloca double, align 8
  %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
  %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
  call void @llvm.lifetime.start.p0(i64 8192, ptr nonnull %4) #7
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %5) #7
  store i32 0, ptr %5, align 4, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %6) #7
  store i32 0, ptr %6, align 4, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %7) #7
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %8) #7
  %35 = getelementptr inbounds i8, ptr %1, i64 624
  %36 = load i32, ptr %35, align 8, !tbaa !9
  %37 = sext i32 %36 to i64
  %38 = shl nsw i64 %37, 3
  %39 = tail call noalias ptr @malloc(i64 noundef %38) #8
  %40 = icmp eq ptr %39, null
  br i1 %40, label %41, label %44

41:                                               ; preds = %3
  %42 = shl i32 %36, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %9) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %9, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.38, i32 noundef %42) #7
  %43 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %43, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1575, ptr noundef nonnull %9) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %9) #7
  br label %221

44:                                               ; preds = %3
  %45 = getelementptr inbounds i8, ptr %1, i64 16
  %46 = load double, ptr %45, align 8, !tbaa !14
  %47 = fmul double %46, 0x3DD7AFE2C34414DC
  %48 = fdiv double 1.000000e-01, %47
  %49 = fadd double %48, -1.000000e+00
  %50 = call i32 @background_tau_of_z(ptr noundef nonnull %1, double noundef %49, ptr noundef nonnull %8) #7
  %51 = icmp eq i32 %50, 1
  br i1 %51, label %52, label %55

52:                                               ; preds = %44
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %10) #7
  %53 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %10, ptr noundef nonnull @.str, ptr noundef nonnull @.str.104, ptr noundef nonnull %53) #7
  %54 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %54, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1584, ptr noundef nonnull %10) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %10) #7
  br label %221

55:                                               ; preds = %44
  %56 = load double, ptr %8, align 8, !tbaa !15
  %57 = getelementptr inbounds i8, ptr %1, i64 844
  %58 = load i16, ptr %57, align 4, !tbaa !16
  %59 = getelementptr inbounds i8, ptr %1, i64 846
  %60 = load i16, ptr %59, align 2, !tbaa !17
  %61 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %56, i16 noundef signext %58, i16 noundef signext %60, ptr noundef nonnull %7, ptr noundef nonnull %39) #7
  %62 = icmp eq i32 %61, 1
  br i1 %62, label %63, label %66

63:                                               ; preds = %55
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %11) #7
  %64 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %11, ptr noundef nonnull @.str, ptr noundef nonnull @.str.105, ptr noundef nonnull %64) #7
  %65 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %65, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1593, ptr noundef nonnull %11) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %11) #7
  br label %221

66:                                               ; preds = %55
  %67 = getelementptr inbounds i8, ptr %1, i64 576
  %68 = load i32, ptr %67, align 8, !tbaa !18
  %69 = sext i32 %68 to i64
  %70 = getelementptr inbounds double, ptr %39, i64 %69
  %71 = load double, ptr %70, align 8, !tbaa !15
  %72 = getelementptr inbounds i8, ptr %1, i64 580
  %73 = load i32, ptr %72, align 4, !tbaa !19
  %74 = sext i32 %73 to i64
  %75 = getelementptr inbounds double, ptr %39, i64 %74
  %76 = load double, ptr %75, align 8, !tbaa !15
  %77 = getelementptr inbounds i8, ptr %1, i64 476
  %78 = load i32, ptr %77, align 4, !tbaa !20
  %79 = sext i32 %78 to i64
  %80 = getelementptr inbounds double, ptr %39, i64 %79
  %81 = load double, ptr %80, align 8, !tbaa !15
  %82 = fneg double %81
  %83 = call double @llvm.fmuladd.f64(double %71, double %76, double %82)
  %84 = fmul double %81, 0x3FCD11DA430439A9
  %85 = fdiv double %83, %84
  call void @free(ptr noundef nonnull %39) #7
  %86 = fadd double %85, -3.046000e+00
  %87 = getelementptr inbounds i8, ptr %0, i64 104
  %88 = call noalias ptr @fopen(ptr noundef nonnull %87, ptr noundef nonnull @.str.106)
  %89 = icmp eq ptr %88, null
  br i1 %89, label %91, label %.preheader170

.preheader170:                                    ; preds = %66
  %90 = call ptr @fgets(ptr noundef nonnull %4, i32 noundef 8191, ptr noundef nonnull %88)
  %.not234 = icmp eq ptr %90, null
  br i1 %.not234, label %._crit_edge, label %.preheader

91:                                               ; preds = %66
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %12) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %12, ptr noundef nonnull @.str.107, ptr noundef nonnull @.str.108, ptr noundef nonnull %87, ptr noundef nonnull @.str.109) #7
  %92 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %92, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1619, ptr noundef nonnull %12) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %12) #7
  br label %221

.preheader:                                       ; preds = %.preheader170, %167
  %.0143241 = phi ptr [ %.1, %167 ], [ null, %.preheader170 ]
  %.0144240 = phi ptr [ %.1145, %167 ], [ null, %.preheader170 ]
  %.0146239 = phi ptr [ %.1147, %167 ], [ null, %.preheader170 ]
  %.0148238 = phi ptr [ %.1149, %167 ], [ null, %.preheader170 ]
  %.0150237 = phi ptr [ %.1151, %167 ], [ null, %.preheader170 ]
  %.0152236 = phi ptr [ %.1153, %167 ], [ null, %.preheader170 ]
  %.0154235 = phi i32 [ %.1155, %167 ], [ 0, %.preheader170 ]
  br label %93

93:                                               ; preds = %93, %.preheader
  %.0142 = phi ptr [ %96, %93 ], [ %4, %.preheader ]
  %94 = load i8, ptr %.0142, align 1, !tbaa !21
  %95 = icmp eq i8 %94, 32
  %96 = getelementptr inbounds i8, ptr %.0142, i64 1
  br i1 %95, label %93, label %97, !llvm.loop !22

97:                                               ; preds = %93
  %98 = icmp sgt i8 %94, 39
  br i1 %98, label %99, label %167

99:                                               ; preds = %97
  %100 = load i32, ptr %5, align 4, !tbaa !5
  %101 = icmp eq i32 %100, 0
  %102 = load i32, ptr %6, align 4
  %103 = icmp eq i32 %102, 0
  %or.cond = select i1 %101, i1 %103, i1 false
  br i1 %or.cond, label %104, label %153

104:                                              ; preds = %99
  %105 = call i32 (ptr, ptr, ...) @__isoc99_sscanf(ptr noundef nonnull %4, ptr noundef nonnull @.str.110, ptr noundef nonnull %5, ptr noundef nonnull %6) #7
  %.not169 = icmp eq i32 %105, 2
  br i1 %.not169, label %108, label %106

106:                                              ; preds = %104
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %13) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %13, ptr noundef nonnull @.str.111, ptr noundef nonnull %87) #7
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %14) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %14, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.112, ptr noundef nonnull %13) #7
  %107 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %107, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1645, ptr noundef nonnull %14) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %14) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %13) #7
  br label %221

108:                                              ; preds = %104
  %109 = load i32, ptr %5, align 4, !tbaa !5
  %110 = sext i32 %109 to i64
  %111 = shl nsw i64 %110, 3
  %112 = call noalias ptr @malloc(i64 noundef %111) #8
  %113 = icmp eq ptr %112, null
  br i1 %113, label %114, label %117

114:                                              ; preds = %108
  %115 = shl i32 %109, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %15) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %15, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.113, i32 noundef %115) #7
  %116 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %116, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1647, ptr noundef nonnull %15) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %15) #7
  br label %221

117:                                              ; preds = %108
  %118 = load i32, ptr %6, align 4, !tbaa !5
  %119 = sext i32 %118 to i64
  %120 = shl nsw i64 %119, 3
  %121 = call noalias ptr @malloc(i64 noundef %120) #8
  %122 = icmp eq ptr %121, null
  br i1 %122, label %123, label %126

123:                                              ; preds = %117
  %124 = shl i32 %118, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %16) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %16, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.114, i32 noundef %124) #7
  %125 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %125, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1648, ptr noundef nonnull %16) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %16) #7
  br label %221

126:                                              ; preds = %117
  %127 = mul nsw i32 %118, %109
  %128 = sext i32 %127 to i64
  %129 = shl nsw i64 %128, 3
  %130 = call noalias ptr @malloc(i64 noundef %129) #8
  %131 = icmp eq ptr %130, null
  br i1 %131, label %132, label %135

132:                                              ; preds = %126
  %133 = shl i32 %127, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %17) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %17, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.115, i32 noundef %133) #7
  %134 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %134, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1649, ptr noundef nonnull %17) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %17) #7
  br label %221

135:                                              ; preds = %126
  %136 = call noalias ptr @malloc(i64 noundef %129) #8
  %137 = icmp eq ptr %136, null
  br i1 %137, label %138, label %141

138:                                              ; preds = %135
  %139 = shl i32 %127, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %18) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %18, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.116, i32 noundef %139) #7
  %140 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %140, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1650, ptr noundef nonnull %18) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %18) #7
  br label %221

141:                                              ; preds = %135
  %142 = call noalias ptr @malloc(i64 noundef %111) #8
  %143 = icmp eq ptr %142, null
  br i1 %143, label %144, label %147

144:                                              ; preds = %141
  %145 = shl i32 %109, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %19) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %19, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.117, i32 noundef %145) #7
  %146 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %146, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1651, ptr noundef nonnull %19) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %19) #7
  br label %221

147:                                              ; preds = %141
  %148 = call noalias ptr @malloc(i64 noundef %111) #8
  %149 = icmp eq ptr %148, null
  br i1 %149, label %150, label %167

150:                                              ; preds = %147
  %151 = shl i32 %109, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %20) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %20, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.118, i32 noundef %151) #7
  %152 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %152, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1652, ptr noundef nonnull %20) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %20) #7
  br label %221

153:                                              ; preds = %99
  %154 = srem i32 %.0154235, %100
  %155 = sext i32 %154 to i64
  %156 = getelementptr inbounds double, ptr %.0143241, i64 %155
  %157 = sdiv i32 %.0154235, %100
  %158 = sext i32 %157 to i64
  %159 = getelementptr inbounds double, ptr %.0144240, i64 %158
  %160 = sext i32 %.0154235 to i64
  %161 = getelementptr inbounds double, ptr %.0146239, i64 %160
  %162 = call i32 (ptr, ptr, ...) @__isoc99_sscanf(ptr noundef nonnull %4, ptr noundef nonnull @.str.119, ptr noundef %156, ptr noundef %159, ptr noundef %161) #7
  %.not168 = icmp eq i32 %162, 3
  br i1 %.not168, label %165, label %163

163:                                              ; preds = %153
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %21) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %21, ptr noundef nonnull @.str.120, ptr noundef nonnull %87) #7
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %22) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %22, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.121, ptr noundef nonnull %21) #7
  %164 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %164, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1665, ptr noundef nonnull %22) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %22) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %21) #7
  br label %221

165:                                              ; preds = %153
  %166 = add nsw i32 %.0154235, 1
  br label %167

167:                                              ; preds = %147, %165, %97
  %.1155 = phi i32 [ %166, %165 ], [ %.0154235, %97 ], [ 0, %147 ]
  %.1153 = phi ptr [ %.0152236, %165 ], [ %.0152236, %97 ], [ %148, %147 ]
  %.1151 = phi ptr [ %.0150237, %165 ], [ %.0150237, %97 ], [ %142, %147 ]
  %.1149 = phi ptr [ %.0148238, %165 ], [ %.0148238, %97 ], [ %136, %147 ]
  %.1147 = phi ptr [ %.0146239, %165 ], [ %.0146239, %97 ], [ %130, %147 ]
  %.1145 = phi ptr [ %.0144240, %165 ], [ %.0144240, %97 ], [ %121, %147 ]
  %.1 = phi ptr [ %.0143241, %165 ], [ %.0143241, %97 ], [ %112, %147 ]
  %168 = call ptr @fgets(ptr noundef nonnull %4, i32 noundef 8191, ptr noundef nonnull %88)
  %.not = icmp eq ptr %168, null
  br i1 %.not, label %._crit_edge, label %.preheader, !llvm.loop !24

._crit_edge:                                      ; preds = %167, %.preheader170
  %.0152.lcssa = phi ptr [ null, %.preheader170 ], [ %.1153, %167 ]
  %.0150.lcssa = phi ptr [ null, %.preheader170 ], [ %.1151, %167 ]
  %.0148.lcssa = phi ptr [ null, %.preheader170 ], [ %.1149, %167 ]
  %.0146.lcssa = phi ptr [ null, %.preheader170 ], [ %.1147, %167 ]
  %.0144.lcssa = phi ptr [ null, %.preheader170 ], [ %.1145, %167 ]
  %.0143.lcssa = phi ptr [ null, %.preheader170 ], [ %.1, %167 ]
  %169 = call i32 @fclose(ptr noundef nonnull %88)
  %170 = load i32, ptr %6, align 4, !tbaa !5
  %171 = load i32, ptr %5, align 4, !tbaa !5
  %172 = getelementptr inbounds i8, ptr %2, i64 606
  %173 = call i32 @array_spline_table_lines(ptr noundef %.0144.lcssa, i32 noundef %170, ptr noundef %.0146.lcssa, i32 noundef %171, ptr noundef %.0148.lcssa, i16 noundef signext 0, ptr noundef nonnull %172) #7
  %174 = icmp eq i32 %173, 1
  br i1 %174, label %175, label %176

175:                                              ; preds = %._crit_edge
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %23) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %23, ptr noundef nonnull @.str, ptr noundef nonnull @.str.122, ptr noundef nonnull %172) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1682, ptr noundef nonnull %23) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %23) #7
  br label %221

176:                                              ; preds = %._crit_edge
  %177 = getelementptr inbounds i8, ptr %1, i64 24
  %178 = load double, ptr %177, align 8, !tbaa !25
  %179 = getelementptr inbounds i8, ptr %1, i64 328
  %180 = load double, ptr %179, align 8, !tbaa !26
  %181 = fmul double %178, %180
  %182 = fmul double %180, %181
  %183 = load double, ptr %.0143.lcssa, align 8, !tbaa !15
  %184 = fcmp olt double %182, %183
  br i1 %184, label %185, label %186

185:                                              ; preds = %176
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %24) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %24, ptr noundef nonnull @.str.123, double noundef %182) #7
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %25) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %25, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.124, ptr noundef nonnull %24) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1690, ptr noundef nonnull %25) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %25) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %24) #7
  call void @free(ptr noundef nonnull %.0143.lcssa) #7
  call void @free(ptr noundef %.0144.lcssa) #7
  call void @free(ptr noundef %.0146.lcssa) #7
  call void @free(ptr noundef %.0148.lcssa) #7
  call void @free(ptr noundef %.0150.lcssa) #7
  call void @free(ptr noundef %.0152.lcssa) #7
  br label %221

186:                                              ; preds = %176
  %187 = load i32, ptr %5, align 4, !tbaa !5
  %188 = add nsw i32 %187, -1
  %189 = sext i32 %188 to i64
  %190 = getelementptr inbounds double, ptr %.0143.lcssa, i64 %189
  %191 = load double, ptr %190, align 8, !tbaa !15
  %192 = fcmp ogt double %182, %191
  br i1 %192, label %193, label %194

193:                                              ; preds = %186
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %26) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %26, ptr noundef nonnull @.str.125, double noundef %182) #7
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %27) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %27, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.126, ptr noundef nonnull %26) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1696, ptr noundef nonnull %27) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %27) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %26) #7
  call void @free(ptr noundef nonnull %.0143.lcssa) #7
  call void @free(ptr noundef %.0144.lcssa) #7
  call void @free(ptr noundef %.0146.lcssa) #7
  call void @free(ptr noundef %.0148.lcssa) #7
  call void @free(ptr noundef %.0150.lcssa) #7
  call void @free(ptr noundef %.0152.lcssa) #7
  br label %221

194:                                              ; preds = %186
  %195 = load double, ptr %.0144.lcssa, align 8, !tbaa !15
  %196 = fcmp olt double %86, %195
  br i1 %196, label %197, label %198

197:                                              ; preds = %194
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %28) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %28, ptr noundef nonnull @.str.127, double noundef %86) #7
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %29) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %29, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.128, ptr noundef nonnull %28) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1702, ptr noundef nonnull %29) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %29) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %28) #7
  call void @free(ptr noundef nonnull %.0143.lcssa) #7
  call void @free(ptr noundef nonnull %.0144.lcssa) #7
  call void @free(ptr noundef %.0146.lcssa) #7
  call void @free(ptr noundef %.0148.lcssa) #7
  call void @free(ptr noundef %.0150.lcssa) #7
  call void @free(ptr noundef %.0152.lcssa) #7
  br label %221

198:                                              ; preds = %194
  %199 = load i32, ptr %6, align 4, !tbaa !5
  %200 = add nsw i32 %199, -1
  %201 = sext i32 %200 to i64
  %202 = getelementptr inbounds double, ptr %.0144.lcssa, i64 %201
  %203 = load double, ptr %202, align 8, !tbaa !15
  %204 = fcmp ogt double %86, %203
  br i1 %204, label %205, label %206

205:                                              ; preds = %198
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %30) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %30, ptr noundef nonnull @.str.129, double noundef %86) #7
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %31) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %31, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.130, ptr noundef nonnull %30) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1708, ptr noundef nonnull %31) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %31) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %30) #7
  call void @free(ptr noundef nonnull %.0143.lcssa) #7
  call void @free(ptr noundef nonnull %.0144.lcssa) #7
  call void @free(ptr noundef %.0146.lcssa) #7
  call void @free(ptr noundef %.0148.lcssa) #7
  call void @free(ptr noundef %.0150.lcssa) #7
  call void @free(ptr noundef %.0152.lcssa) #7
  br label %221

206:                                              ; preds = %198
  %207 = call i32 @array_interpolate_spline(ptr noundef nonnull %.0144.lcssa, i32 noundef %199, ptr noundef %.0146.lcssa, ptr noundef %.0148.lcssa, i32 noundef %187, double noundef %86, ptr noundef nonnull %7, ptr noundef %.0150.lcssa, i32 noundef %187, ptr noundef nonnull %172) #7
  %208 = icmp eq i32 %207, 1
  br i1 %208, label %209, label %210

209:                                              ; preds = %206
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %32) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %32, ptr noundef nonnull @.str, ptr noundef nonnull @.str.131, ptr noundef nonnull %172) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1722, ptr noundef nonnull %32) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %32) #7
  br label %221

210:                                              ; preds = %206
  %211 = load i32, ptr %5, align 4, !tbaa !5
  %212 = call i32 @array_spline_table_lines(ptr noundef nonnull %.0143.lcssa, i32 noundef %211, ptr noundef %.0150.lcssa, i32 noundef 1, ptr noundef %.0152.lcssa, i16 noundef signext 0, ptr noundef nonnull %172) #7
  %213 = icmp eq i32 %212, 1
  br i1 %213, label %214, label %215

214:                                              ; preds = %210
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %33) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %33, ptr noundef nonnull @.str, ptr noundef nonnull @.str.132, ptr noundef nonnull %172) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1733, ptr noundef nonnull %33) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %33) #7
  br label %221

215:                                              ; preds = %210
  %216 = load i32, ptr %5, align 4, !tbaa !5
  %217 = call i32 @array_interpolate_spline(ptr noundef nonnull %.0143.lcssa, i32 noundef %216, ptr noundef %.0150.lcssa, ptr noundef %.0152.lcssa, i32 noundef 1, double noundef %182, ptr noundef nonnull %7, ptr noundef %2, i32 noundef 1, ptr noundef nonnull %172) #7
  %218 = icmp eq i32 %217, 1
  br i1 %218, label %219, label %220

219:                                              ; preds = %215
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %34) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %34, ptr noundef nonnull @.str, ptr noundef nonnull @.str.133, ptr noundef nonnull %172) #7
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %172, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_helium_from_bbn, i32 noundef 1747, ptr noundef nonnull %34) #7
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %34) #7
  br label %221

220:                                              ; preds = %215
  call void @free(ptr noundef nonnull %.0143.lcssa) #7
  call void @free(ptr noundef nonnull %.0144.lcssa) #7
  call void @free(ptr noundef %.0146.lcssa) #7
  call void @free(ptr noundef %.0148.lcssa) #7
  call void @free(ptr noundef %.0150.lcssa) #7
  call void @free(ptr noundef %.0152.lcssa) #7
  br label %221

221:                                              ; preds = %220, %219, %214, %209, %205, %197, %193, %185, %175, %163, %150, %144, %138, %132, %123, %114, %106, %91, %63, %52, %41
  %.0 = phi i32 [ 1, %41 ], [ 1, %52 ], [ 1, %63 ], [ 1, %91 ], [ 1, %106 ], [ 1, %114 ], [ 1, %123 ], [ 1, %132 ], [ 1, %138 ], [ 1, %144 ], [ 1, %150 ], [ 1, %163 ], [ 1, %175 ], [ 1, %185 ], [ 1, %193 ], [ 1, %197 ], [ 1, %205 ], [ 1, %209 ], [ 1, %214 ], [ 1, %219 ], [ 0, %220 ]
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %8) #7
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %7) #7
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %6) #7
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %5) #7
  call void @llvm.lifetime.end.p0(i64 8192, ptr nonnull %4) #7
  ret i32 %.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 #2

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

declare i32 @array_spline_table_lines(ptr noundef, i32 noundef, ptr noundef, i32 noundef, ptr noundef, i16 noundef signext, ptr noundef) local_unnamed_addr #2

; Function Attrs: nofree nounwind
declare noalias noundef ptr @fopen(ptr nocapture noundef readonly, ptr nocapture noundef readonly) local_unnamed_addr #6

; Function Attrs: nofree nounwind
declare noundef ptr @fgets(ptr noundef, i32 noundef, ptr nocapture noundef) local_unnamed_addr #6

; Function Attrs: nofree nounwind
declare noundef i32 @__isoc99_sscanf(ptr nocapture noundef readonly, ptr nocapture noundef readonly, ...) local_unnamed_addr #6

; Function Attrs: nofree nounwind
declare noundef i32 @fclose(ptr nocapture noundef) local_unnamed_addr #6

attributes #0 = { mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #2 = { "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { 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 #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 = { nofree 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 #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 = !{!"int", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !6, i64 624}
!10 = !{!"background", !11, i64 0, !11, i64 8, !11, i64 16, !11, i64 24, !11, i64 32, !11, i64 40, !11, i64 48, !6, i64 56, !11, i64 64, !11, i64 72, !11, i64 80, !11, i64 88, !12, 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, !12, i64 176, !11, i64 184, !11, i64 192, !13, i64 200, !6, i64 208, !6, i64 212, !11, i64 216, !6, i64 224, !13, i64 232, !13, i64 240, !11, i64 248, !13, i64 256, !11, i64 264, !13, i64 272, !11, i64 280, !13, i64 288, !11, i64 296, !13, i64 304, !13, i64 312, !13, i64 320, !11, i64 328, !11, i64 336, !11, i64 344, !11, i64 352, !6, i64 360, !13, 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, !6, i64 464, !6, i64 468, !6, i64 472, !6, i64 476, !6, i64 480, !6, i64 484, !6, i64 488, !6, i64 492, !6, i64 496, !6, i64 500, !6, i64 504, !6, i64 508, !6, i64 512, !6, i64 516, !6, i64 520, !6, i64 524, !6, i64 528, !6, i64 532, !6, i64 536, !6, i64 540, !6, i64 544, !6, i64 548, !6, i64 552, !6, i64 556, !6, i64 560, !6, i64 564, !6, i64 568, !6, i64 572, !6, i64 576, !6, i64 580, !6, i64 584, !6, i64 588, !6, i64 592, !6, i64 596, !6, i64 600, !6, i64 604, !6, i64 608, !6, i64 612, !6, i64 616, !6, i64 620, !6, i64 624, !6, i64 628, !13, i64 632, !13, i64 640, !13, i64 648, !13, i64 656, !13, i64 664, !6, i64 672, !6, i64 676, !6, i64 680, !6, i64 684, !6, i64 688, !6, i64 692, !6, i64 696, !6, i64 700, !6, i64 704, !6, i64 708, !6, i64 712, !6, i64 716, !6, i64 720, !12, i64 724, !12, i64 726, !12, i64 728, !12, i64 730, !12, i64 732, !12, i64 734, !12, i64 736, !12, i64 738, !12, i64 740, !12, i64 742, !12, i64 744, !13, i64 752, !13, i64 760, !13, i64 768, !13, i64 776, !13, i64 784, !13, i64 792, !13, i64 800, !13, i64 808, !13, i64 816, !13, i64 824, !13, i64 832, !12, i64 840, !12, i64 842, !12, i64 844, !12, i64 846, !12, i64 848, !12, i64 850, !7, i64 852, !12, i64 2900, !7, i64 2902}
!11 = !{!"double", !7, i64 0}
!12 = !{!"short", !7, i64 0}
!13 = !{!"any pointer", !7, i64 0}
!14 = !{!10, !11, i64 16}
!15 = !{!11, !11, i64 0}
!16 = !{!10, !12, i64 844}
!17 = !{!10, !12, i64 846}
!18 = !{!10, !6, i64 576}
!19 = !{!10, !6, i64 580}
!20 = !{!10, !6, i64 476}
!21 = !{!7, !7, i64 0}
!22 = distinct !{!22, !23}
!23 = !{!"llvm.loop.mustprogress"}
!24 = distinct !{!24, !23}
!25 = !{!10, !11, i64 24}
!26 = !{!10, !11, i64 328}


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