; LLVM IR from ivan-pi/periodic-lbm :: third-party (dugks_cavity)

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

@ex = external dso_local global [3 x double], align 16
@ey = external dso_local global [3 x double], align 16
@RT = external dso_local global double, align 8
@nu = external dso_local global double, align 8
@dx = external dso_local global double, align 8
@ti = external dso_local global [33 x double], align 16
@yc = external dso_local global [32 x double], align 16
@xc = external dso_local global [32 x double], align 16
@dy = external dso_local global double, align 8
@dt = external dso_local global double, align 8
@tau = external dso_local global double, align 8
@w = external dso_local global double, align 8
@wb = external dso_local global double, align 8
@.str = external hidden unnamed_addr constant [7 x i8], align 1
@ux = external dso_local global [32 x [32 x double]], align 16
@.str.1 = external hidden unnamed_addr constant [13 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [3 x i8], align 1
@.str.3 = external hidden unnamed_addr constant [19 x i8], align 1
@.str.4 = external hidden unnamed_addr constant [20 x i8], align 1
@.str.5 = external hidden unnamed_addr constant [24 x i8], align 1

; Function Attrs: nounwind uwtable
define dso_local void @main() #0 {
  %1 = alloca i32, align 4
  %2 = alloca i32, align 4
  %3 = alloca i32, align 4
  %4 = alloca double, align 8
  %5 = alloca double, align 8
  %6 = alloca double, align 8
  call void @llvm.lifetime.start.p0(i64 4, ptr %1) #5
  call void @llvm.lifetime.start.p0(i64 4, ptr %2) #5
  call void @llvm.lifetime.start.p0(i64 4, ptr %3) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %4) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %5) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %6) #5
  store double 0x3FD5555555555555, ptr @RT, align 8, !tbaa !5
  store i32 0, ptr %1, align 4, !tbaa !9
  br label %7

7:                                                ; preds = %33, %0
  %8 = load i32, ptr %1, align 4, !tbaa !9
  %9 = icmp slt i32 %8, 3
  br i1 %9, label %10, label %36

10:                                               ; preds = %7
  %11 = load double, ptr @RT, align 8, !tbaa !5
  %12 = fmul double 3.000000e+00, %11
  %13 = call double @sqrt(double noundef %12) #5
  %14 = load i32, ptr %1, align 4, !tbaa !9
  %15 = sext i32 %14 to i64
  %16 = getelementptr inbounds [3 x double], ptr @ex, i64 0, i64 %15
  %17 = load double, ptr %16, align 8, !tbaa !5
  %18 = fmul double %13, %17
  %19 = load i32, ptr %1, align 4, !tbaa !9
  %20 = sext i32 %19 to i64
  %21 = getelementptr inbounds [3 x double], ptr @ex, i64 0, i64 %20
  store double %18, ptr %21, align 8, !tbaa !5
  %22 = load double, ptr @RT, align 8, !tbaa !5
  %23 = fmul double 3.000000e+00, %22
  %24 = call double @sqrt(double noundef %23) #5
  %25 = load i32, ptr %1, align 4, !tbaa !9
  %26 = sext i32 %25 to i64
  %27 = getelementptr inbounds [3 x double], ptr @ey, i64 0, i64 %26
  %28 = load double, ptr %27, align 8, !tbaa !5
  %29 = fmul double %24, %28
  %30 = load i32, ptr %1, align 4, !tbaa !9
  %31 = sext i32 %30 to i64
  %32 = getelementptr inbounds [3 x double], ptr @ey, i64 0, i64 %31
  store double %29, ptr %32, align 8, !tbaa !5
  br label %33

33:                                               ; preds = %10
  %34 = load i32, ptr %1, align 4, !tbaa !9
  %35 = add nsw i32 %34, 1
  store i32 %35, ptr %1, align 4, !tbaa !9
  br label %7, !llvm.loop !11

36:                                               ; preds = %7
  store double 1.000000e-04, ptr @nu, align 8, !tbaa !5
  store double 0x3FA1111111111111, ptr @dx, align 8, !tbaa !5
  store i32 1, ptr %1, align 4, !tbaa !9
  br label %37

