void histU(const double *x, const double *y, const int *knots, const double *wx,
           const double *wy, double *hist12, double *hist1, double *hist2,
           int numSamples, int splineOrder, int numBins) {
  int curSample;
  int curBin, curBinX, curBinY;
  for (curBin = 0; curBin < numBins; curBin++) {
    int nbDef = 0;
    for (curSample = 0; curSample < numSamples; curSample++)
      if ((wx[curBin * numSamples + curSample] != 1e+30) &&
          (wy[curBin * numSamples + curSample] != 1e+30))
        nbDef += 1;
    hist1[curBin] = 0.0;
    hist2[curBin] = 0.0;
    for (curSample = 0; curSample < numSamples; curSample++)
      if ((wx[curBin * numSamples + curSample] != 1e+30) &&
          (wy[curBin * numSamples + curSample] != 1e+30)) {
        hist1[curBin] += wx[curBin * numSamples + curSample] / (double)nbDef;
        hist2[curBin] += wy[curBin * numSamples + curSample] / (double)nbDef;
      }
  }
  for (curBinX = 0; curBinX < numBins; curBinX++)
    for (curBinY = 0; curBinY < numBins; curBinY++) {
      int nbDef = 0;
      for (curSample = 0; curSample < numSamples; curSample++)
        if ((wx[curBinX * numSamples + curSample] != 1e+30) &&
            (wy[curBinY * numSamples + curSample] != 1e+30))
          nbDef += 1;
      hist12[curBinX * numBins + curBinY] = 0.0;
      for (curSample = 0; curSample < numSamples; curSample++)
        if ((wx[curBinX * numSamples + curSample] != 1e+30) &&
            (wy[curBinY * numSamples + curSample] != 1e+30))
          hist12[curBinX * numBins + curBinY] +=
              wx[curBinX * numSamples + curSample] *
              wy[curBinY * numSamples + curSample] / (double)nbDef;
    }
}


