; LLVM IR from zhenrong-wang/co2 :: co2_property_calc (co2_prop_calc_th)

; 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.co2_prop = type { double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double }

; Function Attrs: nounwind uwtable
declare dso_local void @co2_prop_calc_tr(ptr noundef, double noundef, double noundef) #0

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

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

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

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

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

; Function Attrs: nounwind uwtable
define dso_local i32 @co2_prop_calc_th(ptr noundef %0, double noundef %1, double noundef %2) #0 {
  %4 = alloca i32, align 4
  %5 = alloca ptr, align 8
  %6 = alloca double, align 8
  %7 = alloca double, align 8
  %8 = alloca double, align 8
  %9 = alloca double, align 8
  %10 = alloca double, align 8
  %11 = alloca double, align 8
  %12 = alloca double, align 8
  %13 = alloca i32, align 4
  %14 = alloca %struct.co2_prop, align 8
  %15 = alloca %struct.co2_prop, align 8
  %16 = alloca i32, align 4
  store ptr %0, ptr %5, align 8, !tbaa !5
  store double %1, ptr %6, align 8, !tbaa !9
  store double %2, ptr %7, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %8) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %9) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %10) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %11) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #3
  call void @llvm.lifetime.start.p0(i64 4, ptr %13) #3
  call void @llvm.lifetime.start.p0(i64 128, ptr %14) #3
  call void @llvm.lifetime.start.p0(i64 128, ptr %15) #3
  %17 = load ptr, ptr %5, align 8, !tbaa !5
  %18 = getelementptr inbounds %struct.co2_prop, ptr %17, i32 0, i32 14
  store double -1.000000e+00, ptr %18, align 8, !tbaa !11
  %19 = load ptr, ptr %5, align 8, !tbaa !5
  %20 = getelementptr inbounds %struct.co2_prop, ptr %19, i32 0, i32 15
  store double -1.000000e+00, ptr %20, align 8, !tbaa !13
  %21 = load double, ptr %6, align 8, !tbaa !9
  %22 = fcmp olt double %21, 0x4073020D1B71758E
  br i1 %22, label %23, label %172

23:                                               ; preds = %3
  %24 = load double, ptr %6, align 8, !tbaa !9
  %25 = call i32 @calc_psat_t(ptr noundef %10, ptr noundef %8, ptr noundef %9, double noundef %24)
  store i32 %25, ptr %13, align 4, !tbaa !14
  %26 = load i32, ptr %13, align 4, !tbaa !14
  %27 = icmp eq i32 %26, 0
  br i1 %27, label %28, label %171

28:                                               ; preds = %23
  %29 = load double, ptr %6, align 8, !tbaa !9
  %30 = load double, ptr %8, align 8, !tbaa !9
  call void @co2_prop_calc_tr(ptr noundef %14, double noundef %29, double noundef %30)
  %31 = load double, ptr %6, align 8, !tbaa !9
  %32 = load double, ptr %9, align 8, !tbaa !9
  call void @co2_prop_calc_tr(ptr noundef %15, double noundef %31, double noundef %32)
  %33 = load double, ptr %7, align 8, !tbaa !9
  %34 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 7
  %35 = load double, ptr %34, align 8, !tbaa !16
  %36 = fcmp ogt double %33, %35
  br i1 %36, label %37, label %170

37:                                               ; preds = %28
  %38 = load double, ptr %7, align 8, !tbaa !9
  %39 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 7
  %40 = load double, ptr %39, align 8, !tbaa !16
  %41 = fcmp olt double %38, %40
  br i1 %41, label %42, label %170

