#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);
void oi_rr_fast(int *ng, int *no, double *xg, double *yg, double *zg,
                double *xo, double *yo, double *zo, double *Dh, double *Dz,
                double *xb, double *vec, double *xa) {
  int ngv = ng[0];
  int nov = no[0];
  double Dhv = Dh[0];
  double Dzv = Dz[0];
  int i, j;
  double g, Dh2, Dz2, hd2, vd2;
  Dh2 = Dhv * Dhv;
  Dz2 = Dzv * Dzv;
  for (i = 0; i < ngv; i++) {
    xa[i] = xb[i];
  }
  for (j = 0; j < nov; j++) {
    for (i = 0; i < ngv; i++) {
      hd2 = ((xg[i] - xo[j]) * (xg[i] - xo[j]) +
             (yg[i] - yo[j]) * (yg[i] - yo[j])) /
            (1000. * 1000.);
      vd2 = (zg[i] - zo[j]) * (zg[i] - zo[j]);
      g = exp(-0.5 * (hd2 / Dh2 + vd2 / Dz2));
      xa[i] = xa[i] + g * vec[j];
    }
  }
  for (i = 0; i < ngv; i++) {
    if (xa[i] < 0)
      xa[i] = 0.;
  }
  return;
}


x86 assembly (clang, -O0):
<oi_rr_fast>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
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    0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11d4 <oi_rr_fast+0xc4>
mov    0x30(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11a0 <oi_rr_fast+0x90>
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    133f <oi_rr_fast+0x22f>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    132c <oi_rr_fast+0x21c>
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x30(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm3
mov    0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xd92(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    0x18(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x70(%rbp),%xmm1
divsd  -0x60(%rbp),%xmm1
movsd  -0x78(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
addsd  %xmm0,%xmm1
movsd  0xd24(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1030 <exp@plt>
movsd  %xmm0,-0x58(%rbp)
mov    0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x58(%rbp),%xmm0
mov    0x38(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11ee <oi_rr_fast+0xde>
jmp    1331 <oi_rr_fast+0x221>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    11db <oi_rr_fast+0xcb>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    138b <oi_rr_fast+0x27b>
mov    0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    1378 <oi_rr_fast+0x268>
mov    0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    137d <oi_rr_fast+0x26d>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1346 <oi_rr_fast+0x236>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<oi_rr_fast>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rcx,%r15
mov    0x90(%rsp),%rbx
mov    0x78(%rsp),%rax
mov    0x70(%rsp),%rbp
mov    (%rdi),%r13d
mov    (%rsi),%esi
mov    %rsi,(%rsp)
movsd  0x0(%rbp),%xmm2
movsd  (%rax),%xmm3
test   %r13d,%r13d
jle    12ca <oi_rr_fast+0x1ba>
mov    0x80(%rsp),%r14
cmp    $0x4,%r13d
jb     1176 <oi_rr_fast+0x66>
lea    (%r14,%r13,8),%rsi
cmp    %rbx,%rsi
jbe    11f0 <oi_rr_fast+0xe0>
lea    (%rbx,%r13,8),%rsi
cmp    %r14,%rsi
jbe    11f0 <oi_rr_fast+0xe0>
xor    %ebp,%ebp
mov    %rbp,%r10
not    %r10
add    %r13,%r10
mov    %r13,%rsi
and    $0x3,%rsi
je     11a5 <oi_rr_fast+0x95>
nopw   0x0(%rax,%rax,1)
movsd  (%r14,%rbp,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
add    $0x1,%rbp
add    $0xffffffffffffffff,%rsi
jne    1190 <oi_rr_fast+0x80>
cmp    $0x3,%r10
jb     12ca <oi_rr_fast+0x1ba>
nop
movsd  (%r14,%rbp,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
movsd  0x8(%r14,%rbp,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rbp,8)
movsd  0x10(%r14,%rbp,8),%xmm0
movsd  %xmm0,0x10(%rbx,%rbp,8)
movsd  0x18(%r14,%rbp,8),%xmm0
movsd  %xmm0,0x18(%rbx,%rbp,8)
add    $0x4,%rbp
cmp    %rbp,%r13
jne    11b0 <oi_rr_fast+0xa0>
jmp    12ca <oi_rr_fast+0x1ba>
mov    %r13d,%ebp
and    $0xfffffffc,%ebp
lea    -0x4(%rbp),%r11
mov    %r11,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%r10d
and    $0x3,%r10d
cmp    $0xc,%r11
jae    1216 <oi_rr_fast+0x106>
xor    %eax,%eax
jmp    1284 <oi_rr_fast+0x174>
and    $0xfffffffffffffffc,%rsi
xor    %eax,%eax
nopl   0x0(%rax)
movups (%r14,%rax,8),%xmm0
movups 0x10(%r14,%rax,8),%xmm1
movups %xmm0,(%rbx,%rax,8)
movups %xmm1,0x10(%rbx,%rax,8)
movups 0x20(%r14,%rax,8),%xmm0
movups 0x30(%r14,%rax,8),%xmm1
movups %xmm0,0x20(%rbx,%rax,8)
movups %xmm1,0x30(%rbx,%rax,8)
movups 0x40(%r14,%rax,8),%xmm0
movups 0x50(%r14,%rax,8),%xmm1
movups %xmm0,0x40(%rbx,%rax,8)
movups %xmm1,0x50(%rbx,%rax,8)
movupd 0x60(%r14,%rax,8),%xmm0
movupd 0x70(%r14,%rax,8),%xmm1
movupd %xmm0,0x60(%rbx,%rax,8)
movupd %xmm1,0x70(%rbx,%rax,8)
add    $0x10,%rax
add    $0xfffffffffffffffc,%rsi
jne    1220 <oi_rr_fast+0x110>
test   %r10,%r10
je     12c1 <oi_rr_fast+0x1b1>
lea    (%rbx,%rax,8),%r11
add    $0x10,%r11
lea    (%r14,%rax,8),%rax
add    $0x10,%rax
shl    $0x5,%r10
xor    %esi,%esi
nop
movupd -0x10(%rax,%rsi,1),%xmm0
movupd (%rax,%rsi,1),%xmm1
movupd %xmm0,-0x10(%r11,%rsi,1)
movupd %xmm1,(%r11,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%r10
jne    12a0 <oi_rr_fast+0x190>
cmp    %r13,%rbp
jne    1178 <oi_rr_fast+0x68>
cmpl   $0x0,(%rsp)
jle    13e2 <oi_rr_fast+0x2d2>
mov    0x68(%rsp),%rax
mov    0x60(%rsp),%rsi
mulsd  %xmm2,%xmm2
mulsd  %xmm3,%xmm3
xor    %edi,%edi
movsd  0xd10(%rip),%xmm4
movsd  0xd10(%rip),%xmm5
mov    %rdx,0x18(%rsp)
movsd  %xmm2,0x10(%rsp)
movsd  %xmm3,0x8(%rsp)
jmp    131e <oi_rr_fast+0x20e>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    (%rsp),%rdi
je     13e2 <oi_rr_fast+0x2d2>
test   %r13d,%r13d
jle    1310 <oi_rr_fast+0x200>
xor    %r12d,%r12d
mov    %rdi,0x20(%rsp)
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%r12,8),%xmm1
subsd  (%r9,%rdi,8),%xmm1
movsd  (%r15,%r12,8),%xmm0
subsd  (%rsi,%rdi,8),%xmm0
mulsd  %xmm0,%xmm0
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
divsd  %xmm4,%xmm1
divsd  %xmm2,%xmm1
movsd  (%r8,%r12,8),%xmm0
subsd  (%rax,%rdi,8),%xmm0
mulsd  %xmm0,%xmm0
divsd  %xmm3,%xmm0
addsd  %xmm1,%xmm0
mulsd  %xmm5,%xmm0
mov    %r9,%rbp
mov    %r8,%r14
call   1030 <exp@plt>
movsd  0xc7f(%rip),%xmm5
movsd  0xc6f(%rip),%xmm4
mov    0x20(%rsp),%rdi
mov    0x60(%rsp),%rsi
mov    0x88(%rsp),%rax
movsd  0x8(%rsp),%xmm3
movsd  0x10(%rsp),%xmm2
mov    %r14,%r8
mov    0x18(%rsp),%rdx
mov    %rbp,%r9
mulsd  (%rax,%rdi,8),%xmm0
mov    0x68(%rsp),%rax
addsd  (%rbx,%r12,8),%xmm0
movsd  %xmm0,(%rbx,%r12,8)
add    $0x1,%r12
cmp    %r12,%r13
jne    1330 <oi_rr_fast+0x220>
jmp    1310 <oi_rr_fast+0x200>
test   %r13d,%r13d
jle    1408 <oi_rr_fast+0x2f8>
cmp    $0x1,%r13d
jne    1417 <oi_rr_fast+0x307>
xor    %eax,%eax
test   $0x1,%r13b
je     1408 <oi_rr_fast+0x2f8>
xorpd  %xmm0,%xmm0
ucomisd (%rbx,%rax,8),%xmm0
jbe    1408 <oi_rr_fast+0x2f8>
movq   $0x0,(%rbx,%rax,8)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %r13d,%ecx
and    $0xfffffffe,%ecx
xor    %eax,%eax
xorpd  %xmm0,%xmm0
jmp    1439 <oi_rr_fast+0x329>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x2,%rax
cmp    %rax,%rcx
je     13ef <oi_rr_fast+0x2df>
ucomisd (%rbx,%rax,8),%xmm0
jbe    1448 <oi_rr_fast+0x338>
movq   $0x0,(%rbx,%rax,8)
ucomisd 0x8(%rbx,%rax,8),%xmm0
jbe    1430 <oi_rr_fast+0x320>
movq   $0x0,0x8(%rbx,%rax,8)
jmp    1430 <oi_rr_fast+0x320>

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

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