; LLVM IR from rdguez-mariano/fast_imas_IPOL :: libUSAC/utils/HomographyFunctions.ll (ZN6HTools17computeDataMatrixEPdjS0_)

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

; Function Attrs: mustprogress nounwind uwtable
define dso_local void @_ZN6HTools17computeDataMatrixEPdjS0_(ptr noundef %0, i32 noundef %1, ptr noundef %2) #0 {
  %4 = alloca ptr, align 8
  %5 = alloca i32, align 4
  %6 = alloca ptr, align 8
  %7 = alloca ptr, align 8
  %8 = alloca ptr, align 8
  %9 = alloca i32, align 4
  %10 = alloca i32, align 4
  store ptr %0, ptr %4, align 8, !tbaa !5
  store i32 %1, ptr %5, align 4, !tbaa !9
  store ptr %2, ptr %6, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %7) #2
  call void @llvm.lifetime.start.p0(i64 8, ptr %8) #2
  call void @llvm.lifetime.start.p0(i64 4, ptr %9) #2
  %11 = load i32, ptr %5, align 4, !tbaa !9
  %12 = mul i32 2, %11
  store i32 %12, ptr %9, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 4, ptr %10) #2
  store i32 0, ptr %10, align 4, !tbaa !9
  br label %13

13:                                               ; preds = %181, %3
  %14 = load i32, ptr %10, align 4, !tbaa !9
  %15 = load i32, ptr %5, align 4, !tbaa !9
  %16 = icmp ult i32 %14, %15
  br i1 %16, label %18, label %17

17:                                               ; preds = %13
  call void @llvm.lifetime.end.p0(i64 4, ptr %10) #2
  br label %184