42:                                               ; preds = %37
  %43 = load double, ptr %7, align 8, !tbaa !9
  %44 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 7
  %45 = load double, ptr %44, align 8, !tbaa !16
  %46 = fsub double %43, %45
  %47 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 7
  %48 = load double, ptr %47, align 8, !tbaa !16
  %49 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 7
  %50 = load double, ptr %49, align 8, !tbaa !16
  %51 = fsub double %48, %50
  %52 = fdiv double %46, %51
  store double %52, ptr %11, align 8, !tbaa !9
  %53 = load double, ptr %10, align 8, !tbaa !9
  %54 = load ptr, ptr %5, align 8, !tbaa !5
  %55 = getelementptr inbounds %struct.co2_prop, ptr %54, i32 0, i32 0
  store double %53, ptr %55, align 8, !tbaa !17
  %56 = load double, ptr %6, align 8, !tbaa !9
  %57 = load ptr, ptr %5, align 8, !tbaa !5
  %58 = getelementptr inbounds %struct.co2_prop, ptr %57, i32 0, i32 1
  store double %56, ptr %58, align 8, !tbaa !18
  %59 = load double, ptr %11, align 8, !tbaa !9
  %60 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 2
  %61 = load double, ptr %60, align 8, !tbaa !19
  %62 = load double, ptr %11, align 8, !tbaa !9
  %63 = fsub double 1.000000e+00, %62
  %64 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 2
  %65 = load double, ptr %64, align 8, !tbaa !19
  %66 = fmul double %63, %65
  %67 = call double @llvm.fmuladd.f64(double %59, double %61, double %66)
  %68 = load ptr, ptr %5, align 8, !tbaa !5
  %69 = getelementptr inbounds %struct.co2_prop, ptr %68, i32 0, i32 2
  store double %67, ptr %69, align 8, !tbaa !19
  %70 = load ptr, ptr %5, align 8, !tbaa !5
  %71 = getelementptr inbounds %struct.co2_prop, ptr %70, i32 0, i32 2
  %72 = load double, ptr %71, align 8, !tbaa !19
  %73 = fdiv double 1.000000e+00, %72
  %74 = load ptr, ptr %5, align 8, !tbaa !5
  %75 = getelementptr inbounds %struct.co2_prop, ptr %74, i32 0, i32 3
  store double %73, ptr %75, align 8, !tbaa !20
  %76 = load double, ptr %11, align 8, !tbaa !9
  %77 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 4
  %78 = load double, ptr %77, align 8, !tbaa !21
  %79 = load double, ptr %11, align 8, !tbaa !9
  %80 = fsub double 1.000000e+00, %79
  %81 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 4
  %82 = load double, ptr %81, align 8, !tbaa !21
  %83 = fmul double %80, %82
  %84 = call double @llvm.fmuladd.f64(double %76, double %78, double %83)
  %85 = load ptr, ptr %5, align 8, !tbaa !5
  %86 = getelementptr inbounds %struct.co2_prop, ptr %85, i32 0, i32 4
  store double %84, ptr %86, align 8, !tbaa !21
  %87 = load double, ptr %11, align 8, !tbaa !9
  %88 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 5
  %89 = load double, ptr %88, align 8, !tbaa !22
  %90 = load double, ptr %11, align 8, !tbaa !9
  %91 = fsub double 1.000000e+00, %90
  %92 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 5
  %93 = load double, ptr %92, align 8, !tbaa !22
  %94 = fmul double %91, %93
  %95 = call double @llvm.fmuladd.f64(double %87, double %89, double %94)
  %96 = load ptr, ptr %5, align 8, !tbaa !5
  %97 = getelementptr inbounds %struct.co2_prop, ptr %96, i32 0, i32 5
  store double %95, ptr %97, align 8, !tbaa !22
  %98 = load double, ptr %11, align 8, !tbaa !9
  %99 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 6
  %100 = load double, ptr %99, align 8, !tbaa !23
  %101 = load double, ptr %11, align 8, !tbaa !9
  %102 = fsub double 1.000000e+00, %101
  %103 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 6
  %104 = load double, ptr %103, align 8, !tbaa !23
  %105 = fmul double %102, %104
  %106 = call double @llvm.fmuladd.f64(double %98, double %100, double %105)
  %107 = load ptr, ptr %5, align 8, !tbaa !5
  %108 = getelementptr inbounds %struct.co2_prop, ptr %107, i32 0, i32 6
  store double %106, ptr %108, align 8, !tbaa !23
  %109 = load double, ptr %7, align 8, !tbaa !9
  %110 = load ptr, ptr %5, align 8, !tbaa !5
  %111 = getelementptr inbounds %struct.co2_prop, ptr %110, i32 0, i32 7
  store double %109, ptr %111, align 8, !tbaa !16
  %112 = load double, ptr %11, align 8, !tbaa !9
  %113 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 8
  %114 = load double, ptr %113, align 8, !tbaa !24
  %115 = load double, ptr %11, align 8, !tbaa !9
  %116 = fsub double 1.000000e+00, %115
  %117 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 8
  %118 = load double, ptr %117, align 8, !tbaa !24
  %119 = fmul double %116, %118
  %120 = call double @llvm.fmuladd.f64(double %112, double %114, double %119)
  %121 = load ptr, ptr %5, align 8, !tbaa !5
  %122 = getelementptr inbounds %struct.co2_prop, ptr %121, i32 0, i32 8
  store double %120, ptr %122, align 8, !tbaa !24
  %123 = load double, ptr %11, align 8, !tbaa !9
  %124 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 9
  %125 = load double, ptr %124, align 8, !tbaa !25
  %126 = load double, ptr %11, align 8, !tbaa !9
  %127 = fsub double 1.000000e+00, %126
  %128 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 9
  %129 = load double, ptr %128, align 8, !tbaa !25
  %130 = fmul double %127, %129
  %131 = call double @llvm.fmuladd.f64(double %123, double %125, double %130)
  %132 = load ptr, ptr %5, align 8, !tbaa !5
  %133 = getelementptr inbounds %struct.co2_prop, ptr %132, i32 0, i32 9
  store double %131, ptr %133, align 8, !tbaa !25
  %134 = load double, ptr %11, align 8, !tbaa !9
  %135 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 10
  %136 = load double, ptr %135, align 8, !tbaa !26
  %137 = load double, ptr %11, align 8, !tbaa !9
  %138 = fsub double 1.000000e+00, %137
  %139 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 10
  %140 = load double, ptr %139, align 8, !tbaa !26
  %141 = fmul double %138, %140
  %142 = call double @llvm.fmuladd.f64(double %134, double %136, double %141)
  %143 = load ptr, ptr %5, align 8, !tbaa !5
  %144 = getelementptr inbounds %struct.co2_prop, ptr %143, i32 0, i32 10
  store double %142, ptr %144, align 8, !tbaa !26
  %145 = load double, ptr %11, align 8, !tbaa !9
  %146 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 11
  %147 = load double, ptr %146, align 8, !tbaa !27
  %148 = load double, ptr %11, align 8, !tbaa !9
  %149 = fsub double 1.000000e+00, %148
  %150 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 11
  %151 = load double, ptr %150, align 8, !tbaa !27
  %152 = fmul double %149, %151
  %153 = call double @llvm.fmuladd.f64(double %145, double %147, double %152)
  %154 = load ptr, ptr %5, align 8, !tbaa !5
  %155 = getelementptr inbounds %struct.co2_prop, ptr %154, i32 0, i32 11
  store double %153, ptr %155, align 8, !tbaa !27
  %156 = load double, ptr %11, align 8, !tbaa !9
  %157 = getelementptr inbounds %struct.co2_prop, ptr %15, i32 0, i32 12
  %158 = load double, ptr %157, align 8, !tbaa !28
  %159 = load double, ptr %11, align 8, !tbaa !9
  %160 = fsub double 1.000000e+00, %159
  %161 = getelementptr inbounds %struct.co2_prop, ptr %14, i32 0, i32 12
  %162 = load double, ptr %161, align 8, !tbaa !28
  %163 = fmul double %160, %162
  %164 = call double @llvm.fmuladd.f64(double %156, double %158, double %163)
  %165 = load ptr, ptr %5, align 8, !tbaa !5
  %166 = getelementptr inbounds %struct.co2_prop, ptr %165, i32 0, i32 12
  store double %164, ptr %166, align 8, !tbaa !28
  %167 = load double, ptr %11, align 8, !tbaa !9
  %168 = load ptr, ptr %5, align 8, !tbaa !5
  %169 = getelementptr inbounds %struct.co2_prop, ptr %168, i32 0, i32 13
  store double %167, ptr %169, align 8, !tbaa !29
  store i32 0, ptr %4, align 4
  store i32 1, ptr %16, align 4
  br label %183

