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

@.str = external hidden unnamed_addr constant [18 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [13 x i8], align 1
@.str.10 = external hidden unnamed_addr constant [27 x i8], align 1
@.str.37 = external hidden unnamed_addr constant [35 x i8], align 1
@.str.212 = external hidden unnamed_addr constant [27 x i8], align 1
@__func__.thermodynamics_recombination_with_recfast = external hidden unnamed_addr constant [42 x i8], align 1
@.str.213 = external hidden unnamed_addr constant [38 x i8], align 1
@.str.214 = external hidden unnamed_addr constant [40 x i8], align 1
@.str.215 = external hidden unnamed_addr constant [70 x i8], align 1
@.str.216 = external hidden unnamed_addr constant [41 x i8], align 1
@.str.217 = external hidden unnamed_addr constant [59 x i8], align 1
@.str.218 = external hidden unnamed_addr constant [113 x i8], align 1
@.str.219 = external hidden unnamed_addr constant [31 x i8], align 1
@.str.220 = external hidden unnamed_addr constant [146 x i8], align 1
@.str.221 = external hidden unnamed_addr constant [74 x i8], align 1
@.str.222 = external hidden unnamed_addr constant [32 x i8], align 1

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

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

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

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

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

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

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

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

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

; Function Attrs: nounwind uwtable
define dso_local i32 @thermodynamics_recombination_with_recfast(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, ptr noundef %4) #5 {
  %6 = alloca i32, align 4
  %7 = alloca ptr, align 8
  %8 = alloca ptr, align 8
  %9 = alloca ptr, align 8
  %10 = alloca ptr, align 8
  %11 = alloca ptr, align 8
  %12 = alloca [3 x double], align 16
  %13 = alloca [3 x double], align 16
  %14 = alloca double, align 8
  %15 = alloca double, align 8
  %16 = alloca double, align 8
  %17 = alloca double, align 8
  %18 = alloca double, align 8
  %19 = alloca double, align 8
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca double, align 8
  %23 = alloca double, align 8
  %24 = alloca double, align 8
  %25 = alloca double, align 8
  %26 = alloca double, align 8
  %27 = alloca double, align 8
  %28 = alloca i32, align 4
  %29 = alloca i32, align 4
  %30 = alloca double, align 8
  %31 = alloca double, align 8
  %32 = alloca double, align 8
  %33 = alloca double, align 8
  %34 = alloca %struct.generic_integrator_workspace, align 8
  %35 = alloca %struct.thermodynamics_parameters_and_workspace, align 8
  %36 = alloca i32, align 4
  %37 = alloca [2048 x i8], align 16
  %38 = alloca i32, align 4
  %39 = alloca [2048 x i8], align 16
  %40 = alloca [2048 x i8], align 16
  %41 = alloca [2048 x i8], align 16
  %42 = alloca [2048 x i8], align 16
  %43 = alloca [2048 x i8], align 16
  %44 = alloca [2048 x i8], align 16
  %45 = alloca [2048 x i8], align 16
  %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
  store ptr %0, ptr %7, align 8, !tbaa !5
  store ptr %1, ptr %8, align 8, !tbaa !5
  store ptr %2, ptr %9, align 8, !tbaa !5
  store ptr %3, ptr %10, align 8, !tbaa !5
  store ptr %4, ptr %11, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 24, ptr %12) #6
  call void @llvm.lifetime.start.p0(i64 24, ptr %13) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %17) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %18) #6
  store double 0.000000e+00, ptr %18, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %19) #6
  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 8, 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
  call void @llvm.lifetime.start.p0(i64 8, ptr %27) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %28) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %29) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %30) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %31) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %32) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %33) #6
  call void @llvm.lifetime.start.p0(i64 2152, ptr %34) #6
  call void @llvm.lifetime.start.p0(i64 32, ptr %35) #6
  %50 = load ptr, ptr %7, align 8, !tbaa !5
  %51 = getelementptr inbounds %struct.precision, ptr %50, i32 0, i32 15
  %52 = load i32, ptr %51, align 8, !tbaa !11
  %53 = load ptr, ptr %10, align 8, !tbaa !5
  %54 = getelementptr inbounds %struct.recombination, ptr %53, i32 0, i32 7
  store i32 %52, ptr %54, align 4, !tbaa !14
  %55 = load ptr, ptr %10, align 8, !tbaa !5
  %56 = getelementptr inbounds %struct.recombination, ptr %55, i32 0, i32 6
  %57 = load i32, ptr %56, align 8, !tbaa !17
  %58 = load ptr, ptr %10, align 8, !tbaa !5
  %59 = getelementptr inbounds %struct.recombination, ptr %58, i32 0, i32 7
  %60 = load i32, ptr %59, align 4, !tbaa !14
  %61 = mul nsw i32 %57, %60
  %62 = sext i32 %61 to i64
  %63 = mul i64 %62, 8
  %64 = call noalias ptr @malloc(i64 noundef %63) #7
  %65 = load ptr, ptr %10, align 8, !tbaa !5
  %66 = getelementptr inbounds %struct.recombination, ptr %65, i32 0, i32 8
  store ptr %64, ptr %66, align 8, !tbaa !18
  %67 = load ptr, ptr %10, align 8, !tbaa !5
  %68 = getelementptr inbounds %struct.recombination, ptr %67, i32 0, i32 8
  %69 = load ptr, ptr %68, align 8, !tbaa !18
  %70 = icmp eq ptr %69, null
  br i1 %70, label %71, label %88

71:                                               ; preds = %5
  call void @llvm.lifetime.start.p0(i64 4, ptr %36) #6
  %72 = load ptr, ptr %10, align 8, !tbaa !5
  %73 = getelementptr inbounds %struct.recombination, ptr %72, i32 0, i32 6
  %74 = load i32, ptr %73, align 8, !tbaa !17
  %75 = load ptr, ptr %10, align 8, !tbaa !5
  %76 = getelementptr inbounds %struct.recombination, ptr %75, i32 0, i32 7
  %77 = load i32, ptr %76, align 4, !tbaa !14
  %78 = mul nsw i32 %74, %77
  %79 = sext i32 %78 to i64
  %80 = mul i64 %79, 8
  %81 = trunc i64 %80 to i32
  store i32 %81, ptr %36, align 4, !tbaa !19
  call void @llvm.lifetime.start.p0(i64 2048, ptr %37) #6
  %82 = getelementptr inbounds [2048 x i8], ptr %37, i64 0, i64 0
  %83 = load i32, ptr %36, align 4, !tbaa !19
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %82, ptr noundef @.str.37, ptr noundef @.str.212, i32 noundef %83)
  %84 = load ptr, ptr %9, align 8, !tbaa !5
  %85 = getelementptr inbounds %struct.thermo, ptr %84, i32 0, i32 93
  %86 = getelementptr inbounds [2048 x i8], ptr %85, i64 0, i64 0
  %87 = getelementptr inbounds [2048 x i8], ptr %37, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %86, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3478, ptr noundef %87)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %37) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %36) #6
  br label %1257

88:                                               ; preds = %5
  %89 = call i32 @initialize_generic_integrator(i32 noundef 3, ptr noundef %34)
  %90 = icmp eq i32 %89, 1
  br i1 %90, label %91, label %99

91:                                               ; preds = %88
  call void @llvm.lifetime.start.p0(i64 2048, ptr %39) #6
  %92 = getelementptr inbounds [2048 x i8], ptr %39, i64 0, i64 0
  %93 = getelementptr inbounds %struct.generic_integrator_workspace, ptr %34, i32 0, i32 13
  %94 = getelementptr inbounds [2048 x i8], ptr %93, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %92, ptr noundef @.str, ptr noundef @.str.213, ptr noundef %94)
  %95 = load ptr, ptr %9, align 8, !tbaa !5
  %96 = getelementptr inbounds %struct.thermo, ptr %95, i32 0, i32 93
  %97 = getelementptr inbounds [2048 x i8], ptr %96, i64 0, i64 0
  %98 = getelementptr inbounds [2048 x i8], ptr %39, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %97, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3483, ptr noundef %98)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %39) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

99:                                               ; preds = %88
  %100 = load ptr, ptr %7, align 8, !tbaa !5
  %101 = getelementptr inbounds %struct.precision, ptr %100, i32 0, i32 15
  %102 = load i32, ptr %101, align 8, !tbaa !11
  store i32 %102, ptr %29, align 4, !tbaa !19
  %103 = load ptr, ptr %8, align 8, !tbaa !5
  %104 = getelementptr inbounds %struct.background, ptr %103, i32 0, i32 0
  %105 = load double, ptr %104, align 8, !tbaa !20
  %106 = fmul double %105, 0x41B1DE784A000000
  %107 = fdiv double %106, 0x449A22FF260E3826
  %108 = load ptr, ptr %10, align 8, !tbaa !5
  %109 = getelementptr inbounds %struct.recombination, ptr %108, i32 0, i32 26
  store double %107, ptr %109, align 8, !tbaa !22
  %110 = load ptr, ptr %8, align 8, !tbaa !5
  %111 = getelementptr inbounds %struct.background, ptr %110, i32 0, i32 3
  %112 = load double, ptr %111, align 8, !tbaa !23
  store double %112, ptr %14, align 8, !tbaa !9
  %113 = load ptr, ptr %9, align 8, !tbaa !5
  %114 = getelementptr inbounds %struct.thermo, ptr %113, i32 0, i32 0
  %115 = load double, ptr %114, align 8, !tbaa !24
  %116 = load ptr, ptr %10, align 8, !tbaa !5
  %117 = getelementptr inbounds %struct.recombination, ptr %116, i32 0, i32 27
  store double %115, ptr %117, align 8, !tbaa !26
  %118 = load ptr, ptr %8, align 8, !tbaa !5
  %119 = getelementptr inbounds %struct.background, ptr %118, i32 0, i32 2
  %120 = load double, ptr %119, align 8, !tbaa !27
  %121 = load ptr, ptr %10, align 8, !tbaa !5
  %122 = getelementptr inbounds %struct.recombination, ptr %121, i32 0, i32 20
  store double %120, ptr %122, align 8, !tbaa !28
  %123 = load ptr, ptr %7, align 8, !tbaa !5
  %124 = getelementptr inbounds %struct.precision, ptr %123, i32 0, i32 14
  %125 = load double, ptr %124, align 8, !tbaa !29
  store double %125, ptr %15, align 8, !tbaa !9
  %126 = load ptr, ptr %7, align 8, !tbaa !5
  %127 = getelementptr inbounds %struct.precision, ptr %126, i32 0, i32 43
  %128 = load double, ptr %127, align 8, !tbaa !30
  %129 = load ptr, ptr %10, align 8, !tbaa !5
  %130 = getelementptr inbounds %struct.recombination, ptr %129, i32 0, i32 19
  store double %128, ptr %130, align 8, !tbaa !31
  %131 = load ptr, ptr %7, align 8, !tbaa !5
  %132 = getelementptr inbounds %struct.precision, ptr %131, i32 0, i32 22
  %133 = load i32, ptr %132, align 8, !tbaa !32
  %134 = icmp ne i32 %133, 1
  br i1 %134, label %135, label %148

135:                                              ; preds = %99
  %136 = load ptr, ptr %7, align 8, !tbaa !5
  %137 = getelementptr inbounds %struct.precision, ptr %136, i32 0, i32 22
  %138 = load i32, ptr %137, align 8, !tbaa !32
  %139 = icmp ne i32 %138, 0
  br i1 %139, label %140, label %148

140:                                              ; preds = %135
  call void @llvm.lifetime.start.p0(i64 2048, ptr %40) #6
  %141 = getelementptr inbounds [2048 x i8], ptr %40, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %141, ptr noundef @.str.214)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %41) #6
  %142 = getelementptr inbounds [2048 x i8], ptr %41, i64 0, i64 0
  %143 = getelementptr inbounds [2048 x i8], ptr %40, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %142, ptr noundef @.str.10, ptr noundef @.str.215, ptr noundef %143)
  %144 = load ptr, ptr %9, align 8, !tbaa !5
  %145 = getelementptr inbounds %struct.thermo, ptr %144, i32 0, i32 93
  %146 = getelementptr inbounds [2048 x i8], ptr %145, i64 0, i64 0
  %147 = getelementptr inbounds [2048 x i8], ptr %41, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %146, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3511, ptr noundef %147)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %41) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %40) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

148:                                              ; preds = %135, %99
  %149 = load ptr, ptr %7, align 8, !tbaa !5
  %150 = getelementptr inbounds %struct.precision, ptr %149, i32 0, i32 23
  %151 = load double, ptr %150, align 8, !tbaa !33
  %152 = load ptr, ptr %10, align 8, !tbaa !5
  %153 = getelementptr inbounds %struct.recombination, ptr %152, i32 0, i32 18
  store double %151, ptr %153, align 8, !tbaa !34
  %154 = load ptr, ptr %7, align 8, !tbaa !5
  %155 = getelementptr inbounds %struct.precision, ptr %154, i32 0, i32 22
  %156 = load i32, ptr %155, align 8, !tbaa !32
  %157 = icmp eq i32 %156, 1
  br i1 %157, label %158, label %166

158:                                              ; preds = %148
  %159 = load ptr, ptr %7, align 8, !tbaa !5
  %160 = getelementptr inbounds %struct.precision, ptr %159, i32 0, i32 24
  %161 = load double, ptr %160, align 8, !tbaa !35
  %162 = load ptr, ptr %10, align 8, !tbaa !5
  %163 = getelementptr inbounds %struct.recombination, ptr %162, i32 0, i32 18
  %164 = load double, ptr %163, align 8, !tbaa !34
  %165 = fadd double %164, %161
  store double %165, ptr %163, align 8, !tbaa !34
  br label %166

166:                                              ; preds = %158, %148
  %167 = load ptr, ptr %7, align 8, !tbaa !5
  %168 = getelementptr inbounds %struct.precision, ptr %167, i32 0, i32 20
  %169 = load i32, ptr %168, align 8, !tbaa !36
  %170 = icmp slt i32 %169, 0
  br i1 %170, label %176, label %171

171:                                              ; preds = %166
  %172 = load ptr, ptr %7, align 8, !tbaa !5
  %173 = getelementptr inbounds %struct.precision, ptr %172, i32 0, i32 20
  %174 = load i32, ptr %173, align 8, !tbaa !36
  %175 = icmp sgt i32 %174, 6
  br i1 %175, label %176, label %184

176:                                              ; preds = %171, %166
  call void @llvm.lifetime.start.p0(i64 2048, ptr %42) #6
  %177 = getelementptr inbounds [2048 x i8], ptr %42, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %177, ptr noundef @.str.216)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %43) #6
  %178 = getelementptr inbounds [2048 x i8], ptr %43, i64 0, i64 0
  %179 = getelementptr inbounds [2048 x i8], ptr %42, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %178, ptr noundef @.str.10, ptr noundef @.str.217, ptr noundef %179)
  %180 = load ptr, ptr %9, align 8, !tbaa !5
  %181 = getelementptr inbounds %struct.thermo, ptr %180, i32 0, i32 93
  %182 = getelementptr inbounds [2048 x i8], ptr %181, i64 0, i64 0
  %183 = getelementptr inbounds [2048 x i8], ptr %43, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %182, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3519, ptr noundef %183)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %43) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %42) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