37:                                               ; preds = %57, %36
  %38 = load i32, ptr %1, align 4, !tbaa !9
  %39 = icmp sle i32 %38, 31
  br i1 %39, label %40, label %60

40:                                               ; preds = %37
  %41 = load i32, ptr %1, align 4, !tbaa !9
  %42 = sub nsw i32 %41, 1
  %43 = sitofp i32 %42 to double
  %44 = load double, ptr @dx, align 8, !tbaa !5
  %45 = fmul double %43, %44
  store double %45, ptr %6, align 8, !tbaa !5
  %46 = load double, ptr %6, align 8, !tbaa !5
  %47 = fsub double %46, 5.000000e-01
  %48 = fmul double 2.000000e+00, %47
  %49 = call double @tanh(double noundef %48) #5
  %50 = call double @tanh(double noundef 1.000000e+00) #5
  %51 = fdiv double %49, %50
  %52 = fadd double %51, 1.000000e+00
  %53 = fdiv double %52, 2.000000e+00
  %54 = load i32, ptr %1, align 4, !tbaa !9
  %55 = sext i32 %54 to i64
  %56 = getelementptr inbounds [33 x double], ptr @ti, i64 0, i64 %55
  store double %53, ptr %56, align 8, !tbaa !5
  br label %57

57:                                               ; preds = %40
  %58 = load i32, ptr %1, align 4, !tbaa !9
  %59 = add nsw i32 %58, 1
  store i32 %59, ptr %1, align 4, !tbaa !9
  br label %37, !llvm.loop !13

60:                                               ; preds = %37
  store i32 1, ptr %1, align 4, !tbaa !9
  br label %61

61:                                               ; preds = %82, %60
  %62 = load i32, ptr %1, align 4, !tbaa !9
  %63 = icmp sle i32 %62, 30
  br i1 %63, label %64, label %85

64:                                               ; preds = %61
  %65 = load i32, ptr %1, align 4, !tbaa !9
  %66 = sext i32 %65 to i64
  %67 = getelementptr inbounds [33 x double], ptr @ti, i64 0, i64 %66
  %68 = load double, ptr %67, align 8, !tbaa !5
  %69 = load i32, ptr %1, align 4, !tbaa !9
  %70 = add nsw i32 %69, 1
  %71 = sext i32 %70 to i64
  %72 = getelementptr inbounds [33 x double], ptr @ti, i64 0, i64 %71
  %73 = load double, ptr %72, align 8, !tbaa !5
  %74 = fadd double %68, %73
  %75 = fmul double 5.000000e-01, %74
  %76 = load i32, ptr %1, align 4, !tbaa !9
  %77 = sext i32 %76 to i64
  %78 = getelementptr inbounds [32 x double], ptr @yc, i64 0, i64 %77
  store double %75, ptr %78, align 8, !tbaa !5
  %79 = load i32, ptr %1, align 4, !tbaa !9
  %80 = sext i32 %79 to i64
  %81 = getelementptr inbounds [32 x double], ptr @xc, i64 0, i64 %80
  store double %75, ptr %81, align 8, !tbaa !5
  br label %82

82:                                               ; preds = %64
  %83 = load i32, ptr %1, align 4, !tbaa !9
  %84 = add nsw i32 %83, 1
  store i32 %84, ptr %1, align 4, !tbaa !9
  br label %61, !llvm.loop !14