x86 assembly (clang, -O0):
<histU>:
push   %rbp
mov    %rsp,%rbp
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
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)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x30(%rbp),%eax
jge    12fd <histU+0x1fd>
movl   $0x0,-0x44(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    11dc <histU+0xdc>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe88(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    118d <histU+0x8d>
jp     118d <histU+0x8d>
jmp    11c9 <histU+0xc9>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe55(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    11c0 <histU+0xc0>
jp     11c0 <histU+0xc0>
jmp    11c9 <histU+0xc9>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11ce <histU+0xce>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    114e <histU+0x4e>
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    12ea <histU+0x1ea>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xdd3(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1242 <histU+0x142>
jp     1242 <histU+0x142>
jmp    12d7 <histU+0x1d7>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xda0(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1275 <histU+0x175>
jp     1275 <histU+0x175>
jmp    12d7 <histU+0x1d7>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x44(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x44(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    12dc <histU+0x1dc>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1203 <histU+0x103>
jmp    12ef <histU+0x1ef>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1134 <histU+0x34>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x30(%rbp),%eax
jge    14da <histU+0x3da>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    0x30(%rbp),%eax
jge    14c7 <histU+0x3c7>
movl   $0x0,-0x48(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    13bf <histU+0x2bf>
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xca5(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1370 <histU+0x270>
jp     1370 <histU+0x270>
jmp    13ac <histU+0x2ac>
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xc72(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    13a3 <histU+0x2a3>
jp     13a3 <histU+0x2a3>
jmp    13ac <histU+0x2ac>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    13b1 <histU+0x2b1>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1331 <histU+0x231>
mov    -0x30(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   0x30(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    14b4 <histU+0x3b4>
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xbf7(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    141e <histU+0x31e>
jp     141e <histU+0x31e>
jmp    14a1 <histU+0x3a1>
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xbc4(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1451 <histU+0x351>
jp     1451 <histU+0x351>
jmp    14a1 <histU+0x3a1>
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   0x30(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    14a6 <histU+0x3a6>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    13df <histU+0x2df>
jmp    14b9 <histU+0x3b9>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1317 <histU+0x217>
jmp    14cc <histU+0x3cc>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1304 <histU+0x204>
pop    %rbp
ret

x86 assembly (clang, -O2):
<histU>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8,-0x38(%rsp)
mov    %rcx,-0x30(%rsp)
mov    0x58(%rsp),%eax
test   %eax,%eax
jle    147c <histU+0x37c>
mov    0x40(%rsp),%r12
mov    0x38(%rsp),%r13
movslq 0x48(%rsp),%rcx
mov    %eax,%r10d
mov    %ecx,%edi
lea    -0x1(%rdi),%rax
mov    %rax,-0x20(%rsp)
mov    %edi,%r11d
and    $0xfffffffe,%r11d
mov    -0x38(%rsp),%r8
lea    0x8(%r8),%rax
mov    %rcx,-0x28(%rsp)
lea    0x0(,%rcx,8),%r15
mov    -0x30(%rsp),%rsi
lea    0x8(%rsi),%rbp
xor    %r14d,%r14d
movsd  0xe92(%rip),%xmm0
jmp    119a <histU+0x9a>
movq   $0x0,0x0(%r13,%r14,8)
movq   $0x0,(%r12,%r14,8)
add    $0x1,%r14
add    %r15,%rax
add    %r15,%rbp
add    %r15,%rsi
add    %r15,%r8
cmp    %r10,%r14
je     12c2 <histU+0x1c2>
cmpl   $0x0,0x48(%rsp)
jle    1170 <histU+0x70>
cmpq   $0x0,-0x20(%rsp)
je     1200 <histU+0x100>
xor    %edx,%edx
xor    %ecx,%ecx
jmp    11b9 <histU+0xb9>
nop
add    $0x2,%rdx
cmp    %rdx,%r11
je     1204 <histU+0x104>
movsd  -0x8(%rbp,%rdx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    11c7 <histU+0xc7>
jnp    11d8 <histU+0xd8>
movsd  -0x8(%rax,%rdx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    11d5 <histU+0xd5>
jnp    11d8 <histU+0xd8>
add    $0x1,%ecx
movsd  0x0(%rbp,%rdx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    11e6 <histU+0xe6>
jnp    11b0 <histU+0xb0>
movsd  (%rax,%rdx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    11f3 <histU+0xf3>
jnp    11b0 <histU+0xb0>
add    $0x1,%ecx
jmp    11b0 <histU+0xb0>
nopl   0x0(%rax,%rax,1)
xor    %edx,%edx
xor    %ecx,%ecx
test   $0x1,%dil
je     1240 <histU+0x140>
mov    %r14,%rbx
imul   -0x28(%rsp),%rbx
add    %rbx,%rdx
mov    -0x30(%rsp),%rbx
movsd  (%rbx,%rdx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1228 <histU+0x128>
jnp    1240 <histU+0x140>
mov    -0x38(%rsp),%rbx
movsd  (%rbx,%rdx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    123a <histU+0x13a>
jnp    1240 <histU+0x140>
add    $0x1,%ecx
nopl   (%rax)
movq   $0x0,0x0(%r13,%r14,8)
movq   $0x0,(%r12,%r14,8)
cmpl   $0x0,0x48(%rsp)
jle    1181 <histU+0x81>
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
xor    %ecx,%ecx
jmp    127d <histU+0x17d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rcx
cmp    %rcx,%rdi
je     1181 <histU+0x81>
movsd  (%rsi,%rcx,8),%xmm2
ucomisd %xmm0,%xmm2
jne    128a <histU+0x18a>
jnp    1270 <histU+0x170>
movsd  (%r8,%rcx,8),%xmm3
ucomisd %xmm0,%xmm3
jne    1298 <histU+0x198>
jnp    1270 <histU+0x170>
divsd  %xmm1,%xmm2
addsd  0x0(%r13,%r14,8),%xmm2
movsd  %xmm2,0x0(%r13,%r14,8)
movsd  (%r8,%rcx,8),%xmm2
divsd  %xmm1,%xmm2
addsd  (%r12,%r14,8),%xmm2
movsd  %xmm2,(%r12,%r14,8)
jmp    1270 <histU+0x170>
cmpl   $0x0,0x58(%rsp)
mov    0x48(%rsp),%r13d
jle    147c <histU+0x37c>
mov    %edi,%ebx
and    $0xfffffffe,%ebx
mov    -0x30(%rsp),%rbp
lea    0x8(%rbp),%rax
mov    -0x38(%rsp),%rcx
add    $0x8,%rcx
mov    %rcx,-0x18(%rsp)
xor    %r12d,%r12d
jmp    1318 <histU+0x218>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    -0x10(%rsp),%r12
add    $0x1,%r12
add    %r15,%rax
add    %r15,%rbp
cmp    %r10,%r12
je     147c <histU+0x37c>
mov    %r12,%rcx
imul   -0x28(%rsp),%rcx
mov    %rcx,-0x8(%rsp)
mov    %r12,-0x10(%rsp)
imul   %r10,%r12
mov    -0x38(%rsp),%rdx
mov    -0x18(%rsp),%r14
xor    %r8d,%r8d
jmp    135b <histU+0x25b>
xchg   %ax,%ax
lea    (%r8,%r12,1),%rcx
movq   $0x0,(%r9,%rcx,8)
add    $0x1,%r8
add    %r15,%r14
add    %r15,%rdx
cmp    %r10,%r8
je     1300 <histU+0x200>
test   %r13d,%r13d
jle    1340 <histU+0x240>
cmpq   $0x0,-0x20(%rsp)
je     13c0 <histU+0x2c0>
xor    %r11d,%r11d
xor    %ecx,%ecx
jmp    1379 <histU+0x279>
nop
add    $0x2,%r11
cmp    %r11,%rbx
je     13c5 <histU+0x2c5>
movsd  -0x8(%rax,%r11,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1388 <histU+0x288>
jnp    139a <histU+0x29a>
movsd  -0x8(%r14,%r11,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1397 <histU+0x297>
jnp    139a <histU+0x29a>
add    $0x1,%ecx
movsd  (%rax,%r11,8),%xmm1
ucomisd %xmm0,%xmm1
jne    13a8 <histU+0x2a8>
jnp    1370 <histU+0x270>
movsd  (%r14,%r11,8),%xmm1
ucomisd %xmm0,%xmm1
jne    13b6 <histU+0x2b6>
jnp    1370 <histU+0x270>
add    $0x1,%ecx
jmp    1370 <histU+0x270>
nopl   0x0(%rax,%rax,1)
xor    %r11d,%r11d
xor    %ecx,%ecx
test   $0x1,%dil
je     140e <histU+0x30e>
mov    -0x8(%rsp),%rsi
add    %r11,%rsi
mov    -0x30(%rsp),%r13
movsd  0x0(%r13,%rsi,8),%xmm1
ucomisd %xmm0,%xmm1
jne    13e7 <histU+0x2e7>
jnp    1409 <histU+0x309>
mov    %r8,%rsi
imul   -0x28(%rsp),%rsi
add    %rsi,%r11
mov    -0x38(%rsp),%rsi
movsd  (%rsi,%r11,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1406 <histU+0x306>
jnp    1409 <histU+0x309>
add    $0x1,%ecx
mov    0x48(%rsp),%r13d
lea    (%r8,%r12,1),%rsi
movq   $0x0,(%r9,%rsi,8)
test   %r13d,%r13d
jle    134c <histU+0x24c>
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
xorpd  %xmm2,%xmm2
xor    %ecx,%ecx
jmp    144d <histU+0x34d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rcx
cmp    %rcx,%rdi
je     134c <histU+0x24c>
movsd  0x0(%rbp,%rcx,8),%xmm3
ucomisd %xmm0,%xmm3
jne    145b <histU+0x35b>
jnp    1440 <histU+0x340>
movsd  (%rdx,%rcx,8),%xmm4
ucomisd %xmm0,%xmm4
jne    1468 <histU+0x368>
jnp    1440 <histU+0x340>
mulsd  %xmm4,%xmm3
divsd  %xmm1,%xmm3
addsd  %xmm3,%xmm2
movsd  %xmm2,(%r9,%rsi,8)
jmp    1440 <histU+0x340>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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