184:                                              ; preds = %171
  %185 = load double, ptr %15, align 8, !tbaa !9
  store double %185, ptr %19, align 8, !tbaa !9
  %186 = load ptr, ptr %10, align 8, !tbaa !5
  %187 = getelementptr inbounds %struct.recombination, ptr %186, i32 0, i32 27
  %188 = load double, ptr %187, align 8, !tbaa !26
  %189 = fsub double 1.000000e+00, %188
  %190 = fdiv double 1.000000e+00, %189
  store double %190, ptr %20, align 8, !tbaa !9
  %191 = load ptr, ptr %10, align 8, !tbaa !5
  %192 = getelementptr inbounds %struct.recombination, ptr %191, i32 0, i32 27
  %193 = load double, ptr %192, align 8, !tbaa !26
  %194 = load ptr, ptr %10, align 8, !tbaa !5
  %195 = getelementptr inbounds %struct.recombination, ptr %194, i32 0, i32 27
  %196 = load double, ptr %195, align 8, !tbaa !26
  %197 = fsub double 1.000000e+00, %196
  %198 = fmul double 3.971500e+00, %197
  %199 = fdiv double %193, %198
  %200 = load ptr, ptr %10, align 8, !tbaa !5
  %201 = getelementptr inbounds %struct.recombination, ptr %200, i32 0, i32 14
  store double %199, ptr %201, align 8, !tbaa !37
  %202 = load ptr, ptr %10, align 8, !tbaa !5
  %203 = getelementptr inbounds %struct.recombination, ptr %202, i32 0, i32 26
  %204 = load double, ptr %203, align 8, !tbaa !22
  %205 = fmul double 3.000000e+00, %204
  %206 = load ptr, ptr %10, align 8, !tbaa !5
  %207 = getelementptr inbounds %struct.recombination, ptr %206, i32 0, i32 26
  %208 = load double, ptr %207, align 8, !tbaa !22
  %209 = fmul double %205, %208
  %210 = load double, ptr %14, align 8, !tbaa !9
  %211 = fmul double %209, %210
  %212 = load double, ptr %20, align 8, !tbaa !9
  %213 = fmul double 0x3E1CD169BDE2F4EC, %212
  %214 = fmul double %213, 0x3A60930446E176B5
  %215 = fdiv double %211, %214
  %216 = load ptr, ptr %10, align 8, !tbaa !5
  %217 = getelementptr inbounds %struct.recombination, ptr %216, i32 0, i32 21
  store double %215, ptr %217, align 8, !tbaa !38
  %218 = load ptr, ptr %10, align 8, !tbaa !5
  %219 = getelementptr inbounds %struct.recombination, ptr %218, i32 0, i32 21
  %220 = load double, ptr %219, align 8, !tbaa !38
  %221 = load ptr, ptr %9, align 8, !tbaa !5
  %222 = getelementptr inbounds %struct.thermo, ptr %221, i32 0, i32 89
  store double %220, ptr %222, align 8, !tbaa !39
  %223 = load ptr, ptr %9, align 8, !tbaa !5
  %224 = getelementptr inbounds %struct.thermo, ptr %223, i32 0, i32 23
  %225 = load double, ptr %224, align 8, !tbaa !40
  %226 = load ptr, ptr %10, align 8, !tbaa !5
  %227 = getelementptr inbounds %struct.recombination, ptr %226, i32 0, i32 28
  store double %225, ptr %227, align 8, !tbaa !41
  %228 = load ptr, ptr %9, align 8, !tbaa !5
  %229 = getelementptr inbounds %struct.thermo, ptr %228, i32 0, i32 24
  %230 = load i16, ptr %229, align 8, !tbaa !42
  %231 = load ptr, ptr %10, align 8, !tbaa !5
  %232 = getelementptr inbounds %struct.recombination, ptr %231, i32 0, i32 29
  store i16 %230, ptr %232, align 8, !tbaa !43
  %233 = load ptr, ptr %9, align 8, !tbaa !5
  %234 = getelementptr inbounds %struct.thermo, ptr %233, i32 0, i32 26
  %235 = load double, ptr %234, align 8, !tbaa !44
  %236 = load ptr, ptr %10, align 8, !tbaa !5
  %237 = getelementptr inbounds %struct.recombination, ptr %236, i32 0, i32 31
  store double %235, ptr %237, align 8, !tbaa !45
  %238 = load ptr, ptr %9, align 8, !tbaa !5
  %239 = getelementptr inbounds %struct.thermo, ptr %238, i32 0, i32 27
  %240 = load double, ptr %239, align 8, !tbaa !46
  %241 = load ptr, ptr %10, align 8, !tbaa !5
  %242 = getelementptr inbounds %struct.recombination, ptr %241, i32 0, i32 32
  store double %240, ptr %242, align 8, !tbaa !47
  %243 = load ptr, ptr %9, align 8, !tbaa !5
  %244 = getelementptr inbounds %struct.thermo, ptr %243, i32 0, i32 28
  %245 = load double, ptr %244, align 8, !tbaa !48
  %246 = load ptr, ptr %10, align 8, !tbaa !5
  %247 = getelementptr inbounds %struct.recombination, ptr %246, i32 0, i32 33
  store double %245, ptr %247, align 8, !tbaa !49
  %248 = load ptr, ptr %9, align 8, !tbaa !5
  %249 = getelementptr inbounds %struct.thermo, ptr %248, i32 0, i32 29
  %250 = load double, ptr %249, align 8, !tbaa !50
  %251 = load ptr, ptr %10, align 8, !tbaa !5
  %252 = getelementptr inbounds %struct.recombination, ptr %251, i32 0, i32 34
  store double %250, ptr %252, align 8, !tbaa !51
  %253 = load ptr, ptr %9, align 8, !tbaa !5
  %254 = getelementptr inbounds %struct.thermo, ptr %253, i32 0, i32 25
  %255 = load double, ptr %254, align 8, !tbaa !52
  %256 = load ptr, ptr %10, align 8, !tbaa !5
  %257 = getelementptr inbounds %struct.recombination, ptr %256, i32 0, i32 30
  store double %255, ptr %257, align 8, !tbaa !53
  %258 = load ptr, ptr %9, align 8, !tbaa !5
  %259 = getelementptr inbounds %struct.thermo, ptr %258, i32 0, i32 30
  %260 = load double, ptr %259, align 8, !tbaa !54
  %261 = load ptr, ptr %10, align 8, !tbaa !5
  %262 = getelementptr inbounds %struct.recombination, ptr %261, i32 0, i32 35
  store double %260, ptr %262, align 8, !tbaa !55
  %263 = load ptr, ptr %9, align 8, !tbaa !5
  %264 = getelementptr inbounds %struct.thermo, ptr %263, i32 0, i32 31
  %265 = load double, ptr %264, align 8, !tbaa !56
  %266 = load ptr, ptr %10, align 8, !tbaa !5
  %267 = getelementptr inbounds %struct.recombination, ptr %266, i32 0, i32 36
  store double %265, ptr %267, align 8, !tbaa !57
  store double 0x3E805102A7EF7C2D, ptr %21, align 8, !tbaa !9
  store double 0x3E6F5F07D0F2CB01, ptr %22, align 8, !tbaa !9
  store double 0x3C241850C63D88D9, ptr %23, align 8, !tbaa !9
  %268 = load double, ptr %23, align 8, !tbaa !9
  %269 = fdiv double %268, 0x3B30B0E7F1528C6D
  %270 = load ptr, ptr %10, align 8, !tbaa !5
  %271 = getelementptr inbounds %struct.recombination, ptr %270, i32 0, i32 9
  store double %269, ptr %271, align 8, !tbaa !58
  store double 0x3C2779EAF6C6E1C1, ptr %24, align 8, !tbaa !9
  %272 = load double, ptr %24, align 8, !tbaa !9
  %273 = fdiv double %272, 0x3B30B0E7F1528C6D
  %274 = load ptr, ptr %10, align 8, !tbaa !5
  %275 = getelementptr inbounds %struct.recombination, ptr %274, i32 0, i32 15
  store double %273, ptr %275, align 8, !tbaa !59
  %276 = load ptr, ptr %10, align 8, !tbaa !5
  %277 = getelementptr inbounds %struct.recombination, ptr %276, i32 0, i32 23
  store double 0x41034358C2F6C513, ptr %277, align 8, !tbaa !60
  %278 = load ptr, ptr %10, align 8, !tbaa !5
  %279 = getelementptr inbounds %struct.recombination, ptr %278, i32 0, i32 24
  store double 0x41116A32737408D4, ptr %279, align 8, !tbaa !61
  %280 = load ptr, ptr %10, align 8, !tbaa !5
  %281 = getelementptr inbounds %struct.recombination, ptr %280, i32 0, i32 25
  store double 0x412345866807EAE4, ptr %281, align 8, !tbaa !62
  %282 = load ptr, ptr %10, align 8, !tbaa !5
  %283 = getelementptr inbounds %struct.recombination, ptr %282, i32 0, i32 10
  store double 0x42E4764E04F59C55, ptr %283, align 8, !tbaa !63
  %284 = load double, ptr %21, align 8, !tbaa !9
  %285 = call double @pow(double noundef %284, double noundef 3.000000e+00) #6
  %286 = fdiv double %285, 0x403921FB54442D18
  %287 = load ptr, ptr %10, align 8, !tbaa !5
  %288 = getelementptr inbounds %struct.recombination, ptr %287, i32 0, i32 11
  store double %286, ptr %288, align 8, !tbaa !64
  %289 = load double, ptr %22, align 8, !tbaa !9
  %290 = call double @pow(double noundef %289, double noundef 3.000000e+00) #6
  %291 = fdiv double %290, 0x403921FB54442D18
  %292 = load ptr, ptr %10, align 8, !tbaa !5
  %293 = getelementptr inbounds %struct.recombination, ptr %292, i32 0, i32 16
  store double %291, ptr %293, align 8, !tbaa !65
  %294 = load double, ptr %21, align 8, !tbaa !9
  %295 = fmul double 0x3B30B0E7F1528C6D, %294
  %296 = fdiv double 0x3ACEBD1F10915C4E, %295
  %297 = load ptr, ptr %10, align 8, !tbaa !5
  %298 = getelementptr inbounds %struct.recombination, ptr %297, i32 0, i32 12
  store double %296, ptr %298, align 8, !tbaa !66
  %299 = load ptr, ptr %10, align 8, !tbaa !5
  %300 = getelementptr inbounds %struct.recombination, ptr %299, i32 0, i32 17
  store double 0x410D341EBBDFE7F9, ptr %300, align 8, !tbaa !67
  %301 = call double @pow(double noundef 0x400921FB54442D18, double noundef 5.000000e+00) #6
  %302 = fmul double 8.000000e+00, %301
  %303 = call double @pow(double noundef 0x3B30B0E7F1528C6D, double noundef 4.000000e+00) #6
  %304 = fmul double %302, %303
  %305 = fdiv double %304, 1.500000e+01
  %306 = call double @pow(double noundef 0x390B860B8B13E1DA, double noundef 3.000000e+00) #6
  %307 = fdiv double %305, %306
  %308 = call double @pow(double noundef 0x41B1DE784A000000, double noundef 3.000000e+00) #6
  %309 = fdiv double %307, %308
  %310 = fmul double 0x3EA5CBF6997CAED0, %309
  %311 = load ptr, ptr %10, align 8, !tbaa !5
  %312 = getelementptr inbounds %struct.recombination, ptr %311, i32 0, i32 13
  store double %310, ptr %312, align 8, !tbaa !68
  %313 = load ptr, ptr %10, align 8, !tbaa !5
  %314 = getelementptr inbounds %struct.recombination, ptr %313, i32 0, i32 22
  store double 0x40BB4CC9D960EA35, ptr %314, align 8, !tbaa !69
  %315 = load ptr, ptr %8, align 8, !tbaa !5
  %316 = getelementptr inbounds %struct.thermodynamics_parameters_and_workspace, ptr %35, i32 0, i32 0
  store ptr %315, ptr %316, align 8, !tbaa !70
  %317 = load ptr, ptr %7, align 8, !tbaa !5
  %318 = getelementptr inbounds %struct.thermodynamics_parameters_and_workspace, ptr %35, i32 0, i32 1
  store ptr %317, ptr %318, align 8, !tbaa !72
  %319 = load ptr, ptr %10, align 8, !tbaa !5
  %320 = getelementptr inbounds %struct.thermodynamics_parameters_and_workspace, ptr %35, i32 0, i32 2
  store ptr %319, ptr %320, align 8, !tbaa !73
  %321 = load ptr, ptr %11, align 8, !tbaa !5
  %322 = getelementptr inbounds %struct.thermodynamics_parameters_and_workspace, ptr %35, i32 0, i32 3
  store ptr %321, ptr %322, align 8, !tbaa !74
  %323 = load double, ptr %15, align 8, !tbaa !9
  %324 = load ptr, ptr %7, align 8, !tbaa !5
  %325 = getelementptr inbounds %struct.precision, ptr %324, i32 0, i32 35
  %326 = load double, ptr %325, align 8, !tbaa !75
  %327 = fcmp olt double %323, %326
  br i1 %327, label %328, label %336

328:                                              ; preds = %184
  call void @llvm.lifetime.start.p0(i64 2048, ptr %44) #6
  %329 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %329, ptr noundef @.str.218)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %45) #6
  %330 = getelementptr inbounds [2048 x i8], ptr %45, i64 0, i64 0
  %331 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %330, ptr noundef @.str.10, ptr noundef @.str.219, ptr noundef %331)
  %332 = load ptr, ptr %9, 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 %45, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %334, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3571, ptr noundef %335)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %45) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %44) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

336:                                              ; preds = %184
  %337 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  store double 1.000000e+00, ptr %337, align 16, !tbaa !9
  %338 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  store double 1.000000e+00, ptr %338, align 8, !tbaa !9
  %339 = load ptr, ptr %10, align 8, !tbaa !5
  %340 = getelementptr inbounds %struct.recombination, ptr %339, i32 0, i32 14
  %341 = load double, ptr %340, align 8, !tbaa !37
  %342 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %341, double 1.000000e+00)
  store double %342, ptr %17, align 8, !tbaa !9
  %343 = load ptr, ptr %10, align 8, !tbaa !5
  %344 = getelementptr inbounds %struct.recombination, ptr %343, i32 0, i32 20
  %345 = load double, ptr %344, align 8, !tbaa !28
  %346 = load double, ptr %19, align 8, !tbaa !9
  %347 = fadd double 1.000000e+00, %346
  %348 = fmul double %345, %347
  %349 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  store double %348, ptr %349, align 16, !tbaa !9
  store i32 0, ptr %28, align 4, !tbaa !19
  br label %350

350:                                              ; preds = %1242, %336
  %351 = load i32, ptr %28, align 4, !tbaa !19
  %352 = load i32, ptr %29, align 4, !tbaa !19
  %353 = icmp slt i32 %351, %352
  br i1 %353, label %354, label %1245

354:                                              ; preds = %350
  %355 = load double, ptr %15, align 8, !tbaa !9
  %356 = load i32, ptr %29, align 4, !tbaa !19
  %357 = load i32, ptr %28, align 4, !tbaa !19
  %358 = sub nsw i32 %356, %357
  %359 = sitofp i32 %358 to double
  %360 = fmul double %355, %359
  %361 = load i32, ptr %29, align 4, !tbaa !19
  %362 = sitofp i32 %361 to double
  %363 = fdiv double %360, %362
  store double %363, ptr %25, align 8, !tbaa !9
  %364 = load double, ptr %15, align 8, !tbaa !9
  %365 = load i32, ptr %29, align 4, !tbaa !19
  %366 = load i32, ptr %28, align 4, !tbaa !19
  %367 = sub nsw i32 %365, %366
  %368 = sub nsw i32 %367, 1
  %369 = sitofp i32 %368 to double
  %370 = fmul double %364, %369
  %371 = load i32, ptr %29, align 4, !tbaa !19
  %372 = sitofp i32 %371 to double
  %373 = fdiv double %370, %372
  store double %373, ptr %26, align 8, !tbaa !9
  %374 = load double, ptr %26, align 8, !tbaa !9
  store double %374, ptr %19, align 8, !tbaa !9
  %375 = load double, ptr %19, align 8, !tbaa !9
  %376 = load ptr, ptr %7, align 8, !tbaa !5
  %377 = getelementptr inbounds %struct.precision, ptr %376, i32 0, i32 31
  %378 = load double, ptr %377, align 8, !tbaa !76
  %379 = load ptr, ptr %7, align 8, !tbaa !5
  %380 = getelementptr inbounds %struct.precision, ptr %379, i32 0, i32 32
  %381 = load double, ptr %380, align 8, !tbaa !77
  %382 = fadd double %378, %381
  %383 = fcmp ogt double %375, %382
  br i1 %383, label %384, label %400

384:                                              ; preds = %354
  store double 1.000000e+00, ptr %18, align 8, !tbaa !9
  store double 1.000000e+00, ptr %16, align 8, !tbaa !9
  %385 = load ptr, ptr %10, align 8, !tbaa !5
  %386 = getelementptr inbounds %struct.recombination, ptr %385, i32 0, i32 14
  %387 = load double, ptr %386, align 8, !tbaa !37
  %388 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %387, double 1.000000e+00)
  store double %388, ptr %17, align 8, !tbaa !9
  %389 = load double, ptr %18, align 8, !tbaa !9
  %390 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  store double %389, ptr %390, align 16, !tbaa !9
  %391 = load double, ptr %16, align 8, !tbaa !9
  %392 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  store double %391, ptr %392, align 8, !tbaa !9
  %393 = load ptr, ptr %10, align 8, !tbaa !5
  %394 = getelementptr inbounds %struct.recombination, ptr %393, i32 0, i32 20
  %395 = load double, ptr %394, align 8, !tbaa !28
  %396 = load double, ptr %19, align 8, !tbaa !9
  %397 = fadd double 1.000000e+00, %396
  %398 = fmul double %395, %397
  %399 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  store double %398, ptr %399, align 16, !tbaa !9
  br label %1054

400:                                              ; preds = %354
  %401 = load double, ptr %19, align 8, !tbaa !9
  %402 = load ptr, ptr %7, align 8, !tbaa !5
  %403 = getelementptr inbounds %struct.precision, ptr %402, i32 0, i32 33
  %404 = load double, ptr %403, align 8, !tbaa !78
  %405 = load ptr, ptr %7, align 8, !tbaa !5
  %406 = getelementptr inbounds %struct.precision, ptr %405, i32 0, i32 34
  %407 = load double, ptr %406, align 8, !tbaa !79
  %408 = fadd double %404, %407
  %409 = fcmp ogt double %401, %408
  br i1 %409, label %410, label %536

410:                                              ; preds = %400
  store double 1.000000e+00, ptr %18, align 8, !tbaa !9
  store double 1.000000e+00, ptr %16, align 8, !tbaa !9
  %411 = load ptr, ptr %10, align 8, !tbaa !5
  %412 = getelementptr inbounds %struct.recombination, ptr %411, i32 0, i32 10
  %413 = load double, ptr %412, align 8, !tbaa !63
  %414 = load ptr, ptr %10, align 8, !tbaa !5
  %415 = getelementptr inbounds %struct.recombination, ptr %414, i32 0, i32 20
  %416 = load double, ptr %415, align 8, !tbaa !28
  %417 = fmul double %413, %416
  %418 = load double, ptr %19, align 8, !tbaa !9
  %419 = fadd double 1.000000e+00, %418
  %420 = fdiv double %417, %419
  %421 = call double @log(double noundef %420) #6
  %422 = load ptr, ptr %10, align 8, !tbaa !5
  %423 = getelementptr inbounds %struct.recombination, ptr %422, i32 0, i32 25
  %424 = load double, ptr %423, align 8, !tbaa !62
  %425 = load ptr, ptr %10, align 8, !tbaa !5
  %426 = getelementptr inbounds %struct.recombination, ptr %425, i32 0, i32 20
  %427 = load double, ptr %426, align 8, !tbaa !28
  %428 = load double, ptr %19, align 8, !tbaa !9
  %429 = fadd double 1.000000e+00, %428
  %430 = fmul double %427, %429
  %431 = fdiv double %424, %430
  %432 = fneg double %431
  %433 = call double @llvm.fmuladd.f64(double 1.500000e+00, double %421, double %432)
  %434 = call double @exp(double noundef %433) #6
  %435 = load ptr, ptr %10, align 8, !tbaa !5
  %436 = getelementptr inbounds %struct.recombination, ptr %435, i32 0, i32 21
  %437 = load double, ptr %436, align 8, !tbaa !38
  %438 = fdiv double %434, %437
  store double %438, ptr %27, align 8, !tbaa !9
  %439 = load double, ptr %19, align 8, !tbaa !9
  %440 = load ptr, ptr %7, align 8, !tbaa !5
  %441 = getelementptr inbounds %struct.precision, ptr %440, i32 0, i32 31
  %442 = load double, ptr %441, align 8, !tbaa !76
  %443 = load ptr, ptr %7, align 8, !tbaa !5
  %444 = getelementptr inbounds %struct.precision, ptr %443, i32 0, i32 32
  %445 = load double, ptr %444, align 8, !tbaa !77
  %446 = fsub double %442, %445
  %447 = fcmp ogt double %439, %446
  br i1 %447, label %448, label %500

448:                                              ; preds = %410
  %449 = load ptr, ptr %10, align 8, !tbaa !5
  %450 = getelementptr inbounds %struct.recombination, ptr %449, i32 0, i32 14
  %451 = load double, ptr %450, align 8, !tbaa !37
  %452 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %451, double 1.000000e+00)
  store double %452, ptr %30, align 8, !tbaa !9
  %453 = load double, ptr %27, align 8, !tbaa !9
  %454 = fsub double %453, 1.000000e+00
  %455 = load ptr, ptr %10, align 8, !tbaa !5
  %456 = getelementptr inbounds %struct.recombination, ptr %455, i32 0, i32 14
  %457 = load double, ptr %456, align 8, !tbaa !37
  %458 = fsub double %454, %457
  %459 = call double @pow(double noundef %458, double noundef 2.000000e+00) #6
  %460 = load ptr, ptr %10, align 8, !tbaa !5
  %461 = getelementptr inbounds %struct.recombination, ptr %460, i32 0, i32 14
  %462 = load double, ptr %461, align 8, !tbaa !37
  %463 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %462, double 1.000000e+00)
  %464 = fmul double 4.000000e+00, %463
  %465 = load double, ptr %27, align 8, !tbaa !9
  %466 = call double @llvm.fmuladd.f64(double %464, double %465, double %459)
  %467 = call double @sqrt(double noundef %466) #6
  %468 = load double, ptr %27, align 8, !tbaa !9
  %469 = fsub double %468, 1.000000e+00
  %470 = load ptr, ptr %10, align 8, !tbaa !5
  %471 = getelementptr inbounds %struct.recombination, ptr %470, i32 0, i32 14
  %472 = load double, ptr %471, align 8, !tbaa !37
  %473 = fsub double %469, %472
  %474 = fsub double %467, %473
  %475 = fmul double 5.000000e-01, %474
  store double %475, ptr %31, align 8, !tbaa !9
  %476 = load ptr, ptr %7, align 8, !tbaa !5
  %477 = getelementptr inbounds %struct.precision, ptr %476, i32 0, i32 31
  %478 = load double, ptr %477, align 8, !tbaa !76
  %479 = load double, ptr %19, align 8, !tbaa !9
  %480 = fsub double %478, %479
  %481 = load ptr, ptr %7, align 8, !tbaa !5
  %482 = getelementptr inbounds %struct.precision, ptr %481, i32 0, i32 32
  %483 = load double, ptr %482, align 8, !tbaa !77
  %484 = fdiv double %480, %483
  store double %484, ptr %32, align 8, !tbaa !9
  %485 = load double, ptr %32, align 8, !tbaa !9
  %486 = fmul double -7.500000e-01, %485
  %487 = load double, ptr %32, align 8, !tbaa !9
  %488 = load double, ptr %32, align 8, !tbaa !9
  %489 = fmul double %487, %488
  %490 = fdiv double %489, 3.000000e+00
  %491 = fsub double %490, 1.000000e+00
  %492 = call double @llvm.fmuladd.f64(double %486, double %491, double 5.000000e-01)
  store double %492, ptr %33, align 8, !tbaa !9
  %493 = load double, ptr %33, align 8, !tbaa !9
  %494 = load double, ptr %31, align 8, !tbaa !9
  %495 = load double, ptr %33, align 8, !tbaa !9
  %496 = fsub double 1.000000e+00, %495
  %497 = load double, ptr %30, align 8, !tbaa !9
  %498 = fmul double %496, %497
  %499 = call double @llvm.fmuladd.f64(double %493, double %494, double %498)
  store double %499, ptr %17, align 8, !tbaa !9
  br label %524

