void inv_hat(int n, double hat_hat[], double vec[], double vec1[],
             double eigs[], double lam0, double lam) {
  int i, j, k;
  double tmp[n][n], sum;
  for (i = 0; i < n; i++) {
    for (j = i; j < n; j++) {
      for (sum = 0, k = 0; k < n; k++)
        sum += hat_hat[k + i * n] * hat_hat[k + j * n] / (eigs[k] + lam - lam0);
      tmp[i][j] = sum;
    }
  }
  for (i = 0; i < n; i++) {
    for (j = i; j < n; j++)
      hat_hat[j + i * n] = hat_hat[i + j * n] = tmp[i][j];
    for (sum = 0, j = 0; j < n; j++)
      sum += hat_hat[j + i * n] * vec[j];
    vec1[i] = sum;
  }
}


x86 assembly (clang, -O0):
<inv_hat>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,%ecx
mov    %rcx,%rax
mov    %rax,-0x78(%rbp)
mov    %rsp,%rdx
mov    %rdx,-0x50(%rbp)
mov    %rcx,%rdx
imul   %rdx,%rdx
lea    0xf(,%rdx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x70(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x58(%rbp)
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    124c <inv_hat+0x14c>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1239 <inv_hat+0x139>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1206 <inv_hat+0x106>
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x40(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
addsd  -0x38(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
divsd  %xmm1,%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1197 <inv_hat+0x97>
mov    -0x70(%rbp),%rax
mov    -0x78(%rbp),%rdx
movsd  -0x68(%rbp),%xmm0
movslq -0x3c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    117c <inv_hat+0x7c>
jmp    123e <inv_hat+0x13e>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    116a <inv_hat+0x6a>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    134f <inv_hat+0x24f>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12cf <inv_hat+0x1cf>
mov    -0x70(%rbp),%rax
mov    -0x78(%rbp),%rdx
movslq -0x3c(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x40(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1265 <inv_hat+0x165>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    132f <inv_hat+0x22f>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    12de <inv_hat+0x1de>
movsd  -0x68(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1253 <inv_hat+0x153>
mov    -0x50(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<inv_hat>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rcx,-0x48(%rbp)
mov    %rsi,-0x30(%rbp)
mov    %edi,%r12d
mov    %r12,%rax
imul   %r12,%rax
mov    %rsp,%r9
lea    0xf(,%rax,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%r9
mov    %r9,%rsp
test   %edi,%edi
jle    130c <inv_hat+0x20c>
lea    0x0(,%r12,8),%r10
xor    %r11d,%r11d
mov    -0x30(%rbp),%rcx
nop
mov    %r11,%rax
imul   %r12,%rax
lea    (%r9,%rax,8),%r14
mov    %rcx,%rbx
mov    %r11,%r15
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorpd  %xmm2,%xmm2
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
movsd  (%rcx,%rax,8),%xmm3
mulsd  (%rbx,%rax,8),%xmm3
movsd  (%r8,%rax,8),%xmm4
addsd  %xmm1,%xmm4
subsd  %xmm0,%xmm4
divsd  %xmm4,%xmm3
addsd  %xmm3,%xmm2
add    $0x1,%rax
cmp    %rax,%r12
jne    1180 <inv_hat+0x80>
movsd  %xmm2,(%r14,%r15,8)
add    $0x1,%r15
add    %r10,%rbx
cmp    %r12,%r15
jne    1170 <inv_hat+0x70>
add    $0x1,%r11
add    %r10,%rcx
cmp    %r12,%r11
jne    1150 <inv_hat+0x50>
test   %edi,%edi
jle    130c <inv_hat+0x20c>
lea    -0x1(%r12),%rax
mov    %rax,-0x40(%rbp)
mov    %r12d,%r14d
and    $0x3,%r14d
mov    %r12d,%r11d
and    $0xfffffffc,%r11d
lea    0x8(,%r12,8),%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rdi
lea    0x18(%rdi),%r15
xor    %esi,%esi
mov    %r12,%r13
jmp    123b <inv_hat+0x13b>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    -0x48(%rbp),%rax
movsd  %xmm0,(%rax,%rsi,8)
add    $0x1,%rsi
mov    -0x38(%rbp),%rax
add    %rax,%rdi
add    %rax,%r9
add    $0xffffffffffffffff,%r13
add    %r10,%r15
add    %r10,-0x30(%rbp)
cmp    %r12,%rsi
je     130c <inv_hat+0x20c>
xor    %eax,%eax
xor    %ecx,%ecx
nop
movsd  (%r9,%rcx,8),%xmm0
movsd  %xmm0,(%rdi,%rax,8)
movsd  %xmm0,(%rdi,%rcx,8)
add    $0x1,%rcx
add    %r12,%rax
cmp    %rcx,%r13
jne    1240 <inv_hat+0x140>
cmpq   $0x3,-0x40(%rbp)
jae    1270 <inv_hat+0x170>
xorpd  %xmm0,%xmm0
xor    %r8d,%r8d
jmp    12cf <inv_hat+0x1cf>
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
xor    %r8d,%r8d
nopw   0x0(%rax,%rax,1)
movsd  -0x18(%r15,%r8,8),%xmm1
movsd  -0x10(%r15,%r8,8),%xmm2
mulsd  (%rdx,%r8,8),%xmm1
mulsd  0x8(%rdx,%r8,8),%xmm2
addsd  %xmm0,%xmm1
movsd  -0x8(%r15,%r8,8),%xmm3
mulsd  0x10(%rdx,%r8,8),%xmm3
addsd  %xmm1,%xmm2
movsd  (%r15,%r8,8),%xmm0
mulsd  0x18(%rdx,%r8,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%r8
cmp    %r8,%r11
jne    1280 <inv_hat+0x180>
test   %r14,%r14
je     1210 <inv_hat+0x110>
lea    (%rdx,%r8,8),%rcx
mov    -0x30(%rbp),%rax
lea    (%rax,%r8,8),%rax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax,%rbx,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rbx
cmp    %rbx,%r14
jne    12f0 <inv_hat+0x1f0>
jmp    1210 <inv_hat+0x110>
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @inv_hat(i32 %0, double* nocapture %1, double* nocapture %2, double* nocapture writeonly %3, double* nocapture %4, double %5, double %6) {
%8 = zext i32 %0 to i64
%9 = mul nuw i64 %8, %8
%10 = alloca double, i64 %9
%11 = icmp sgt i32 %0, 0
br i1 %11, label %12, label %129
12: ; preds = %7
%13 = zext i32 %0 to i64
br label %14
14: ; preds = %52, %12
%15 = phi i64 [ 0, %12 ], [ %53, %52 ]
%16 = mul nsw i64 %15, %13
%17 = mul nuw nsw i64 %15, %8
%18 = getelementptr inbounds double, double* %10, i64 %17
br label %27
19: ; preds = %52
br i1 %11, label %20, label %129
20: ; preds = %19
%21 = zext i32 %0 to i64
%22 = add nsw i64 %8, -1
%23 = and i64 %8, 3
%24 = icmp ult i64 %22, 3
%25 = and i64 %8, 4294967292
%26 = icmp eq i64 %23, 0
br label %55
27: ; preds = %48, %14
%28 = phi i64 [ %15, %14 ], [ %50, %48 ]
%29 = mul nsw i64 %28, %13
br label %30
30: ; preds = %27, %30
%31 = phi i64 [ 0, %27 ], [ %46, %30 ]
%32 = phi double [ 0.000000e+00, %27 ], [ %45, %30 ]
%33 = add nuw nsw i64 %31, %16
%34 = getelementptr inbounds double, double* %1, i64 %33
%35 = load double, double* %34
%36 = add nuw nsw i64 %31, %29
%37 = getelementptr inbounds double, double* %1, i64 %36
%38 = load double, double* %37
%39 = fmul double %35, %38
%40 = getelementptr inbounds double, double* %4, i64 %31
%41 = load double, double* %40
%42 = fadd double %41, %6
%43 = fsub double %42, %5
%44 = fdiv double %39, %43
%45 = fadd double %32, %44
%46 = add nuw nsw i64 %31, 1
%47 = icmp eq i64 %46, %8
br i1 %47, label %48, label %30
48: ; preds = %30
%49 = getelementptr inbounds double, double* %18, i64 %28
store double %45, double* %49
%50 = add nuw nsw i64 %28, 1
%51 = icmp eq i64 %50, %8
br i1 %51, label %52, label %27
52: ; preds = %48
%53 = add nuw nsw i64 %15, 1
%54 = icmp eq i64 %53, %8
br i1 %54, label %19, label %14
55: ; preds = %124, %20
%56 = phi i64 [ 0, %20 ], [ %127, %124 ]
%57 = mul nuw nsw i64 %56, %8
%58 = getelementptr inbounds double, double* %10, i64 %57
%59 = mul nsw i64 %56, %21
br label %62
60: ; preds = %62
%61 = mul nsw i64 %56, %21
br i1 %24, label %107, label %73
62: ; preds = %55, %62
%63 = phi i64 [ %56, %55 ], [ %71, %62 ]
%64 = getelementptr inbounds double, double* %58, i64 %63
%65 = load double, double* %64
%66 = mul nsw i64 %63, %21
%67 = add nuw nsw i64 %66, %56
%68 = getelementptr inbounds double, double* %1, i64 %67
store double %65, double* %68
%69 = add nuw nsw i64 %63, %59
%70 = getelementptr inbounds double, double* %1, i64 %69
store double %65, double* %70
%71 = add nuw nsw i64 %63, 1
%72 = icmp eq i64 %71, %8
br i1 %72, label %60, label %62
73: ; preds = %60, %73
%74 = phi i64 [ %104, %73 ], [ 0, %60 ]
%75 = phi double [ %103, %73 ], [ 0.000000e+00, %60 ]
%76 = phi i64 [ %105, %73 ], [ 0, %60 ]
%77 = add nuw nsw i64 %74, %61
%78 = getelementptr inbounds double, double* %1, i64 %77
%79 = load double, double* %78
%80 = getelementptr inbounds double, double* %2, i64 %74
%81 = load double, double* %80
%82 = tail call double @llvm.fmuladd.f64(double %79, double %81, double %75)
%83 = or i64 %74, 1
%84 = add nuw nsw i64 %83, %61
%85 = getelementptr inbounds double, double* %1, i64 %84
%86 = load double, double* %85
%87 = getelementptr inbounds double, double* %2, i64 %83
%88 = load double, double* %87
%89 = tail call double @llvm.fmuladd.f64(double %86, double %88, double %82)
%90 = or i64 %74, 2
%91 = add nuw nsw i64 %90, %61
%92 = getelementptr inbounds double, double* %1, i64 %91
%93 = load double, double* %92
%94 = getelementptr inbounds double, double* %2, i64 %90
%95 = load double, double* %94
%96 = tail call double @llvm.fmuladd.f64(double %93, double %95, double %89)
%97 = or i64 %74, 3
%98 = add nuw nsw i64 %97, %61
%99 = getelementptr inbounds double, double* %1, i64 %98
%100 = load double, double* %99
%101 = getelementptr inbounds double, double* %2, i64 %97
%102 = load double, double* %101
%103 = tail call double @llvm.fmuladd.f64(double %100, double %102, double %96)
%104 = add nuw nsw i64 %74, 4
%105 = add i64 %76, 4
%106 = icmp eq i64 %105, %25
br i1 %106, label %107, label %73
107: ; preds = %73, %60
%108 = phi double [ undef, %60 ], [ %103, %73 ]
%109 = phi i64 [ 0, %60 ], [ %104, %73 ]
%110 = phi double [ 0.000000e+00, %60 ], [ %103, %73 ]
br i1 %26, label %124, label %111
111: ; preds = %107, %111
%112 = phi i64 [ %121, %111 ], [ %109, %107 ]
%113 = phi double [ %120, %111 ], [ %110, %107 ]
%114 = phi i64 [ %122, %111 ], [ 0, %107 ]
%115 = add nuw nsw i64 %112, %61
%116 = getelementptr inbounds double, double* %1, i64 %115
%117 = load double, double* %116
%118 = getelementptr inbounds double, double* %2, i64 %112
%119 = load double, double* %118
%120 = tail call double @llvm.fmuladd.f64(double %117, double %119, double %113)
%121 = add nuw nsw i64 %112, 1
%122 = add i64 %114, 1
%123 = icmp eq i64 %122, %23
br i1 %123, label %124, label %111
124: ; preds = %111, %107
%125 = phi double [ %108, %107 ], [ %120, %111 ]
%126 = getelementptr inbounds double, double* %3, i64 %56
store double %125, double* %126
%127 = add nuw nsw i64 %56, 1
%128 = icmp eq i64 %127, %8
br i1 %128, label %129, label %55
129: ; preds = %124, %7, %19
ret void
}