85:                                               ; preds = %61
  %86 = load double, ptr getelementptr inbounds ([32 x double], ptr @xc, i64 0, i64 1), align 8, !tbaa !5
  %87 = fneg double %86
  store double %87, ptr @xc, align 16, !tbaa !5
  %88 = load double, ptr getelementptr inbounds ([32 x double], ptr @xc, i64 0, i64 30), align 16, !tbaa !5
  %89 = fsub double 2.000000e+00, %88
  store double %89, ptr getelementptr inbounds ([32 x double], ptr @xc, i64 0, i64 31), align 8, !tbaa !5
  %90 = load double, ptr getelementptr inbounds ([32 x double], ptr @yc, i64 0, i64 1), align 8, !tbaa !5
  %91 = fneg double %90
  store double %91, ptr @yc, align 16, !tbaa !5
  %92 = load double, ptr getelementptr inbounds ([32 x double], ptr @yc, i64 0, i64 30), align 16, !tbaa !5
  %93 = fsub double 2.000000e+00, %92
  store double %93, ptr getelementptr inbounds ([32 x double], ptr @yc, i64 0, i64 31), align 8, !tbaa !5
  %94 = load double, ptr getelementptr inbounds ([32 x double], ptr @xc, i64 0, i64 1), align 8, !tbaa !5
  %95 = load double, ptr @xc, align 16, !tbaa !5
  %96 = fsub double %94, %95
  store double %96, ptr @dx, align 8, !tbaa !5
  store double %96, ptr @dy, align 8, !tbaa !5
  %97 = load double, ptr @dx, align 8, !tbaa !5
  %98 = fmul double 1.000000e-01, %97
  %99 = load double, ptr @RT, align 8, !tbaa !5
  %100 = fmul double 6.000000e+00, %99
  %101 = call double @sqrt(double noundef %100) #5
  %102 = fdiv double %98, %101
  store double %102, ptr @dt, align 8, !tbaa !5
  %103 = load double, ptr @nu, align 8, !tbaa !5
  %104 = load double, ptr @RT, align 8, !tbaa !5
  %105 = fdiv double %103, %104
  store double %105, ptr @tau, align 8, !tbaa !5
  %106 = load double, ptr @dt, align 8, !tbaa !5
  %107 = fmul double 1.500000e+00, %106
  %108 = load double, ptr @tau, align 8, !tbaa !5
  %109 = load double, ptr @dt, align 8, !tbaa !5
  %110 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %108, double %109)
  %111 = fdiv double %107, %110
  store double %111, ptr @w, align 8, !tbaa !5
  %112 = load double, ptr @dt, align 8, !tbaa !5
  %113 = load double, ptr @tau, align 8, !tbaa !5
  %114 = load double, ptr @dt, align 8, !tbaa !5
  %115 = fmul double 5.000000e-01, %114
  %116 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %113, double %115)
  %117 = fdiv double %112, %116
  store double %117, ptr @wb, align 8, !tbaa !5
  %118 = load double, ptr @w, align 8, !tbaa !5
  %119 = call i32 (ptr, ...) @printf(ptr noundef @.str, double noundef %118)
  call void @gks_ini()
  %120 = load double, ptr getelementptr inbounds ([32 x double], ptr getelementptr inbounds ([32 x [32 x double]], ptr @ux, i64 0, i64 16), i64 0, i64 16), align 16, !tbaa !5
  store double %120, ptr %5, align 8, !tbaa !5
  store i32 0, ptr %1, align 4, !tbaa !9
  br label %121

121:                                              ; preds = %162, %85
  %122 = call i32 (ptr, ...) @printf(ptr noundef @.str.1)
  %123 = call i32 (ptr, ...) @__isoc99_scanf(ptr noundef @.str.2, ptr noundef %3)
  %124 = load i32, ptr %1, align 4, !tbaa !9
  %125 = load i32, ptr %3, align 4, !tbaa !9
  %126 = add nsw i32 %125, %124
  store i32 %126, ptr %3, align 4, !tbaa !9
  %127 = load double, ptr @dt, align 8, !tbaa !5
  %128 = load i32, ptr %3, align 4, !tbaa !9
  %129 = call i32 (ptr, ...) @printf(ptr noundef @.str.3, double noundef %127, i32 noundef %128)
  store double 1.000000e+00, ptr %4, align 8, !tbaa !5
  br label %130

130:                                              ; preds = %156, %121
  %131 = load i32, ptr %1, align 4, !tbaa !9
  %132 = load i32, ptr %3, align 4, !tbaa !9
  %133 = icmp slt i32 %131, %132
  br i1 %133, label %134, label %137

134:                                              ; preds = %130
  %135 = load double, ptr %4, align 8, !tbaa !5
  %136 = fcmp ogt double %135, 0x3EB0C6F7A0B5ED8D
  br label %137

137:                                              ; preds = %134, %130
  %138 = phi i1 [ false, %130 ], [ %136, %134 ]
  br i1 %138, label %139, label %157