500:                                              ; preds = %410
  %501 = load double, ptr %27, align 8, !tbaa !9
  %502 = fsub double %501, 1.000000e+00
  %503 = load ptr, ptr %10, align 8, !tbaa !5
  %504 = getelementptr inbounds %struct.recombination, ptr %503, i32 0, i32 14
  %505 = load double, ptr %504, align 8, !tbaa !37
  %506 = fsub double %502, %505
  %507 = call double @pow(double noundef %506, double noundef 2.000000e+00) #6
  %508 = load ptr, ptr %10, align 8, !tbaa !5
  %509 = getelementptr inbounds %struct.recombination, ptr %508, i32 0, i32 14
  %510 = load double, ptr %509, align 8, !tbaa !37
  %511 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %510, double 1.000000e+00)
  %512 = fmul double 4.000000e+00, %511
  %513 = load double, ptr %27, align 8, !tbaa !9
  %514 = call double @llvm.fmuladd.f64(double %512, double %513, double %507)
  %515 = call double @sqrt(double noundef %514) #6
  %516 = load double, ptr %27, align 8, !tbaa !9
  %517 = fsub double %516, 1.000000e+00
  %518 = load ptr, ptr %10, align 8, !tbaa !5
  %519 = getelementptr inbounds %struct.recombination, ptr %518, i32 0, i32 14
  %520 = load double, ptr %519, align 8, !tbaa !37
  %521 = fsub double %517, %520
  %522 = fsub double %515, %521
  %523 = fmul double 5.000000e-01, %522
  store double %523, ptr %17, align 8, !tbaa !9
  br label %524

524:                                              ; preds = %500, %448
  %525 = load double, ptr %18, align 8, !tbaa !9
  %526 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  store double %525, ptr %526, align 16, !tbaa !9
  %527 = load double, ptr %16, align 8, !tbaa !9
  %528 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  store double %527, ptr %528, align 8, !tbaa !9
  %529 = load ptr, ptr %10, align 8, !tbaa !5
  %530 = getelementptr inbounds %struct.recombination, ptr %529, i32 0, i32 20
  %531 = load double, ptr %530, align 8, !tbaa !28
  %532 = load double, ptr %19, align 8, !tbaa !9
  %533 = fadd double 1.000000e+00, %532
  %534 = fmul double %531, %533
  %535 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  store double %534, ptr %535, align 16, !tbaa !9
  br label %1053

536:                                              ; preds = %400
  %537 = load double, ptr %19, align 8, !tbaa !9
  %538 = load ptr, ptr %7, align 8, !tbaa !5
  %539 = getelementptr inbounds %struct.precision, ptr %538, i32 0, i32 35
  %540 = load double, ptr %539, align 8, !tbaa !75
  %541 = load ptr, ptr %7, align 8, !tbaa !5
  %542 = getelementptr inbounds %struct.precision, ptr %541, i32 0, i32 36
  %543 = load double, ptr %542, align 8, !tbaa !80
  %544 = fadd double %540, %543
  %545 = fcmp ogt double %537, %544
  br i1 %545, label %546, label %653

546:                                              ; preds = %536
  store double 1.000000e+00, ptr %18, align 8, !tbaa !9
  store double 1.000000e+00, ptr %16, align 8, !tbaa !9
  %547 = load double, ptr %19, align 8, !tbaa !9
  %548 = load ptr, ptr %7, align 8, !tbaa !5
  %549 = getelementptr inbounds %struct.precision, ptr %548, i32 0, i32 33
  %550 = load double, ptr %549, align 8, !tbaa !78
  %551 = load ptr, ptr %7, align 8, !tbaa !5
  %552 = getelementptr inbounds %struct.precision, ptr %551, i32 0, i32 34
  %553 = load double, ptr %552, align 8, !tbaa !79
  %554 = fsub double %550, %553
  %555 = fcmp ogt double %547, %554
  br i1 %555, label %556, label %636

556:                                              ; preds = %546
  %557 = load ptr, ptr %10, align 8, !tbaa !5
  %558 = getelementptr inbounds %struct.recombination, ptr %557, i32 0, i32 10
  %559 = load double, ptr %558, align 8, !tbaa !63
  %560 = load ptr, ptr %10, align 8, !tbaa !5
  %561 = getelementptr inbounds %struct.recombination, ptr %560, i32 0, i32 20
  %562 = load double, ptr %561, align 8, !tbaa !28
  %563 = fmul double %559, %562
  %564 = load double, ptr %19, align 8, !tbaa !9
  %565 = fadd double 1.000000e+00, %564
  %566 = fdiv double %563, %565
  %567 = call double @log(double noundef %566) #6
  %568 = load ptr, ptr %10, align 8, !tbaa !5
  %569 = getelementptr inbounds %struct.recombination, ptr %568, i32 0, i32 25
  %570 = load double, ptr %569, align 8, !tbaa !62
  %571 = load ptr, ptr %10, align 8, !tbaa !5
  %572 = getelementptr inbounds %struct.recombination, ptr %571, i32 0, i32 20
  %573 = load double, ptr %572, align 8, !tbaa !28
  %574 = load double, ptr %19, align 8, !tbaa !9
  %575 = fadd double 1.000000e+00, %574
  %576 = fmul double %573, %575
  %577 = fdiv double %570, %576
  %578 = fneg double %577
  %579 = call double @llvm.fmuladd.f64(double 1.500000e+00, double %567, double %578)
  %580 = call double @exp(double noundef %579) #6
  %581 = load ptr, ptr %10, align 8, !tbaa !5
  %582 = getelementptr inbounds %struct.recombination, ptr %581, i32 0, i32 21
  %583 = load double, ptr %582, align 8, !tbaa !38
  %584 = fdiv double %580, %583
  store double %584, ptr %27, align 8, !tbaa !9
  %585 = load double, ptr %27, align 8, !tbaa !9
  %586 = fsub double %585, 1.000000e+00
  %587 = load ptr, ptr %10, align 8, !tbaa !5
  %588 = getelementptr inbounds %struct.recombination, ptr %587, i32 0, i32 14
  %589 = load double, ptr %588, align 8, !tbaa !37
  %590 = fsub double %586, %589
  %591 = call double @pow(double noundef %590, double noundef 2.000000e+00) #6
  %592 = load ptr, ptr %10, align 8, !tbaa !5
  %593 = getelementptr inbounds %struct.recombination, ptr %592, i32 0, i32 14
  %594 = load double, ptr %593, align 8, !tbaa !37
  %595 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %594, double 1.000000e+00)
  %596 = fmul double 4.000000e+00, %595
  %597 = load double, ptr %27, align 8, !tbaa !9
  %598 = call double @llvm.fmuladd.f64(double %596, double %597, double %591)
  %599 = call double @sqrt(double noundef %598) #6
  %600 = load double, ptr %27, align 8, !tbaa !9
  %601 = fsub double %600, 1.000000e+00
  %602 = load ptr, ptr %10, align 8, !tbaa !5
  %603 = getelementptr inbounds %struct.recombination, ptr %602, i32 0, i32 14
  %604 = load double, ptr %603, align 8, !tbaa !37
  %605 = fsub double %601, %604
  %606 = fsub double %599, %605
  %607 = fmul double 5.000000e-01, %606
  store double %607, ptr %30, align 8, !tbaa !9
  %608 = load ptr, ptr %10, align 8, !tbaa !5
  %609 = getelementptr inbounds %struct.recombination, ptr %608, i32 0, i32 14
  %610 = load double, ptr %609, align 8, !tbaa !37
  %611 = fadd double 1.000000e+00, %610
  store double %611, ptr %31, align 8, !tbaa !9
  %612 = load ptr, ptr %7, align 8, !tbaa !5
  %613 = getelementptr inbounds %struct.precision, ptr %612, i32 0, i32 33
  %614 = load double, ptr %613, align 8, !tbaa !78
  %615 = load double, ptr %19, align 8, !tbaa !9
  %616 = fsub double %614, %615
  %617 = load ptr, ptr %7, align 8, !tbaa !5
  %618 = getelementptr inbounds %struct.precision, ptr %617, i32 0, i32 34
  %619 = load double, ptr %618, align 8, !tbaa !79
  %620 = fdiv double %616, %619
  store double %620, ptr %32, align 8, !tbaa !9
  %621 = load double, ptr %32, align 8, !tbaa !9
  %622 = fmul double -7.500000e-01, %621
  %623 = load double, ptr %32, align 8, !tbaa !9
  %624 = load double, ptr %32, align 8, !tbaa !9
  %625 = fmul double %623, %624
  %626 = fdiv double %625, 3.000000e+00
  %627 = fsub double %626, 1.000000e+00
  %628 = call double @llvm.fmuladd.f64(double %622, double %627, double 5.000000e-01)
  store double %628, ptr %33, align 8, !tbaa !9
  %629 = load double, ptr %33, align 8, !tbaa !9
  %630 = load double, ptr %31, align 8, !tbaa !9
  %631 = load double, ptr %33, align 8, !tbaa !9
  %632 = fsub double 1.000000e+00, %631
  %633 = load double, ptr %30, align 8, !tbaa !9
  %634 = fmul double %632, %633
  %635 = call double @llvm.fmuladd.f64(double %629, double %630, double %634)
  store double %635, ptr %17, align 8, !tbaa !9
  br label %641

636:                                              ; preds = %546
  %637 = load ptr, ptr %10, align 8, !tbaa !5
  %638 = getelementptr inbounds %struct.recombination, ptr %637, i32 0, i32 14
  %639 = load double, ptr %638, align 8, !tbaa !37
  %640 = fadd double 1.000000e+00, %639
  store double %640, ptr %17, align 8, !tbaa !9
  br label %641

641:                                              ; preds = %636, %556
  %642 = load double, ptr %18, align 8, !tbaa !9
  %643 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  store double %642, ptr %643, align 16, !tbaa !9
  %644 = load double, ptr %16, align 8, !tbaa !9
  %645 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  store double %644, ptr %645, align 8, !tbaa !9
  %646 = load ptr, ptr %10, align 8, !tbaa !5
  %647 = getelementptr inbounds %struct.recombination, ptr %646, i32 0, i32 20
  %648 = load double, ptr %647, align 8, !tbaa !28
  %649 = load double, ptr %19, align 8, !tbaa !9
  %650 = fadd double 1.000000e+00, %649
  %651 = fmul double %648, %650
  %652 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  store double %651, ptr %652, align 16, !tbaa !9
  br label %1052

653:                                              ; preds = %536
  %654 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  %655 = load double, ptr %654, align 8, !tbaa !9
  %656 = load ptr, ptr %7, align 8, !tbaa !5
  %657 = getelementptr inbounds %struct.precision, ptr %656, i32 0, i32 37
  %658 = load double, ptr %657, align 8, !tbaa !81
  %659 = fcmp ogt double %655, %658
  br i1 %659, label %660, label %764

660:                                              ; preds = %653
  store double 1.000000e+00, ptr %18, align 8, !tbaa !9
  %661 = load ptr, ptr %10, align 8, !tbaa !5
  %662 = getelementptr inbounds %struct.recombination, ptr %661, i32 0, i32 10
  %663 = load double, ptr %662, align 8, !tbaa !63
  %664 = load ptr, ptr %10, align 8, !tbaa !5
  %665 = getelementptr inbounds %struct.recombination, ptr %664, i32 0, i32 20
  %666 = load double, ptr %665, align 8, !tbaa !28
  %667 = fmul double %663, %666
  %668 = load double, ptr %19, align 8, !tbaa !9
  %669 = fadd double 1.000000e+00, %668
  %670 = fdiv double %667, %669
  %671 = call double @log(double noundef %670) #6
  %672 = load ptr, ptr %10, align 8, !tbaa !5
  %673 = getelementptr inbounds %struct.recombination, ptr %672, i32 0, i32 24
  %674 = load double, ptr %673, align 8, !tbaa !61
  %675 = load ptr, ptr %10, align 8, !tbaa !5
  %676 = getelementptr inbounds %struct.recombination, ptr %675, i32 0, i32 20
  %677 = load double, ptr %676, align 8, !tbaa !28
  %678 = load double, ptr %19, align 8, !tbaa !9
  %679 = fadd double 1.000000e+00, %678
  %680 = fmul double %677, %679
  %681 = fdiv double %674, %680
  %682 = fneg double %681
  %683 = call double @llvm.fmuladd.f64(double 1.500000e+00, double %671, double %682)
  %684 = call double @exp(double noundef %683) #6
  %685 = fmul double 4.000000e+00, %684
  %686 = load ptr, ptr %10, align 8, !tbaa !5
  %687 = getelementptr inbounds %struct.recombination, ptr %686, i32 0, i32 21
  %688 = load double, ptr %687, align 8, !tbaa !38
  %689 = fdiv double %685, %688
  store double %689, ptr %27, align 8, !tbaa !9
  %690 = load double, ptr %27, align 8, !tbaa !9
  %691 = fsub double %690, 1.000000e+00
  %692 = call double @pow(double noundef %691, double noundef 2.000000e+00) #6
  %693 = load ptr, ptr %10, align 8, !tbaa !5
  %694 = getelementptr inbounds %struct.recombination, ptr %693, i32 0, i32 14
  %695 = load double, ptr %694, align 8, !tbaa !37
  %696 = fadd double 1.000000e+00, %695
  %697 = fmul double 4.000000e+00, %696
  %698 = load double, ptr %27, align 8, !tbaa !9
  %699 = call double @llvm.fmuladd.f64(double %697, double %698, double %692)
  %700 = call double @sqrt(double noundef %699) #6
  %701 = load double, ptr %27, align 8, !tbaa !9
  %702 = fsub double %701, 1.000000e+00
  %703 = fsub double %700, %702
  %704 = fmul double 5.000000e-01, %703
  store double %704, ptr %16, align 8, !tbaa !9
  %705 = load double, ptr %19, align 8, !tbaa !9
  %706 = load ptr, ptr %7, align 8, !tbaa !5
  %707 = getelementptr inbounds %struct.precision, ptr %706, i32 0, i32 35
  %708 = load double, ptr %707, align 8, !tbaa !75
  %709 = load ptr, ptr %7, align 8, !tbaa !5
  %710 = getelementptr inbounds %struct.precision, ptr %709, i32 0, i32 36
  %711 = load double, ptr %710, align 8, !tbaa !80
  %712 = fsub double %708, %711
  %713 = fcmp ogt double %705, %712
  br i1 %713, label %714, label %744

714:                                              ; preds = %660
  %715 = load ptr, ptr %10, align 8, !tbaa !5
  %716 = getelementptr inbounds %struct.recombination, ptr %715, i32 0, i32 14
  %717 = load double, ptr %716, align 8, !tbaa !37
  %718 = fadd double 1.000000e+00, %717
  store double %718, ptr %30, align 8, !tbaa !9
  %719 = load double, ptr %16, align 8, !tbaa !9
  store double %719, ptr %31, align 8, !tbaa !9
  %720 = load ptr, ptr %7, align 8, !tbaa !5
  %721 = getelementptr inbounds %struct.precision, ptr %720, i32 0, i32 35
  %722 = load double, ptr %721, align 8, !tbaa !75
  %723 = load double, ptr %19, align 8, !tbaa !9
  %724 = fsub double %722, %723
  %725 = load ptr, ptr %7, align 8, !tbaa !5
  %726 = getelementptr inbounds %struct.precision, ptr %725, i32 0, i32 36
  %727 = load double, ptr %726, align 8, !tbaa !80
  %728 = fdiv double %724, %727
  store double %728, ptr %32, align 8, !tbaa !9
  %729 = load double, ptr %32, align 8, !tbaa !9
  %730 = fmul double -7.500000e-01, %729
  %731 = load double, ptr %32, align 8, !tbaa !9
  %732 = load double, ptr %32, align 8, !tbaa !9
  %733 = fmul double %731, %732
  %734 = fdiv double %733, 3.000000e+00
  %735 = fsub double %734, 1.000000e+00
  %736 = call double @llvm.fmuladd.f64(double %730, double %735, double 5.000000e-01)
  store double %736, ptr %33, align 8, !tbaa !9
  %737 = load double, ptr %33, align 8, !tbaa !9
  %738 = load double, ptr %31, align 8, !tbaa !9
  %739 = load double, ptr %33, align 8, !tbaa !9
  %740 = fsub double 1.000000e+00, %739
  %741 = load double, ptr %30, align 8, !tbaa !9
  %742 = fmul double %740, %741
  %743 = call double @llvm.fmuladd.f64(double %737, double %738, double %742)
  store double %743, ptr %17, align 8, !tbaa !9
  br label %746

744:                                              ; preds = %660
  %745 = load double, ptr %16, align 8, !tbaa !9
  store double %745, ptr %17, align 8, !tbaa !9
  br label %746

746:                                              ; preds = %744, %714
  %747 = load double, ptr %17, align 8, !tbaa !9
  %748 = fsub double %747, 1.000000e+00
  %749 = load ptr, ptr %10, align 8, !tbaa !5
  %750 = getelementptr inbounds %struct.recombination, ptr %749, i32 0, i32 14
  %751 = load double, ptr %750, align 8, !tbaa !37
  %752 = fdiv double %748, %751
  store double %752, ptr %16, align 8, !tbaa !9
  %753 = load double, ptr %18, align 8, !tbaa !9
  %754 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  store double %753, ptr %754, align 16, !tbaa !9
  %755 = load double, ptr %16, align 8, !tbaa !9
  %756 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  store double %755, ptr %756, align 8, !tbaa !9
  %757 = load ptr, ptr %10, align 8, !tbaa !5
  %758 = getelementptr inbounds %struct.recombination, ptr %757, i32 0, i32 20
  %759 = load double, ptr %758, align 8, !tbaa !28
  %760 = load double, ptr %19, align 8, !tbaa !9
  %761 = fadd double 1.000000e+00, %760
  %762 = fmul double %759, %761
  %763 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  store double %762, ptr %763, align 16, !tbaa !9
  br label %1051

764:                                              ; preds = %653
  %765 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %766 = load double, ptr %765, align 16, !tbaa !9
  %767 = load ptr, ptr %7, align 8, !tbaa !5
  %768 = getelementptr inbounds %struct.precision, ptr %767, i32 0, i32 40
  %769 = load double, ptr %768, align 8, !tbaa !82
  %770 = fcmp ogt double %766, %769
  br i1 %770, label %771, label %928

