#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double log(double);
double _forward(double *alpha, const double *A, const double *pobs,
                const double *pi, int N, int T) {
  int i, j, t;
  double sum, logprob, scaling;
  scaling = 0.0;
  for (i = 0; i < N; i++) {
    alpha[i] = pi[i] * pobs[i];
    scaling += alpha[i];
  }
  logprob = log(scaling);
  if (scaling != 0)
    for (i = 0; i < N; i++)
      alpha[i] /= scaling;
  for (t = 0; t < T - 1; t++) {
    scaling = 0.0;
    for (j = 0; j < N; j++) {
      sum = 0.0;
      for (i = 0; i < N; i++) {
        sum += alpha[t * N + i] * A[i * N + j];
      }
      alpha[(t + 1) * N + j] = sum * pobs[(t + 1) * N + j];
      scaling += alpha[(t + 1) * N + j];
    }
    if (scaling != 0)
      for (j = 0; j < N; j++)
        alpha[(t + 1) * N + j] /= scaling;
    logprob += log(scaling);
  }
  return logprob;
}


x86 assembly (clang, -O0):
<_forward>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1197 <_forward+0x87>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    113f <_forward+0x2f>
movsd  -0x50(%rbp),%xmm0
call   1030 <log@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11c3 <_forward+0xb3>
jp     11c3 <_forward+0xb3>
jmp    1204 <_forward+0xf4>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11ff <_forward+0xef>
movsd  -0x50(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ca <_forward+0xba>
jmp    1204 <_forward+0xf4>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1392 <_forward+0x282>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1306 <_forward+0x1f6>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    129e <_forward+0x18e>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1246 <_forward+0x136>
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
imul   -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
imul   -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
imul   -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    122b <_forward+0x11b>
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1323 <_forward+0x213>
jp     1323 <_forward+0x213>
jmp    1370 <_forward+0x260>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    136b <_forward+0x25b>
movsd  -0x50(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
imul   -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    132a <_forward+0x21a>
jmp    1370 <_forward+0x260>
movsd  -0x50(%rbp),%xmm0
call   1030 <log@plt>
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    120b <_forward+0xfb>
movsd  -0x48(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<_forward>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r9d,%r14d
mov    %rdx,%rbp
mov    %rsi,0x18(%rsp)
mov    %rdi,%rbx
mov    %r8d,%r15d
test   %r8d,%r8d
mov    %r8d,0xc(%rsp)
jle    1147 <_forward+0x37>
cmp    $0x1,%r8d
jne    115e <_forward+0x4e>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    11a3 <_forward+0x93>
xorpd  %xmm0,%xmm0
call   1030 <log@plt>
mov    0xc(%rsp),%r11d
movapd %xmm0,%xmm4
jmp    1287 <_forward+0x177>
mov    %r15d,%edx
and    $0xfffffffe,%edx
xorpd  %xmm0,%xmm0
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movsd  (%rcx,%rax,8),%xmm1
mulsd  0x0(%rbp,%rax,8),%xmm1
movsd  %xmm1,(%rbx,%rax,8)
addsd  %xmm0,%xmm1
movsd  0x8(%rcx,%rax,8),%xmm0
mulsd  0x8(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rax,8)
addsd  %xmm1,%xmm0
add    $0x2,%rax
cmp    %rax,%rdx
jne    1170 <_forward+0x60>
test   $0x1,%r15b
je     11bd <_forward+0xad>
movsd  (%rcx,%rax,8),%xmm1
mulsd  0x0(%rbp,%rax,8),%xmm1
movsd  %xmm1,(%rbx,%rax,8)
addsd  %xmm1,%xmm0
movapd %xmm0,0x20(%rsp)
call   1030 <log@plt>
movapd 0x20(%rsp),%xmm3
movapd %xmm0,%xmm4
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm3
jne    11de <_forward+0xce>
jnp    11f6 <_forward+0xe6>
mov    0xc(%rsp),%r11d
test   %r11d,%r11d
jle    1287 <_forward+0x177>
cmp    $0x2,%r11d
jae    1200 <_forward+0xf0>
xor    %eax,%eax
jmp    1270 <_forward+0x160>
mov    0xc(%rsp),%r11d
jmp    1287 <_forward+0x177>
mov    %r15d,%eax
and    $0xfffffffe,%eax
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
lea    -0x2(%rax),%rdx
mov    %rdx,%rcx
shr    %rcx
add    $0x1,%rcx
test   %rdx,%rdx
je     150f <_forward+0x3ff>
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %edx,%edx
xchg   %ax,%ax
movupd (%rbx,%rdx,8),%xmm1
movupd 0x10(%rbx,%rdx,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,(%rbx,%rdx,8)
divpd  %xmm0,%xmm2
movupd %xmm2,0x10(%rbx,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1230 <_forward+0x120>
test   $0x1,%cl
je     126b <_forward+0x15b>
movupd (%rbx,%rdx,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%rbx,%rdx,8)
cmp    %r15,%rax
je     1287 <_forward+0x177>
movsd  (%rbx,%rax,8),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    1270 <_forward+0x160>
mov    %r14,%rax
cmp    $0x2,%eax
jl     14fc <_forward+0x3ec>
mov    %rbp,0x20(%rsp)
add    $0xffffffff,%eax
movslq %r11d,%r12
mov    %r15d,%ebp
and    $0xfffffffe,%ebp
lea    -0x2(%rbp),%rcx
mov    %rcx,0x48(%rsp)
mov    %rcx,%rdx
shr    %rdx
add    $0x1,%rdx
mov    %rdx,%rcx
mov    %rdx,0x40(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x38(%rsp)
lea    0x8(%rbx),%rcx
mov    %rcx,0x50(%rsp)
mov    %r12,%r14
shl    $0x4,%r14
mov    %rbx,%r13
lea    0x10(%rbx),%rcx
mov    %rcx,0x30(%rsp)
xor    %r8d,%r8d
mov    %rax,%r10
mov    %rax,0x58(%rsp)
jmp    132c <_forward+0x21c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
mov    %r8,%rbx
call   1030 <log@plt>
movsd  0x60(%rsp),%xmm4
addsd  %xmm0,%xmm4
mov    %rbx,%r8
mov    0x58(%rsp),%r10
cmp    %r10,%rbx
mov    0xc(%rsp),%r11d
je     14fc <_forward+0x3ec>
xorpd  %xmm0,%xmm0
test   %r11d,%r11d
movsd  %xmm4,0x60(%rsp)
jle    1300 <_forward+0x1f0>
mov    %r8,%r9
imul   %r12,%r9
lea    0x1(%r8),%rax
mov    %rax,0x10(%rsp)
imul   %r12,%rax
mov    0x50(%rsp),%rcx
lea    (%rcx,%r9,8),%rsi
xorpd  %xmm0,%xmm0
mov    0x18(%rsp),%rdx
xor    %ecx,%ecx
jmp    139a <_forward+0x28a>
cs nopw 0x0(%rax,%rax,1)
nop
lea    (%rcx,%rax,1),%rdi
mov    0x20(%rsp),%rbx
mulsd  (%rbx,%rdi,8),%xmm1
movsd  %xmm1,0x0(%r13,%rdi,8)
addsd  %xmm1,%xmm0
add    $0x1,%rcx
add    $0x8,%rdx
cmp    %r15,%rcx
je     1420 <_forward+0x310>
cmp    $0x1,%r11d
jne    13b0 <_forward+0x2a0>
xor    %ebx,%ebx
xorpd  %xmm1,%xmm1
jmp    13ed <_forward+0x2dd>
nopl   0x0(%rax,%rax,1)
xorpd  %xmm2,%xmm2
mov    %rdx,%rdi
xor    %ebx,%ebx
nopl   0x0(%rax)
movsd  -0x8(%rsi,%rbx,8),%xmm3
mulsd  (%rdi),%xmm3
movsd  (%rsi,%rbx,8),%xmm1
addsd  %xmm2,%xmm3
mulsd  (%rdi,%r12,8),%xmm1
addsd  %xmm3,%xmm1
add    $0x2,%rbx
add    %r14,%rdi
movapd %xmm1,%xmm2
cmp    %rbx,%rbp
jne    13c0 <_forward+0x2b0>
test   $0x1,%r15b
je     1370 <_forward+0x260>
lea    (%rbx,%r9,1),%rdi
movsd  0x0(%r13,%rdi,8),%xmm2
imul   %r12,%rbx
add    %rcx,%rbx
mov    0x18(%rsp),%rdi
mulsd  (%rdi,%rbx,8),%xmm2
addsd  %xmm2,%xmm1
jmp    1370 <_forward+0x260>
nopl   0x0(%rax)
ucomisd 0xbd8(%rip),%xmm0
jne    1430 <_forward+0x320>
jnp    1300 <_forward+0x1f0>
test   %r11d,%r11d
jle    1300 <_forward+0x1f0>
cmp    $0x2,%r11d
jae    1443 <_forward+0x333>
xor    %ecx,%ecx
jmp    14bf <_forward+0x3af>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
cmpq   $0x0,0x48(%rsp)
je     14f1 <_forward+0x3e1>
mov    0x30(%rsp),%rcx
lea    (%rcx,%rax,8),%rdx
mov    0x38(%rsp),%rsi
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movupd -0x10(%rdx,%rcx,8),%xmm2
movupd (%rdx,%rcx,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rdx,%rcx,8)
divpd  %xmm1,%xmm3
movupd %xmm3,(%rdx,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1470 <_forward+0x360>
testb  $0x1,0x40(%rsp)
je     14b7 <_forward+0x3a7>
add    %rax,%rcx
mov    %r13,%rdx
movupd 0x0(%r13,%rcx,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,0x0(%r13,%rcx,8)
mov    %rbp,%rcx
cmp    %r15,%rbp
je     14e7 <_forward+0x3d7>
lea    0x0(,%rax,8),%rax
add    %r13,%rax
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r15
jne    14d0 <_forward+0x3c0>
mov    0x10(%rsp),%rbx
jmp    1307 <_forward+0x1f7>
xor    %ecx,%ecx
testb  $0x1,0x40(%rsp)
jne    149f <_forward+0x38f>
jmp    14b7 <_forward+0x3a7>
movapd %xmm4,%xmm0
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%cl
jne    125d <_forward+0x14d>
jmp    126b <_forward+0x15b>

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

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