170:                                              ; preds = %37, %28
  br label %171

171:                                              ; preds = %170, %23
  br label %172

172:                                              ; preds = %171, %3
  %173 = load double, ptr %6, align 8, !tbaa !9
  %174 = load double, ptr %7, align 8, !tbaa !9
  %175 = call i32 @dens_th(ptr noundef %12, double noundef %173, double noundef %174)
  store i32 %175, ptr %13, align 4, !tbaa !14
  %176 = load i32, ptr %13, align 4, !tbaa !14
  %177 = icmp eq i32 %176, -1
  br i1 %177, label %178, label %179

178:                                              ; preds = %172
  store i32 -1, ptr %4, align 4
  store i32 1, ptr %16, align 4
  br label %183

179:                                              ; preds = %172
  %180 = load ptr, ptr %5, align 8, !tbaa !5
  %181 = load double, ptr %6, align 8, !tbaa !9
  %182 = load double, ptr %12, align 8, !tbaa !9
  call void @co2_prop_calc_tr(ptr noundef %180, double noundef %181, double noundef %182)
  store i32 0, ptr %4, align 4
  store i32 1, ptr %16, align 4
  br label %183

183:                                              ; preds = %179, %178, %42
  call void @llvm.lifetime.end.p0(i64 128, ptr %15) #3
  call void @llvm.lifetime.end.p0(i64 128, ptr %14) #3
  call void @llvm.lifetime.end.p0(i64 4, ptr %13) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %11) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %10) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %9) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %8) #3
  %184 = load i32, ptr %4, align 4
  ret i32 %184
}

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

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

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