771:                                              ; preds = %764
  %772 = load ptr, ptr %10, align 8, !tbaa !5
  %773 = getelementptr inbounds %struct.recombination, ptr %772, i32 0, i32 10
  %774 = load double, ptr %773, align 8, !tbaa !63
  %775 = load ptr, ptr %10, align 8, !tbaa !5
  %776 = getelementptr inbounds %struct.recombination, ptr %775, i32 0, i32 20
  %777 = load double, ptr %776, align 8, !tbaa !28
  %778 = fmul double %774, %777
  %779 = load double, ptr %19, align 8, !tbaa !9
  %780 = fadd double 1.000000e+00, %779
  %781 = fdiv double %778, %780
  %782 = call double @log(double noundef %781) #6
  %783 = load ptr, ptr %10, align 8, !tbaa !5
  %784 = getelementptr inbounds %struct.recombination, ptr %783, i32 0, i32 23
  %785 = load double, ptr %784, align 8, !tbaa !60
  %786 = load ptr, ptr %10, align 8, !tbaa !5
  %787 = getelementptr inbounds %struct.recombination, ptr %786, i32 0, i32 20
  %788 = load double, ptr %787, align 8, !tbaa !28
  %789 = load double, ptr %19, align 8, !tbaa !9
  %790 = fadd double 1.000000e+00, %789
  %791 = fmul double %788, %790
  %792 = fdiv double %785, %791
  %793 = fneg double %792
  %794 = call double @llvm.fmuladd.f64(double 1.500000e+00, double %782, double %793)
  %795 = call double @exp(double noundef %794) #6
  %796 = load ptr, ptr %10, align 8, !tbaa !5
  %797 = getelementptr inbounds %struct.recombination, ptr %796, i32 0, i32 21
  %798 = load double, ptr %797, align 8, !tbaa !38
  %799 = fdiv double %795, %798
  store double %799, ptr %27, align 8, !tbaa !9
  %800 = load double, ptr %27, align 8, !tbaa !9
  %801 = call double @pow(double noundef %800, double noundef 2.000000e+00) #6
  %802 = load double, ptr %27, align 8, !tbaa !9
  %803 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %802, double %801)
  %804 = call double @sqrt(double noundef %803) #6
  %805 = load double, ptr %27, align 8, !tbaa !9
  %806 = fsub double %804, %805
  %807 = fmul double 5.000000e-01, %806
  store double %807, ptr %18, align 8, !tbaa !9
  %808 = load double, ptr %25, align 8, !tbaa !9
  %809 = load double, ptr %26, align 8, !tbaa !9
  %810 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %811 = load ptr, ptr %7, align 8, !tbaa !5
  %812 = getelementptr inbounds %struct.precision, ptr %811, i32 0, i32 19
  %813 = load double, ptr %812, align 8, !tbaa !83
  %814 = load ptr, ptr %7, align 8, !tbaa !5
  %815 = getelementptr inbounds %struct.precision, ptr %814, i32 0, i32 175
  %816 = load double, ptr %815, align 8, !tbaa !84
  %817 = call i32 @generic_integrator(ptr noundef @thermodynamics_derivs_with_recfast, double noundef %808, double noundef %809, ptr noundef %810, ptr noundef %35, double noundef %813, double noundef %816, ptr noundef %34)
  %818 = icmp eq i32 %817, 1
  br i1 %818, label %819, label %827

819:                                              ; preds = %771
  call void @llvm.lifetime.start.p0(i64 2048, ptr %46) #6
  %820 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  %821 = getelementptr inbounds %struct.generic_integrator_workspace, ptr %34, i32 0, i32 13
  %822 = getelementptr inbounds [2048 x i8], ptr %821, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %820, ptr noundef @.str, ptr noundef @.str.220, ptr noundef %822)
  %823 = load ptr, ptr %9, align 8, !tbaa !5
  %824 = getelementptr inbounds %struct.thermo, ptr %823, i32 0, i32 93
  %825 = getelementptr inbounds [2048 x i8], ptr %824, i64 0, i64 0
  %826 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %825, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3707, ptr noundef %826)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %46) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

827:                                              ; preds = %771
  %828 = load double, ptr %18, align 8, !tbaa !9
  %829 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  store double %828, ptr %829, align 16, !tbaa !9
  %830 = load ptr, ptr %7, align 8, !tbaa !5
  %831 = getelementptr inbounds %struct.precision, ptr %830, i32 0, i32 37
  %832 = load double, ptr %831, align 8, !tbaa !81
  %833 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  %834 = load double, ptr %833, align 8, !tbaa !9
  %835 = fsub double %832, %834
  %836 = load ptr, ptr %7, align 8, !tbaa !5
  %837 = getelementptr inbounds %struct.precision, ptr %836, i32 0, i32 39
  %838 = load double, ptr %837, align 8, !tbaa !85
  %839 = fcmp olt double %835, %838
  br i1 %839, label %840, label %918

840:                                              ; preds = %827
  %841 = load ptr, ptr %10, align 8, !tbaa !5
  %842 = getelementptr inbounds %struct.recombination, ptr %841, i32 0, i32 10
  %843 = load double, ptr %842, align 8, !tbaa !63
  %844 = load ptr, ptr %10, align 8, !tbaa !5
  %845 = getelementptr inbounds %struct.recombination, ptr %844, i32 0, i32 20
  %846 = load double, ptr %845, align 8, !tbaa !28
  %847 = fmul double %843, %846
  %848 = load double, ptr %19, align 8, !tbaa !9
  %849 = fadd double 1.000000e+00, %848
  %850 = fdiv double %847, %849
  %851 = call double @log(double noundef %850) #6
  %852 = load ptr, ptr %10, align 8, !tbaa !5
  %853 = getelementptr inbounds %struct.recombination, ptr %852, i32 0, i32 24
  %854 = load double, ptr %853, align 8, !tbaa !61
  %855 = load ptr, ptr %10, align 8, !tbaa !5
  %856 = getelementptr inbounds %struct.recombination, ptr %855, i32 0, i32 20
  %857 = load double, ptr %856, align 8, !tbaa !28
  %858 = load double, ptr %19, align 8, !tbaa !9
  %859 = fadd double 1.000000e+00, %858
  %860 = fmul double %857, %859
  %861 = fdiv double %854, %860
  %862 = fneg double %861
  %863 = call double @llvm.fmuladd.f64(double 1.500000e+00, double %851, double %862)
  %864 = call double @exp(double noundef %863) #6
  %865 = fmul double 4.000000e+00, %864
  %866 = load ptr, ptr %10, align 8, !tbaa !5
  %867 = getelementptr inbounds %struct.recombination, ptr %866, i32 0, i32 21
  %868 = load double, ptr %867, align 8, !tbaa !38
  %869 = fdiv double %865, %868
  store double %869, ptr %27, align 8, !tbaa !9
  %870 = load double, ptr %27, align 8, !tbaa !9
  %871 = fsub double %870, 1.000000e+00
  %872 = call double @pow(double noundef %871, double noundef 2.000000e+00) #6
  %873 = load ptr, ptr %10, align 8, !tbaa !5
  %874 = getelementptr inbounds %struct.recombination, ptr %873, i32 0, i32 14
  %875 = load double, ptr %874, align 8, !tbaa !37
  %876 = fadd double 1.000000e+00, %875
  %877 = fmul double 4.000000e+00, %876
  %878 = load double, ptr %27, align 8, !tbaa !9
  %879 = call double @llvm.fmuladd.f64(double %877, double %878, double %872)
  %880 = call double @sqrt(double noundef %879) #6
  %881 = load double, ptr %27, align 8, !tbaa !9
  %882 = fsub double %881, 1.000000e+00
  %883 = fsub double %880, %882
  %884 = fmul double 5.000000e-01, %883
  store double %884, ptr %30, align 8, !tbaa !9
  %885 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %886 = load double, ptr %885, align 16, !tbaa !9
  %887 = load ptr, ptr %10, align 8, !tbaa !5
  %888 = getelementptr inbounds %struct.recombination, ptr %887, i32 0, i32 14
  %889 = load double, ptr %888, align 8, !tbaa !37
  %890 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  %891 = load double, ptr %890, align 8, !tbaa !9
  %892 = call double @llvm.fmuladd.f64(double %889, double %891, double %886)
  store double %892, ptr %31, align 8, !tbaa !9
  %893 = load ptr, ptr %7, align 8, !tbaa !5
  %894 = getelementptr inbounds %struct.precision, ptr %893, i32 0, i32 37
  %895 = load double, ptr %894, align 8, !tbaa !81
  %896 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  %897 = load double, ptr %896, align 8, !tbaa !9
  %898 = fsub double %895, %897
  %899 = load ptr, ptr %7, align 8, !tbaa !5
  %900 = getelementptr inbounds %struct.precision, ptr %899, i32 0, i32 39
  %901 = load double, ptr %900, align 8, !tbaa !85
  %902 = fdiv double %898, %901
  store double %902, ptr %32, align 8, !tbaa !9
  %903 = load double, ptr %32, align 8, !tbaa !9
  %904 = load double, ptr %32, align 8, !tbaa !9
  %905 = fmul double %903, %904
  %906 = load double, ptr %32, align 8, !tbaa !9
  %907 = fdiv double %906, 3.000000e+00
  %908 = fsub double 5.000000e-01, %907
  %909 = fmul double %905, %908
  %910 = fmul double %909, 6.000000e+00
  store double %910, ptr %33, align 8, !tbaa !9
  %911 = load double, ptr %33, align 8, !tbaa !9
  %912 = load double, ptr %31, align 8, !tbaa !9
  %913 = load double, ptr %33, align 8, !tbaa !9
  %914 = fsub double 1.000000e+00, %913
  %915 = load double, ptr %30, align 8, !tbaa !9
  %916 = fmul double %914, %915
  %917 = call double @llvm.fmuladd.f64(double %911, double %912, double %916)
  store double %917, ptr %17, align 8, !tbaa !9
  br label %927

918:                                              ; preds = %827
  %919 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %920 = load double, ptr %919, align 16, !tbaa !9
  %921 = load ptr, ptr %10, align 8, !tbaa !5
  %922 = getelementptr inbounds %struct.recombination, ptr %921, i32 0, i32 14
  %923 = load double, ptr %922, align 8, !tbaa !37
  %924 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  %925 = load double, ptr %924, align 8, !tbaa !9
  %926 = call double @llvm.fmuladd.f64(double %923, double %925, double %920)
  store double %926, ptr %17, align 8, !tbaa !9
  br label %927

927:                                              ; preds = %918, %840
  br label %1050

928:                                              ; preds = %764
  %929 = load ptr, ptr %7, align 8, !tbaa !5
  %930 = getelementptr inbounds %struct.precision, ptr %929, i32 0, i32 40
  %931 = load double, ptr %930, align 8, !tbaa !82
  %932 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %933 = load double, ptr %932, align 16, !tbaa !9
  %934 = fsub double %931, %933
  %935 = load ptr, ptr %7, align 8, !tbaa !5
  %936 = getelementptr inbounds %struct.precision, ptr %935, i32 0, i32 42
  %937 = load double, ptr %936, align 8, !tbaa !86
  %938 = fcmp olt double %934, %937
  br i1 %938, label %939, label %976

939:                                              ; preds = %928
  %940 = load ptr, ptr %10, align 8, !tbaa !5
  %941 = getelementptr inbounds %struct.recombination, ptr %940, i32 0, i32 10
  %942 = load double, ptr %941, align 8, !tbaa !63
  %943 = load ptr, ptr %10, align 8, !tbaa !5
  %944 = getelementptr inbounds %struct.recombination, ptr %943, i32 0, i32 20
  %945 = load double, ptr %944, align 8, !tbaa !28
  %946 = fmul double %942, %945
  %947 = load double, ptr %19, align 8, !tbaa !9
  %948 = fadd double 1.000000e+00, %947
  %949 = fdiv double %946, %948
  %950 = call double @log(double noundef %949) #6
  %951 = load ptr, ptr %10, align 8, !tbaa !5
  %952 = getelementptr inbounds %struct.recombination, ptr %951, i32 0, i32 23
  %953 = load double, ptr %952, align 8, !tbaa !60
  %954 = load ptr, ptr %10, align 8, !tbaa !5
  %955 = getelementptr inbounds %struct.recombination, ptr %954, i32 0, i32 20
  %956 = load double, ptr %955, align 8, !tbaa !28
  %957 = load double, ptr %19, align 8, !tbaa !9
  %958 = fadd double 1.000000e+00, %957
  %959 = fmul double %956, %958
  %960 = fdiv double %953, %959
  %961 = fneg double %960
  %962 = call double @llvm.fmuladd.f64(double 1.500000e+00, double %950, double %961)
  %963 = call double @exp(double noundef %962) #6
  %964 = load ptr, ptr %10, align 8, !tbaa !5
  %965 = getelementptr inbounds %struct.recombination, ptr %964, i32 0, i32 21
  %966 = load double, ptr %965, align 8, !tbaa !38
  %967 = fdiv double %963, %966
  store double %967, ptr %27, align 8, !tbaa !9
  %968 = load double, ptr %27, align 8, !tbaa !9
  %969 = call double @pow(double noundef %968, double noundef 2.000000e+00) #6
  %970 = load double, ptr %27, align 8, !tbaa !9
  %971 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %970, double %969)
  %972 = call double @sqrt(double noundef %971) #6
  %973 = load double, ptr %27, align 8, !tbaa !9
  %974 = fsub double %972, %973
  %975 = fmul double 5.000000e-01, %974
  store double %975, ptr %18, align 8, !tbaa !9
  br label %976

976:                                              ; preds = %939, %928
  %977 = load double, ptr %25, align 8, !tbaa !9
  %978 = load double, ptr %26, align 8, !tbaa !9
  %979 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %980 = load ptr, ptr %7, align 8, !tbaa !5
  %981 = getelementptr inbounds %struct.precision, ptr %980, i32 0, i32 19
  %982 = load double, ptr %981, align 8, !tbaa !83
  %983 = load ptr, ptr %7, align 8, !tbaa !5
  %984 = getelementptr inbounds %struct.precision, ptr %983, i32 0, i32 175
  %985 = load double, ptr %984, align 8, !tbaa !84
  %986 = call i32 @generic_integrator(ptr noundef @thermodynamics_derivs_with_recfast, double noundef %977, double noundef %978, ptr noundef %979, ptr noundef %35, double noundef %982, double noundef %985, ptr noundef %34)
  %987 = icmp eq i32 %986, 1
  br i1 %987, label %988, label %996

988:                                              ; preds = %976
  call void @llvm.lifetime.start.p0(i64 2048, ptr %47) #6
  %989 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  %990 = getelementptr inbounds %struct.generic_integrator_workspace, ptr %34, i32 0, i32 13
  %991 = getelementptr inbounds [2048 x i8], ptr %990, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %989, ptr noundef @.str, ptr noundef @.str.220, ptr noundef %991)
  %992 = load ptr, ptr %9, align 8, !tbaa !5
  %993 = getelementptr inbounds %struct.thermo, ptr %992, i32 0, i32 93
  %994 = getelementptr inbounds [2048 x i8], ptr %993, i64 0, i64 0
  %995 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %994, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3749, ptr noundef %995)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %47) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

996:                                              ; preds = %976
  %997 = load ptr, ptr %7, align 8, !tbaa !5
  %998 = getelementptr inbounds %struct.precision, ptr %997, i32 0, i32 40
  %999 = load double, ptr %998, align 8, !tbaa !82
  %1000 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %1001 = load double, ptr %1000, align 16, !tbaa !9
  %1002 = fsub double %999, %1001
  %1003 = load ptr, ptr %7, align 8, !tbaa !5
  %1004 = getelementptr inbounds %struct.precision, ptr %1003, i32 0, i32 42
  %1005 = load double, ptr %1004, align 8, !tbaa !86
  %1006 = fcmp olt double %1002, %1005
  br i1 %1006, label %1007, label %1040

1007:                                             ; preds = %996
  %1008 = load ptr, ptr %7, align 8, !tbaa !5
  %1009 = getelementptr inbounds %struct.precision, ptr %1008, i32 0, i32 40
  %1010 = load double, ptr %1009, align 8, !tbaa !82
  %1011 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %1012 = load double, ptr %1011, align 16, !tbaa !9
  %1013 = fsub double %1010, %1012
  %1014 = load ptr, ptr %7, align 8, !tbaa !5
  %1015 = getelementptr inbounds %struct.precision, ptr %1014, i32 0, i32 42
  %1016 = load double, ptr %1015, align 8, !tbaa !86
  %1017 = fdiv double %1013, %1016
  store double %1017, ptr %32, align 8, !tbaa !9
  %1018 = load double, ptr %32, align 8, !tbaa !9
  %1019 = load double, ptr %32, align 8, !tbaa !9
  %1020 = fmul double %1018, %1019
  %1021 = load double, ptr %32, align 8, !tbaa !9
  %1022 = fdiv double %1021, 3.000000e+00
  %1023 = fsub double 5.000000e-01, %1022
  %1024 = fmul double %1020, %1023
  %1025 = fmul double %1024, 6.000000e+00
  store double %1025, ptr %33, align 8, !tbaa !9
  %1026 = load double, ptr %33, align 8, !tbaa !9
  %1027 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %1028 = load double, ptr %1027, align 16, !tbaa !9
  %1029 = load double, ptr %33, align 8, !tbaa !9
  %1030 = fsub double 1.000000e+00, %1029
  %1031 = load double, ptr %18, align 8, !tbaa !9
  %1032 = fmul double %1030, %1031
  %1033 = call double @llvm.fmuladd.f64(double %1026, double %1028, double %1032)
  %1034 = load ptr, ptr %10, align 8, !tbaa !5
  %1035 = getelementptr inbounds %struct.recombination, ptr %1034, i32 0, i32 14
  %1036 = load double, ptr %1035, align 8, !tbaa !37
  %1037 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  %1038 = load double, ptr %1037, align 8, !tbaa !9
  %1039 = call double @llvm.fmuladd.f64(double %1036, double %1038, double %1033)
  store double %1039, ptr %17, align 8, !tbaa !9
  br label %1049

1040:                                             ; preds = %996
  %1041 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %1042 = load double, ptr %1041, align 16, !tbaa !9
  %1043 = load ptr, ptr %10, align 8, !tbaa !5
  %1044 = getelementptr inbounds %struct.recombination, ptr %1043, i32 0, i32 14
  %1045 = load double, ptr %1044, align 8, !tbaa !37
  %1046 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 1
  %1047 = load double, ptr %1046, align 8, !tbaa !9
  %1048 = call double @llvm.fmuladd.f64(double %1045, double %1047, double %1042)
  store double %1048, ptr %17, align 8, !tbaa !9
  br label %1049

1049:                                             ; preds = %1040, %1007
  br label %1050

1050:                                             ; preds = %1049, %927
  br label %1051

1051:                                             ; preds = %1050, %746
  br label %1052

1052:                                             ; preds = %1051, %641
  br label %1053

1053:                                             ; preds = %1052, %524
  br label %1054

