; LLVM IR from jmbr/libpme6 :: linear-algebra.ll (ZN5ewald14matrix_inverseEPKdPd)

; 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: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.start.p0(i64 immarg, ptr nocapture) #0

; Function Attrs: nocallback nofree nounwind willreturn memory(argmem: write)
declare void @llvm.memset.p0.i64(ptr nocapture writeonly, i8, i64, i1 immarg) #1

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

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

; Function Attrs: mustprogress nounwind uwtable
declare dso_local noundef double @_ZN5ewald11determinantEPKd(ptr noundef) #3

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

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nocallback nofree nounwind willreturn memory(argmem: write) }
attributes #2 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #3 = { 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 #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 = !{!"any pointer", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"double", !7, i64 0}


; 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 nocallback nofree nounwind willreturn memory(argmem: write)
declare void @llvm.memset.p0.i64(ptr nocapture writeonly, i8, i64, i1 immarg) #0

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

; Function Attrs: mustprogress nounwind uwtable
declare dso_local noundef double @_ZN5ewald11determinantEPKd(ptr noundef) local_unnamed_addr #2

; Function Attrs: mustprogress nounwind uwtable
define dso_local void @_ZN5ewald14matrix_inverseEPKdPd(ptr noundef %0, ptr nocapture noundef writeonly %1) local_unnamed_addr #2 {
  %3 = tail call noundef double @_ZN5ewald11determinantEPKd(ptr noundef %0)
  %4 = fdiv double 1.000000e+00, %3
  tail call void @llvm.memset.p0.i64(ptr noundef nonnull align 8 dereferenceable(72) %1, i8 0, i64 72, i1 false)
  %5 = getelementptr inbounds i8, ptr %0, i64 32
  %6 = load double, ptr %5, align 8, !tbaa !5
  %7 = getelementptr inbounds i8, ptr %0, i64 64
  %8 = load double, ptr %7, align 8, !tbaa !5
  %9 = getelementptr inbounds i8, ptr %0, i64 40
  %10 = load double, ptr %9, align 8, !tbaa !5
  %11 = getelementptr inbounds i8, ptr %0, i64 56
  %12 = load double, ptr %11, align 8, !tbaa !5
  %13 = fneg double %10
  %14 = fmul double %12, %13
  %15 = tail call double @llvm.fmuladd.f64(double %6, double %8, double %14)
  %16 = fmul double %4, %15
  store double %16, ptr %1, align 8, !tbaa !5
  %17 = getelementptr inbounds i8, ptr %0, i64 8
  %18 = load double, ptr %17, align 8, !tbaa !5
  %19 = fneg double %18
  %20 = load double, ptr %7, align 8, !tbaa !5
  %21 = getelementptr inbounds i8, ptr %0, i64 16
  %22 = load double, ptr %21, align 8, !tbaa !5
  %23 = load double, ptr %11, align 8, !tbaa !5
  %24 = fmul double %22, %23
  %25 = tail call double @llvm.fmuladd.f64(double %19, double %20, double %24)
  %26 = fmul double %4, %25
  %27 = getelementptr inbounds i8, ptr %1, i64 8
  store double %26, ptr %27, align 8, !tbaa !5
  %28 = load double, ptr %17, align 8, !tbaa !5
  %29 = load double, ptr %9, align 8, !tbaa !5
  %30 = load double, ptr %21, align 8, !tbaa !5
  %31 = load double, ptr %5, align 8, !tbaa !5
  %32 = fneg double %30
  %33 = fmul double %31, %32
  %34 = tail call double @llvm.fmuladd.f64(double %28, double %29, double %33)
  %35 = fmul double %4, %34
  %36 = getelementptr inbounds i8, ptr %1, i64 16
  store double %35, ptr %36, align 8, !tbaa !5
  %37 = getelementptr inbounds i8, ptr %0, i64 24
  %38 = load double, ptr %37, align 8, !tbaa !5
  %39 = fneg double %38
  %40 = load double, ptr %7, align 8, !tbaa !5
  %41 = load double, ptr %9, align 8, !tbaa !5
  %42 = getelementptr inbounds i8, ptr %0, i64 48
  %43 = load double, ptr %42, align 8, !tbaa !5
  %44 = fmul double %41, %43
  %45 = tail call double @llvm.fmuladd.f64(double %39, double %40, double %44)
  %46 = fmul double %4, %45
  %47 = getelementptr inbounds i8, ptr %1, i64 24
  store double %46, ptr %47, align 8, !tbaa !5
  %48 = load double, ptr %0, align 8, !tbaa !5
  %49 = load double, ptr %7, align 8, !tbaa !5
  %50 = load double, ptr %21, align 8, !tbaa !5
  %51 = load double, ptr %42, align 8, !tbaa !5
  %52 = fneg double %50
  %53 = fmul double %51, %52
  %54 = tail call double @llvm.fmuladd.f64(double %48, double %49, double %53)
  %55 = fmul double %4, %54
  %56 = getelementptr inbounds i8, ptr %1, i64 32
  store double %55, ptr %56, align 8, !tbaa !5
  %57 = load double, ptr %0, align 8, !tbaa !5
  %58 = fneg double %57
  %59 = load double, ptr %9, align 8, !tbaa !5
  %60 = load double, ptr %21, align 8, !tbaa !5
  %61 = load double, ptr %37, align 8, !tbaa !5
  %62 = fmul double %60, %61
  %63 = tail call double @llvm.fmuladd.f64(double %58, double %59, double %62)
  %64 = fmul double %4, %63
  %65 = getelementptr inbounds i8, ptr %1, i64 40
  store double %64, ptr %65, align 8, !tbaa !5
  %66 = load double, ptr %37, align 8, !tbaa !5
  %67 = load double, ptr %11, align 8, !tbaa !5
  %68 = load double, ptr %5, align 8, !tbaa !5
  %69 = load double, ptr %42, align 8, !tbaa !5
  %70 = fneg double %68
  %71 = fmul double %69, %70
  %72 = tail call double @llvm.fmuladd.f64(double %66, double %67, double %71)
  %73 = fmul double %4, %72
  %74 = getelementptr inbounds i8, ptr %1, i64 48
  store double %73, ptr %74, align 8, !tbaa !5
  %75 = load double, ptr %0, align 8, !tbaa !5
  %76 = fneg double %75
  %77 = load double, ptr %11, align 8, !tbaa !5
  %78 = load double, ptr %17, align 8, !tbaa !5
  %79 = load double, ptr %42, align 8, !tbaa !5
  %80 = fmul double %78, %79
  %81 = tail call double @llvm.fmuladd.f64(double %76, double %77, double %80)
  %82 = fmul double %4, %81
  %83 = getelementptr inbounds i8, ptr %1, i64 56
  store double %82, ptr %83, align 8, !tbaa !5
  %84 = load double, ptr %0, align 8, !tbaa !5
  %85 = load double, ptr %5, align 8, !tbaa !5
  %86 = load double, ptr %17, align 8, !tbaa !5
  %87 = load double, ptr %37, align 8, !tbaa !5
  %88 = fneg double %86
  %89 = fmul double %87, %88
  %90 = tail call double @llvm.fmuladd.f64(double %84, double %85, double %89)
  %91 = fmul double %4, %90
  %92 = getelementptr inbounds i8, ptr %1, i64 64
  store double %91, ptr %92, align 8, !tbaa !5
  ret void
}

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

!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"}


; Effective O3 passes: EarlyCSEPass,GlobalOptPass,InferFunctionAttrsPass,InstCombinePass,PostOrderFunctionAttrsPass,ReassociatePass,SROAPass,TailCallElimPass