139:                                              ; preds = %137
  %140 = load i32, ptr %1, align 4, !tbaa !9
  %141 = add nsw i32 %140, 1
  store i32 %141, ptr %1, align 4, !tbaa !9
  call void @Evol()
  %142 = load i32, ptr %1, align 4, !tbaa !9
  %143 = srem i32 %142, 1000
  %144 = icmp eq i32 %143, 0
  br i1 %144, label %145, label %156

145:                                              ; preds = %139
  %146 = load double, ptr getelementptr inbounds ([32 x double], ptr getelementptr inbounds ([32 x [32 x double]], ptr @ux, i64 0, i64 16), i64 0, i64 16), align 16, !tbaa !5
  %147 = load double, ptr %5, align 8, !tbaa !5
  %148 = fsub double %146, %147
  %149 = call double @llvm.fabs.f64(double %148)
  %150 = fdiv double %149, 1.000000e-01
  store double %150, ptr %4, align 8, !tbaa !5
  %151 = load double, ptr getelementptr inbounds ([32 x double], ptr getelementptr inbounds ([32 x [32 x double]], ptr @ux, i64 0, i64 16), i64 0, i64 16), align 16, !tbaa !5
  store double %151, ptr %5, align 8, !tbaa !5
  %152 = load double, ptr %4, align 8, !tbaa !5
  %153 = load double, ptr %5, align 8, !tbaa !5
  %154 = load i32, ptr %1, align 4, !tbaa !9
  %155 = call i32 (ptr, ...) @printf(ptr noundef @.str.4, double noundef %152, double noundef %153, i32 noundef %154)
  br label %156

156:                                              ; preds = %145, %139
  br label %130, !llvm.loop !15

157:                                              ; preds = %137
  call void @datadeal()
  %158 = call i32 (ptr, ...) @printf(ptr noundef @.str.5)
  %159 = call i32 (ptr, ...) @__isoc99_scanf(ptr noundef @.str.2, ptr noundef %2)
  %160 = load i32, ptr %2, align 4, !tbaa !9
  %161 = icmp ne i32 %160, 0
  br i1 %161, label %162, label %163

162:                                              ; preds = %157
  br label %121

163:                                              ; preds = %157
  call void @llvm.lifetime.end.p0(i64 8, ptr %6) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %5) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %4) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %3) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %2) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %1) #5
  ret void
}

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

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

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

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

declare i32 @printf(ptr noundef, ...) #4

; Function Attrs: nounwind uwtable
declare dso_local void @gks_ini() #0

declare i32 @__isoc99_scanf(ptr noundef, ...) #4

; Function Attrs: nounwind uwtable
declare dso_local void @Evol() #0

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

; Function Attrs: nounwind uwtable
declare dso_local void @datadeal() #0

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

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 = { 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 #3 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #4 = { "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #5 = { nounwind }

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"double", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"int", !7, i64 0}
!11 = distinct !{!11, !12}
!12 = !{!"llvm.loop.mustprogress"}
!13 = distinct !{!13, !12}
!14 = distinct !{!14, !12}
!15 = distinct !{!15, !12}


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

@RT = external dso_local local_unnamed_addr global double, align 8
@nu = external dso_local local_unnamed_addr global double, align 8
@dx = external dso_local local_unnamed_addr global double, align 8
@ti = external dso_local local_unnamed_addr global [33 x double], align 16
@yc = external dso_local local_unnamed_addr global [32 x double], align 16
@xc = external dso_local local_unnamed_addr global [32 x double], align 16
@dy = external dso_local local_unnamed_addr global double, align 8
@dt = external dso_local local_unnamed_addr global double, align 8
@tau = external dso_local local_unnamed_addr global double, align 8
@w = external dso_local local_unnamed_addr global double, align 8
@wb = external dso_local local_unnamed_addr global double, align 8
@.str = external hidden unnamed_addr constant [7 x i8], align 1
@ux = external dso_local local_unnamed_addr global [32 x [32 x double]], align 16
@.str.1 = external hidden unnamed_addr constant [13 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [3 x i8], align 1
@.str.3 = external hidden unnamed_addr constant [19 x i8], align 1
@.str.4 = external hidden unnamed_addr constant [20 x i8], align 1
@.str.5 = external hidden unnamed_addr constant [24 x i8], align 1

; Function Attrs: nounwind uwtable
define dso_local void @main() local_unnamed_addr #0 {
  %1 = alloca i32, align 4
  %2 = alloca i32, align 4
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %1) #5
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %2) #5
  store double 0x3FD5555555555555, ptr @RT, align 8, !tbaa !5
  store double 1.000000e-04, ptr @nu, align 8, !tbaa !5
  store double 0x3FA1111111111111, ptr @dx, align 8, !tbaa !5
  br label %3

