#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 exp(double);
double pow(int, int);
void c_model2_sr(int *nin, int *itot, int *nsin, double *x, double *x1,
                 double *titre, double *titrepred, double *dd, int *ntheta,
                 double *theta) {
  int n = nin[0];
  int total_inf = itot[0];
  int nsamp = nsin[0];
  double T_1 = theta[1];
  double T_2 = theta[2];
  double mu = theta[0];
  int t_sample = n;
  double yrTitre[n * nsamp];
  int maskedInfectionHistory[n];
  double cumInfectionHistory[n];
  int k;
  int i;
  int j;
  int m;
  double xx2;
  for (k = 0; k < nsamp; k++) {
    for (j = 0; j < n; j++) {
      xx2 = 0;
      for (m = 0; m < n; m++) {
        if (m <= j) {
          maskedInfectionHistory[m] = x[m];
        } else {
          maskedInfectionHistory[m] = 0;
        }
      }
      cumInfectionHistory[0] = maskedInfectionHistory[0];
      for (m = 1; m < n; m++) {
        cumInfectionHistory[m] =
            cumInfectionHistory[m - 1] + maskedInfectionHistory[m];
      }
      for (i = 0; i < n; i++) {
        x1[i] = mu * dd[k * n + i] * maskedInfectionHistory[i] *
                pow(1 + T_1, (total_inf - cumInfectionHistory[i])) *
                exp(-1.0 * T_2 * (cumInfectionHistory[i] - 1.0));
      }
      for (i = 0; i < n; i++) {
        xx2 = xx2 + x1[i];
      }
      yrTitre[k + j * nsamp] = xx2;
    }
  }
  for (k = 0; k < nsamp; k++) {
    titrepred[k] = yrTitre[k + (t_sample - 1) * nsamp];
  }
}