; 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.co2_prop = type { double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double }

; Function Attrs: nounwind uwtable
declare dso_local void @co2_prop_calc_tr(ptr noundef, double noundef, double noundef) local_unnamed_addr #0

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

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

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

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

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

; Function Attrs: nounwind uwtable
define dso_local range(i32 -1, 1) i32 @co2_prop_calc_th(ptr noundef %0, double noundef %1, double noundef %2) local_unnamed_addr #0 {
  %4 = alloca double, align 8
  %5 = alloca double, align 8
  %6 = alloca double, align 8
  %7 = alloca double, align 8
  %8 = alloca %struct.co2_prop, align 8
  %9 = alloca %struct.co2_prop, align 8
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %4) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %5) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %6) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %7) #4
  call void @llvm.lifetime.start.p0(i64 128, ptr nonnull %8) #4
  call void @llvm.lifetime.start.p0(i64 128, ptr nonnull %9) #4
  %10 = getelementptr inbounds i8, ptr %0, i64 112
  store <2 x double> <double -1.000000e+00, double -1.000000e+00>, ptr %10, align 8, !tbaa !5
  %11 = fcmp olt double %1, 0x4073020D1B71758E
  br i1 %11, label %12, label %82

12:                                               ; preds = %3
  %13 = call i32 @calc_psat_t(ptr noundef nonnull %6, ptr noundef nonnull %4, ptr noundef nonnull %5, double noundef %1)
  %14 = icmp eq i32 %13, 0
  br i1 %14, label %15, label %82

15:                                               ; preds = %12
  %16 = load double, ptr %4, align 8, !tbaa !5
  call void @co2_prop_calc_tr(ptr noundef nonnull %8, double noundef %1, double noundef %16)
  %17 = load double, ptr %5, align 8, !tbaa !5
  call void @co2_prop_calc_tr(ptr noundef nonnull %9, double noundef %1, double noundef %17)
  %18 = getelementptr inbounds i8, ptr %8, i64 56
  %19 = load double, ptr %18, align 8, !tbaa !9
  %20 = fcmp olt double %19, %2
  br i1 %20, label %21, label %82

21:                                               ; preds = %15
  %22 = getelementptr inbounds i8, ptr %9, i64 56
  %23 = load double, ptr %22, align 8, !tbaa !9
  %24 = fcmp ogt double %23, %2
  br i1 %24, label %25, label %82