1054:                                             ; preds = %1053, %384
  %1055 = load double, ptr %26, align 8, !tbaa !9
  %1056 = load ptr, ptr %10, align 8, !tbaa !5
  %1057 = getelementptr inbounds %struct.recombination, ptr %1056, i32 0, i32 8
  %1058 = load ptr, ptr %1057, align 8, !tbaa !18
  %1059 = load i32, ptr %29, align 4, !tbaa !19
  %1060 = load i32, ptr %28, align 4, !tbaa !19
  %1061 = sub nsw i32 %1059, %1060
  %1062 = sub nsw i32 %1061, 1
  %1063 = load ptr, ptr %10, align 8, !tbaa !5
  %1064 = getelementptr inbounds %struct.recombination, ptr %1063, i32 0, i32 6
  %1065 = load i32, ptr %1064, align 8, !tbaa !17
  %1066 = mul nsw i32 %1062, %1065
  %1067 = sext i32 %1066 to i64
  %1068 = getelementptr inbounds double, ptr %1058, i64 %1067
  %1069 = load ptr, ptr %10, align 8, !tbaa !5
  %1070 = getelementptr inbounds %struct.recombination, ptr %1069, i32 0, i32 0
  %1071 = load i32, ptr %1070, align 8, !tbaa !87
  %1072 = sext i32 %1071 to i64
  %1073 = getelementptr inbounds double, ptr %1068, i64 %1072
  store double %1055, ptr %1073, align 8, !tbaa !9
  %1074 = load double, ptr %17, align 8, !tbaa !9
  %1075 = load ptr, ptr %10, align 8, !tbaa !5
  %1076 = getelementptr inbounds %struct.recombination, ptr %1075, i32 0, i32 8
  %1077 = load ptr, ptr %1076, align 8, !tbaa !18
  %1078 = load i32, ptr %29, align 4, !tbaa !19
  %1079 = load i32, ptr %28, align 4, !tbaa !19
  %1080 = sub nsw i32 %1078, %1079
  %1081 = sub nsw i32 %1080, 1
  %1082 = load ptr, ptr %10, align 8, !tbaa !5
  %1083 = getelementptr inbounds %struct.recombination, ptr %1082, i32 0, i32 6
  %1084 = load i32, ptr %1083, align 8, !tbaa !17
  %1085 = mul nsw i32 %1081, %1084
  %1086 = sext i32 %1085 to i64
  %1087 = getelementptr inbounds double, ptr %1077, i64 %1086
  %1088 = load ptr, ptr %10, align 8, !tbaa !5
  %1089 = getelementptr inbounds %struct.recombination, ptr %1088, i32 0, i32 1
  %1090 = load i32, ptr %1089, align 4, !tbaa !88
  %1091 = sext i32 %1090 to i64
  %1092 = getelementptr inbounds double, ptr %1087, i64 %1091
  store double %1074, ptr %1092, align 8, !tbaa !9
  %1093 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  %1094 = load double, ptr %1093, align 16, !tbaa !9
  %1095 = load ptr, ptr %10, align 8, !tbaa !5
  %1096 = getelementptr inbounds %struct.recombination, ptr %1095, i32 0, i32 8
  %1097 = load ptr, ptr %1096, align 8, !tbaa !18
  %1098 = load i32, ptr %29, align 4, !tbaa !19
  %1099 = load i32, ptr %28, align 4, !tbaa !19
  %1100 = sub nsw i32 %1098, %1099
  %1101 = sub nsw i32 %1100, 1
  %1102 = load ptr, ptr %10, align 8, !tbaa !5
  %1103 = getelementptr inbounds %struct.recombination, ptr %1102, i32 0, i32 6
  %1104 = load i32, ptr %1103, align 8, !tbaa !17
  %1105 = mul nsw i32 %1101, %1104
  %1106 = sext i32 %1105 to i64
  %1107 = getelementptr inbounds double, ptr %1097, i64 %1106
  %1108 = load ptr, ptr %10, align 8, !tbaa !5
  %1109 = getelementptr inbounds %struct.recombination, ptr %1108, i32 0, i32 2
  %1110 = load i32, ptr %1109, align 8, !tbaa !89
  %1111 = sext i32 %1110 to i64
  %1112 = getelementptr inbounds double, ptr %1107, i64 %1111
  store double %1094, ptr %1112, align 8, !tbaa !9
  %1113 = load double, ptr %26, align 8, !tbaa !9
  %1114 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 0
  %1115 = getelementptr inbounds [3 x double], ptr %13, i64 0, i64 0
  %1116 = load ptr, ptr %9, align 8, !tbaa !5
  %1117 = getelementptr inbounds %struct.thermo, ptr %1116, i32 0, i32 93
  %1118 = getelementptr inbounds [2048 x i8], ptr %1117, i64 0, i64 0
  %1119 = call i32 @thermodynamics_derivs_with_recfast(double noundef %1113, ptr noundef %1114, ptr noundef %1115, ptr noundef %35, ptr noundef %1118)
  %1120 = icmp eq i32 %1119, 1
  br i1 %1120, label %1121, label %1130

1121:                                             ; preds = %1054
  call void @llvm.lifetime.start.p0(i64 2048, ptr %48) #6
  %1122 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  %1123 = load ptr, ptr %9, align 8, !tbaa !5
  %1124 = getelementptr inbounds %struct.thermo, ptr %1123, i32 0, i32 93
  %1125 = getelementptr inbounds [2048 x i8], ptr %1124, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1122, ptr noundef @.str, ptr noundef @.str.221, ptr noundef %1125)
  %1126 = load ptr, ptr %9, align 8, !tbaa !5
  %1127 = getelementptr inbounds %struct.thermo, ptr %1126, i32 0, i32 93
  %1128 = getelementptr inbounds [2048 x i8], ptr %1127, i64 0, i64 0
  %1129 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1128, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3782, ptr noundef %1129)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %48) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

1130:                                             ; preds = %1054
  %1131 = load ptr, ptr %10, align 8, !tbaa !5
  %1132 = getelementptr inbounds %struct.recombination, ptr %1131, i32 0, i32 27
  %1133 = load double, ptr %1132, align 8, !tbaa !26
  %1134 = call double @llvm.fmuladd.f64(double 0xBFE7F14DA42541E8, double %1133, double 1.000000e+00)
  %1135 = load double, ptr %17, align 8, !tbaa !9
  %1136 = load ptr, ptr %10, align 8, !tbaa !5
  %1137 = getelementptr inbounds %struct.recombination, ptr %1136, i32 0, i32 27
  %1138 = load double, ptr %1137, align 8, !tbaa !26
  %1139 = fsub double 1.000000e+00, %1138
  %1140 = call double @llvm.fmuladd.f64(double %1135, double %1139, double %1134)
  %1141 = fmul double 0x3D39D63187BB40C4, %1140
  %1142 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  %1143 = load double, ptr %1142, align 16, !tbaa !9
  %1144 = fmul double %1141, %1143
  %1145 = load ptr, ptr %10, align 8, !tbaa !5
  %1146 = getelementptr inbounds %struct.recombination, ptr %1145, i32 0, i32 8
  %1147 = load ptr, ptr %1146, align 8, !tbaa !18
  %1148 = load i32, ptr %29, align 4, !tbaa !19
  %1149 = load i32, ptr %28, align 4, !tbaa !19
  %1150 = sub nsw i32 %1148, %1149
  %1151 = sub nsw i32 %1150, 1
  %1152 = load ptr, ptr %10, align 8, !tbaa !5
  %1153 = getelementptr inbounds %struct.recombination, ptr %1152, i32 0, i32 6
  %1154 = load i32, ptr %1153, align 8, !tbaa !17
  %1155 = mul nsw i32 %1151, %1154
  %1156 = sext i32 %1155 to i64
  %1157 = getelementptr inbounds double, ptr %1147, i64 %1156
  %1158 = load ptr, ptr %10, align 8, !tbaa !5
  %1159 = getelementptr inbounds %struct.recombination, ptr %1158, i32 0, i32 3
  %1160 = load i32, ptr %1159, align 4, !tbaa !90
  %1161 = sext i32 %1160 to i64
  %1162 = getelementptr inbounds double, ptr %1157, i64 %1161
  store double %1144, ptr %1162, align 8, !tbaa !9
  %1163 = load ptr, ptr %10, align 8, !tbaa !5
  %1164 = getelementptr inbounds %struct.recombination, ptr %1163, i32 0, i32 8
  %1165 = load ptr, ptr %1164, align 8, !tbaa !18
  %1166 = load i32, ptr %29, align 4, !tbaa !19
  %1167 = load i32, ptr %28, align 4, !tbaa !19
  %1168 = sub nsw i32 %1166, %1167
  %1169 = sub nsw i32 %1168, 1
  %1170 = load ptr, ptr %10, align 8, !tbaa !5
  %1171 = getelementptr inbounds %struct.recombination, ptr %1170, i32 0, i32 6
  %1172 = load i32, ptr %1171, align 8, !tbaa !17
  %1173 = mul nsw i32 %1169, %1172
  %1174 = sext i32 %1173 to i64
  %1175 = getelementptr inbounds double, ptr %1165, i64 %1174
  %1176 = load ptr, ptr %10, align 8, !tbaa !5
  %1177 = getelementptr inbounds %struct.recombination, ptr %1176, i32 0, i32 3
  %1178 = load i32, ptr %1177, align 4, !tbaa !90
  %1179 = sext i32 %1178 to i64
  %1180 = getelementptr inbounds double, ptr %1175, i64 %1179
  %1181 = load double, ptr %1180, align 8, !tbaa !9
  %1182 = load double, ptr %26, align 8, !tbaa !9
  %1183 = fadd double 1.000000e+00, %1182
  %1184 = getelementptr inbounds [3 x double], ptr %13, i64 0, i64 2
  %1185 = load double, ptr %1184, align 16, !tbaa !9
  %1186 = fmul double %1183, %1185
  %1187 = getelementptr inbounds [3 x double], ptr %12, i64 0, i64 2
  %1188 = load double, ptr %1187, align 16, !tbaa !9
  %1189 = fdiv double %1186, %1188
  %1190 = fdiv double %1189, 3.000000e+00
  %1191 = fadd double 1.000000e+00, %1190
  %1192 = fmul double %1181, %1191
  %1193 = load ptr, ptr %10, align 8, !tbaa !5
  %1194 = getelementptr inbounds %struct.recombination, ptr %1193, i32 0, i32 8
  %1195 = load ptr, ptr %1194, align 8, !tbaa !18
  %1196 = load i32, ptr %29, align 4, !tbaa !19
  %1197 = load i32, ptr %28, align 4, !tbaa !19
  %1198 = sub nsw i32 %1196, %1197
  %1199 = sub nsw i32 %1198, 1
  %1200 = load ptr, ptr %10, align 8, !tbaa !5
  %1201 = getelementptr inbounds %struct.recombination, ptr %1200, i32 0, i32 6
  %1202 = load i32, ptr %1201, align 8, !tbaa !17
  %1203 = mul nsw i32 %1199, %1202
  %1204 = sext i32 %1203 to i64
  %1205 = getelementptr inbounds double, ptr %1195, i64 %1204
  %1206 = load ptr, ptr %10, align 8, !tbaa !5
  %1207 = getelementptr inbounds %struct.recombination, ptr %1206, i32 0, i32 4
  %1208 = load i32, ptr %1207, align 8, !tbaa !91
  %1209 = sext i32 %1208 to i64
  %1210 = getelementptr inbounds double, ptr %1205, i64 %1209
  store double %1192, ptr %1210, align 8, !tbaa !9
  %1211 = load double, ptr %26, align 8, !tbaa !9
  %1212 = fadd double 1.000000e+00, %1211
  %1213 = load double, ptr %26, align 8, !tbaa !9
  %1214 = fadd double 1.000000e+00, %1213
  %1215 = fmul double %1212, %1214
  %1216 = load ptr, ptr %10, align 8, !tbaa !5
  %1217 = getelementptr inbounds %struct.recombination, ptr %1216, i32 0, i32 21
  %1218 = load double, ptr %1217, align 8, !tbaa !38
  %1219 = fmul double %1215, %1218
  %1220 = load double, ptr %17, align 8, !tbaa !9
  %1221 = fmul double %1219, %1220
  %1222 = fmul double %1221, 0x3A15151E37FBF236
  %1223 = fmul double %1222, 0x449A22FF260E3826
  %1224 = load ptr, ptr %10, align 8, !tbaa !5
  %1225 = getelementptr inbounds %struct.recombination, ptr %1224, i32 0, i32 8
  %1226 = load ptr, ptr %1225, align 8, !tbaa !18
  %1227 = load i32, ptr %29, align 4, !tbaa !19
  %1228 = load i32, ptr %28, align 4, !tbaa !19
  %1229 = sub nsw i32 %1227, %1228
  %1230 = sub nsw i32 %1229, 1
  %1231 = load ptr, ptr %10, align 8, !tbaa !5
  %1232 = getelementptr inbounds %struct.recombination, ptr %1231, i32 0, i32 6
  %1233 = load i32, ptr %1232, align 8, !tbaa !17
  %1234 = mul nsw i32 %1230, %1233
  %1235 = sext i32 %1234 to i64
  %1236 = getelementptr inbounds double, ptr %1226, i64 %1235
  %1237 = load ptr, ptr %10, align 8, !tbaa !5
  %1238 = getelementptr inbounds %struct.recombination, ptr %1237, i32 0, i32 5
  %1239 = load i32, ptr %1238, align 4, !tbaa !92
  %1240 = sext i32 %1239 to i64
  %1241 = getelementptr inbounds double, ptr %1236, i64 %1240
  store double %1223, ptr %1241, align 8, !tbaa !9
  br label %1242

1242:                                             ; preds = %1130
  %1243 = load i32, ptr %28, align 4, !tbaa !19
  %1244 = add nsw i32 %1243, 1
  store i32 %1244, ptr %28, align 4, !tbaa !19
  br label %350, !llvm.loop !93

1245:                                             ; preds = %350
  %1246 = call i32 @cleanup_generic_integrator(ptr noundef %34)
  %1247 = icmp eq i32 %1246, 1
  br i1 %1247, label %1248, label %1256

1248:                                             ; preds = %1245
  call void @llvm.lifetime.start.p0(i64 2048, ptr %49) #6
  %1249 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  %1250 = getelementptr inbounds %struct.generic_integrator_workspace, ptr %34, i32 0, i32 13
  %1251 = getelementptr inbounds [2048 x i8], ptr %1250, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1249, ptr noundef @.str, ptr noundef @.str.222, ptr noundef %1251)
  %1252 = load ptr, ptr %9, align 8, !tbaa !5
  %1253 = getelementptr inbounds %struct.thermo, ptr %1252, i32 0, i32 93
  %1254 = getelementptr inbounds [2048 x i8], ptr %1253, i64 0, i64 0
  %1255 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1254, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3812, ptr noundef %1255)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %49) #6
  store i32 1, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

1256:                                             ; preds = %1245
  store i32 0, ptr %6, align 4
  store i32 1, ptr %38, align 4
  br label %1257

1257:                                             ; preds = %1256, %1248, %1121, %988, %819, %328, %176, %140, %91, %71
  call void @llvm.lifetime.end.p0(i64 32, ptr %35) #6
  call void @llvm.lifetime.end.p0(i64 2152, ptr %34) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %33) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %32) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %31) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %30) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %29) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %28) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %27) #6
  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 8, 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 8, 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 24, ptr %13) #6
  call void @llvm.lifetime.end.p0(i64 24, ptr %12) #6
  %1258 = load i32, ptr %6, align 4
  ret i32 %1258
}

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

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

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