x86 assembly (clang, -O0):
<c_model2_sr>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x3c(%rbp)
mov    0x28(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    0x28(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x3c(%rbp),%ecx
imul   %ecx,%eax
mov    %rsp,%rcx
mov    %rcx,-0x68(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0xb0(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x70(%rbp)
mov    -0x34(%rbp),%eax
lea    0xf(,%rax,4),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0xa8(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x78(%rbp)
mov    -0x34(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0xa0(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x80(%rbp)
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    14fa <c_model2_sr+0x3da>
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    14e1 <c_model2_sr+0x3c1>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12cf <c_model2_sr+0x1af>
mov    -0x90(%rbp),%eax
cmp    -0x8c(%rbp),%eax
jg     12a1 <c_model2_sr+0x181>
mov    -0xa8(%rbp),%rax
mov    -0x20(%rbp),%rcx
movslq -0x90(%rbp),%rdx
cvttsd2si (%rcx,%rdx,8),%edx
movslq -0x90(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    12b6 <c_model2_sr+0x196>
mov    -0xa8(%rbp),%rax
movslq -0x90(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
jmp    12bb <c_model2_sr+0x19b>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    125a <c_model2_sr+0x13a>
mov    -0xa0(%rbp),%rax
mov    -0xa8(%rbp),%rcx
cvtsi2sdl (%rcx),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    134d <c_model2_sr+0x22d>
mov    -0xa0(%rbp),%rax
mov    -0xa8(%rbp),%rcx
mov    -0x90(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
movsd  (%rax,%rdx,8),%xmm0
movslq -0x90(%rbp),%rdx
cvtsi2sdl (%rcx,%rdx,4),%xmm1
addsd  %xmm1,%xmm0
movslq -0x90(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    12ef <c_model2_sr+0x1cf>
movl   $0x0,-0x88(%rbp)
mov    -0x88(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1457 <c_model2_sr+0x337>
mov    -0xa0(%rbp),%rax
mov    -0xa8(%rbp),%rcx
movsd  -0x58(%rbp),%xmm0
mov    0x18(%rbp),%rdx
mov    -0x84(%rbp),%esi
imul   -0x34(%rbp),%esi
add    -0x88(%rbp),%esi
movslq %esi,%rsi
mulsd  (%rdx,%rsi,8),%xmm0
movslq -0x88(%rbp),%rdx
cvtsi2sdl (%rcx,%rdx,4),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  0xc4b(%rip),%xmm0
addsd  -0x48(%rbp),%xmm0
cvttsd2si %xmm0,%edi
cvtsi2sdl -0x38(%rbp),%xmm0
movslq -0x88(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
cvttsd2si %xmm0,%esi
call   1040 <pow@plt>
mov    -0xa0(%rbp),%rax
movaps %xmm0,%xmm1
movsd  -0xc0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0xc0a(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
movslq -0x88(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xbe9(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
call   1030 <exp@plt>
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x88(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x88(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x88(%rbp)
jmp    1357 <c_model2_sr+0x237>
movl   $0x0,-0x88(%rbp)
mov    -0x88(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    14a4 <c_model2_sr+0x384>
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x28(%rbp),%rax
movslq -0x88(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x88(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x88(%rbp)
jmp    1461 <c_model2_sr+0x341>
mov    -0xb0(%rbp),%rax
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x84(%rbp),%ecx
mov    -0x8c(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    1236 <c_model2_sr+0x116>
jmp    14e6 <c_model2_sr+0x3c6>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    121d <c_model2_sr+0xfd>
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1558 <c_model2_sr+0x438>
mov    -0xb0(%rbp),%rax
mov    -0x84(%rbp),%ecx
mov    -0x5c(%rbp),%edx
sub    $0x1,%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    1504 <c_model2_sr+0x3e4>
mov    -0x68(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<c_model2_sr>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %r8,%rbx
mov    (%rdi),%r9d
mov    (%rdx),%edx
mov    %edx,%eax
imul   %r9d,%eax
lea    0xf(,%rax,8),%r8
and    $0xfffffffffffffff0,%r8
mov    %rsp,%rax
sub    %r8,%rax
mov    %rax,-0x80(%rbp)
mov    %rax,%rsp
mov    %rsp,%r14
lea    0xf(,%r9,4),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%r14
mov    %r14,%rsp
mov    %rsp,%r12
lea    0xf(,%r9,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%r12
mov    %r12,%rsp
test   %edx,%edx
jle    1547 <c_model2_sr+0x417>
mov    0x28(%rbp),%r8
mov    0x18(%rbp),%rax
mov    %rax,-0x88(%rbp)
movsd  0x8(%r8),%xmm0
addsd  0xe46(%rip),%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x44(%rbp)
xorps  %xmm0,%xmm0
cvtsi2sdl (%rsi),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0x10(%r8),%xmm0
xorps  0xe33(%rip),%xmm0
movaps %xmm0,-0xd0(%rbp)
movsd  (%r8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movslq %r9d,%rax
lea    -0x1(%r9),%r8
lea    -0x2(%r9),%rsi
mov    %rsi,-0xa8(%rbp)
mov    %eax,%r11d
and    $0xfffffffe,%r11d
mov    %r8d,%r10d
and    $0x3,%r10d
mov    %r8,-0x38(%rbp)
and    $0xfffffffffffffffc,%r8
mov    %r8,-0x98(%rbp)
mov    %eax,%edi
and    $0x7,%edi
mov    %eax,%r8d
and    $0xfffffff8,%r8d
shl    $0x3,%rax
mov    %rax,-0x90(%rbp)
xor    %eax,%eax
mov    %rax,-0x30(%rbp)
mov    %rdx,-0x78(%rbp)
mov    %r9,-0x70(%rbp)
mov    %rcx,-0x68(%rbp)
mov    %r11,-0x60(%rbp)
mov    %r10,-0x58(%rbp)
mov    %rdi,-0x50(%rbp)
mov    %r8,-0xa0(%rbp)
jmp    1286 <c_model2_sr+0x156>
xchg   %ax,%ax
mov    -0x30(%rbp),%rsi
add    $0x1,%rsi
mov    -0x90(%rbp),%rax
add    %rax,-0x88(%rbp)
mov    %rsi,%rax
mov    %rsi,-0x30(%rbp)
cmp    %rdx,%rsi
je     1523 <c_model2_sr+0x3f3>
test   %r9d,%r9d
jle    1260 <c_model2_sr+0x130>
xor    %r13d,%r13d
jmp    12ad <c_model2_sr+0x17d>
mov    %r13,%rax
imul   %rdx,%rax
add    -0x30(%rbp),%rax
mov    -0x80(%rbp),%rsi
movsd  %xmm0,(%rsi,%rax,8)
add    $0x1,%r13
cmp    %r9,%r13
je     1260 <c_model2_sr+0x130>
cmpq   $0x0,-0x38(%rbp)
je     12f0 <c_model2_sr+0x1c0>
xor    %eax,%eax
jmp    12ce <c_model2_sr+0x19e>
nopl   0x0(%rax,%rax,1)
mov    %esi,0x4(%r14,%rax,4)
add    $0x2,%rax
cmp    %rax,%r11
je     12f2 <c_model2_sr+0x1c2>
xor    %esi,%esi
mov    $0x0,%edx
cmp    %r13,%rax
ja     12df <c_model2_sr+0x1af>
cvttsd2si (%rcx,%rax,8),%edx
mov    %edx,(%r14,%rax,4)
cmp    %r13,%rax
jae    12c0 <c_model2_sr+0x190>
cvttsd2si 0x8(%rcx,%rax,8),%esi
jmp    12c0 <c_model2_sr+0x190>
xor    %eax,%eax
test   $0x1,%r9b
je     1308 <c_model2_sr+0x1d8>
xor    %edx,%edx
cmp    %r13,%rax
ja     1304 <c_model2_sr+0x1d4>
cvttsd2si (%rcx,%rax,8),%edx
mov    %edx,(%r14,%rax,4)
xorps  %xmm2,%xmm2
cvtsi2sdl (%r14),%xmm2
movsd  %xmm2,(%r12)
cmp    $0x2,%r9d
jl     13ca <c_model2_sr+0x29a>
mov    $0x1,%eax
movapd %xmm2,%xmm0
cmpq   $0x3,-0xa8(%rbp)
03
jb     13a1 <c_model2_sr+0x271>
movapd %xmm2,%xmm0
xor    %eax,%eax
mov    -0x98(%rbp),%rdx
xorps  %xmm1,%xmm1
cvtsi2sdl 0x4(%r14,%rax,4),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%r12,%rax,8)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x8(%r14,%rax,4),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x10(%r12,%rax,8)
xorps  %xmm1,%xmm1
cvtsi2sdl 0xc(%r14,%rax,4),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%r12,%rax,8)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x10(%r14,%rax,4),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%r12,%rax,8)
add    $0x4,%rax
cmp    %rax,%rdx
jne    1340 <c_model2_sr+0x210>
add    $0x1,%rax
test   %r10,%r10
je     13ca <c_model2_sr+0x29a>
lea    (%r12,%rax,8),%rsi
lea    (%r14,%rax,4),%rax
xor    %edx,%edx
xorps  %xmm1,%xmm1
cvtsi2sdl (%rax,%rdx,4),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r10
jne    13b0 <c_model2_sr+0x280>
xor    %r15d,%r15d
nopl   (%rax)
movsd  %xmm2,-0xc0(%rbp)
ff
mov    -0x88(%rbp),%rax
movsd  (%rax,%r15,8),%xmm0
mulsd  -0xb0(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
cvtsi2sdl (%r14,%r15,4),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x40(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
subsd  %xmm2,%xmm0
cvttsd2si %xmm0,%esi
mov    -0x44(%rbp),%edi
call   1050 <pow@plt>
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0xc0(%rbp),%xmm0
ff
addsd  0xbd7(%rip),%xmm0
mulsd  -0xd0(%rbp),%xmm0
ff
call   1040 <exp@plt>
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,(%rbx,%r15,8)
cmp    %r15,-0x38(%rbp)
je     1460 <c_model2_sr+0x330>
movsd  0x8(%r12,%r15,8),%xmm2
add    $0x1,%r15
jmp    13d0 <c_model2_sr+0x2a0>
nop
cmpq   $0x7,-0x38(%rbp)
jae    1490 <c_model2_sr+0x360>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
mov    -0x78(%rbp),%rdx
mov    -0x70(%rbp),%r9
mov    -0x68(%rbp),%rcx
mov    -0x60(%rbp),%r11
mov    -0x58(%rbp),%r10
mov    -0x50(%rbp),%rdi
jmp    14f8 <c_model2_sr+0x3c8>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
xor    %eax,%eax
mov    -0x78(%rbp),%rdx
mov    -0x70(%rbp),%r9
mov    -0x68(%rbp),%rcx
mov    -0x60(%rbp),%r11
mov    -0x58(%rbp),%r10
mov    -0x50(%rbp),%rdi
mov    -0xa0(%rbp),%r8
cs nopw 0x0(%rax,%rax,1)
nop
addsd  (%rbx,%rax,8),%xmm0
addsd  0x8(%rbx,%rax,8),%xmm0
addsd  0x10(%rbx,%rax,8),%xmm0
addsd  0x18(%rbx,%rax,8),%xmm0
addsd  0x20(%rbx,%rax,8),%xmm0
addsd  0x28(%rbx,%rax,8),%xmm0
addsd  0x30(%rbx,%rax,8),%xmm0
addsd  0x38(%rbx,%rax,8),%xmm0
add    $0x8,%rax
cmp    %rax,%r8
jne    14c0 <c_model2_sr+0x390>
test   %rdi,%rdi
je     1290 <c_model2_sr+0x160>
lea    (%rbx,%rax,8),%rax
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
addsd  (%rax,%rsi,8),%xmm0
add    $0x1,%rsi
cmp    %rsi,%rdi
jne    1510 <c_model2_sr+0x3e0>
jmp    1290 <c_model2_sr+0x160>
test   %edx,%edx
jle    1547 <c_model2_sr+0x417>
add    $0xffffffff,%r9d
imul   %edx,%r9d
movslq %r9d,%rax
mov    -0x80(%rbp),%rcx
lea    (%rcx,%rax,8),%rsi
shl    $0x3,%rdx
mov    0x10(%rbp),%rdi
call   1030 <memcpy@plt>
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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