3:                                                ; preds = %0, %3
  %indvars.iv = phi i64 [ 1, %0 ], [ %indvars.iv.next, %3 ]
  %4 = trunc i64 %indvars.iv to i32
  %5 = add i32 %4, -1
  %6 = sitofp i32 %5 to double
  %7 = load double, ptr @dx, align 8, !tbaa !5
  %8 = fmul double %7, %6
  %9 = fadd double %8, -5.000000e-01
  %10 = fmul double %9, 2.000000e+00
  %11 = tail call double @tanh(double noundef %10) #5
  %12 = tail call double @tanh(double noundef 1.000000e+00) #5
  %13 = fdiv double %11, 0x3FE85EFAB514F394
  %14 = fadd double %13, 1.000000e+00
  %15 = fmul double %14, 5.000000e-01
  %16 = getelementptr inbounds [33 x double], ptr @ti, i64 0, i64 %indvars.iv
  store double %15, ptr %16, align 8, !tbaa !5
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next, 32
  br i1 %exitcond.not, label %.preheader.preheader, label %3, !llvm.loop !9

.preheader.preheader:                             ; preds = %3
  %17 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 8), align 8, !tbaa !5
  %18 = shufflevector <2 x double> %17, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %19 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 24), align 8, !tbaa !5
  %20 = shufflevector <2 x double> %18, <2 x double> %19, <2 x i32> <i32 0, i32 2>
  %21 = fadd <2 x double> %17, %20
  %22 = fmul <2 x double> %21, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %22, ptr getelementptr inbounds (i8, ptr @yc, i64 8), align 8, !tbaa !5
  store <2 x double> %22, ptr getelementptr inbounds (i8, ptr @xc, i64 8), align 8, !tbaa !5
  %23 = shufflevector <2 x double> %19, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %24 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 40), align 8, !tbaa !5
  %25 = shufflevector <2 x double> %23, <2 x double> %24, <2 x i32> <i32 0, i32 2>
  %26 = fadd <2 x double> %19, %25
  %27 = fmul <2 x double> %26, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %27, ptr getelementptr inbounds (i8, ptr @yc, i64 24), align 8, !tbaa !5
  store <2 x double> %27, ptr getelementptr inbounds (i8, ptr @xc, i64 24), align 8, !tbaa !5
  %28 = shufflevector <2 x double> %24, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %29 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 56), align 8, !tbaa !5
  %30 = shufflevector <2 x double> %28, <2 x double> %29, <2 x i32> <i32 0, i32 2>
  %31 = fadd <2 x double> %24, %30
  %32 = fmul <2 x double> %31, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %32, ptr getelementptr inbounds (i8, ptr @yc, i64 40), align 8, !tbaa !5
  store <2 x double> %32, ptr getelementptr inbounds (i8, ptr @xc, i64 40), align 8, !tbaa !5
  %33 = shufflevector <2 x double> %29, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %34 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 72), align 8, !tbaa !5
  %35 = shufflevector <2 x double> %33, <2 x double> %34, <2 x i32> <i32 0, i32 2>
  %36 = fadd <2 x double> %29, %35
  %37 = fmul <2 x double> %36, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %37, ptr getelementptr inbounds (i8, ptr @yc, i64 56), align 8, !tbaa !5
  store <2 x double> %37, ptr getelementptr inbounds (i8, ptr @xc, i64 56), align 8, !tbaa !5
  %38 = shufflevector <2 x double> %34, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %39 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 88), align 8, !tbaa !5
  %40 = shufflevector <2 x double> %38, <2 x double> %39, <2 x i32> <i32 0, i32 2>
  %41 = fadd <2 x double> %34, %40
  %42 = fmul <2 x double> %41, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %42, ptr getelementptr inbounds (i8, ptr @yc, i64 72), align 8, !tbaa !5
  store <2 x double> %42, ptr getelementptr inbounds (i8, ptr @xc, i64 72), align 8, !tbaa !5
  %43 = shufflevector <2 x double> %39, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %44 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 104), align 8, !tbaa !5
  %45 = shufflevector <2 x double> %43, <2 x double> %44, <2 x i32> <i32 0, i32 2>
  %46 = fadd <2 x double> %39, %45
  %47 = fmul <2 x double> %46, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %47, ptr getelementptr inbounds (i8, ptr @yc, i64 88), align 8, !tbaa !5
  store <2 x double> %47, ptr getelementptr inbounds (i8, ptr @xc, i64 88), align 8, !tbaa !5
  %48 = shufflevector <2 x double> %44, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %49 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 120), align 8, !tbaa !5
  %50 = shufflevector <2 x double> %48, <2 x double> %49, <2 x i32> <i32 0, i32 2>
  %51 = fadd <2 x double> %44, %50
  %52 = fmul <2 x double> %51, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %52, ptr getelementptr inbounds (i8, ptr @yc, i64 104), align 8, !tbaa !5
  store <2 x double> %52, ptr getelementptr inbounds (i8, ptr @xc, i64 104), align 8, !tbaa !5
  %53 = shufflevector <2 x double> %49, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %54 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 136), align 8, !tbaa !5
  %55 = shufflevector <2 x double> %53, <2 x double> %54, <2 x i32> <i32 0, i32 2>
  %56 = fadd <2 x double> %49, %55
  %57 = fmul <2 x double> %56, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %57, ptr getelementptr inbounds (i8, ptr @yc, i64 120), align 8, !tbaa !5
  store <2 x double> %57, ptr getelementptr inbounds (i8, ptr @xc, i64 120), align 8, !tbaa !5
  %58 = shufflevector <2 x double> %54, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %59 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 152), align 8, !tbaa !5
  %60 = shufflevector <2 x double> %58, <2 x double> %59, <2 x i32> <i32 0, i32 2>
  %61 = fadd <2 x double> %54, %60
  %62 = fmul <2 x double> %61, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %62, ptr getelementptr inbounds (i8, ptr @yc, i64 136), align 8, !tbaa !5
  store <2 x double> %62, ptr getelementptr inbounds (i8, ptr @xc, i64 136), align 8, !tbaa !5
  %63 = shufflevector <2 x double> %59, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %64 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 168), align 8, !tbaa !5
  %65 = shufflevector <2 x double> %63, <2 x double> %64, <2 x i32> <i32 0, i32 2>
  %66 = fadd <2 x double> %59, %65
  %67 = fmul <2 x double> %66, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %67, ptr getelementptr inbounds (i8, ptr @yc, i64 152), align 8, !tbaa !5
  store <2 x double> %67, ptr getelementptr inbounds (i8, ptr @xc, i64 152), align 8, !tbaa !5
  %68 = shufflevector <2 x double> %64, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %69 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 184), align 8, !tbaa !5
  %70 = shufflevector <2 x double> %68, <2 x double> %69, <2 x i32> <i32 0, i32 2>
  %71 = fadd <2 x double> %64, %70
  %72 = fmul <2 x double> %71, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %72, ptr getelementptr inbounds (i8, ptr @yc, i64 168), align 8, !tbaa !5
  store <2 x double> %72, ptr getelementptr inbounds (i8, ptr @xc, i64 168), align 8, !tbaa !5
  %73 = shufflevector <2 x double> %69, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %74 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 200), align 8, !tbaa !5
  %75 = shufflevector <2 x double> %73, <2 x double> %74, <2 x i32> <i32 0, i32 2>
  %76 = fadd <2 x double> %69, %75
  %77 = fmul <2 x double> %76, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %77, ptr getelementptr inbounds (i8, ptr @yc, i64 184), align 8, !tbaa !5
  store <2 x double> %77, ptr getelementptr inbounds (i8, ptr @xc, i64 184), align 8, !tbaa !5
  %78 = shufflevector <2 x double> %74, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %79 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 216), align 8, !tbaa !5
  %80 = shufflevector <2 x double> %78, <2 x double> %79, <2 x i32> <i32 0, i32 2>
  %81 = fadd <2 x double> %74, %80
  %82 = fmul <2 x double> %81, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %82, ptr getelementptr inbounds (i8, ptr @yc, i64 200), align 8, !tbaa !5
  store <2 x double> %82, ptr getelementptr inbounds (i8, ptr @xc, i64 200), align 8, !tbaa !5
  %83 = shufflevector <2 x double> %79, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %84 = load double, ptr getelementptr inbounds (i8, ptr @ti, i64 248), align 8, !tbaa !5
  %85 = load <2 x double>, ptr getelementptr inbounds (i8, ptr @ti, i64 232), align 8, !tbaa !5
  %86 = shufflevector <2 x double> %83, <2 x double> %85, <2 x i32> <i32 0, i32 2>
  %87 = fadd <2 x double> %79, %86
  %88 = fmul <2 x double> %87, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %88, ptr getelementptr inbounds (i8, ptr @yc, i64 216), align 8, !tbaa !5
  store <2 x double> %88, ptr getelementptr inbounds (i8, ptr @xc, i64 216), align 8, !tbaa !5
  %89 = shufflevector <2 x double> %85, <2 x double> poison, <2 x i32> <i32 1, i32 poison>
  %90 = insertelement <2 x double> %89, double %84, i64 1
  %91 = fadd <2 x double> %85, %90
  %92 = fmul <2 x double> %91, <double 5.000000e-01, double 5.000000e-01>
  store <2 x double> %92, ptr getelementptr inbounds (i8, ptr @yc, i64 232), align 8, !tbaa !5
  store <2 x double> %92, ptr getelementptr inbounds (i8, ptr @xc, i64 232), align 8, !tbaa !5
  %93 = load double, ptr getelementptr inbounds (i8, ptr @xc, i64 8), align 8, !tbaa !5
  %94 = fneg double %93
  store double %94, ptr @xc, align 16, !tbaa !5
  %95 = extractelement <2 x double> %92, i64 1
  %96 = fsub double 2.000000e+00, %95
  store double %96, ptr getelementptr inbounds (i8, ptr @xc, i64 248), align 8, !tbaa !5
  %97 = load double, ptr getelementptr inbounds (i8, ptr @yc, i64 8), align 8, !tbaa !5
  %98 = fneg double %97
  store double %98, ptr @yc, align 16, !tbaa !5
  store double %96, ptr getelementptr inbounds (i8, ptr @yc, i64 248), align 8, !tbaa !5
  %99 = fadd double %93, %93
  store double %99, ptr @dx, align 8, !tbaa !5
  store double %99, ptr @dy, align 8, !tbaa !5
  %100 = fmul double %99, 1.000000e-01
  %101 = load double, ptr @RT, align 8, !tbaa !5
  %102 = fmul double %101, 6.000000e+00
  %103 = tail call double @sqrt(double noundef %102) #5
  %104 = fdiv double %100, %103
  store double %104, ptr @dt, align 8, !tbaa !5
  %105 = load double, ptr @nu, align 8, !tbaa !5
  %106 = load double, ptr @RT, align 8, !tbaa !5
  %107 = fdiv double %105, %106
  store double %107, ptr @tau, align 8, !tbaa !5
  %108 = fmul double %104, 1.500000e+00
  %109 = tail call double @llvm.fmuladd.f64(double %107, double 2.000000e+00, double %104)
  %110 = fdiv double %108, %109
  store double %110, ptr @w, align 8, !tbaa !5
  %111 = fmul double %104, 5.000000e-01
  %112 = tail call double @llvm.fmuladd.f64(double %107, double 2.000000e+00, double %111)
  %113 = fdiv double %104, %112
  store double %113, ptr @wb, align 8, !tbaa !5
  %114 = tail call i32 (ptr, ...) @printf(ptr noundef nonnull dereferenceable(1) @.str, double noundef %110)
  tail call void @gks_ini()
  %115 = load double, ptr getelementptr inbounds (i8, ptr @ux, i64 4224), align 16, !tbaa !5
  br label %116