18:                                               ; preds = %13
  %19 = load ptr, ptr %6, align 8, !tbaa !5
  %20 = load i32, ptr %10, align 4, !tbaa !9
  %21 = mul i32 6, %20
  %22 = zext i32 %21 to i64
  %23 = getelementptr inbounds double, ptr %19, i64 %22
  store ptr %23, ptr %7, align 8, !tbaa !5
  %24 = load ptr, ptr %4, align 8, !tbaa !5
  %25 = load i32, ptr %10, align 4, !tbaa !9
  %26 = mul i32 2, %25
  %27 = zext i32 %26 to i64
  %28 = getelementptr inbounds double, ptr %24, i64 %27
  store ptr %28, ptr %8, align 8, !tbaa !5
  %29 = load ptr, ptr %8, align 8, !tbaa !5
  store double 0.000000e+00, ptr %29, align 8, !tbaa !11
  %30 = load ptr, ptr %8, align 8, !tbaa !5
  %31 = load i32, ptr %9, align 4, !tbaa !9
  %32 = zext i32 %31 to i64
  %33 = getelementptr inbounds double, ptr %30, i64 %32
  store double 0.000000e+00, ptr %33, align 8, !tbaa !11
  %34 = load ptr, ptr %8, align 8, !tbaa !5
  %35 = load i32, ptr %9, align 4, !tbaa !9
  %36 = mul i32 2, %35
  %37 = zext i32 %36 to i64
  %38 = getelementptr inbounds double, ptr %34, i64 %37
  store double 0.000000e+00, ptr %38, align 8, !tbaa !11
  %39 = load ptr, ptr %7, align 8, !tbaa !5
  %40 = getelementptr inbounds double, ptr %39, i64 0
  %41 = load double, ptr %40, align 8, !tbaa !11
  %42 = fneg double %41
  %43 = load ptr, ptr %8, align 8, !tbaa !5
  %44 = load i32, ptr %9, align 4, !tbaa !9
  %45 = mul i32 3, %44
  %46 = zext i32 %45 to i64
  %47 = getelementptr inbounds double, ptr %43, i64 %46
  store double %42, ptr %47, align 8, !tbaa !11
  %48 = load ptr, ptr %7, align 8, !tbaa !5
  %49 = getelementptr inbounds double, ptr %48, i64 1
  %50 = load double, ptr %49, align 8, !tbaa !11
  %51 = fneg double %50
  %52 = load ptr, ptr %8, align 8, !tbaa !5
  %53 = load i32, ptr %9, align 4, !tbaa !9
  %54 = mul i32 4, %53
  %55 = zext i32 %54 to i64
  %56 = getelementptr inbounds double, ptr %52, i64 %55
  store double %51, ptr %56, align 8, !tbaa !11
  %57 = load ptr, ptr %7, align 8, !tbaa !5
  %58 = getelementptr inbounds double, ptr %57, i64 2
  %59 = load double, ptr %58, align 8, !tbaa !11
  %60 = fneg double %59
  %61 = load ptr, ptr %8, align 8, !tbaa !5
  %62 = load i32, ptr %9, align 4, !tbaa !9
  %63 = mul i32 5, %62
  %64 = zext i32 %63 to i64
  %65 = getelementptr inbounds double, ptr %61, i64 %64
  store double %60, ptr %65, align 8, !tbaa !11
  %66 = load ptr, ptr %7, align 8, !tbaa !5
  %67 = getelementptr inbounds double, ptr %66, i64 4
  %68 = load double, ptr %67, align 8, !tbaa !11
  %69 = load ptr, ptr %7, align 8, !tbaa !5
  %70 = getelementptr inbounds double, ptr %69, i64 0
  %71 = load double, ptr %70, align 8, !tbaa !11
  %72 = fmul double %68, %71
  %73 = load ptr, ptr %8, align 8, !tbaa !5
  %74 = load i32, ptr %9, align 4, !tbaa !9
  %75 = mul i32 6, %74
  %76 = zext i32 %75 to i64
  %77 = getelementptr inbounds double, ptr %73, i64 %76
  store double %72, ptr %77, align 8, !tbaa !11
  %78 = load ptr, ptr %7, align 8, !tbaa !5
  %79 = getelementptr inbounds double, ptr %78, i64 4
  %80 = load double, ptr %79, align 8, !tbaa !11
  %81 = load ptr, ptr %7, align 8, !tbaa !5
  %82 = getelementptr inbounds double, ptr %81, i64 1
  %83 = load double, ptr %82, align 8, !tbaa !11
  %84 = fmul double %80, %83
  %85 = load ptr, ptr %8, align 8, !tbaa !5
  %86 = load i32, ptr %9, align 4, !tbaa !9
  %87 = mul i32 7, %86
  %88 = zext i32 %87 to i64
  %89 = getelementptr inbounds double, ptr %85, i64 %88
  store double %84, ptr %89, align 8, !tbaa !11
  %90 = load ptr, ptr %7, align 8, !tbaa !5
  %91 = getelementptr inbounds double, ptr %90, i64 4
  %92 = load double, ptr %91, align 8, !tbaa !11
  %93 = load ptr, ptr %7, align 8, !tbaa !5
  %94 = getelementptr inbounds double, ptr %93, i64 2
  %95 = load double, ptr %94, align 8, !tbaa !11
  %96 = fmul double %92, %95
  %97 = load ptr, ptr %8, align 8, !tbaa !5
  %98 = load i32, ptr %9, align 4, !tbaa !9
  %99 = mul i32 8, %98
  %100 = zext i32 %99 to i64
  %101 = getelementptr inbounds double, ptr %97, i64 %100
  store double %96, ptr %101, align 8, !tbaa !11
  %102 = load ptr, ptr %4, align 8, !tbaa !5
  %103 = load i32, ptr %10, align 4, !tbaa !9
  %104 = mul i32 2, %103
  %105 = zext i32 %104 to i64
  %106 = getelementptr inbounds double, ptr %102, i64 %105
  %107 = getelementptr inbounds double, ptr %106, i64 1
  store ptr %107, ptr %8, align 8, !tbaa !5
  %108 = load ptr, ptr %7, align 8, !tbaa !5
  %109 = getelementptr inbounds double, ptr %108, i64 0
  %110 = load double, ptr %109, align 8, !tbaa !11
  %111 = load ptr, ptr %8, align 8, !tbaa !5
  store double %110, ptr %111, align 8, !tbaa !11
  %112 = load ptr, ptr %7, align 8, !tbaa !5
  %113 = getelementptr inbounds double, ptr %112, i64 1
  %114 = load double, ptr %113, align 8, !tbaa !11
  %115 = load ptr, ptr %8, align 8, !tbaa !5
  %116 = load i32, ptr %9, align 4, !tbaa !9
  %117 = zext i32 %116 to i64
  %118 = getelementptr inbounds double, ptr %115, i64 %117
  store double %114, ptr %118, align 8, !tbaa !11
  %119 = load ptr, ptr %7, align 8, !tbaa !5
  %120 = getelementptr inbounds double, ptr %119, i64 2
  %121 = load double, ptr %120, align 8, !tbaa !11
  %122 = load ptr, ptr %8, align 8, !tbaa !5
  %123 = load i32, ptr %9, align 4, !tbaa !9
  %124 = mul i32 2, %123
  %125 = zext i32 %124 to i64
  %126 = getelementptr inbounds double, ptr %122, i64 %125
  store double %121, ptr %126, align 8, !tbaa !11
  %127 = load ptr, ptr %8, align 8, !tbaa !5
  %128 = load i32, ptr %9, align 4, !tbaa !9
  %129 = mul i32 3, %128
  %130 = zext i32 %129 to i64
  %131 = getelementptr inbounds double, ptr %127, i64 %130
  store double 0.000000e+00, ptr %131, align 8, !tbaa !11
  %132 = load ptr, ptr %8, align 8, !tbaa !5
  %133 = load i32, ptr %9, align 4, !tbaa !9
  %134 = mul i32 4, %133
  %135 = zext i32 %134 to i64
  %136 = getelementptr inbounds double, ptr %132, i64 %135
  store double 0.000000e+00, ptr %136, align 8, !tbaa !11
  %137 = load ptr, ptr %8, align 8, !tbaa !5
  %138 = load i32, ptr %9, align 4, !tbaa !9
  %139 = mul i32 5, %138
  %140 = zext i32 %139 to i64
  %141 = getelementptr inbounds double, ptr %137, i64 %140
  store double 0.000000e+00, ptr %141, align 8, !tbaa !11
  %142 = load ptr, ptr %7, align 8, !tbaa !5
  %143 = getelementptr inbounds double, ptr %142, i64 3
  %144 = load double, ptr %143, align 8, !tbaa !11
  %145 = fneg double %144
  %146 = load ptr, ptr %7, align 8, !tbaa !5
  %147 = getelementptr inbounds double, ptr %146, i64 0
  %148 = load double, ptr %147, align 8, !tbaa !11
  %149 = fmul double %145, %148
  %150 = load ptr, ptr %8, align 8, !tbaa !5
  %151 = load i32, ptr %9, align 4, !tbaa !9
  %152 = mul i32 6, %151
  %153 = zext i32 %152 to i64
  %154 = getelementptr inbounds double, ptr %150, i64 %153
  store double %149, ptr %154, align 8, !tbaa !11
  %155 = load ptr, ptr %7, align 8, !tbaa !5
  %156 = getelementptr inbounds double, ptr %155, i64 3
  %157 = load double, ptr %156, align 8, !tbaa !11
  %158 = fneg double %157
  %159 = load ptr, ptr %7, align 8, !tbaa !5
  %160 = getelementptr inbounds double, ptr %159, i64 1
  %161 = load double, ptr %160, align 8, !tbaa !11
  %162 = fmul double %158, %161
  %163 = load ptr, ptr %8, align 8, !tbaa !5
  %164 = load i32, ptr %9, align 4, !tbaa !9
  %165 = mul i32 7, %164
  %166 = zext i32 %165 to i64
  %167 = getelementptr inbounds double, ptr %163, i64 %166
  store double %162, ptr %167, align 8, !tbaa !11
  %168 = load ptr, ptr %7, align 8, !tbaa !5
  %169 = getelementptr inbounds double, ptr %168, i64 3
  %170 = load double, ptr %169, align 8, !tbaa !11
  %171 = fneg double %170
  %172 = load ptr, ptr %7, align 8, !tbaa !5
  %173 = getelementptr inbounds double, ptr %172, i64 2
  %174 = load double, ptr %173, align 8, !tbaa !11
  %175 = fmul double %171, %174
  %176 = load ptr, ptr %8, align 8, !tbaa !5
  %177 = load i32, ptr %9, align 4, !tbaa !9
  %178 = mul i32 8, %177
  %179 = zext i32 %178 to i64
  %180 = getelementptr inbounds double, ptr %176, i64 %179
  store double %175, ptr %180, align 8, !tbaa !11
  br label %181

