void _backward(double *beta, const double *A, const double *pobs, int N,
               int T) {
  int i, j, t;
  double sum, scaling;
  scaling = 0.0;
  for (i = 0; i < N; i++) {
    beta[(T - 1) * N + i] = 1.0;
    scaling += beta[(T - 1) * N + i];
  }
  for (i = 0; i < N; i++)
    beta[(T - 1) * N + i] /= scaling;
  for (t = T - 2; t >= 0; t--) {
    scaling = 0.0;
    for (i = 0; i < N; i++) {
      sum = 0.0;
      for (j = 0; j < N; j++) {
        sum += A[i * N + j] * pobs[(t + 1) * N + j] * beta[(t + 1) * N + j];
      }
      beta[t * N + i] = sum;
      scaling += sum;
    }
    if (scaling != 0)
      for (j = 0; j < N; j++)
        beta[t * N + j] /= scaling;
  }
}


x86 assembly (clang, -O0):
<_backward>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1184 <_backward+0x84>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xeb2(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1126 <_backward+0x26>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11cc <_backward+0xcc>
movsd  -0x40(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    118b <_backward+0x8b>
mov    -0x20(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jl     132f <_backward+0x22f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12b5 <_backward+0x1b5>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    127d <_backward+0x17d>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1209 <_backward+0x109>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11ee <_backward+0xee>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12d2 <_backward+0x1d2>
jp     12d2 <_backward+0x1d2>
jmp    131c <_backward+0x21c>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1317 <_backward+0x217>
movsd  -0x40(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12d9 <_backward+0x1d9>
jmp    131c <_backward+0x21c>
jmp    1321 <_backward+0x221>
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    11d5 <_backward+0xd5>
pop    %rbp
ret

x86 assembly (clang, -O2):
<_backward>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8d,%r9d
mov    %rdx,-0x48(%rsp)
mov    %rsi,-0x50(%rsp)
mov    %ecx,%r11d
test   %ecx,%ecx
jle    1287 <_backward+0x187>
mov    %r9,%r10
lea    -0x1(%r9),%eax
imul   %ecx,%eax
movslq %eax,%r8
lea    -0x1(%r11),%rax
mov    %r11d,%edx
and    $0x3,%edx
movabs $0x3ff0000000000000,%rbp
cmp    $0x3,%rax
jae    1151 <_backward+0x51>
xorpd  %xmm0,%xmm0
xor    %esi,%esi
jmp    119c <_backward+0x9c>
mov    %r11d,%ebx
and    $0xfffffffc,%ebx
lea    (%rdi,%r8,8),%rax
add    $0x18,%rax
xorpd  %xmm0,%xmm0
xor    %esi,%esi
movsd  0xe93(%rip),%xmm1
nopl   (%rax)
mov    %rbp,-0x18(%rax,%rsi,8)
addsd  %xmm1,%xmm0
mov    %rbp,-0x10(%rax,%rsi,8)
addsd  %xmm1,%xmm0
mov    %rbp,-0x8(%rax,%rsi,8)
addsd  %xmm1,%xmm0
mov    %rbp,(%rax,%rsi,8)
addsd  %xmm1,%xmm0
add    $0x4,%rsi
cmp    %rsi,%rbx
jne    1170 <_backward+0x70>
test   %rdx,%rdx
je     11d1 <_backward+0xd1>
add    %r8,%rsi
lea    (%rdi,%rsi,8),%rax
xor    %esi,%esi
movsd  0xe4e(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %rbp,(%rax,%rsi,8)
addsd  %xmm1,%xmm0
add    $0x1,%rsi
cmp    %rsi,%rdx
jne    11c0 <_backward+0xc0>
test   %ecx,%ecx
mov    %r10,%r9
jle    1287 <_backward+0x187>
cmp    $0x1,%ecx
jne    11e8 <_backward+0xe8>
xor    %edx,%edx
jmp    1267 <_backward+0x167>
mov    %r11d,%edx
and    $0xfffffffe,%edx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rdx),%rax
mov    %rax,%rsi
shr    %rsi
add    $0x1,%rsi
test   %rax,%rax
je     14f4 <_backward+0x3f4>
mov    %rsi,%rax
and    $0xfffffffffffffffe,%rax
lea    (%rdi,%r8,8),%rbx
add    $0x10,%rbx
xor    %ebp,%ebp
xchg   %ax,%ax
movupd -0x10(%rbx,%rbp,8),%xmm2
movupd (%rbx,%rbp,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rbx,%rbp,8)
divpd  %xmm1,%xmm3
movupd %xmm3,(%rbx,%rbp,8)
add    $0x4,%rbp
add    $0xfffffffffffffffe,%rax
jne    1220 <_backward+0x120>
test   $0x1,%sil
je     125f <_backward+0x15f>
add    %r8,%rbp
movupd (%rdi,%rbp,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdi,%rbp,8)
cmp    %r11,%rdx
mov    %r10,%r9
je     1287 <_backward+0x187>
lea    (%rdi,%r8,8),%rax
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r11
jne    1270 <_backward+0x170>
cmp    $0x2,%r9d
jl     14e9 <_backward+0x3e9>
add    $0xfffffffe,%r9d
movslq %ecx,%rbp
mov    %r11d,%r10d
and    $0xfffffffe,%r10d
lea    -0x2(%r10),%rax
mov    %rax,-0x28(%rsp)
mov    %rax,%rdx
shr    %rdx
add    $0x1,%rdx
mov    %rdx,%rax
mov    %rdx,-0x30(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x38(%rsp)
mov    -0x50(%rsp),%rax
add    $0x8,%rax
mov    %rax,-0x20(%rsp)
lea    0x0(,%rbp,8),%r13
lea    0x1(%r9),%rdx
imul   %rbp,%rdx
lea    (%rdi,%rdx,8),%rbx
add    $0x8,%rbx
mov    %r13,%rax
neg    %rax
mov    %rax,-0x18(%rsp)
mov    -0x48(%rsp),%rax
lea    (%rax,%rdx,8),%r14
add    $0x8,%r14
mov    %rbp,-0x40(%rsp)
mov    %rbp,%rdx
imul   %r9,%rdx
lea    (%rdi,%rdx,8),%rsi
add    $0x10,%rsi
lea    (%rdi,%rdx,8),%r15
xorpd  %xmm0,%xmm0
jmp    1341 <_backward+0x241>
nop
lea    -0x1(%r9),%rdx
mov    -0x18(%rsp),%rax
add    %rax,%rbx
add    %rax,%r14
add    %rax,%rsi
add    %rax,%r15
test   %r9,%r9
mov    %rdx,%r9
jle    14e9 <_backward+0x3e9>
test   %ecx,%ecx
jle    1320 <_backward+0x220>
mov    %rsi,-0x10(%rsp)
lea    0x1(%r9),%r12
mov    -0x40(%rsp),%rax
imul   %rax,%r12
mov    %r9,-0x8(%rsp)
mov    %r9,%r8
imul   %rax,%r8
xorpd  %xmm1,%xmm1
mov    -0x20(%rsp),%rbp
xor    %esi,%esi
jmp    138d <_backward+0x28d>
lea    (%rsi,%r8,1),%rdx
movsd  %xmm2,(%rdi,%rdx,8)
addsd  %xmm2,%xmm1
add    $0x1,%rsi
add    %r13,%rbp
cmp    %r11,%rsi
je     1430 <_backward+0x330>
cmp    $0x1,%ecx
jne    13a0 <_backward+0x2a0>
xor    %edx,%edx
xorpd  %xmm2,%xmm2
jmp    13e9 <_backward+0x2e9>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm3,%xmm3
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%rbp,%rdx,8),%xmm4
movsd  0x0(%rbp,%rdx,8),%xmm2
mulsd  -0x8(%r14,%rdx,8),%xmm4
mulsd  -0x8(%rbx,%rdx,8),%xmm4
mulsd  (%r14,%rdx,8),%xmm2
mulsd  (%rbx,%rdx,8),%xmm2
addsd  %xmm3,%xmm4
addsd  %xmm4,%xmm2
add    $0x2,%rdx
movapd %xmm2,%xmm3
cmp    %rdx,%r10
jne    13b0 <_backward+0x2b0>
test   $0x1,%r11b
je     1370 <_backward+0x270>
mov    %rsi,%r9
imul   -0x40(%rsp),%r9
add    %rdx,%r9
mov    -0x50(%rsp),%rax
movsd  (%rax,%r9,8),%xmm3
add    %r12,%rdx
mov    -0x48(%rsp),%rax
mulsd  (%rax,%rdx,8),%xmm3
mulsd  (%rdi,%rdx,8),%xmm3
addsd  %xmm3,%xmm2
jmp    1370 <_backward+0x270>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
ucomisd %xmm0,%xmm1
mov    -0x8(%rsp),%r9
mov    -0x10(%rsp),%rsi
jne    1446 <_backward+0x346>
jnp    1320 <_backward+0x220>
test   %ecx,%ecx
jle    1320 <_backward+0x220>
cmp    $0x2,%ecx
jae    1457 <_backward+0x357>
xor    %edx,%edx
jmp    14c0 <_backward+0x3c0>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
cmpq   $0x0,-0x28(%rsp)
je     14de <_backward+0x3de>
mov    -0x38(%rsp),%rbp
xor    %edx,%edx
xchg   %ax,%ax
movupd -0x10(%rsi,%rdx,8),%xmm3
movupd (%rsi,%rdx,8),%xmm4
divpd  %xmm2,%xmm3
movupd %xmm3,-0x10(%rsi,%rdx,8)
divpd  %xmm2,%xmm4
movupd %xmm4,(%rsi,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rbp
jne    1470 <_backward+0x370>
testb  $0x1,-0x30(%rsp)
je     14b0 <_backward+0x3b0>
add    %r8,%rdx
movupd (%rdi,%rdx,8),%xmm3
divpd  %xmm2,%xmm3
movupd %xmm3,(%rdi,%rdx,8)
mov    %r10,%rdx
cmp    %r11,%r10
je     1320 <_backward+0x220>
nopl   0x0(%rax)
movsd  (%r15,%rdx,8),%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%r15,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r11
jne    14c0 <_backward+0x3c0>
jmp    1320 <_backward+0x220>
xor    %edx,%edx
testb  $0x1,-0x30(%rsp)
jne    149f <_backward+0x39f>
jmp    14b0 <_backward+0x3b0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ebp,%ebp
test   $0x1,%sil
jne    124e <_backward+0x14e>
jmp    125f <_backward+0x15f>

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

LLVM IR (-O2):
define void @_backward(double* nocapture %0, double* nocapture %1, double* nocapture %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %3, 0
br i1 %6, label %7, label %114
7: ; preds = %5
%8 = add nsw i32 %4, -1
%9 = mul nsw i32 %8, %3
%10 = sext i32 %9 to i64
%11 = zext i32 %3 to i64
%12 = add nsw i64 %11, -1
%13 = and i64 %11, 3
%14 = icmp ult i64 %12, 3
br i1 %14, label %17, label %15
15: ; preds = %7
%16 = and i64 %11, 4294967292
br label %84
17: ; preds = %84, %7
%18 = phi double [ undef, %7 ], [ %102, %84 ]
%19 = phi i64 [ 0, %7 ], [ %103, %84 ]
%20 = phi double [ 0.000000e+00, %7 ], [ %102, %84 ]
%21 = icmp eq i64 %13, 0
br i1 %21, label %32, label %22
22: ; preds = %17, %22
%23 = phi i64 [ %29, %22 ], [ %19, %17 ]
%24 = phi double [ %28, %22 ], [ %20, %17 ]
%25 = phi i64 [ %30, %22 ], [ 0, %17 ]
%26 = add nsw i64 %23, %10
%27 = getelementptr inbounds double, double* %0, i64 %26
store double 1.000000e+00, double* %27
%28 = fadd double %24, 1.000000e+00
%29 = add nuw nsw i64 %23, 1
%30 = add i64 %25, 1
%31 = icmp eq i64 %30, %13
br i1 %31, label %32, label %22
32: ; preds = %22, %17
%33 = phi double [ %18, %17 ], [ %28, %22 ]
br i1 %6, label %34, label %114
34: ; preds = %32
%35 = add nsw i32 %4, -1
%36 = mul nsw i32 %35, %3
%37 = sext i32 %36 to i64
%38 = zext i32 %3 to i64
%39 = icmp eq i32 %3, 1
br i1 %39, label %82, label %40
40: ; preds = %34
%41 = and i64 %11, 4294967294
%42 = insertelement <2 x double> poison, double %33, i64 0
%43 = shufflevector <2 x double> %42, <2 x double> poison, <2 x i32> zeroinitializer
%44 = add nsw i64 %41, -2
%45 = lshr exact i64 %44, 1
%46 = add nuw i64 %45, 1
%47 = and i64 %46, 1
%48 = icmp eq i64 %44, 0
br i1 %48, label %70, label %49
49: ; preds = %40
%50 = and i64 %46, -2
br label %51
51: ; preds = %51, %49
%52 = phi i64 [ 0, %49 ], [ %67, %51 ]
%53 = phi i64 [ 0, %49 ], [ %68, %51 ]
%54 = add nsw i64 %52, %37
%55 = getelementptr inbounds double, double* %0, i64 %54
%56 = bitcast double* %55 to <2 x double>*
%57 = load <2 x double>, <2 x double>* %56
%58 = fdiv <2 x double> %57, %43
%59 = bitcast double* %55 to <2 x double>*
store <2 x double> %58, <2 x double>* %59
%60 = or i64 %52, 2
%61 = add nsw i64 %60, %37
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = fdiv <2 x double> %64, %43
%66 = bitcast double* %62 to <2 x double>*
store <2 x double> %65, <2 x double>* %66
%67 = add nuw i64 %52, 4
%68 = add nuw i64 %53, 2
%69 = icmp eq i64 %68, %50
br i1 %69, label %70, label %51
70: ; preds = %51, %40
%71 = phi i64 [ 0, %40 ], [ %67, %51 ]
%72 = icmp eq i64 %47, 0
br i1 %72, label %80, label %73
73: ; preds = %70
%74 = add nsw i64 %71, %37
%75 = getelementptr inbounds double, double* %0, i64 %74
%76 = bitcast double* %75 to <2 x double>*
%77 = load <2 x double>, <2 x double>* %76
%78 = fdiv <2 x double> %77, %43
%79 = bitcast double* %75 to <2 x double>*
store <2 x double> %78, <2 x double>* %79
br label %80
80: ; preds = %70, %73
%81 = icmp eq i64 %41, %11
br i1 %81, label %114, label %82
82: ; preds = %34, %80
%83 = phi i64 [ 0, %34 ], [ %41, %80 ]
br label %106
84: ; preds = %84, %15
%85 = phi i64 [ 0, %15 ], [ %103, %84 ]
%86 = phi double [ 0.000000e+00, %15 ], [ %102, %84 ]
%87 = phi i64 [ 0, %15 ], [ %104, %84 ]
%88 = add nsw i64 %85, %10
%89 = getelementptr inbounds double, double* %0, i64 %88
store double 1.000000e+00, double* %89
%90 = fadd double %86, 1.000000e+00
%91 = or i64 %85, 1
%92 = add nsw i64 %91, %10
%93 = getelementptr inbounds double, double* %0, i64 %92
store double 1.000000e+00, double* %93
%94 = fadd double %90, 1.000000e+00
%95 = or i64 %85, 2
%96 = add nsw i64 %95, %10
%97 = getelementptr inbounds double, double* %0, i64 %96
store double 1.000000e+00, double* %97
%98 = fadd double %94, 1.000000e+00
%99 = or i64 %85, 3
%100 = add nsw i64 %99, %10
%101 = getelementptr inbounds double, double* %0, i64 %100
store double 1.000000e+00, double* %101
%102 = fadd double %98, 1.000000e+00
%103 = add nuw nsw i64 %85, 4
%104 = add i64 %87, 4
%105 = icmp eq i64 %104, %16
br i1 %105, label %17, label %84
106: ; preds = %82, %106
%107 = phi i64 [ %112, %106 ], [ %83, %82 ]
%108 = add nsw i64 %107, %37
%109 = getelementptr inbounds double, double* %0, i64 %108
%110 = load double, double* %109
%111 = fdiv double %110, %33
store double %111, double* %109
%112 = add nuw nsw i64 %107, 1
%113 = icmp eq i64 %112, %38
br i1 %113, label %114, label %106
114: ; preds = %106, %80, %5, %32
%115 = icmp sgt i32 %4, 1
br i1 %115, label %116, label %248
116: ; preds = %114
%117 = add nsw i32 %4, -2
%118 = icmp slt i32 %3, 1
%119 = sext i32 %3 to i64
%120 = zext i32 %117 to i64
%121 = zext i32 %3 to i64
%122 = zext i32 %3 to i64
%123 = and i64 %121, 4294967294
%124 = add nsw i64 %123, -2
%125 = lshr exact i64 %124, 1
%126 = add nuw i64 %125, 1
%127 = and i64 %121, 1
%128 = icmp eq i32 %3, 1
%129 = and i64 %121, 4294967294
%130 = icmp eq i64 %127, 0
%131 = icmp ult i32 %3, 2
%132 = and i64 %121, 4294967294
%133 = and i64 %126, 1
%134 = icmp eq i64 %124, 0
%135 = and i64 %126, -2
%136 = icmp eq i64 %133, 0
%137 = icmp eq i64 %132, %121
br label %138
138: ; preds = %116, %245
%139 = phi i64 [ %120, %116 ], [ %246, %245 ]
br i1 %6, label %140, label %245
140: ; preds = %138
%141 = add nuw nsw i64 %139, 1
%142 = mul nsw i64 %141, %119
%143 = mul nsw i64 %139, %119
br label %144
144: ; preds = %191, %140
%145 = phi i64 [ 0, %140 ], [ %196, %191 ]
%146 = phi double [ 0.000000e+00, %140 ], [ %195, %191 ]
%147 = mul nsw i64 %145, %119
br i1 %128, label %176, label %148
148: ; preds = %144, %148
%149 = phi i64 [ %173, %148 ], [ 0, %144 ]
%150 = phi double [ %172, %148 ], [ 0.000000e+00, %144 ]
%151 = phi i64 [ %174, %148 ], [ 0, %144 ]
%152 = add nsw i64 %149, %147
%153 = getelementptr inbounds double, double* %1, i64 %152
%154 = load double, double* %153
%155 = add nsw i64 %149, %142
%156 = getelementptr inbounds double, double* %2, i64 %155
%157 = load double, double* %156
%158 = fmul double %154, %157
%159 = getelementptr inbounds double, double* %0, i64 %155
%160 = load double, double* %159
%161 = tail call double @llvm.fmuladd.f64(double %158, double %160, double %150)
%162 = or i64 %149, 1
%163 = add nsw i64 %162, %147
%164 = getelementptr inbounds double, double* %1, i64 %163
%165 = load double, double* %164
%166 = add nsw i64 %162, %142
%167 = getelementptr inbounds double, double* %2, i64 %166
%168 = load double, double* %167
%169 = fmul double %165, %168
%170 = getelementptr inbounds double, double* %0, i64 %166
%171 = load double, double* %170
%172 = tail call double @llvm.fmuladd.f64(double %169, double %171, double %161)
%173 = add nuw nsw i64 %149, 2
%174 = add i64 %151, 2
%175 = icmp eq i64 %174, %129
br i1 %175, label %176, label %148
176: ; preds = %148, %144
%177 = phi double [ undef, %144 ], [ %172, %148 ]
%178 = phi i64 [ 0, %144 ], [ %173, %148 ]
%179 = phi double [ 0.000000e+00, %144 ], [ %172, %148 ]
br i1 %130, label %191, label %180
180: ; preds = %176
%181 = add nsw i64 %178, %147
%182 = getelementptr inbounds double, double* %1, i64 %181
%183 = load double, double* %182
%184 = add nsw i64 %178, %142
%185 = getelementptr inbounds double, double* %2, i64 %184
%186 = load double, double* %185
%187 = fmul double %183, %186
%188 = getelementptr inbounds double, double* %0, i64 %184
%189 = load double, double* %188
%190 = tail call double @llvm.fmuladd.f64(double %187, double %189, double %179)
br label %191
191: ; preds = %176, %180
%192 = phi double [ %177, %176 ], [ %190, %180 ]
%193 = add nsw i64 %145, %143
%194 = getelementptr inbounds double, double* %0, i64 %193
store double %192, double* %194
%195 = fadd double %146, %192
%196 = add nuw nsw i64 %145, 1
%197 = icmp eq i64 %196, %121
br i1 %197, label %198, label %144
198: ; preds = %191
%199 = fcmp oeq double %195, 0.000000e+00
%200 = or i1 %199, %118
br i1 %200, label %245, label %201
201: ; preds = %198
%202 = mul nsw i64 %139, %119
br i1 %131, label %235, label %203
203: ; preds = %201
%204 = insertelement <2 x double> poison, double %195, i64 0
%205 = shufflevector <2 x double> %204, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %134, label %225, label %206
206: ; preds = %203, %206
%207 = phi i64 [ %222, %206 ], [ 0, %203 ]
%208 = phi i64 [ %223, %206 ], [ 0, %203 ]
%209 = add nsw i64 %207, %202
%210 = getelementptr inbounds double, double* %0, i64 %209
%211 = bitcast double* %210 to <2 x double>*
%212 = load <2 x double>, <2 x double>* %211
%213 = fdiv <2 x double> %212, %205
%214 = bitcast double* %210 to <2 x double>*
store <2 x double> %213, <2 x double>* %214
%215 = or i64 %207, 2
%216 = add nsw i64 %215, %202
%217 = getelementptr inbounds double, double* %0, i64 %216
%218 = bitcast double* %217 to <2 x double>*
%219 = load <2 x double>, <2 x double>* %218
%220 = fdiv <2 x double> %219, %205
%221 = bitcast double* %217 to <2 x double>*
store <2 x double> %220, <2 x double>* %221
%222 = add nuw i64 %207, 4
%223 = add nuw i64 %208, 2
%224 = icmp eq i64 %223, %135
br i1 %224, label %225, label %206
225: ; preds = %206, %203
%226 = phi i64 [ 0, %203 ], [ %222, %206 ]
br i1 %136, label %234, label %227
227: ; preds = %225
%228 = add nsw i64 %226, %202
%229 = getelementptr inbounds double, double* %0, i64 %228
%230 = bitcast double* %229 to <2 x double>*
%231 = load <2 x double>, <2 x double>* %230
%232 = fdiv <2 x double> %231, %205
%233 = bitcast double* %229 to <2 x double>*
store <2 x double> %232, <2 x double>* %233
br label %234
234: ; preds = %225, %227
br i1 %137, label %245, label %235
235: ; preds = %201, %234
%236 = phi i64 [ 0, %201 ], [ %132, %234 ]
br label %237
237: ; preds = %235, %237
%238 = phi i64 [ %243, %237 ], [ %236, %235 ]
%239 = add nsw i64 %238, %202
%240 = getelementptr inbounds double, double* %0, i64 %239
%241 = load double, double* %240
%242 = fdiv double %241, %195
store double %242, double* %240
%243 = add nuw nsw i64 %238, 1
%244 = icmp eq i64 %243, %122
br i1 %244, label %245, label %237
245: ; preds = %237, %234, %138, %198
%246 = add nsw i64 %139, -1
%247 = icmp sgt i64 %139, 0
br i1 %247, label %138, label %248
248: ; preds = %245, %114
ret void
}