116:                                              ; preds = %._crit_edge, %.preheader.preheader
  %.025 = phi double [ %115, %.preheader.preheader ], [ %.126.lcssa, %._crit_edge ]
  %.3 = phi i32 [ 0, %.preheader.preheader ], [ %.4.lcssa, %._crit_edge ]
  %117 = call i32 (ptr, ...) @printf(ptr noundef nonnull dereferenceable(1) @.str.1)
  %118 = call i32 (ptr, ...) @__isoc99_scanf(ptr noundef nonnull @.str.2, ptr noundef nonnull %2)
  %119 = load i32, ptr %2, align 4, !tbaa !11
  %120 = add nsw i32 %119, %.3
  store i32 %120, ptr %2, align 4, !tbaa !11
  %121 = load double, ptr @dt, align 8, !tbaa !5
  %122 = call i32 (ptr, ...) @printf(ptr noundef nonnull dereferenceable(1) @.str.3, double noundef %121, i32 noundef %120)
  %123 = load i32, ptr %2, align 4, !tbaa !11
  %124 = icmp slt i32 %.3, %123
  br i1 %124, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %116, %134
  %.435 = phi i32 [ %125, %134 ], [ %.3, %116 ]
  %.12634 = phi double [ %.227, %134 ], [ %.025, %116 ]
  %.02833 = phi double [ %.129, %134 ], [ 1.000000e+00, %116 ]
  %125 = add nsw i32 %.435, 1
  call void @Evol()
  %126 = srem i32 %125, 1000
  %127 = icmp eq i32 %126, 0
  br i1 %127, label %128, label %134