181:                                              ; preds = %18
  %182 = load i32, ptr %10, align 4, !tbaa !9
  %183 = add i32 %182, 1
  store i32 %183, ptr %10, align 4, !tbaa !9
  br label %13, !llvm.loop !13

184:                                              ; preds = %17
  call void @llvm.lifetime.end.p0(i64 4, ptr %9) #2
  call void @llvm.lifetime.end.p0(i64 8, ptr %8) #2
  call void @llvm.lifetime.end.p0(i64 8, ptr %7) #2
  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: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.end.p0(i64 immarg, ptr nocapture) #1

attributes #0 = { mustprogress 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 }

!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++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"int", !7, i64 0}
!11 = !{!12, !12, i64 0}
!12 = !{!"double", !7, i64 0}
!13 = distinct !{!13, !14}
!14 = !{!"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"

; Function Attrs: mustprogress nofree norecurse nosync nounwind memory(argmem: readwrite) uwtable
define dso_local void @_ZN6HTools17computeDataMatrixEPdjS0_(ptr nocapture noundef writeonly %0, i32 noundef %1, ptr nocapture noundef readonly %2) local_unnamed_addr #0 {
  %.not = icmp eq i32 %1, 0
  br i1 %.not, label %._crit_edge, label %.lr.ph

.lr.ph:                                           ; preds = %3
  %4 = shl i32 %1, 1
  %5 = zext i32 %4 to i64
  %6 = shl i32 %1, 2
  %7 = zext i32 %6 to i64
  %8 = mul i32 %1, 6
  %9 = zext i32 %8 to i64
  %10 = shl i32 %1, 3
  %11 = zext i32 %10 to i64
  %12 = mul i32 %1, 10
  %13 = zext i32 %12 to i64
  %14 = mul i32 %1, 12
  %15 = zext i32 %14 to i64
  %16 = mul i32 %1, 14
  %17 = zext i32 %16 to i64
  %18 = shl i32 %1, 4
  %19 = zext i32 %18 to i64
  %wide.trip.count = zext i32 %1 to i64
  br label %20

._crit_edge:                                      ; preds = %20, %3
  ret void

20:                                               ; preds = %.lr.ph, %20
  %indvars.iv = phi i64 [ 0, %.lr.ph ], [ %indvars.iv.next, %20 ]
  %21 = mul i64 %indvars.iv, 6
  %22 = and i64 %21, 4294967294
  %23 = getelementptr inbounds double, ptr %2, i64 %22
  %24 = shl nuw i64 %indvars.iv, 1
  %25 = and i64 %24, 4294967294
  %26 = getelementptr inbounds double, ptr %0, i64 %25
  store double 0.000000e+00, ptr %26, align 8, !tbaa !5
  %27 = getelementptr inbounds double, ptr %26, i64 %5
  store double 0.000000e+00, ptr %27, align 8, !tbaa !5
  %28 = getelementptr inbounds double, ptr %26, i64 %7
  store double 0.000000e+00, ptr %28, align 8, !tbaa !5
  %29 = load double, ptr %23, align 8, !tbaa !5
  %30 = fneg double %29
  %31 = getelementptr inbounds double, ptr %26, i64 %9
  store double %30, ptr %31, align 8, !tbaa !5
  %32 = getelementptr inbounds i8, ptr %23, i64 8
  %33 = load double, ptr %32, align 8, !tbaa !5
  %34 = fneg double %33
  %35 = getelementptr inbounds double, ptr %26, i64 %11
  store double %34, ptr %35, align 8, !tbaa !5
  %36 = getelementptr inbounds i8, ptr %23, i64 16
  %37 = load double, ptr %36, align 8, !tbaa !5
  %38 = fneg double %37
  %39 = getelementptr inbounds double, ptr %26, i64 %13
  store double %38, ptr %39, align 8, !tbaa !5
  %40 = getelementptr inbounds i8, ptr %23, i64 32
  %41 = load double, ptr %40, align 8, !tbaa !5
  %42 = load double, ptr %23, align 8, !tbaa !5
  %43 = fmul double %41, %42
  %44 = getelementptr inbounds double, ptr %26, i64 %15
  store double %43, ptr %44, align 8, !tbaa !5
  %45 = load double, ptr %40, align 8, !tbaa !5
  %46 = load double, ptr %32, align 8, !tbaa !5
  %47 = fmul double %45, %46
  %48 = getelementptr inbounds double, ptr %26, i64 %17
  store double %47, ptr %48, align 8, !tbaa !5
  %49 = load double, ptr %40, align 8, !tbaa !5
  %50 = load double, ptr %36, align 8, !tbaa !5
  %51 = fmul double %49, %50
  %52 = getelementptr inbounds double, ptr %26, i64 %19
  store double %51, ptr %52, align 8, !tbaa !5
  %53 = getelementptr inbounds i8, ptr %26, i64 8
  %54 = load double, ptr %23, align 8, !tbaa !5
  store double %54, ptr %53, align 8, !tbaa !5
  %55 = load double, ptr %32, align 8, !tbaa !5
  %56 = getelementptr inbounds double, ptr %53, i64 %5
  store double %55, ptr %56, align 8, !tbaa !5
  %57 = load double, ptr %36, align 8, !tbaa !5
  %58 = getelementptr inbounds double, ptr %53, i64 %7
  store double %57, ptr %58, align 8, !tbaa !5
  %59 = getelementptr inbounds double, ptr %53, i64 %9
  store double 0.000000e+00, ptr %59, align 8, !tbaa !5
  %60 = getelementptr inbounds double, ptr %53, i64 %11
  store double 0.000000e+00, ptr %60, align 8, !tbaa !5
  %61 = getelementptr inbounds double, ptr %53, i64 %13
  store double 0.000000e+00, ptr %61, align 8, !tbaa !5
  %62 = getelementptr inbounds i8, ptr %23, i64 24
  %63 = load double, ptr %62, align 8, !tbaa !5
  %64 = fneg double %63
  %65 = load double, ptr %23, align 8, !tbaa !5
  %66 = fmul double %65, %64
  %67 = getelementptr inbounds double, ptr %53, i64 %15
  store double %66, ptr %67, align 8, !tbaa !5
  %68 = load double, ptr %62, align 8, !tbaa !5
  %69 = fneg double %68
  %70 = load double, ptr %32, align 8, !tbaa !5
  %71 = fmul double %70, %69
  %72 = getelementptr inbounds double, ptr %53, i64 %17
  store double %71, ptr %72, align 8, !tbaa !5
  %73 = load double, ptr %62, align 8, !tbaa !5
  %74 = fneg double %73
  %75 = load double, ptr %36, align 8, !tbaa !5
  %76 = fmul double %75, %74
  %77 = getelementptr inbounds double, ptr %53, i64 %19
  store double %76, ptr %77, align 8, !tbaa !5
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next, %wide.trip.count
  br i1 %exitcond.not, label %._crit_edge, label %20, !llvm.loop !9
}

attributes #0 = { mustprogress nofree norecurse nosync nounwind memory(argmem: readwrite) 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" }

!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++ TBAA"}
!9 = distinct !{!9, !10}
!10 = !{!"llvm.loop.mustprogress"}


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