declare i32 @cleanup_generic_integrator(ptr noundef) #3

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nounwind "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #3 = { "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #4 = { nounwind allocsize(0) "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #5 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #6 = { 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 = !{!"double", !7, i64 0}
!11 = !{!12, !13, i64 368}
!12 = !{!"precision", !10, i64 0, !10, i64 8, !10, i64 16, !10, i64 24, !10, i64 32, !10, i64 40, !10, i64 48, !10, i64 56, !10, i64 64, !10, i64 72, !10, i64 80, !10, i64 88, !10, i64 96, !7, i64 104, !10, i64 360, !13, i64 368, !10, i64 376, !13, i64 384, !13, i64 388, !10, i64 392, !13, i64 400, !10, i64 408, !13, i64 416, !10, i64 424, !10, i64 432, !10, i64 440, !10, i64 448, !10, i64 456, !10, i64 464, !10, i64 472, !10, i64 480, !10, i64 488, !10, i64 496, !10, i64 504, !10, i64 512, !10, i64 520, !10, i64 528, !10, i64 536, !10, i64 544, !10, i64 552, !10, i64 560, !10, i64 568, !10, i64 576, !10, i64 584, !10, i64 592, !10, i64 600, !10, i64 608, !10, i64 616, !13, i64 624, !7, i64 628, !7, i64 884, !7, i64 1140, !10, i64 1400, !10, i64 1408, !10, i64 1416, !10, i64 1424, !10, i64 1432, !10, i64 1440, !10, i64 1448, !10, i64 1456, !10, i64 1464, !10, i64 1472, !10, i64 1480, !10, i64 1488, !10, i64 1496, !10, i64 1504, !10, i64 1512, !10, i64 1520, !13, i64 1528, !10, i64 1536, !10, i64 1544, !13, i64 1552, !13, i64 1556, !13, i64 1560, !13, i64 1564, !13, i64 1568, !13, i64 1572, !13, i64 1576, !13, i64 1580, !10, i64 1584, !10, i64 1592, !10, i64 1600, !10, i64 1608, !10, i64 1616, !10, i64 1624, !10, i64 1632, !10, i64 1640, !13, i64 1648, !10, i64 1656, !10, i64 1664, !13, i64 1672, !10, i64 1680, !10, i64 1688, !13, i64 1696, !10, i64 1704, !13, i64 1712, !10, i64 1720, !10, i64 1728, !13, i64 1736, !10, i64 1744, !10, i64 1752, !10, i64 1760, !13, i64 1768, !10, i64 1776, !10, i64 1784, !10, i64 1792, !10, i64 1800, !10, i64 1808, !13, i64 1816, !10, i64 1824, !10, i64 1832, !10, i64 1840, !10, i64 1848, !10, i64 1856, !10, i64 1864, !10, i64 1872, !13, i64 1880, !10, i64 1888, !10, i64 1896, !10, i64 1904, !10, i64 1912, !10, i64 1920, !10, i64 1928, !10, i64 1936, !10, i64 1944, !10, i64 1952, !10, i64 1960, !10, i64 1968, !10, i64 1976, !10, i64 1984, !10, i64 1992, !10, i64 2000, !10, i64 2008, !10, i64 2016, !10, i64 2024, !10, i64 2032, !10, i64 2040, !10, i64 2048, !10, i64 2056, !10, i64 2064, !10, i64 2072, !10, i64 2080, !10, i64 2088, !10, i64 2096, !10, i64 2104, !10, i64 2112, !10, i64 2120, !10, i64 2128, !10, i64 2136, !10, i64 2144, !10, i64 2152, !10, i64 2160, !10, i64 2168, !10, i64 2176, !10, i64 2184, !10, i64 2192, !10, i64 2200, !10, i64 2208, !10, i64 2216, !10, i64 2224, !10, i64 2232, !10, i64 2240, !10, i64 2248, !13, i64 2256, !10, i64 2264, !10, i64 2272, !10, i64 2280, !10, i64 2288, !13, i64 2296, !10, i64 2304, !10, i64 2312, !13, i64 2320, !13, i64 2324, !13, i64 2328, !10, i64 2336, !10, i64 2344, !7, i64 2352}
!13 = !{!"int", !7, i64 0}
!14 = !{!15, !13, i64 28}
!15 = !{!"recombination", !13, i64 0, !13, i64 4, !13, i64 8, !13, i64 12, !13, i64 16, !13, i64 20, !13, i64 24, !13, i64 28, !6, i64 32, !10, i64 40, !10, i64 48, !10, i64 56, !10, i64 64, !10, i64 72, !10, i64 80, !10, i64 88, !10, i64 96, !10, i64 104, !10, i64 112, !10, i64 120, !10, i64 128, !10, i64 136, !10, i64 144, !10, i64 152, !10, i64 160, !10, i64 168, !10, i64 176, !10, i64 184, !10, i64 192, !16, i64 200, !10, i64 208, !10, i64 216, !10, i64 224, !10, i64 232, !10, i64 240, !10, i64 248, !10, i64 256}
!16 = !{!"short", !7, i64 0}
!17 = !{!15, !13, i64 24}
!18 = !{!15, !6, i64 32}
!19 = !{!13, !13, i64 0}
!20 = !{!21, !10, i64 0}
!21 = !{!"background", !10, i64 0, !10, i64 8, !10, i64 16, !10, i64 24, !10, i64 32, !10, i64 40, !10, i64 48, !13, i64 56, !10, i64 64, !10, i64 72, !10, i64 80, !10, i64 88, !16, i64 96, !10, i64 104, !10, i64 112, !10, i64 120, !10, i64 128, !10, i64 136, !10, i64 144, !10, i64 152, !10, i64 160, !10, i64 168, !16, i64 176, !10, i64 184, !10, i64 192, !6, i64 200, !13, i64 208, !13, i64 212, !10, i64 216, !13, i64 224, !6, i64 232, !6, i64 240, !10, i64 248, !6, i64 256, !10, i64 264, !6, i64 272, !10, i64 280, !6, i64 288, !10, i64 296, !6, i64 304, !6, i64 312, !6, i64 320, !10, i64 328, !10, i64 336, !10, i64 344, !10, i64 352, !13, i64 360, !6, i64 368, !10, i64 376, !10, i64 384, !10, i64 392, !10, i64 400, !10, i64 408, !10, i64 416, !10, i64 424, !10, i64 432, !10, i64 440, !10, i64 448, !10, i64 456, !13, i64 464, !13, i64 468, !13, i64 472, !13, i64 476, !13, i64 480, !13, i64 484, !13, i64 488, !13, i64 492, !13, i64 496, !13, i64 500, !13, i64 504, !13, i64 508, !13, i64 512, !13, i64 516, !13, i64 520, !13, i64 524, !13, i64 528, !13, i64 532, !13, i64 536, !13, i64 540, !13, i64 544, !13, i64 548, !13, i64 552, !13, i64 556, !13, i64 560, !13, i64 564, !13, i64 568, !13, i64 572, !13, i64 576, !13, i64 580, !13, i64 584, !13, i64 588, !13, i64 592, !13, i64 596, !13, i64 600, !13, i64 604, !13, i64 608, !13, i64 612, !13, i64 616, !13, i64 620, !13, i64 624, !13, i64 628, !6, i64 632, !6, i64 640, !6, i64 648, !6, i64 656, !6, i64 664, !13, i64 672, !13, i64 676, !13, i64 680, !13, i64 684, !13, i64 688, !13, i64 692, !13, i64 696, !13, i64 700, !13, i64 704, !13, i64 708, !13, i64 712, !13, i64 716, !13, i64 720, !16, i64 724, !16, i64 726, !16, i64 728, !16, i64 730, !16, i64 732, !16, i64 734, !16, i64 736, !16, i64 738, !16, i64 740, !16, i64 742, !16, 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, !16, i64 840, !16, i64 842, !16, i64 844, !16, i64 846, !16, i64 848, !16, i64 850, !7, i64 852, !16, i64 2900, !7, i64 2902}
!22 = !{!15, !10, i64 176}
!23 = !{!21, !10, i64 24}
!24 = !{!25, !10, i64 0}
!25 = !{!"thermo", !10, i64 0, !13, i64 8, !13, i64 12, !13, i64 16, !10, i64 24, !10, i64 32, !16, i64 40, !16, i64 42, !10, i64 48, !10, i64 56, !10, i64 64, !10, i64 72, !13, i64 80, !6, i64 88, !6, i64 96, !10, i64 104, !13, i64 112, !6, i64 120, !6, i64 128, !10, i64 136, !13, i64 144, !6, i64 152, !6, i64 160, !10, i64 168, !16, i64 176, !10, i64 184, !10, i64 192, !10, i64 200, !10, i64 208, !10, i64 216, !10, i64 224, !10, i64 232, !10, i64 240, !10, i64 248, !10, i64 256, !10, i64 264, !13, i64 272, !13, i64 276, !13, i64 280, !13, i64 284, !13, i64 288, !13, i64 292, !13, i64 296, !13, i64 300, !13, i64 304, !13, i64 308, !13, i64 312, !13, i64 316, !13, i64 320, !13, i64 324, !13, i64 328, !13, i64 332, !13, i64 336, !13, i64 340, !13, i64 344, !13, i64 348, !13, i64 352, !13, i64 356, !13, i64 360, !13, i64 364, !13, i64 368, !13, i64 372, !13, i64 376, !6, i64 384, !6, i64 392, !6, i64 400, !10, i64 408, !10, i64 416, !10, i64 424, !10, i64 432, !10, i64 440, !10, i64 448, !10, i64 456, !10, i64 464, !10, i64 472, !10, i64 480, !10, i64 488, !10, i64 496, !10, i64 504, !10, i64 512, !10, i64 520, !10, i64 528, !10, i64 536, !10, i64 544, !10, i64 552, !10, i64 560, !10, i64 568, !10, i64 576, !10, i64 584, !10, i64 592, !16, i64 600, !16, i64 602, !16, i64 604, !7, i64 606}
!26 = !{!15, !10, i64 184}
!27 = !{!21, !10, i64 16}
!28 = !{!15, !10, i64 128}
!29 = !{!12, !10, i64 360}
!30 = !{!12, !10, i64 584}
!31 = !{!15, !10, i64 120}
!32 = !{!12, !13, i64 416}
!33 = !{!12, !10, i64 424}
!34 = !{!15, !10, i64 112}
!35 = !{!12, !10, i64 432}
!36 = !{!12, !13, i64 400}
!37 = !{!15, !10, i64 80}
!38 = !{!15, !10, i64 136}
!39 = !{!25, !10, i64 592}
!40 = !{!25, !10, i64 168}
!41 = !{!15, !10, i64 192}
!42 = !{!25, !16, i64 176}
!43 = !{!15, !16, i64 200}
!44 = !{!25, !10, i64 192}
!45 = !{!15, !10, i64 216}
!46 = !{!25, !10, i64 200}
!47 = !{!15, !10, i64 224}
!48 = !{!25, !10, i64 208}
!49 = !{!15, !10, i64 232}
!50 = !{!25, !10, i64 216}
!51 = !{!15, !10, i64 240}
!52 = !{!25, !10, i64 184}
!53 = !{!15, !10, i64 208}
!54 = !{!25, !10, i64 224}
!55 = !{!15, !10, i64 248}
!56 = !{!25, !10, i64 232}
!57 = !{!15, !10, i64 256}
!58 = !{!15, !10, i64 40}
!59 = !{!15, !10, i64 88}
!60 = !{!15, !10, i64 152}
!61 = !{!15, !10, i64 160}
!62 = !{!15, !10, i64 168}
!63 = !{!15, !10, i64 48}
!64 = !{!15, !10, i64 56}
!65 = !{!15, !10, i64 96}
!66 = !{!15, !10, i64 64}
!67 = !{!15, !10, i64 104}
!68 = !{!15, !10, i64 72}
!69 = !{!15, !10, i64 144}
!70 = !{!71, !6, i64 0}
!71 = !{!"thermodynamics_parameters_and_workspace", !6, i64 0, !6, i64 8, !6, i64 16, !6, i64 24}
!72 = !{!71, !6, i64 8}
!73 = !{!71, !6, i64 16}
!74 = !{!71, !6, i64 24}
!75 = !{!12, !10, i64 520}
!76 = !{!12, !10, i64 488}
!77 = !{!12, !10, i64 496}
!78 = !{!12, !10, i64 504}
!79 = !{!12, !10, i64 512}
!80 = !{!12, !10, i64 528}
!81 = !{!12, !10, i64 536}
!82 = !{!12, !10, i64 560}
!83 = !{!12, !10, i64 392}
!84 = !{!12, !10, i64 2344}
!85 = !{!12, !10, i64 552}
!86 = !{!12, !10, i64 576}
!87 = !{!15, !13, i64 0}
!88 = !{!15, !13, i64 4}
!89 = !{!15, !13, i64 8}
!90 = !{!15, !13, i64 12}
!91 = !{!15, !13, i64 16}
!92 = !{!15, !13, i64 20}
!93 = distinct !{!93, !94}
!94 = !{!"llvm.loop.mustprogress"}


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

%struct.generic_integrator_workspace = type { i32, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, double, [2048 x i8] }
%struct.thermodynamics_parameters_and_workspace = type { ptr, ptr, ptr, ptr }

@.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.212 = external hidden unnamed_addr constant [27 x i8], align 1
@__func__.thermodynamics_recombination_with_recfast = external hidden unnamed_addr constant [42 x i8], align 1
@.str.213 = external hidden unnamed_addr constant [38 x i8], align 1
@.str.214 = external hidden unnamed_addr constant [40 x i8], align 1
@.str.215 = external hidden unnamed_addr constant [70 x i8], align 1
@.str.216 = external hidden unnamed_addr constant [41 x i8], align 1
@.str.217 = external hidden unnamed_addr constant [59 x i8], align 1
@.str.218 = external hidden unnamed_addr constant [113 x i8], align 1
@.str.219 = external hidden unnamed_addr constant [31 x i8], align 1
@.str.220 = external hidden unnamed_addr constant [146 x i8], align 1
@.str.221 = external hidden unnamed_addr constant [74 x i8], align 1
@.str.222 = external hidden unnamed_addr constant [32 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

; Function Attrs: mustprogress nofree nounwind willreturn allockind("alloc,uninitialized") allocsize(0) memory(inaccessiblemem: readwrite)
declare noalias noundef ptr @malloc(i64 noundef) local_unnamed_addr #3

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

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

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

; Function Attrs: nounwind uwtable
define dso_local range(i32 0, 2) i32 @thermodynamics_recombination_with_recfast(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, ptr noundef %4) local_unnamed_addr #5 {
  %6 = alloca [3 x double], align 16
  %7 = alloca [3 x double], align 16
  %8 = alloca %struct.generic_integrator_workspace, align 8
  %9 = alloca %struct.thermodynamics_parameters_and_workspace, align 8
  %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
  call void @llvm.lifetime.start.p0(i64 24, ptr nonnull %6) #6
  call void @llvm.lifetime.start.p0(i64 24, ptr nonnull %7) #6
  call void @llvm.lifetime.start.p0(i64 2152, ptr nonnull %8) #6
  call void @llvm.lifetime.start.p0(i64 32, ptr nonnull %9) #6
  %22 = getelementptr inbounds i8, ptr %0, i64 368
  %23 = load i32, ptr %22, align 8, !tbaa !5
  %24 = getelementptr inbounds i8, ptr %3, i64 28
  store i32 %23, ptr %24, align 4, !tbaa !11
  %25 = getelementptr inbounds i8, ptr %3, i64 24
  %26 = load i32, ptr %25, align 8, !tbaa !15
  %27 = mul nsw i32 %26, %23
  %28 = sext i32 %27 to i64
  %29 = shl nsw i64 %28, 3
  %30 = tail call noalias ptr @malloc(i64 noundef %29) #7
  %31 = getelementptr inbounds i8, ptr %3, i64 32
  store ptr %30, ptr %31, align 8, !tbaa !16
  %32 = icmp eq ptr %30, null
  br i1 %32, label %33, label %36

33:                                               ; preds = %5
  %34 = shl i32 %27, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %10) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %10, ptr noundef nonnull @.str.37, ptr noundef nonnull @.str.212, i32 noundef %34) #6
  %35 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %35, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3478, ptr noundef nonnull %10) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %10) #6
  br label %534

36:                                               ; preds = %5
  %37 = call i32 @initialize_generic_integrator(i32 noundef 3, ptr noundef nonnull %8) #6
  %38 = icmp eq i32 %37, 1
  br i1 %38, label %39, label %42

39:                                               ; preds = %36
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %11) #6
  %40 = getelementptr inbounds i8, ptr %8, i64 104
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %11, ptr noundef nonnull @.str, ptr noundef nonnull @.str.213, ptr noundef nonnull %40) #6
  %41 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %41, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3483, ptr noundef nonnull %11) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %11) #6
  br label %534

42:                                               ; preds = %36
  %43 = load i32, ptr %22, align 8, !tbaa !5
  %44 = load double, ptr %1, align 8, !tbaa !17
  %45 = fmul double %44, 0x41B1DE784A000000
  %46 = fdiv double %45, 0x449A22FF260E3826
  %47 = getelementptr inbounds i8, ptr %3, i64 176
  store double %46, ptr %47, align 8, !tbaa !19
  %48 = getelementptr inbounds i8, ptr %1, i64 24
  %49 = load double, ptr %48, align 8, !tbaa !20
  %50 = load double, ptr %2, align 8, !tbaa !21
  %51 = getelementptr inbounds i8, ptr %3, i64 184
  store double %50, ptr %51, align 8, !tbaa !23
  %52 = getelementptr inbounds i8, ptr %1, i64 16
  %53 = load double, ptr %52, align 8, !tbaa !24
  %54 = getelementptr inbounds i8, ptr %3, i64 128
  store double %53, ptr %54, align 8, !tbaa !25
  %55 = getelementptr inbounds i8, ptr %0, i64 360
  %56 = load double, ptr %55, align 8, !tbaa !26
  %57 = getelementptr inbounds i8, ptr %0, i64 584
  %58 = load double, ptr %57, align 8, !tbaa !27
  %59 = getelementptr inbounds i8, ptr %3, i64 120
  store double %58, ptr %59, align 8, !tbaa !28
  %60 = getelementptr inbounds i8, ptr %0, i64 416
  %61 = load i32, ptr %60, align 8, !tbaa !29
  %switch = icmp ult i32 %61, 2
  br i1 %switch, label %64, label %62

62:                                               ; preds = %42
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %12) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %12, ptr noundef nonnull @.str.214) #6
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %13) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %13, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.215, ptr noundef nonnull %12) #6
  %63 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %63, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3511, ptr noundef nonnull %13) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %13) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %12) #6
  br label %534

64:                                               ; preds = %42
  %65 = getelementptr inbounds i8, ptr %0, i64 424
  %66 = load double, ptr %65, align 8, !tbaa !30
  %67 = getelementptr inbounds i8, ptr %3, i64 112
  store double %66, ptr %67, align 8, !tbaa !31
  %68 = icmp eq i32 %61, 1
  br i1 %68, label %69, label %73

69:                                               ; preds = %64
  %70 = getelementptr inbounds i8, ptr %0, i64 432
  %71 = load double, ptr %70, align 8, !tbaa !32
  %72 = fadd double %66, %71
  store double %72, ptr %67, align 8, !tbaa !31
  br label %73

73:                                               ; preds = %69, %64
  %74 = getelementptr inbounds i8, ptr %0, i64 400
  %75 = load i32, ptr %74, align 8, !tbaa !33
  %or.cond = icmp ugt i32 %75, 6
  br i1 %or.cond, label %76, label %78

76:                                               ; preds = %73
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %14) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %14, ptr noundef nonnull @.str.216) #6
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %15) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %15, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.217, ptr noundef nonnull %14) #6
  %77 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %77, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3519, ptr noundef nonnull %15) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %15) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %14) #6
  br label %534

78:                                               ; preds = %73
  %79 = fsub double 1.000000e+00, %50
  %80 = fdiv double 1.000000e+00, %79
  %81 = fmul double %79, 3.971500e+00
  %82 = fdiv double %50, %81
  %83 = getelementptr inbounds i8, ptr %3, i64 80
  store double %82, ptr %83, align 8, !tbaa !34
  %84 = fmul double %46, 3.000000e+00
  %85 = fmul double %46, %84
  %86 = fmul double %49, %85
  %87 = fmul double %80, 0x3E1CD169BDE2F4EC
  %88 = fmul double %87, 0x3A60930446E176B5
  %89 = fdiv double %86, %88
  %90 = getelementptr inbounds i8, ptr %3, i64 136
  store double %89, ptr %90, align 8, !tbaa !35
  %91 = getelementptr inbounds i8, ptr %2, i64 592
  store double %89, ptr %91, align 8, !tbaa !36
  %92 = getelementptr inbounds i8, ptr %2, i64 168
  %93 = load double, ptr %92, align 8, !tbaa !37
  %94 = getelementptr inbounds i8, ptr %3, i64 192
  store double %93, ptr %94, align 8, !tbaa !38
  %95 = getelementptr inbounds i8, ptr %2, i64 176
  %96 = load i16, ptr %95, align 8, !tbaa !39
  %97 = getelementptr inbounds i8, ptr %3, i64 200
  store i16 %96, ptr %97, align 8, !tbaa !40
  %98 = getelementptr inbounds i8, ptr %2, i64 200
  %99 = getelementptr inbounds i8, ptr %3, i64 224
  %100 = load <2 x double>, ptr %98, align 8, !tbaa !41
  store <2 x double> %100, ptr %99, align 8, !tbaa !41
  %101 = getelementptr inbounds i8, ptr %2, i64 216
  %102 = load double, ptr %101, align 8, !tbaa !42
  %103 = getelementptr inbounds i8, ptr %3, i64 240
  store double %102, ptr %103, align 8, !tbaa !43
  %104 = getelementptr inbounds i8, ptr %2, i64 184
  %105 = getelementptr inbounds i8, ptr %3, i64 208
  %106 = load <2 x double>, ptr %104, align 8, !tbaa !41
  store <2 x double> %106, ptr %105, align 8, !tbaa !41
  %107 = getelementptr inbounds i8, ptr %2, i64 224
  %108 = getelementptr inbounds i8, ptr %3, i64 248
  %109 = load <2 x double>, ptr %107, align 8, !tbaa !41
  store <2 x double> %109, ptr %108, align 8, !tbaa !41
  %110 = getelementptr inbounds i8, ptr %3, i64 40
  %111 = getelementptr inbounds i8, ptr %3, i64 88
  %112 = getelementptr inbounds i8, ptr %3, i64 152
  %113 = getelementptr inbounds i8, ptr %3, i64 160
  %114 = getelementptr inbounds i8, ptr %3, i64 168
  store <2 x double> <double 0x41116A32737408D4, double 0x412345866807EAE4>, ptr %113, align 8, !tbaa !41
  %115 = getelementptr inbounds i8, ptr %3, i64 48
  store <2 x double> <double 0x40E34354E1E59D09, double 0x42E4764E04F59C55>, ptr %110, align 8, !tbaa !41
  %116 = getelementptr inbounds i8, ptr %3, i64 56
  store <2 x double> <double 0x40E68118AC20A6BD, double 0x3B2331C2ECCE4D3B>, ptr %111, align 8, !tbaa !41
  store <2 x double> <double 0x3B559AC95028F4FC, double 0x40FCE50714FABBA2>, ptr %116, align 8, !tbaa !41
  %117 = getelementptr inbounds i8, ptr %3, i64 104
  store double 0x410D341EBBDFE7F9, ptr %117, align 8, !tbaa !44
  %118 = getelementptr inbounds i8, ptr %3, i64 72
  store double 0x3B82912DC103EAAB, ptr %118, align 8, !tbaa !45
  %119 = getelementptr inbounds i8, ptr %3, i64 144
  store <2 x double> <double 0x40BB4CC9D960EA35, double 0x41034358C2F6C513>, ptr %119, align 8, !tbaa !41
  store ptr %1, ptr %9, align 8, !tbaa !46
  %120 = getelementptr inbounds i8, ptr %9, i64 8
  store ptr %0, ptr %120, align 8, !tbaa !48
  %121 = getelementptr inbounds i8, ptr %9, i64 16
  store ptr %3, ptr %121, align 8, !tbaa !49
  %122 = getelementptr inbounds i8, ptr %9, i64 24
  store ptr %4, ptr %122, align 8, !tbaa !50
  %123 = getelementptr inbounds i8, ptr %0, i64 520
  %124 = load double, ptr %123, align 8, !tbaa !51
  %125 = fcmp olt double %56, %124
  br i1 %125, label %126, label %128

126:                                              ; preds = %78
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %16) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %16, ptr noundef nonnull @.str.218) #6
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %17) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %17, ptr noundef nonnull @.str.10, ptr noundef nonnull @.str.219, ptr noundef nonnull %16) #6
  %127 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %127, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3571, ptr noundef nonnull %17) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %17) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %16) #6
  br label %534