128:                                              ; preds = %.lr.ph
  %129 = load double, ptr getelementptr inbounds (i8, ptr @ux, i64 4224), align 16, !tbaa !5
  %130 = fsub double %129, %.12634
  %131 = call double @llvm.fabs.f64(double %130)
  %132 = fdiv double %131, 1.000000e-01
  %133 = call i32 (ptr, ...) @printf(ptr noundef nonnull dereferenceable(1) @.str.4, double noundef %132, double noundef %129, i32 noundef %125)
  br label %134

134:                                              ; preds = %128, %.lr.ph
  %.129 = phi double [ %132, %128 ], [ %.02833, %.lr.ph ]
  %.227 = phi double [ %129, %128 ], [ %.12634, %.lr.ph ]
  %135 = load i32, ptr %2, align 4, !tbaa !11
  %136 = icmp slt i32 %125, %135
  %137 = fcmp ogt double %.129, 0x3EB0C6F7A0B5ED8D
  %138 = select i1 %136, i1 %137, i1 false
  br i1 %138, label %.lr.ph, label %._crit_edge, !llvm.loop !13

._crit_edge:                                      ; preds = %134, %116
  %.126.lcssa = phi double [ %.025, %116 ], [ %.227, %134 ]
  %.4.lcssa = phi i32 [ %.3, %116 ], [ %125, %134 ]
  call void @datadeal()
  %139 = call i32 (ptr, ...) @printf(ptr noundef nonnull dereferenceable(1) @.str.5)
  %140 = call i32 (ptr, ...) @__isoc99_scanf(ptr noundef nonnull @.str.2, ptr noundef nonnull %1)
  %141 = load i32, ptr %1, align 4, !tbaa !11
  %.not = icmp eq i32 %141, 0
  br i1 %.not, label %142, label %116

142:                                              ; preds = %._crit_edge
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %2) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %1) #5
  ret void
}

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

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

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

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

; Function Attrs: nounwind uwtable
declare dso_local void @gks_ini() local_unnamed_addr #0

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

; Function Attrs: nounwind uwtable
declare dso_local void @Evol() local_unnamed_addr #0

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

; Function Attrs: nounwind uwtable
declare dso_local void @datadeal() local_unnamed_addr #0

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

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 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 #3 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #4 = { nofree nounwind "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #5 = { nounwind }

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"double", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = distinct !{!9, !10}
!10 = !{!"llvm.loop.mustprogress"}
!11 = !{!12, !12, i64 0}
!12 = !{!"int", !7, i64 0}
!13 = distinct !{!13, !10}


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