void trback(const int n, double *__restrict__ a, double *__restrict__ z,
            const int ldz) {
  register int ik, k, j, i;
  int l, iz;
  double s, h;
  double *zp;
  for (i = 1; i < n; i++) {
    l = i - 1;
    iz = (i * (l) / 2) + i;
    ik = iz + i;
    h = a[ik];
    if (h == 0.0)
      continue;
    for (j = 0; j < n; j++) {
      s = 0.0;
      ik = iz;
      zp = z + j * ldz;
      for (k = 0; k <= l; ++ik, k++) {
        s += a[ik] * (*zp);
        ++zp;
      }
      s = (s / h) / h;
      ik = iz;
      zp = z + j * ldz;
      for (k = 0; k <= l; ++ik, k++) {
        (*zp) -= a[ik] * s;
        ++zp;
      }
    }
  }
}


x86 assembly (clang, -O0):
<trback>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12de <trback+0x1de>
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
imul   -0x30(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    117b <trback+0x7b>
jp     117b <trback+0x7b>
jmp    12d0 <trback+0x1d0>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12cb <trback+0x1cb>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     1212 <trback+0x112>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x50(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x50(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11bc <trback+0xbc>
movsd  -0x40(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     12b8 <trback+0x1b8>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x50(%rbp),%rax
movsd  (%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x50(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    124c <trback+0x14c>
jmp    12bd <trback+0x1bd>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1182 <trback+0x82>
jmp    12d0 <trback+0x1d0>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1119 <trback+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<trback>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    $0x2,%edi
jl     13e6 <trback+0x2e6>
movslq %ecx,%r14
mov    %edi,%r13d
lea    0x10(%rdx),%rax
mov    %rax,-0x28(%rsp)
lea    0x0(,%r14,8),%rax
mov    %rax,-0x10(%rsp)
lea    0x18(%rsi),%rax
mov    %rax,-0x38(%rsp)
mov    $0x1,%edi
xor    %ebx,%ebx
mov    $0x1,%eax
mov    %rax,-0x58(%rsp)
mov    $0x2,%eax
mov    %rax,-0x40(%rsp)
xorpd  %xmm0,%xmm0
jmp    1191 <trback+0x91>
nopl   0x0(%rax)
mov    -0x50(%rsp),%rdi
add    $0x1,%rdi
mov    -0x58(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x58(%rsp)
addq   $0x1,-0x40(%rsp)
addq   $0x8,-0x38(%rsp)
mov    -0x5c(%rsp),%eax
mov    %eax,%ebx
cmp    %r13,%rdi
je     13e6 <trback+0x2e6>
lea    0x1(%rbx),%eax
mov    %eax,-0x5c(%rsp)
lea    -0x1(%rdi),%eax
imul   %edi,%eax
mov    %eax,%ecx
shr    $0x1f,%ecx
add    %eax,%ecx
sar    %ecx
mov    %rdi,-0x50(%rsp)
lea    (%rcx,%rdi,2),%eax
cltq
movsd  (%rsi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
jne    11be <trback+0xbe>
jnp    1160 <trback+0x60>
mov    -0x50(%rsp),%rax
mov    %rax,%r8
and    $0xfffffffffffffffc,%r8
mov    -0x58(%rsp),%rbp
mov    %ebp,%edi
and    $0x3,%edi
and    $0xfffffffc,%ebp
neg    %ebp
movslq %ecx,%r10
lea    (%rax,%r10,1),%r9
mov    %rax,%rcx
movabs $0x7ffffffffffffffc,%rax
ff ff 7f
and    %rax,%rcx
mov    %r9,-0x48(%rsp)
mov    %rcx,-0x30(%rsp)
lea    (%r9,%rcx,1),%rax
mov    %rax,-0x18(%rsp)
mov    -0x40(%rsp),%rax
add    %r10,%rax
lea    (%rsi,%rax,8),%rax
mov    %rax,-0x20(%rsp)
mov    -0x38(%rsp),%rax
lea    (%rax,%r10,8),%r15
mov    -0x28(%rsp),%r10
xor    %r9d,%r9d
mov    %rbx,-0x8(%rsp)
jmp    1247 <trback+0x147>
nopl   0x0(%rax)
add    $0x1,%r9
add    -0x10(%rsp),%r10
cmp    %r13,%r9
mov    -0x8(%rsp),%rbx
je     1160 <trback+0x60>
mov    %r9,%rcx
imul   %r14,%rcx
lea    (%rdx,%rcx,8),%r12
xorpd  %xmm2,%xmm2
mov    -0x48(%rsp),%rcx
mov    %r12,%r11
cmp    $0x3,%ebx
jb     12dd <trback+0x1dd>
mov    %r13,%rbx
mov    %r14,%rax
xorpd  %xmm2,%xmm2
mov    -0x20(%rsp),%r14
xor    %r13d,%r13d
mov    %r12,%r11
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%r14),%xmm3
movsd  (%r14),%xmm4
mulsd  (%r11),%xmm3
addsd  %xmm2,%xmm3
mulsd  0x8(%r11),%xmm4
addsd  %xmm3,%xmm4
movsd  0x8(%r14),%xmm3
mulsd  0x10(%r11),%xmm3
addsd  %xmm4,%xmm3
movsd  0x10(%r14),%xmm2
mulsd  0x18(%r11),%xmm2
addsd  %xmm3,%xmm2
add    $0x20,%r11
add    $0xfffffffffffffffc,%r13
add    $0x20,%r14
cmp    %r13d,%ebp
jne    1280 <trback+0x180>
mov    -0x48(%rsp),%rcx
sub    %r13,%rcx
mov    %rax,%r14
mov    %rbx,%r13
testb  $0x3,-0x5c(%rsp)
je     1307 <trback+0x207>
lea    (%rsi,%rcx,8),%rbx
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movsd  (%rbx,%rcx,8),%xmm3
mulsd  (%r11,%rcx,8),%xmm3
addsd  %xmm3,%xmm2
add    $0x1,%rcx
cmp    %ecx,%edi
jne    12f0 <trback+0x1f0>
xor    %r11d,%r11d
divsd  %xmm1,%xmm2
divsd  %xmm1,%xmm2
cmpq   $0x4,-0x50(%rsp)
jae    1330 <trback+0x230>
mov    -0x48(%rsp),%rcx
jmp    13a8 <trback+0x2a8>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    -0x30(%rsp),%rax
lea    (%r12,%rax,8),%r12
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x10(%r15,%rbx,8),%xmm4
movupd (%r15,%rbx,8),%xmm5
movupd -0x10(%r10,%rbx,8),%xmm6
movupd (%r10,%rbx,8),%xmm7
mulpd  %xmm3,%xmm4
subpd  %xmm4,%xmm6
mulpd  %xmm3,%xmm5
subpd  %xmm5,%xmm7
movupd %xmm6,-0x10(%r10,%rbx,8)
movupd %xmm7,(%r10,%rbx,8)
add    $0x4,%rbx
cmp    %rbx,%r8
jne    1350 <trback+0x250>
mov    -0x18(%rsp),%rcx
mov    -0x30(%rsp),%rax
mov    %eax,%r11d
cmp    %rax,-0x50(%rsp)
je     1230 <trback+0x130>
lea    (%rsi,%rcx,8),%rbx
mov    -0x58(%rsp),%rcx
sub    %r11d,%ecx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
movsd  (%r12,%rax,8),%xmm3
movsd  (%rbx,%rax,8),%xmm4
mulsd  %xmm2,%xmm4
subsd  %xmm4,%xmm3
movsd  %xmm3,(%r12,%rax,8)
add    $0x1,%rax
cmp    %eax,%ecx
jne    13c0 <trback+0x2c0>
jmp    1230 <trback+0x130>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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