128:                                              ; preds = %78
  %129 = getelementptr inbounds i8, ptr %6, i64 8
  store <2 x double> <double 1.000000e+00, double 1.000000e+00>, ptr %6, align 16, !tbaa !41
  %130 = fadd double %56, 1.000000e+00
  %131 = fmul double %53, %130
  %132 = getelementptr inbounds i8, ptr %6, i64 16
  store double %131, ptr %132, align 16, !tbaa !41
  %133 = icmp sgt i32 %43, 0
  br i1 %133, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %128
  %134 = uitofp nneg i32 %43 to double
  %135 = getelementptr inbounds i8, ptr %0, i64 488
  %136 = getelementptr inbounds i8, ptr %0, i64 496
  %137 = getelementptr inbounds i8, ptr %0, i64 504
  %138 = getelementptr inbounds i8, ptr %0, i64 512
  %139 = getelementptr inbounds i8, ptr %0, i64 528
  %140 = getelementptr inbounds i8, ptr %0, i64 536
  %141 = getelementptr inbounds i8, ptr %0, i64 560
  %142 = getelementptr inbounds i8, ptr %0, i64 576
  %143 = getelementptr inbounds i8, ptr %0, i64 392
  %144 = getelementptr inbounds i8, ptr %0, i64 2344
  %145 = getelementptr inbounds i8, ptr %0, i64 552
  %146 = getelementptr inbounds i8, ptr %3, i64 4
  %147 = getelementptr inbounds i8, ptr %3, i64 8
  %148 = getelementptr inbounds i8, ptr %2, i64 606
  %149 = getelementptr inbounds i8, ptr %3, i64 12
  %150 = getelementptr inbounds i8, ptr %7, i64 16
  %151 = getelementptr inbounds i8, ptr %3, i64 16
  %152 = getelementptr inbounds i8, ptr %3, i64 20
  br label %153

153:                                              ; preds = %.lr.ph, %492
  %.0339362 = phi i32 [ 0, %.lr.ph ], [ %528, %492 ]
  %.0341361 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2343, %492 ]
  %154 = sub nsw i32 %43, %.0339362
  %155 = sitofp i32 %154 to double
  %156 = fmul double %56, %155
  %157 = fdiv double %156, %134
  %158 = add nsw i32 %154, -1
  %159 = sitofp i32 %158 to double
  %160 = fmul double %56, %159
  %161 = fdiv double %160, %134
  %162 = load double, ptr %135, align 8, !tbaa !52
  %163 = load double, ptr %136, align 8, !tbaa !53
  %164 = fadd double %162, %163
  %165 = fcmp ogt double %161, %164
  br i1 %165, label %166, label %172

166:                                              ; preds = %153
  %167 = load double, ptr %83, align 8, !tbaa !34
  %168 = call double @llvm.fmuladd.f64(double %167, double 2.000000e+00, double 1.000000e+00)
  store <2 x double> <double 1.000000e+00, double 1.000000e+00>, ptr %6, align 16, !tbaa !41
  %169 = load double, ptr %54, align 8, !tbaa !25
  %170 = fadd double %161, 1.000000e+00
  %171 = fmul double %170, %169
  store double %171, ptr %132, align 16, !tbaa !41
  br label %473

172:                                              ; preds = %153
  %173 = load double, ptr %137, align 8, !tbaa !54
  %174 = load double, ptr %138, align 8, !tbaa !55
  %175 = fadd double %173, %174
  %176 = fcmp ogt double %161, %175
  br i1 %176, label %177, label %236

177:                                              ; preds = %172
  %178 = load double, ptr %115, align 8, !tbaa !56
  %179 = load double, ptr %54, align 8, !tbaa !25
  %180 = fmul double %178, %179
  %181 = fadd double %161, 1.000000e+00
  %182 = fdiv double %180, %181
  %183 = call double @log(double noundef %182) #6
  %184 = load double, ptr %114, align 8, !tbaa !57
  %185 = load double, ptr %54, align 8, !tbaa !25
  %186 = fmul double %181, %185
  %187 = fneg double %184
  %188 = fdiv double %187, %186
  %189 = call double @llvm.fmuladd.f64(double %183, double 1.500000e+00, double %188)
  %190 = call double @exp(double noundef %189) #6
  %191 = load double, ptr %90, align 8, !tbaa !35
  %192 = fdiv double %190, %191
  %193 = load double, ptr %135, align 8, !tbaa !52
  %194 = load double, ptr %136, align 8, !tbaa !53
  %195 = fsub double %193, %194
  %196 = fcmp ogt double %161, %195
  br i1 %196, label %197, label %221

197:                                              ; preds = %177
  %198 = load double, ptr %83, align 8, !tbaa !34
  %199 = call double @llvm.fmuladd.f64(double %198, double 2.000000e+00, double 1.000000e+00)
  %200 = fadd double %192, -1.000000e+00
  %201 = fsub double %200, %198
  %square360 = fmul double %201, %201
  %202 = fmul double %199, 4.000000e+00
  %203 = call double @llvm.fmuladd.f64(double %202, double %192, double %square360)
  %204 = call double @sqrt(double noundef %203) #6
  %205 = load double, ptr %83, align 8, !tbaa !34
  %206 = fsub double %200, %205
  %207 = fsub double %204, %206
  %208 = fmul double %207, 5.000000e-01
  %209 = load double, ptr %135, align 8, !tbaa !52
  %210 = fsub double %209, %161
  %211 = load double, ptr %136, align 8, !tbaa !53
  %212 = fdiv double %210, %211
  %213 = fmul double %212, -7.500000e-01
  %214 = fmul double %212, %212
  %215 = fdiv double %214, 3.000000e+00
  %216 = fadd double %215, -1.000000e+00
  %217 = call double @llvm.fmuladd.f64(double %213, double %216, double 5.000000e-01)
  %218 = fsub double 1.000000e+00, %217
  %219 = fmul double %199, %218
  %220 = call double @llvm.fmuladd.f64(double %217, double %208, double %219)
  br label %233

221:                                              ; preds = %177
  %222 = fadd double %192, -1.000000e+00
  %223 = load double, ptr %83, align 8, !tbaa !34
  %224 = fsub double %222, %223
  %square359 = fmul double %224, %224
  %225 = call double @llvm.fmuladd.f64(double %223, double 2.000000e+00, double 1.000000e+00)
  %226 = fmul double %225, 4.000000e+00
  %227 = call double @llvm.fmuladd.f64(double %226, double %192, double %square359)
  %228 = call double @sqrt(double noundef %227) #6
  %229 = load double, ptr %83, align 8, !tbaa !34
  %230 = fsub double %222, %229
  %231 = fsub double %228, %230
  %232 = fmul double %231, 5.000000e-01
  br label %233

233:                                              ; preds = %221, %197
  %.0340 = phi double [ %220, %197 ], [ %232, %221 ]
  store <2 x double> <double 1.000000e+00, double 1.000000e+00>, ptr %6, align 16, !tbaa !41
  %234 = load double, ptr %54, align 8, !tbaa !25
  %235 = fmul double %181, %234
  store double %235, ptr %132, align 16, !tbaa !41
  br label %473

236:                                              ; preds = %172
  %237 = load double, ptr %123, align 8, !tbaa !51
  %238 = load double, ptr %139, align 8, !tbaa !58
  %239 = fadd double %237, %238
  %240 = fcmp ogt double %161, %239
  br i1 %240, label %241, label %290

241:                                              ; preds = %236
  %242 = fsub double %173, %174
  %243 = fcmp ogt double %161, %242
  br i1 %243, label %244, label %284

244:                                              ; preds = %241
  %245 = load double, ptr %115, align 8, !tbaa !56
  %246 = load double, ptr %54, align 8, !tbaa !25
  %247 = fmul double %245, %246
  %248 = fadd double %161, 1.000000e+00
  %249 = fdiv double %247, %248
  %250 = call double @log(double noundef %249) #6
  %251 = load double, ptr %114, align 8, !tbaa !57
  %252 = load double, ptr %54, align 8, !tbaa !25
  %253 = fmul double %248, %252
  %254 = fneg double %251
  %255 = fdiv double %254, %253
  %256 = call double @llvm.fmuladd.f64(double %250, double 1.500000e+00, double %255)
  %257 = call double @exp(double noundef %256) #6
  %258 = load double, ptr %90, align 8, !tbaa !35
  %259 = fdiv double %257, %258
  %260 = fadd double %259, -1.000000e+00
  %261 = load double, ptr %83, align 8, !tbaa !34
  %262 = fsub double %260, %261
  %square358 = fmul double %262, %262
  %263 = call double @llvm.fmuladd.f64(double %261, double 2.000000e+00, double 1.000000e+00)
  %264 = fmul double %263, 4.000000e+00
  %265 = call double @llvm.fmuladd.f64(double %264, double %259, double %square358)
  %266 = call double @sqrt(double noundef %265) #6
  %267 = load double, ptr %83, align 8, !tbaa !34
  %268 = fsub double %260, %267
  %269 = fsub double %266, %268
  %270 = fmul double %269, 5.000000e-01
  %271 = fadd double %267, 1.000000e+00
  %272 = load double, ptr %137, align 8, !tbaa !54
  %273 = fsub double %272, %161
  %274 = load double, ptr %138, align 8, !tbaa !55
  %275 = fdiv double %273, %274
  %276 = fmul double %275, -7.500000e-01
  %277 = fmul double %275, %275
  %278 = fdiv double %277, 3.000000e+00
  %279 = fadd double %278, -1.000000e+00
  %280 = call double @llvm.fmuladd.f64(double %276, double %279, double 5.000000e-01)
  %281 = fsub double 1.000000e+00, %280
  %282 = fmul double %270, %281
  %283 = call double @llvm.fmuladd.f64(double %280, double %271, double %282)
  br label %287

284:                                              ; preds = %241
  %285 = load double, ptr %83, align 8, !tbaa !34
  %286 = fadd double %285, 1.000000e+00
  %.pre363 = fadd double %161, 1.000000e+00
  br label %287

287:                                              ; preds = %284, %244
  %.pre-phi = phi double [ %.pre363, %284 ], [ %248, %244 ]
  %.1 = phi double [ %286, %284 ], [ %283, %244 ]
  store <2 x double> <double 1.000000e+00, double 1.000000e+00>, ptr %6, align 16, !tbaa !41
  %288 = load double, ptr %54, align 8, !tbaa !25
  %289 = fmul double %.pre-phi, %288
  store double %289, ptr %132, align 16, !tbaa !41
  br label %473

290:                                              ; preds = %236
  %291 = load double, ptr %129, align 8, !tbaa !41
  %292 = load double, ptr %140, align 8, !tbaa !59
  %293 = fcmp ogt double %291, %292
  br i1 %293, label %294, label %340

294:                                              ; preds = %290
  %295 = load double, ptr %115, align 8, !tbaa !56
  %296 = load double, ptr %54, align 8, !tbaa !25
  %297 = fmul double %295, %296
  %298 = fadd double %161, 1.000000e+00
  %299 = fdiv double %297, %298
  %300 = call double @log(double noundef %299) #6
  %301 = load double, ptr %113, align 8, !tbaa !60
  %302 = load double, ptr %54, align 8, !tbaa !25
  %303 = fmul double %298, %302
  %304 = fneg double %301
  %305 = fdiv double %304, %303
  %306 = call double @llvm.fmuladd.f64(double %300, double 1.500000e+00, double %305)
  %307 = call double @exp(double noundef %306) #6
  %308 = fmul double %307, 4.000000e+00
  %309 = load double, ptr %90, align 8, !tbaa !35
  %310 = fdiv double %308, %309
  %311 = fadd double %310, -1.000000e+00
  %square357 = fmul double %311, %311
  %312 = load double, ptr %83, align 8, !tbaa !34
  %313 = fadd double %312, 1.000000e+00
  %314 = fmul double %313, 4.000000e+00
  %315 = call double @llvm.fmuladd.f64(double %314, double %310, double %square357)
  %316 = call double @sqrt(double noundef %315) #6
  %317 = fsub double %316, %311
  %318 = fmul double %317, 5.000000e-01
  %319 = load double, ptr %123, align 8, !tbaa !51
  %320 = load double, ptr %139, align 8, !tbaa !58
  %321 = fsub double %319, %320
  %322 = fcmp ogt double %161, %321
  %.pre = load double, ptr %83, align 8, !tbaa !34
  br i1 %322, label %323, label %335

323:                                              ; preds = %294
  %324 = fadd double %.pre, 1.000000e+00
  %325 = fsub double %319, %161
  %326 = fdiv double %325, %320
  %327 = fmul double %326, -7.500000e-01
  %328 = fmul double %326, %326
  %329 = fdiv double %328, 3.000000e+00
  %330 = fadd double %329, -1.000000e+00
  %331 = call double @llvm.fmuladd.f64(double %327, double %330, double 5.000000e-01)
  %332 = fsub double 1.000000e+00, %331
  %333 = fmul double %332, %324
  %334 = call double @llvm.fmuladd.f64(double %331, double %318, double %333)
  br label %335

335:                                              ; preds = %294, %323
  %.2 = phi double [ %334, %323 ], [ %318, %294 ]
  %336 = fadd double %.2, -1.000000e+00
  %337 = fdiv double %336, %.pre
  store double 1.000000e+00, ptr %6, align 16, !tbaa !41
  store double %337, ptr %129, align 8, !tbaa !41
  %338 = load double, ptr %54, align 8, !tbaa !25
  %339 = fmul double %298, %338
  store double %339, ptr %132, align 16, !tbaa !41
  br label %473

340:                                              ; preds = %290
  %341 = load double, ptr %6, align 16, !tbaa !41
  %342 = load double, ptr %141, align 8, !tbaa !61
  %343 = fcmp ogt double %341, %342
  br i1 %343, label %344, label %419

344:                                              ; preds = %340
  %345 = load double, ptr %115, align 8, !tbaa !56
  %346 = load double, ptr %54, align 8, !tbaa !25
  %347 = fmul double %345, %346
  %348 = fadd double %161, 1.000000e+00
  %349 = fdiv double %347, %348
  %350 = call double @log(double noundef %349) #6
  %351 = load double, ptr %112, align 8, !tbaa !62
  %352 = load double, ptr %54, align 8, !tbaa !25
  %353 = fmul double %348, %352
  %354 = fneg double %351
  %355 = fdiv double %354, %353
  %356 = call double @llvm.fmuladd.f64(double %350, double 1.500000e+00, double %355)
  %357 = call double @exp(double noundef %356) #6
  %358 = load double, ptr %90, align 8, !tbaa !35
  %359 = fdiv double %357, %358
  %square355 = fmul double %359, %359
  %360 = call double @llvm.fmuladd.f64(double %359, double 4.000000e+00, double %square355)
  %361 = call double @sqrt(double noundef %360) #6
  %362 = fsub double %361, %359
  %363 = fmul double %362, 5.000000e-01
  %364 = load double, ptr %143, align 8, !tbaa !63
  %365 = load double, ptr %144, align 8, !tbaa !64
  %366 = call i32 @generic_integrator(ptr noundef nonnull @thermodynamics_derivs_with_recfast, double noundef %157, double noundef %161, ptr noundef nonnull %6, ptr noundef nonnull %9, double noundef %364, double noundef %365, ptr noundef nonnull %8) #6
  %367 = icmp eq i32 %366, 1
  br i1 %367, label %368, label %370

368:                                              ; preds = %344
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %18) #6
  %369 = getelementptr inbounds i8, ptr %8, i64 104
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %18, ptr noundef nonnull @.str, ptr noundef nonnull @.str.220, ptr noundef nonnull %369) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %148, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3707, ptr noundef nonnull %18) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %18) #6
  br label %534

370:                                              ; preds = %344
  store double %363, ptr %6, align 16, !tbaa !41
  %371 = load double, ptr %140, align 8, !tbaa !59
  %372 = load double, ptr %129, align 8, !tbaa !41
  %373 = fsub double %371, %372
  %374 = load double, ptr %145, align 8, !tbaa !65
  %375 = fcmp olt double %373, %374
  br i1 %375, label %376, label %416

376:                                              ; preds = %370
  %377 = load double, ptr %115, align 8, !tbaa !56
  %378 = load double, ptr %54, align 8, !tbaa !25
  %379 = fmul double %377, %378
  %380 = fdiv double %379, %348
  %381 = call double @log(double noundef %380) #6
  %382 = load double, ptr %113, align 8, !tbaa !60
  %383 = load double, ptr %54, align 8, !tbaa !25
  %384 = fmul double %348, %383
  %385 = fneg double %382
  %386 = fdiv double %385, %384
  %387 = call double @llvm.fmuladd.f64(double %381, double 1.500000e+00, double %386)
  %388 = call double @exp(double noundef %387) #6
  %389 = fmul double %388, 4.000000e+00
  %390 = load double, ptr %90, align 8, !tbaa !35
  %391 = fdiv double %389, %390
  %392 = fadd double %391, -1.000000e+00
  %square356 = fmul double %392, %392
  %393 = load double, ptr %83, align 8, !tbaa !34
  %394 = fadd double %393, 1.000000e+00
  %395 = fmul double %394, 4.000000e+00
  %396 = call double @llvm.fmuladd.f64(double %395, double %391, double %square356)
  %397 = call double @sqrt(double noundef %396) #6
  %398 = fsub double %397, %392
  %399 = fmul double %398, 5.000000e-01
  %400 = load double, ptr %6, align 16, !tbaa !41
  %401 = load double, ptr %83, align 8, !tbaa !34
  %402 = load double, ptr %129, align 8, !tbaa !41
  %403 = call double @llvm.fmuladd.f64(double %401, double %402, double %400)
  %404 = load double, ptr %140, align 8, !tbaa !59
  %405 = fsub double %404, %402
  %406 = load double, ptr %145, align 8, !tbaa !65
  %407 = fdiv double %405, %406
  %408 = fmul double %407, %407
  %409 = fdiv double %407, 3.000000e+00
  %410 = fsub double 5.000000e-01, %409
  %411 = fmul double %408, %410
  %412 = fmul double %411, 6.000000e+00
  %413 = fsub double 1.000000e+00, %412
  %414 = fmul double %399, %413
  %415 = call double @llvm.fmuladd.f64(double %412, double %403, double %414)
  br label %473

416:                                              ; preds = %370
  %417 = load double, ptr %83, align 8, !tbaa !34
  %418 = call double @llvm.fmuladd.f64(double %417, double %372, double %363)
  br label %473

419:                                              ; preds = %340
  %420 = fsub double %342, %341
  %421 = load double, ptr %142, align 8, !tbaa !66
  %422 = fcmp olt double %420, %421
  br i1 %422, label %423, label %443

423:                                              ; preds = %419
  %424 = load double, ptr %115, align 8, !tbaa !56
  %425 = load double, ptr %54, align 8, !tbaa !25
  %426 = fmul double %424, %425
  %427 = fadd double %161, 1.000000e+00
  %428 = fdiv double %426, %427
  %429 = call double @log(double noundef %428) #6
  %430 = load double, ptr %112, align 8, !tbaa !62
  %431 = load double, ptr %54, align 8, !tbaa !25
  %432 = fmul double %427, %431
  %433 = fneg double %430
  %434 = fdiv double %433, %432
  %435 = call double @llvm.fmuladd.f64(double %429, double 1.500000e+00, double %434)
  %436 = call double @exp(double noundef %435) #6
  %437 = load double, ptr %90, align 8, !tbaa !35
  %438 = fdiv double %436, %437
  %square = fmul double %438, %438
  %439 = call double @llvm.fmuladd.f64(double %438, double 4.000000e+00, double %square)
  %440 = call double @sqrt(double noundef %439) #6
  %441 = fsub double %440, %438
  %442 = fmul double %441, 5.000000e-01
  br label %443