25:                                               ; preds = %21
  %26 = fsub double %2, %19
  %27 = fsub double %23, %19
  %28 = fdiv double %26, %27
  %29 = load double, ptr %6, align 8, !tbaa !5
  store double %29, ptr %0, align 8, !tbaa !11
  %30 = getelementptr inbounds i8, ptr %0, i64 8
  store double %1, ptr %30, align 8, !tbaa !12
  %31 = getelementptr inbounds i8, ptr %9, i64 16
  %32 = load double, ptr %31, align 8, !tbaa !13
  %33 = fsub double 1.000000e+00, %28
  %34 = getelementptr inbounds i8, ptr %8, i64 16
  %35 = load double, ptr %34, align 8, !tbaa !13
  %36 = fmul double %33, %35
  %37 = call double @llvm.fmuladd.f64(double %28, double %32, double %36)
  %38 = getelementptr inbounds i8, ptr %0, i64 16
  store double %37, ptr %38, align 8, !tbaa !13
  %39 = fdiv double 1.000000e+00, %37
  %40 = getelementptr inbounds i8, ptr %0, i64 24
  store double %39, ptr %40, align 8, !tbaa !14
  %41 = getelementptr inbounds i8, ptr %9, i64 32
  %42 = getelementptr inbounds i8, ptr %8, i64 32
  %43 = getelementptr inbounds i8, ptr %0, i64 32
  %44 = load <2 x double>, ptr %41, align 8, !tbaa !5
  %45 = load <2 x double>, ptr %42, align 8, !tbaa !5
  %46 = insertelement <2 x double> poison, double %33, i64 0
  %47 = shufflevector <2 x double> %46, <2 x double> poison, <2 x i32> zeroinitializer
  %48 = fmul <2 x double> %47, %45
  %49 = insertelement <2 x double> poison, double %28, i64 0
  %50 = shufflevector <2 x double> %49, <2 x double> poison, <2 x i32> zeroinitializer
  %51 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %50, <2 x double> %44, <2 x double> %48)
  store <2 x double> %51, ptr %43, align 8, !tbaa !5
  %52 = getelementptr inbounds i8, ptr %9, i64 48
  %53 = load double, ptr %52, align 8, !tbaa !15
  %54 = getelementptr inbounds i8, ptr %8, i64 48
  %55 = load double, ptr %54, align 8, !tbaa !15
  %56 = fmul double %33, %55
  %57 = call double @llvm.fmuladd.f64(double %28, double %53, double %56)
  %58 = getelementptr inbounds i8, ptr %0, i64 48
  store double %57, ptr %58, align 8, !tbaa !15
  %59 = getelementptr inbounds i8, ptr %0, i64 56
  store double %2, ptr %59, align 8, !tbaa !9
  %60 = getelementptr inbounds i8, ptr %9, i64 64
  %61 = getelementptr inbounds i8, ptr %8, i64 64
  %62 = getelementptr inbounds i8, ptr %0, i64 64
  %63 = load <2 x double>, ptr %60, align 8, !tbaa !5
  %64 = load <2 x double>, ptr %61, align 8, !tbaa !5
  %65 = fmul <2 x double> %47, %64
  %66 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %50, <2 x double> %63, <2 x double> %65)
  store <2 x double> %66, ptr %62, align 8, !tbaa !5
  %67 = getelementptr inbounds i8, ptr %9, i64 80
  %68 = getelementptr inbounds i8, ptr %8, i64 80
  %69 = getelementptr inbounds i8, ptr %0, i64 80
  %70 = load <2 x double>, ptr %67, align 8, !tbaa !5
  %71 = load <2 x double>, ptr %68, align 8, !tbaa !5
  %72 = fmul <2 x double> %47, %71
  %73 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %50, <2 x double> %70, <2 x double> %72)
  store <2 x double> %73, ptr %69, align 8, !tbaa !5
  %74 = getelementptr inbounds i8, ptr %9, i64 96
  %75 = load double, ptr %74, align 8, !tbaa !16
  %76 = getelementptr inbounds i8, ptr %8, i64 96
  %77 = load double, ptr %76, align 8, !tbaa !16
  %78 = fmul double %33, %77
  %79 = call double @llvm.fmuladd.f64(double %28, double %75, double %78)
  %80 = getelementptr inbounds i8, ptr %0, i64 96
  store double %79, ptr %80, align 8, !tbaa !16
  %81 = getelementptr inbounds i8, ptr %0, i64 104
  store double %28, ptr %81, align 8, !tbaa !17
  br label %87

82:                                               ; preds = %12, %21, %15, %3
  %83 = call i32 @dens_th(ptr noundef nonnull %7, double noundef %1, double noundef %2)
  %84 = icmp eq i32 %83, -1
  br i1 %84, label %87, label %85

85:                                               ; preds = %82
  %86 = load double, ptr %7, align 8, !tbaa !5
  call void @co2_prop_calc_tr(ptr noundef nonnull %0, double noundef %1, double noundef %86)
  br label %87

87:                                               ; preds = %82, %85, %25
  %.0 = phi i32 [ 0, %25 ], [ 0, %85 ], [ -1, %82 ]
  call void @llvm.lifetime.end.p0(i64 128, ptr nonnull %9) #4
  call void @llvm.lifetime.end.p0(i64 128, ptr nonnull %8) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %7) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %6) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %5) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %4) #4
  ret i32 %.0
}

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

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

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"double", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !6, i64 56}
!10 = !{!"", !6, i64 0, !6, i64 8, !6, i64 16, !6, i64 24, !6, i64 32, !6, i64 40, !6, i64 48, !6, i64 56, !6, i64 64, !6, i64 72, !6, i64 80, !6, i64 88, !6, i64 96, !6, i64 104, !6, i64 112, !6, i64 120}
!11 = !{!10, !6, i64 0}
!12 = !{!10, !6, i64 8}
!13 = !{!10, !6, i64 16}
!14 = !{!10, !6, i64 24}
!15 = !{!10, !6, i64 48}
!16 = !{!10, !6, i64 96}
!17 = !{!10, !6, i64 104}


; Effective O3 passes: CorrelatedValuePropagationPass,EarlyCSEPass,GlobalOptPass,IPSCCPPass,InferFunctionAttrsPass,InstCombinePass,SLPVectorizerPass,SROAPass,SimplifyCFGPass