443:                                              ; preds = %423, %419
  %.1342 = phi double [ %442, %423 ], [ %.0341361, %419 ]
  %444 = load double, ptr %143, align 8, !tbaa !63
  %445 = load double, ptr %144, align 8, !tbaa !64
  %446 = call i32 @generic_integrator(ptr noundef nonnull @thermodynamics_derivs_with_recfast, double noundef %157, double noundef %161, ptr noundef nonnull %6, ptr noundef nonnull %9, double noundef %444, double noundef %445, ptr noundef nonnull %8) #6
  %447 = icmp eq i32 %446, 1
  br i1 %447, label %448, label %450

448:                                              ; preds = %443
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %19) #6
  %449 = getelementptr inbounds i8, ptr %8, i64 104
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %19, ptr noundef nonnull @.str, ptr noundef nonnull @.str.220, ptr noundef nonnull %449) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %148, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3749, ptr noundef nonnull %19) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %19) #6
  br label %534

450:                                              ; preds = %443
  %451 = load double, ptr %141, align 8, !tbaa !61
  %452 = load double, ptr %6, align 16, !tbaa !41
  %453 = fsub double %451, %452
  %454 = load double, ptr %142, align 8, !tbaa !66
  %455 = fcmp olt double %453, %454
  br i1 %455, label %456, label %469

456:                                              ; preds = %450
  %457 = fdiv double %453, %454
  %458 = fmul double %457, %457
  %459 = fdiv double %457, 3.000000e+00
  %460 = fsub double 5.000000e-01, %459
  %461 = fmul double %458, %460
  %462 = fmul double %461, 6.000000e+00
  %463 = fsub double 1.000000e+00, %462
  %464 = fmul double %.1342, %463
  %465 = call double @llvm.fmuladd.f64(double %462, double %452, double %464)
  %466 = load double, ptr %83, align 8, !tbaa !34
  %467 = load double, ptr %129, align 8, !tbaa !41
  %468 = call double @llvm.fmuladd.f64(double %466, double %467, double %465)
  br label %473

469:                                              ; preds = %450
  %470 = load double, ptr %83, align 8, !tbaa !34
  %471 = load double, ptr %129, align 8, !tbaa !41
  %472 = call double @llvm.fmuladd.f64(double %470, double %471, double %452)
  br label %473

473:                                              ; preds = %233, %335, %456, %469, %376, %416, %287, %166
  %.2343 = phi double [ 1.000000e+00, %166 ], [ 1.000000e+00, %233 ], [ 1.000000e+00, %287 ], [ 1.000000e+00, %335 ], [ %363, %376 ], [ %363, %416 ], [ %.1342, %456 ], [ %.1342, %469 ]
  %.3 = phi double [ %168, %166 ], [ %.0340, %233 ], [ %.1, %287 ], [ %.2, %335 ], [ %415, %376 ], [ %418, %416 ], [ %468, %456 ], [ %472, %469 ]
  %474 = load ptr, ptr %31, align 8, !tbaa !16
  %475 = load i32, ptr %25, align 8, !tbaa !15
  %476 = mul nsw i32 %475, %158
  %477 = sext i32 %476 to i64
  %478 = getelementptr inbounds double, ptr %474, i64 %477
  %479 = load i32, ptr %3, align 8, !tbaa !67
  %480 = sext i32 %479 to i64
  %481 = getelementptr inbounds double, ptr %478, i64 %480
  store double %161, ptr %481, align 8, !tbaa !41
  %482 = load i32, ptr %146, align 4, !tbaa !68
  %483 = sext i32 %482 to i64
  %484 = getelementptr inbounds double, ptr %478, i64 %483
  store double %.3, ptr %484, align 8, !tbaa !41
  %485 = load double, ptr %132, align 16, !tbaa !41
  %486 = load i32, ptr %147, align 8, !tbaa !69
  %487 = sext i32 %486 to i64
  %488 = getelementptr inbounds double, ptr %478, i64 %487
  store double %485, ptr %488, align 8, !tbaa !41
  %489 = call i32 @thermodynamics_derivs_with_recfast(double noundef %161, ptr noundef nonnull %6, ptr noundef nonnull %7, ptr noundef nonnull %9, ptr noundef nonnull %148)
  %490 = icmp eq i32 %489, 1
  br i1 %490, label %491, label %492

491:                                              ; preds = %473
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %20) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %20, ptr noundef nonnull @.str, ptr noundef nonnull @.str.221, ptr noundef nonnull %148) #6
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %148, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3782, ptr noundef nonnull %20) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %20) #6
  br label %534

492:                                              ; preds = %473
  %493 = load double, ptr %51, align 8, !tbaa !23
  %494 = call double @llvm.fmuladd.f64(double %493, double 0xBFE7F14DA42541E8, double 1.000000e+00)
  %495 = fsub double 1.000000e+00, %493
  %496 = call double @llvm.fmuladd.f64(double %.3, double %495, double %494)
  %497 = fmul double %496, 0x3D39D63187BB40C4
  %498 = load double, ptr %132, align 16, !tbaa !41
  %499 = fmul double %498, %497
  %500 = load ptr, ptr %31, align 8, !tbaa !16
  %501 = load i32, ptr %25, align 8, !tbaa !15
  %502 = mul nsw i32 %501, %158
  %503 = sext i32 %502 to i64
  %504 = getelementptr inbounds double, ptr %500, i64 %503
  %505 = load i32, ptr %149, align 4, !tbaa !70
  %506 = sext i32 %505 to i64
  %507 = getelementptr inbounds double, ptr %504, i64 %506
  store double %499, ptr %507, align 8, !tbaa !41
  %508 = fadd double %161, 1.000000e+00
  %509 = load double, ptr %150, align 16, !tbaa !41
  %510 = fmul double %508, %509
  %511 = load double, ptr %132, align 16, !tbaa !41
  %512 = fdiv double %510, %511
  %513 = fdiv double %512, 3.000000e+00
  %514 = fadd double %513, 1.000000e+00
  %515 = fmul double %499, %514
  %516 = load i32, ptr %151, align 8, !tbaa !71
  %517 = sext i32 %516 to i64
  %518 = getelementptr inbounds double, ptr %504, i64 %517
  store double %515, ptr %518, align 8, !tbaa !41
  %519 = fmul double %508, %508
  %520 = load double, ptr %90, align 8, !tbaa !35
  %521 = fmul double %519, %520
  %522 = fmul double %.3, %521
  %523 = fmul double %522, 0x3A15151E37FBF236
  %524 = fmul double %523, 0x449A22FF260E3826
  %525 = load i32, ptr %152, align 4, !tbaa !72
  %526 = sext i32 %525 to i64
  %527 = getelementptr inbounds double, ptr %504, i64 %526
  store double %524, ptr %527, align 8, !tbaa !41
  %528 = add nuw nsw i32 %.0339362, 1
  %exitcond.not = icmp eq i32 %528, %43
  br i1 %exitcond.not, label %._crit_edge, label %153, !llvm.loop !73

._crit_edge:                                      ; preds = %492, %128
  %529 = call i32 @cleanup_generic_integrator(ptr noundef nonnull %8) #6
  %530 = icmp eq i32 %529, 1
  br i1 %530, label %531, label %534

531:                                              ; preds = %._crit_edge
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %21) #6
  %532 = getelementptr inbounds i8, ptr %8, i64 104
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %21, ptr noundef nonnull @.str, ptr noundef nonnull @.str.222, ptr noundef nonnull %532) #6
  %533 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %533, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_recombination_with_recfast, i32 noundef 3812, ptr noundef nonnull %21) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %21) #6
  br label %534

534:                                              ; preds = %._crit_edge, %531, %491, %448, %368, %126, %76, %62, %39, %33
  %.0 = phi i32 [ 1, %33 ], [ 1, %39 ], [ 1, %62 ], [ 1, %76 ], [ 1, %126 ], [ 1, %491 ], [ 1, %368 ], [ 1, %448 ], [ 1, %531 ], [ 0, %._crit_edge ]
  call void @llvm.lifetime.end.p0(i64 32, ptr nonnull %9) #6
  call void @llvm.lifetime.end.p0(i64 2152, ptr nonnull %8) #6
  call void @llvm.lifetime.end.p0(i64 24, ptr nonnull %7) #6
  call void @llvm.lifetime.end.p0(i64 24, ptr nonnull %6) #6
  ret i32 %.0
}

declare i32 @initialize_generic_integrator(i32 noundef, ptr noundef) local_unnamed_addr #2

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

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

declare i32 @cleanup_generic_integrator(ptr noundef) local_unnamed_addr #2

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 = { 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 #4 = { mustprogress nofree nounwind willreturn memory(write) "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #5 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #6 = { 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, !10, i64 368}
!6 = !{!"precision", !7, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !7, i64 40, !7, i64 48, !7, i64 56, !7, i64 64, !7, i64 72, !7, i64 80, !7, i64 88, !7, i64 96, !8, i64 104, !7, i64 360, !10, i64 368, !7, i64 376, !10, i64 384, !10, i64 388, !7, i64 392, !10, i64 400, !7, i64 408, !10, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !7, i64 464, !7, i64 472, !7, i64 480, !7, i64 488, !7, i64 496, !7, i64 504, !7, i64 512, !7, i64 520, !7, i64 528, !7, i64 536, !7, i64 544, !7, i64 552, !7, i64 560, !7, i64 568, !7, i64 576, !7, i64 584, !7, i64 592, !7, i64 600, !7, i64 608, !7, i64 616, !10, i64 624, !8, i64 628, !8, i64 884, !8, i64 1140, !7, i64 1400, !7, i64 1408, !7, i64 1416, !7, i64 1424, !7, i64 1432, !7, i64 1440, !7, i64 1448, !7, i64 1456, !7, i64 1464, !7, i64 1472, !7, i64 1480, !7, i64 1488, !7, i64 1496, !7, i64 1504, !7, i64 1512, !7, i64 1520, !10, i64 1528, !7, i64 1536, !7, i64 1544, !10, i64 1552, !10, i64 1556, !10, i64 1560, !10, i64 1564, !10, i64 1568, !10, i64 1572, !10, i64 1576, !10, i64 1580, !7, i64 1584, !7, i64 1592, !7, i64 1600, !7, i64 1608, !7, i64 1616, !7, i64 1624, !7, i64 1632, !7, i64 1640, !10, i64 1648, !7, i64 1656, !7, i64 1664, !10, i64 1672, !7, i64 1680, !7, i64 1688, !10, i64 1696, !7, i64 1704, !10, i64 1712, !7, i64 1720, !7, i64 1728, !10, i64 1736, !7, i64 1744, !7, i64 1752, !7, i64 1760, !10, i64 1768, !7, i64 1776, !7, i64 1784, !7, i64 1792, !7, i64 1800, !7, i64 1808, !10, i64 1816, !7, i64 1824, !7, i64 1832, !7, i64 1840, !7, i64 1848, !7, i64 1856, !7, i64 1864, !7, i64 1872, !10, i64 1880, !7, i64 1888, !7, i64 1896, !7, i64 1904, !7, i64 1912, !7, i64 1920, !7, i64 1928, !7, i64 1936, !7, i64 1944, !7, i64 1952, !7, i64 1960, !7, i64 1968, !7, i64 1976, !7, i64 1984, !7, i64 1992, !7, i64 2000, !7, i64 2008, !7, i64 2016, !7, i64 2024, !7, i64 2032, !7, i64 2040, !7, i64 2048, !7, i64 2056, !7, i64 2064, !7, i64 2072, !7, i64 2080, !7, i64 2088, !7, i64 2096, !7, i64 2104, !7, i64 2112, !7, i64 2120, !7, i64 2128, !7, i64 2136, !7, i64 2144, !7, i64 2152, !7, i64 2160, !7, i64 2168, !7, i64 2176, !7, i64 2184, !7, i64 2192, !7, i64 2200, !7, i64 2208, !7, i64 2216, !7, i64 2224, !7, i64 2232, !7, i64 2240, !7, i64 2248, !10, i64 2256, !7, i64 2264, !7, i64 2272, !7, i64 2280, !7, i64 2288, !10, i64 2296, !7, i64 2304, !7, i64 2312, !10, i64 2320, !10, i64 2324, !10, i64 2328, !7, i64 2336, !7, i64 2344, !8, i64 2352}
!7 = !{!"double", !8, i64 0}
!8 = !{!"omnipotent char", !9, i64 0}
!9 = !{!"Simple C/C++ TBAA"}
!10 = !{!"int", !8, i64 0}
!11 = !{!12, !10, i64 28}
!12 = !{!"recombination", !10, i64 0, !10, i64 4, !10, i64 8, !10, i64 12, !10, i64 16, !10, i64 20, !10, i64 24, !10, i64 28, !13, i64 32, !7, i64 40, !7, i64 48, !7, i64 56, !7, i64 64, !7, i64 72, !7, i64 80, !7, i64 88, !7, i64 96, !7, i64 104, !7, i64 112, !7, i64 120, !7, i64 128, !7, i64 136, !7, i64 144, !7, i64 152, !7, i64 160, !7, i64 168, !7, i64 176, !7, i64 184, !7, i64 192, !14, i64 200, !7, i64 208, !7, i64 216, !7, i64 224, !7, i64 232, !7, i64 240, !7, i64 248, !7, i64 256}
!13 = !{!"any pointer", !8, i64 0}
!14 = !{!"short", !8, i64 0}
!15 = !{!12, !10, i64 24}
!16 = !{!12, !13, i64 32}
!17 = !{!18, !7, i64 0}
!18 = !{!"background", !7, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !7, i64 40, !7, i64 48, !10, i64 56, !7, i64 64, !7, i64 72, !7, i64 80, !7, i64 88, !14, i64 96, !7, i64 104, !7, i64 112, !7, i64 120, !7, i64 128, !7, i64 136, !7, i64 144, !7, i64 152, !7, i64 160, !7, i64 168, !14, i64 176, !7, i64 184, !7, i64 192, !13, i64 200, !10, i64 208, !10, i64 212, !7, i64 216, !10, i64 224, !13, i64 232, !13, i64 240, !7, i64 248, !13, i64 256, !7, i64 264, !13, i64 272, !7, i64 280, !13, i64 288, !7, i64 296, !13, i64 304, !13, i64 312, !13, i64 320, !7, i64 328, !7, i64 336, !7, i64 344, !7, i64 352, !10, i64 360, !13, i64 368, !7, i64 376, !7, i64 384, !7, i64 392, !7, i64 400, !7, i64 408, !7, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !10, i64 464, !10, i64 468, !10, i64 472, !10, i64 476, !10, i64 480, !10, i64 484, !10, i64 488, !10, i64 492, !10, i64 496, !10, i64 500, !10, i64 504, !10, i64 508, !10, i64 512, !10, i64 516, !10, i64 520, !10, i64 524, !10, i64 528, !10, i64 532, !10, i64 536, !10, i64 540, !10, i64 544, !10, i64 548, !10, i64 552, !10, i64 556, !10, i64 560, !10, i64 564, !10, i64 568, !10, i64 572, !10, i64 576, !10, i64 580, !10, i64 584, !10, i64 588, !10, i64 592, !10, i64 596, !10, i64 600, !10, i64 604, !10, i64 608, !10, i64 612, !10, i64 616, !10, i64 620, !10, i64 624, !10, i64 628, !13, i64 632, !13, i64 640, !13, i64 648, !13, i64 656, !13, 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, !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, !14, i64 840, !14, i64 842, !14, i64 844, !14, i64 846, !14, i64 848, !14, i64 850, !8, i64 852, !14, i64 2900, !8, i64 2902}
!19 = !{!12, !7, i64 176}
!20 = !{!18, !7, i64 24}
!21 = !{!22, !7, i64 0}
!22 = !{!"thermo", !7, i64 0, !10, i64 8, !10, i64 12, !10, i64 16, !7, i64 24, !7, i64 32, !14, i64 40, !14, i64 42, !7, i64 48, !7, i64 56, !7, i64 64, !7, i64 72, !10, i64 80, !13, i64 88, !13, i64 96, !7, i64 104, !10, i64 112, !13, i64 120, !13, i64 128, !7, i64 136, !10, i64 144, !13, i64 152, !13, i64 160, !7, i64 168, !14, i64 176, !7, i64 184, !7, i64 192, !7, i64 200, !7, i64 208, !7, i64 216, !7, i64 224, !7, i64 232, !7, i64 240, !7, i64 248, !7, i64 256, !7, i64 264, !10, i64 272, !10, i64 276, !10, i64 280, !10, i64 284, !10, i64 288, !10, i64 292, !10, i64 296, !10, i64 300, !10, i64 304, !10, i64 308, !10, i64 312, !10, i64 316, !10, i64 320, !10, i64 324, !10, i64 328, !10, i64 332, !10, i64 336, !10, i64 340, !10, i64 344, !10, i64 348, !10, i64 352, !10, i64 356, !10, i64 360, !10, i64 364, !10, i64 368, !10, i64 372, !10, i64 376, !13, i64 384, !13, i64 392, !13, i64 400, !7, i64 408, !7, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !7, i64 464, !7, i64 472, !7, i64 480, !7, i64 488, !7, i64 496, !7, i64 504, !7, i64 512, !7, i64 520, !7, i64 528, !7, i64 536, !7, i64 544, !7, i64 552, !7, i64 560, !7, i64 568, !7, i64 576, !7, i64 584, !7, i64 592, !14, i64 600, !14, i64 602, !14, i64 604, !8, i64 606}
!23 = !{!12, !7, i64 184}
!24 = !{!18, !7, i64 16}
!25 = !{!12, !7, i64 128}
!26 = !{!6, !7, i64 360}
!27 = !{!6, !7, i64 584}
!28 = !{!12, !7, i64 120}
!29 = !{!6, !10, i64 416}
!30 = !{!6, !7, i64 424}
!31 = !{!12, !7, i64 112}
!32 = !{!6, !7, i64 432}
!33 = !{!6, !10, i64 400}
!34 = !{!12, !7, i64 80}
!35 = !{!12, !7, i64 136}
!36 = !{!22, !7, i64 592}
!37 = !{!22, !7, i64 168}
!38 = !{!12, !7, i64 192}
!39 = !{!22, !14, i64 176}
!40 = !{!12, !14, i64 200}
!41 = !{!7, !7, i64 0}
!42 = !{!22, !7, i64 216}
!43 = !{!12, !7, i64 240}
!44 = !{!12, !7, i64 104}
!45 = !{!12, !7, i64 72}
!46 = !{!47, !13, i64 0}
!47 = !{!"thermodynamics_parameters_and_workspace", !13, i64 0, !13, i64 8, !13, i64 16, !13, i64 24}
!48 = !{!47, !13, i64 8}
!49 = !{!47, !13, i64 16}
!50 = !{!47, !13, i64 24}
!51 = !{!6, !7, i64 520}
!52 = !{!6, !7, i64 488}
!53 = !{!6, !7, i64 496}
!54 = !{!6, !7, i64 504}
!55 = !{!6, !7, i64 512}
!56 = !{!12, !7, i64 48}
!57 = !{!12, !7, i64 168}
!58 = !{!6, !7, i64 528}
!59 = !{!6, !7, i64 536}
!60 = !{!12, !7, i64 160}
!61 = !{!6, !7, i64 560}
!62 = !{!12, !7, i64 152}
!63 = !{!6, !7, i64 392}
!64 = !{!6, !7, i64 2344}
!65 = !{!6, !7, i64 552}
!66 = !{!6, !7, i64 576}
!67 = !{!12, !10, i64 0}
!68 = !{!12, !10, i64 4}
!69 = !{!12, !10, i64 8}
!70 = !{!12, !10, i64 12}
!71 = !{!12, !10, i64 16}
!72 = !{!12, !10, i64 20}
!73 = distinct !{!73, !74}
!74 = !{!"llvm.loop.mustprogress"}


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