void ComputeSecondMoment4D(float *ima, float *atlas, float *means,
                           float *Ameans, float *variance, float *Avariance,
                           const int *dims, int f) {
  int indice = 0;
  int i, j, k, t;
  int ii, jj, kk;
  int ni, nj, nk;
  float var = 0.0;
  float Avar = 0.0;
  for (k = 0; k < dims[2]; k++) {
    for (i = 0; i < dims[1]; i++) {
      for (j = 0; j < dims[0]; j++) {
        Avar = 0;
        indice = 0;
        for (ii = -f; ii <= f; ii++) {
          for (jj = -f; jj <= f; jj++) {
            for (kk = -f; kk <= f; kk++) {
              ni = i + ii;
              nj = j + jj;
              nk = k + kk;
              if (ni >= 0 && nj >= 0 && nk > 0 && ni < dims[1] &&
                  nj < dims[0] && nk < dims[2]) {
                Avar = Avar +
                       ((atlas[nk * (dims[0] * dims[1]) + (ni * dims[0]) + nj] -
                         Ameans[k * (dims[0] * dims[1]) + (i * dims[0]) + j]) *
                        (atlas[nk * (dims[0] * dims[1]) + (ni * dims[0]) + nj] -
                         Ameans[k * (dims[0] * dims[1]) + (i * dims[0]) + j]));
                indice = indice + 1;
              }
            }
          }
        }
        Avar = Avar / (indice - 1);
        Avariance[k * (dims[0] * dims[1]) + (i * dims[0]) + j] = Avar;
      }
    }
  }
  for (t = 0; t < dims[3]; t++) {
    for (k = 0; k < dims[2]; k++) {
      for (i = 0; i < dims[1]; i++) {
        for (j = 0; j < dims[0]; j++) {
          var = 0;
          indice = 0;
          for (ii = -f; ii <= f; ii++) {
            for (jj = -f; jj <= f; jj++) {
              for (kk = -f; kk <= f; kk++) {
                ni = i + ii;
                nj = j + jj;
                nk = k + kk;
                if (ni >= 0 && nj >= 0 && nk > 0 && ni < dims[1] &&
                    nj < dims[0] && nk < dims[2]) {
                  var = var +
                        ((ima[t * (dims[0] * dims[1] * dims[2]) +
                              nk * (dims[0] * dims[1]) + (ni * dims[0]) + nj] -
                          means[t * (dims[0] * dims[1] * dims[2]) +
                                k * (dims[0] * dims[1]) + (i * dims[0]) + j]) *
                         (ima[t * (dims[0] * dims[1] * dims[2]) +
                              nk * (dims[0] * dims[1]) + (ni * dims[0]) + nj] -
                          means[t * (dims[0] * dims[1] * dims[2]) +
                                k * (dims[0] * dims[1]) + (i * dims[0]) + j]));
                  indice = indice + 1;
                }
              }
            }
          }
          var = var / (indice - 1);
          variance[t * (dims[0] * dims[1] * dims[2]) + k * (dims[0] * dims[1]) +
                   (i * dims[0]) + j] = var;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<ComputeSecondMoment4D>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%eax
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,-0x34(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    13c1 <ComputeSecondMoment4D+0x2c1>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    13ae <ComputeSecondMoment4D+0x2ae>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    139b <ComputeSecondMoment4D+0x29b>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movl   $0x0,-0x34(%rbp)
xor    %eax,%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     1341 <ComputeSecondMoment4D+0x241>
xor    %eax,%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     132e <ComputeSecondMoment4D+0x22e>
xor    %eax,%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     131b <ComputeSecondMoment4D+0x21b>
mov    -0x38(%rbp),%eax
add    -0x48(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x40(%rbp),%eax
add    -0x50(%rbp),%eax
mov    %eax,-0x5c(%rbp)
cmpl   $0x0,-0x54(%rbp)
jl     1308 <ComputeSecondMoment4D+0x208>
cmpl   $0x0,-0x58(%rbp)
jl     1308 <ComputeSecondMoment4D+0x208>
cmpl   $0x0,-0x5c(%rbp)
jle    1308 <ComputeSecondMoment4D+0x208>
mov    -0x54(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    1308 <ComputeSecondMoment4D+0x208>
mov    -0x58(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    1308 <ComputeSecondMoment4D+0x208>
mov    -0x5c(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    1308 <ComputeSecondMoment4D+0x208>
movss  -0x64(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
imul   %edx,%ecx
mov    -0x54(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x58(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
imul   %edx,%ecx
mov    -0x38(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
imul   %edx,%ecx
mov    -0x54(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x58(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
imul   %edx,%ecx
mov    -0x38(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x64(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    130d <ComputeSecondMoment4D+0x20d>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    11bd <ComputeSecondMoment4D+0xbd>
jmp    1320 <ComputeSecondMoment4D+0x220>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11a9 <ComputeSecondMoment4D+0xa9>
jmp    1333 <ComputeSecondMoment4D+0x233>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1195 <ComputeSecondMoment4D+0x95>
movss  -0x64(%rbp),%xmm0
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x64(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
imul   %edx,%ecx
mov    -0x38(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    116f <ComputeSecondMoment4D+0x6f>
jmp    13a0 <ComputeSecondMoment4D+0x2a0>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1158 <ComputeSecondMoment4D+0x58>
jmp    13b3 <ComputeSecondMoment4D+0x2b3>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1141 <ComputeSecondMoment4D+0x41>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0xc(%rcx),%eax
jge    1708 <ComputeSecondMoment4D+0x608>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    16f5 <ComputeSecondMoment4D+0x5f5>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    16e2 <ComputeSecondMoment4D+0x5e2>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    16cf <ComputeSecondMoment4D+0x5cf>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x60(%rbp)
movl   $0x0,-0x34(%rbp)
xor    %eax,%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     1657 <ComputeSecondMoment4D+0x557>
xor    %eax,%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     1644 <ComputeSecondMoment4D+0x544>
xor    %eax,%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    0x18(%rbp),%eax
jg     1631 <ComputeSecondMoment4D+0x531>
mov    -0x38(%rbp),%eax
add    -0x48(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x40(%rbp),%eax
add    -0x50(%rbp),%eax
mov    %eax,-0x5c(%rbp)
cmpl   $0x0,-0x54(%rbp)
jl     161e <ComputeSecondMoment4D+0x51e>
cmpl   $0x0,-0x58(%rbp)
jl     161e <ComputeSecondMoment4D+0x51e>
cmpl   $0x0,-0x5c(%rbp)
jle    161e <ComputeSecondMoment4D+0x51e>
mov    -0x54(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    161e <ComputeSecondMoment4D+0x51e>
mov    -0x58(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    161e <ComputeSecondMoment4D+0x51e>
mov    -0x5c(%rbp),%eax
mov    0x10(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    161e <ComputeSecondMoment4D+0x51e>
movss  -0x60(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
mov    0x10(%rbp),%rsi
imul   0x8(%rsi),%edx
imul   %edx,%ecx
mov    -0x5c(%rbp),%edx
mov    0x10(%rbp),%rsi
mov    (%rsi),%esi
mov    0x10(%rbp),%rdi
imul   0x4(%rdi),%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x54(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x58(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
mov    0x10(%rbp),%rsi
imul   0x8(%rsi),%edx
imul   %edx,%ecx
mov    -0x40(%rbp),%edx
mov    0x10(%rbp),%rsi
mov    (%rsi),%esi
mov    0x10(%rbp),%rdi
imul   0x4(%rdi),%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x38(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
mov    0x10(%rbp),%rsi
imul   0x8(%rsi),%edx
imul   %edx,%ecx
mov    -0x5c(%rbp),%edx
mov    0x10(%rbp),%rsi
mov    (%rsi),%esi
mov    0x10(%rbp),%rdi
imul   0x4(%rdi),%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x54(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x58(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
mov    0x10(%rbp),%rsi
imul   0x8(%rsi),%edx
imul   %edx,%ecx
mov    -0x40(%rbp),%edx
mov    0x10(%rbp),%rsi
mov    (%rsi),%esi
mov    0x10(%rbp),%rdi
imul   0x4(%rdi),%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x38(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x60(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1623 <ComputeSecondMoment4D+0x523>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    145b <ComputeSecondMoment4D+0x35b>
jmp    1636 <ComputeSecondMoment4D+0x536>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1447 <ComputeSecondMoment4D+0x347>
jmp    1649 <ComputeSecondMoment4D+0x549>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1433 <ComputeSecondMoment4D+0x333>
movss  -0x60(%rbp),%xmm0
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x60(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    0x10(%rbp),%rdx
mov    (%rdx),%edx
mov    0x10(%rbp),%rsi
imul   0x4(%rsi),%edx
mov    0x10(%rbp),%rsi
imul   0x8(%rsi),%edx
imul   %edx,%ecx
mov    -0x40(%rbp),%edx
mov    0x10(%rbp),%rsi
mov    (%rsi),%esi
mov    0x10(%rbp),%rdi
imul   0x4(%rdi),%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x38(%rbp),%edx
mov    0x10(%rbp),%rsi
imul   (%rsi),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    140d <ComputeSecondMoment4D+0x30d>
jmp    16d4 <ComputeSecondMoment4D+0x5d4>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    13f6 <ComputeSecondMoment4D+0x2f6>
jmp    16e7 <ComputeSecondMoment4D+0x5e7>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    13df <ComputeSecondMoment4D+0x2df>
jmp    16fa <ComputeSecondMoment4D+0x5fa>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    13c8 <ComputeSecondMoment4D+0x2c8>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ComputeSecondMoment4D>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r9,0x8(%rsp)
mov    %r8,0x40(%rsp)
mov    %rdx,0x50(%rsp)
mov    %rsi,-0x18(%rsp)
mov    %rdi,0x48(%rsp)
mov    %rcx,%rdi
mov    0x90(%rsp),%rax
movslq 0x8(%rax),%rax
mov    %eax,%ecx
mov    %rcx,-0x10(%rsp)
mov    %rax,0x28(%rsp)
test   %rax,%rax
jle    13bc <ComputeSecondMoment4D+0x2bc>
mov    0x90(%rsp),%rax
mov    0x4(%rax),%ebp
mov    0x98(%rsp),%eax
mov    %eax,%ecx
neg    %ecx
mov    %ecx,-0x80(%rsp)
movslq %ecx,%rcx
mov    %rcx,-0x20(%rsp)
cmovs  %eax,%ecx
mov    %ebp,%edx
imul   %eax,%edx
neg    %edx
mov    %edx,-0x38(%rsp)
mov    %rcx,-0x48(%rsp)
add    %ecx,%eax
add    $0x1,%eax
mov    %eax,0x20(%rsp)
xor    %eax,%eax
mov    %rax,-0x50(%rsp)
jmp    11c6 <ComputeSecondMoment4D+0xc6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x50(%rsp),%rcx
add    $0x1,%rcx
add    %ebp,-0x38(%rsp)
addq   $0x1,-0x20(%rsp)
mov    %rcx,%rax
mov    %rcx,-0x50(%rsp)
cmp    -0x10(%rsp),%rcx
je     13bc <ComputeSecondMoment4D+0x2bc>
test   %ebp,%ebp
jle    11a0 <ComputeSecondMoment4D+0xa0>
mov    0x90(%rsp),%rax
mov    (%rax),%r13d
mov    %r13d,%r11d
imul   %ebp,%r11d
mov    %r11d,%ecx
mov    -0x50(%rsp),%rax
imul   %eax,%ecx
mov    %ecx,(%rsp)
mov    %ebp,%ecx
imul   %eax,%ecx
mov    %rcx,-0x30(%rsp)
mov    %r13d,%eax
imul   -0x38(%rsp),%eax
mov    %r13d,%edx
not    %edx
mov    0x98(%rsp),%ecx
imul   %ecx,%edx
add    %eax,%edx
mov    %edx,-0x78(%rsp)
xor    %eax,%eax
mov    %rax,-0x70(%rsp)
jmp    123d <ComputeSecondMoment4D+0x13d>
nopw   0x0(%rax,%rax,1)
mov    -0x70(%rsp),%rcx
add    $0x1,%ecx
add    %r13d,-0x78(%rsp)
mov    %rcx,%rax
mov    %rcx,-0x70(%rsp)
cmp    %ebp,%ecx
je     11a0 <ComputeSecondMoment4D+0xa0>
test   %r13d,%r13d
jle    1220 <ComputeSecondMoment4D+0x120>
mov    -0x30(%rsp),%rax
mov    -0x70(%rsp),%rcx
add    %ecx,%eax
imul   %r13d,%eax
mov    %rax,0x10(%rsp)
mov    %r13d,%eax
imul   %ecx,%eax
mov    %rax,-0x28(%rsp)
mov    -0x78(%rsp),%esi
xor    %eax,%eax
mov    %rax,-0x60(%rsp)
jmp    12b7 <ComputeSecondMoment4D+0x1b7>
nop
xor    %r9d,%r9d
xorps  %xmm0,%xmm0
add    $0xffffffff,%r9d
xorps  %xmm1,%xmm1
cvtsi2ss %r9d,%xmm1
divss  %xmm1,%xmm0
mov    0x10(%rsp),%rax
mov    -0x60(%rsp),%rdx
add    %edx,%eax
cltq
mov    0x8(%rsp),%rcx
movss  %xmm0,(%rcx,%rax,4)
add    $0x1,%rdx
mov    0x18(%rsp),%esi
add    $0x1,%esi
mov    %rdx,-0x60(%rsp)
cmp    %r13,%rdx
je     1220 <ComputeSecondMoment4D+0x120>
mov    0x98(%rsp),%eax
test   %eax,%eax
mov    %esi,0x18(%rsp)
js     1270 <ComputeSecondMoment4D+0x170>
mov    -0x28(%rsp),%rax
mov    -0x60(%rsp),%rdx
add    %edx,%eax
add    (%rsp),%eax
movslq %eax,%r14
xor    %r9d,%r9d
xorps  %xmm0,%xmm0
mov    -0x80(%rsp),%eax
mov    %eax,%ecx
jmp    130b <ComputeSecondMoment4D+0x20b>
cs nopw 0x0(%rax,%rax,1)
mov    -0x68(%rsp),%rcx
lea    0x1(%rcx),%eax
mov    -0x7c(%rsp),%esi
add    %r13d,%esi
cmp    -0x48(%rsp),%ecx
mov    %eax,%ecx
je     1276 <ComputeSecondMoment4D+0x176>
mov    -0x70(%rsp),%rax
mov    %rcx,-0x68(%rsp)
lea    (%rcx,%rax,1),%r15d
mov    %esi,-0x7c(%rsp)
mov    -0x80(%rsp),%eax
mov    %eax,%ecx
jmp    1343 <ComputeSecondMoment4D+0x243>
cs nopw 0x0(%rax,%rax,1)
nop
mov    -0x58(%rsp),%rcx
lea    0x1(%rcx),%eax
add    $0x1,%esi
cmp    -0x48(%rsp),%ecx
mov    %eax,%ecx
je     12f0 <ComputeSecondMoment4D+0x1f0>
mov    -0x60(%rsp),%rax
mov    %rcx,-0x58(%rsp)
lea    (%rcx,%rax,1),%ebx
mov    %r15d,%eax
or     %ebx,%eax
setns  %r8b
mov    0x20(%rsp),%r12d
mov    -0x20(%rsp),%r10
mov    %esi,%ecx
jmp    137d <ComputeSecondMoment4D+0x27d>
nopw   0x0(%rax,%rax,1)
add    %r11d,%ecx
add    $0x1,%r10
add    $0xffffffff,%r12d
je     1330 <ComputeSecondMoment4D+0x230>
test   %r10,%r10
setg   %al
and    %r8b,%al
cmp    $0x1,%al
jne    1370 <ComputeSecondMoment4D+0x270>
cmp    %ebp,%r15d
jge    1370 <ComputeSecondMoment4D+0x270>
cmp    %r13d,%ebx
jge    1370 <ComputeSecondMoment4D+0x270>
cmp    -0x10(%rsp),%r10
jge    1370 <ComputeSecondMoment4D+0x270>
movslq %ecx,%rax
mov    -0x18(%rsp),%rdx
movss  (%rdx,%rax,4),%xmm1
subss  (%rdi,%r14,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%r9d
jmp    1370 <ComputeSecondMoment4D+0x270>
mov    0x90(%rsp),%rax
mov    0xc(%rax),%eax
mov    %eax,0x34(%rsp)
test   %eax,%eax
jle    169c <ComputeSecondMoment4D+0x59c>
mov    0x98(%rsp),%eax
mov    %eax,%ecx
neg    %ecx
movslq %ecx,%rdx
mov    %rdx,0x38(%rsp)
cmovs  %eax,%edx
mov    0x28(%rsp),%esi
mov    %rsi,-0x50(%rsp)
mov    %rdx,-0x18(%rsp)
add    %edx,%eax
add    $0x1,%eax
mov    %eax,-0x60(%rsp)
mov    -0x10(%rsp),%rax
xor    %edx,%edx
mov    %rdx,-0x8(%rsp)
mov    %ecx,-0x7c(%rsp)
mov    %ecx,-0x3c(%rsp)
mov    0x48(%rsp),%rdi
jmp    1440 <ComputeSecondMoment4D+0x340>
nopl   0x0(%rax,%rax,1)
mov    -0x8(%rsp),%rcx
add    $0x1,%ecx
mov    -0x10(%rsp),%rax
add    %eax,-0x3c(%rsp)
mov    %rcx,-0x8(%rsp)
cmp    0x34(%rsp),%ecx
je     169c <ComputeSecondMoment4D+0x59c>
test   %eax,%eax
jle    1420 <ComputeSecondMoment4D+0x320>
mov    0x90(%rsp),%rax
mov    0x4(%rax),%r9d
mov    -0x8(%rsp),%rax
imul   0x28(%rsp),%eax
mov    %rax,-0x38(%rsp)
mov    %r9d,%eax
imul   -0x3c(%rsp),%eax
mov    %eax,-0x30(%rsp)
mov    0x38(%rsp),%rax
mov    %rax,-0x48(%rsp)
xor    %eax,%eax
mov    %rax,-0x78(%rsp)
jmp    14a7 <ComputeSecondMoment4D+0x3a7>
xchg   %ax,%ax
mov    -0x78(%rsp),%rcx
add    $0x1,%rcx
add    %r9d,-0x30(%rsp)
addq   $0x1,-0x48(%rsp)
mov    %rcx,%rax
mov    %rcx,-0x78(%rsp)
cmp    -0x50(%rsp),%rcx
je     1420 <ComputeSecondMoment4D+0x320>
test   %r9d,%r9d
jle    1480 <ComputeSecondMoment4D+0x380>
mov    0x90(%rsp),%rax
mov    (%rax),%r15d
mov    %r15d,%r12d
imul   %r9d,%r12d
mov    -0x38(%rsp),%rax
mov    %eax,%ecx
imul   %r12d,%ecx
mov    %ecx,0x10(%rsp)
mov    %r12d,%edx
mov    -0x78(%rsp),%rcx
imul   %ecx,%edx
mov    %edx,0x8(%rsp)
add    %ecx,%eax
imul   %r9d,%eax
mov    %rax,(%rsp)
mov    %r15d,%eax
imul   -0x30(%rsp),%eax
mov    %r15d,%edx
not    %edx
mov    0x98(%rsp),%ecx
imul   %ecx,%edx
add    %eax,%edx
mov    %edx,-0x28(%rsp)
xor    %eax,%eax
mov    %rax,-0x68(%rsp)
jmp    152e <ComputeSecondMoment4D+0x42e>
nopl   0x0(%rax)
mov    -0x68(%rsp),%rcx
add    $0x1,%ecx
add    %r15d,-0x28(%rsp)
mov    %rcx,%rax
mov    %rcx,-0x68(%rsp)
cmp    %r9d,%ecx
je     1480 <ComputeSecondMoment4D+0x380>
test   %r15d,%r15d
jle    1510 <ComputeSecondMoment4D+0x410>
mov    (%rsp),%rax
mov    -0x68(%rsp),%rcx
add    %ecx,%eax
imul   %r15d,%eax
mov    %rax,-0x70(%rsp)
mov    %r15d,%eax
imul   %ecx,%eax
mov    %rax,0x18(%rsp)
mov    -0x28(%rsp),%ecx
xor    %eax,%eax
mov    %rax,-0x58(%rsp)
jmp    15a7 <ComputeSecondMoment4D+0x4a7>
nop
xor    %r10d,%r10d
xorps  %xmm0,%xmm0
add    $0xffffffff,%r10d
xorps  %xmm1,%xmm1
cvtsi2ss %r10d,%xmm1
divss  %xmm1,%xmm0
mov    -0x70(%rsp),%rax
mov    -0x58(%rsp),%rdx
add    %edx,%eax
cltq
mov    0x40(%rsp),%rcx
movss  %xmm0,(%rcx,%rax,4)
add    $0x1,%rdx
mov    -0x80(%rsp),%ecx
add    $0x1,%ecx
mov    %rdx,-0x58(%rsp)
cmp    %r15,%rdx
je     1510 <ComputeSecondMoment4D+0x410>
mov    0x98(%rsp),%eax
test   %eax,%eax
mov    %ecx,-0x80(%rsp)
js     1560 <ComputeSecondMoment4D+0x460>
mov    0x18(%rsp),%rax
mov    -0x58(%rsp),%rdx
add    %edx,%eax
add    0x8(%rsp),%eax
add    0x10(%rsp),%eax
movslq %eax,%r14
xor    %r10d,%r10d
xorps  %xmm0,%xmm0
mov    -0x7c(%rsp),%eax
mov    %eax,%edx
jmp    15fb <ComputeSecondMoment4D+0x4fb>
nopl   0x0(%rax,%rax,1)
mov    0x20(%rsp),%rdx
lea    0x1(%rdx),%eax
mov    -0x20(%rsp),%ecx
add    %r15d,%ecx
cmp    -0x18(%rsp),%edx
mov    %eax,%edx
je     1566 <ComputeSecondMoment4D+0x466>
mov    -0x68(%rsp),%rax
mov    %rdx,0x20(%rsp)
lea    (%rdx,%rax,1),%r11d
mov    %ecx,-0x20(%rsp)
mov    %ecx,%eax
mov    -0x7c(%rsp),%ecx
mov    %ecx,%ebx
jmp    162e <ComputeSecondMoment4D+0x52e>
nopw   0x0(%rax,%rax,1)
lea    0x1(%rbx),%ecx
add    $0x1,%eax
cmp    -0x18(%rsp),%ebx
mov    %ecx,%ebx
je     15e0 <ComputeSecondMoment4D+0x4e0>
mov    -0x58(%rsp),%rcx
lea    (%rbx,%rcx,1),%ebp
mov    %r11d,%ecx
or     %ebp,%ecx
setns  %cl
mov    -0x60(%rsp),%r13d
mov    -0x48(%rsp),%r8
mov    %eax,%esi
jmp    165d <ComputeSecondMoment4D+0x55d>
nopl   0x0(%rax)
add    %r12d,%esi
add    $0x1,%r8
add    $0xffffffff,%r13d
je     1620 <ComputeSecondMoment4D+0x520>
test   %r8,%r8
setg   %dl
and    %cl,%dl
cmp    $0x1,%dl
jne    1650 <ComputeSecondMoment4D+0x550>
cmp    %r9d,%r11d
jge    1650 <ComputeSecondMoment4D+0x550>
cmp    %r15d,%ebp
jge    1650 <ComputeSecondMoment4D+0x550>
cmp    0x28(%rsp),%r8
jge    1650 <ComputeSecondMoment4D+0x550>
movslq %esi,%rdx
movss  (%rdi,%rdx,4),%xmm1
mov    0x50(%rsp),%rdx
subss  (%rdx,%r14,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%r10d
jmp    1650 <ComputeSecondMoment4D+0x550>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @ComputeSecondMoment4D(float* %0, float* %1, float* %2, float* %3, float* %4, float* %5, i32* %6, i32 %7) {
%9 = alloca float*
%10 = alloca float*
%11 = alloca float*
%12 = alloca float*
%13 = alloca float*
%14 = alloca float*
%15 = alloca i32*
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%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 float
%29 = alloca float
store float* %0, float** %9
store float* %1, float** %10
store float* %2, float** %11
store float* %3, float** %12
store float* %4, float** %13
store float* %5, float** %14
store i32* %6, i32** %15
store i32 %7, i32* %16
store i32 0, i32* %17
store float 0.000000e+00, float* %28
store float 0.000000e+00, float* %29
store i32 0, i32* %20
br label %30
30: ; preds = %245, %8
%31 = load i32, i32* %20
%32 = load i32*, i32** %15
%33 = getelementptr inbounds i32, i32* %32, i64 2
%34 = load i32, i32* %33
%35 = icmp slt i32 %31, %34
br i1 %35, label %36, label %248
36: ; preds = %30
store i32 0, i32* %18
br label %37
37: ; preds = %241, %36
%38 = load i32, i32* %18
%39 = load i32*, i32** %15
%40 = getelementptr inbounds i32, i32* %39, i64 1
%41 = load i32, i32* %40
%42 = icmp slt i32 %38, %41
br i1 %42, label %43, label %244
43: ; preds = %37
store i32 0, i32* %19
br label %44
44: ; preds = %237, %43
%45 = load i32, i32* %19
%46 = load i32*, i32** %15
%47 = getelementptr inbounds i32, i32* %46, i64 0
%48 = load i32, i32* %47
%49 = icmp slt i32 %45, %48
br i1 %49, label %50, label %240
50: ; preds = %44
store float 0.000000e+00, float* %29
store i32 0, i32* %17
%51 = load i32, i32* %16
%52 = sub nsw i32 0, %51
store i32 %52, i32* %22
br label %53
53: ; preds = %207, %50
%54 = load i32, i32* %22
%55 = load i32, i32* %16
%56 = icmp sle i32 %54, %55
br i1 %56, label %57, label %210
57: ; preds = %53
%58 = load i32, i32* %16
%59 = sub nsw i32 0, %58
store i32 %59, i32* %23
br label %60
60: ; preds = %203, %57
%61 = load i32, i32* %23
%62 = load i32, i32* %16
%63 = icmp sle i32 %61, %62
br i1 %63, label %64, label %206
64: ; preds = %60
%65 = load i32, i32* %16
%66 = sub nsw i32 0, %65
store i32 %66, i32* %24
br label %67
67: ; preds = %199, %64
%68 = load i32, i32* %24
%69 = load i32, i32* %16
%70 = icmp sle i32 %68, %69
br i1 %70, label %71, label %202
71: ; preds = %67
%72 = load i32, i32* %18
%73 = load i32, i32* %22
%74 = add nsw i32 %72, %73
store i32 %74, i32* %25
%75 = load i32, i32* %19
%76 = load i32, i32* %23
%77 = add nsw i32 %75, %76
store i32 %77, i32* %26
%78 = load i32, i32* %20
%79 = load i32, i32* %24
%80 = add nsw i32 %78, %79
store i32 %80, i32* %27
%81 = load i32, i32* %25
%82 = icmp sge i32 %81, 0
br i1 %82, label %83, label %198
83: ; preds = %71
%84 = load i32, i32* %26
%85 = icmp sge i32 %84, 0
br i1 %85, label %86, label %198
86: ; preds = %83
%87 = load i32, i32* %27
%88 = icmp sgt i32 %87, 0
br i1 %88, label %89, label %198
89: ; preds = %86
%90 = load i32, i32* %25
%91 = load i32*, i32** %15
%92 = getelementptr inbounds i32, i32* %91, i64 1
%93 = load i32, i32* %92
%94 = icmp slt i32 %90, %93
br i1 %94, label %95, label %198
95: ; preds = %89
%96 = load i32, i32* %26
%97 = load i32*, i32** %15
%98 = getelementptr inbounds i32, i32* %97, i64 0
%99 = load i32, i32* %98
%100 = icmp slt i32 %96, %99
br i1 %100, label %101, label %198
101: ; preds = %95
%102 = load i32, i32* %27
%103 = load i32*, i32** %15
%104 = getelementptr inbounds i32, i32* %103, i64 2
%105 = load i32, i32* %104
%106 = icmp slt i32 %102, %105
br i1 %106, label %107, label %198
107: ; preds = %101
%108 = load float, float* %29
%109 = load float*, float** %10
%110 = load i32, i32* %27
%111 = load i32*, i32** %15
%112 = getelementptr inbounds i32, i32* %111, i64 0
%113 = load i32, i32* %112
%114 = load i32*, i32** %15
%115 = getelementptr inbounds i32, i32* %114, i64 1
%116 = load i32, i32* %115
%117 = mul nsw i32 %113, %116
%118 = mul nsw i32 %110, %117
%119 = load i32, i32* %25
%120 = load i32*, i32** %15
%121 = getelementptr inbounds i32, i32* %120, i64 0
%122 = load i32, i32* %121
%123 = mul nsw i32 %119, %122
%124 = add nsw i32 %118, %123
%125 = load i32, i32* %26
%126 = add nsw i32 %124, %125
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds float, float* %109, i64 %127
%129 = load float, float* %128
%130 = load float*, float** %12
%131 = load i32, i32* %20
%132 = load i32*, i32** %15
%133 = getelementptr inbounds i32, i32* %132, i64 0
%134 = load i32, i32* %133
%135 = load i32*, i32** %15
%136 = getelementptr inbounds i32, i32* %135, i64 1
%137 = load i32, i32* %136
%138 = mul nsw i32 %134, %137
%139 = mul nsw i32 %131, %138
%140 = load i32, i32* %18
%141 = load i32*, i32** %15
%142 = getelementptr inbounds i32, i32* %141, i64 0
%143 = load i32, i32* %142
%144 = mul nsw i32 %140, %143
%145 = add nsw i32 %139, %144
%146 = load i32, i32* %19
%147 = add nsw i32 %145, %146
%148 = sext i32 %147 to i64
%149 = getelementptr inbounds float, float* %130, i64 %148
%150 = load float, float* %149
%151 = fsub float %129, %150
%152 = load float*, float** %10
%153 = load i32, i32* %27
%154 = load i32*, i32** %15
%155 = getelementptr inbounds i32, i32* %154, i64 0
%156 = load i32, i32* %155
%157 = load i32*, i32** %15
%158 = getelementptr inbounds i32, i32* %157, i64 1
%159 = load i32, i32* %158
%160 = mul nsw i32 %156, %159
%161 = mul nsw i32 %153, %160
%162 = load i32, i32* %25
%163 = load i32*, i32** %15
%164 = getelementptr inbounds i32, i32* %163, i64 0
%165 = load i32, i32* %164
%166 = mul nsw i32 %162, %165
%167 = add nsw i32 %161, %166
%168 = load i32, i32* %26
%169 = add nsw i32 %167, %168
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds float, float* %152, i64 %170
%172 = load float, float* %171
%173 = load float*, float** %12
%174 = load i32, i32* %20
%175 = load i32*, i32** %15
%176 = getelementptr inbounds i32, i32* %175, i64 0
%177 = load i32, i32* %176
%178 = load i32*, i32** %15
%179 = getelementptr inbounds i32, i32* %178, i64 1
%180 = load i32, i32* %179
%181 = mul nsw i32 %177, %180
%182 = mul nsw i32 %174, %181
%183 = load i32, i32* %18
%184 = load i32*, i32** %15
%185 = getelementptr inbounds i32, i32* %184, i64 0
%186 = load i32, i32* %185
%187 = mul nsw i32 %183, %186
%188 = add nsw i32 %182, %187
%189 = load i32, i32* %19
%190 = add nsw i32 %188, %189
%191 = sext i32 %190 to i64
%192 = getelementptr inbounds float, float* %173, i64 %191
%193 = load float, float* %192
%194 = fsub float %172, %193
%195 = call float @llvm.fmuladd.f32(float %151, float %194, float %108)
store float %195, float* %29
%196 = load i32, i32* %17
%197 = add nsw i32 %196, 1
store i32 %197, i32* %17
br label %198
198: ; preds = %107, %101, %95, %89, %86, %83, %71
br label %199
199: ; preds = %198
%200 = load i32, i32* %24
%201 = add nsw i32 %200, 1
store i32 %201, i32* %24
br label %67
202: ; preds = %67
br label %203
203: ; preds = %202
%204 = load i32, i32* %23
%205 = add nsw i32 %204, 1
store i32 %205, i32* %23
br label %60
206: ; preds = %60
br label %207
207: ; preds = %206
%208 = load i32, i32* %22
%209 = add nsw i32 %208, 1
store i32 %209, i32* %22
br label %53
210: ; preds = %53
%211 = load float, float* %29
%212 = load i32, i32* %17
%213 = sub nsw i32 %212, 1
%214 = sitofp i32 %213 to float
%215 = fdiv float %211, %214
store float %215, float* %29
%216 = load float, float* %29
%217 = load float*, float** %14
%218 = load i32, i32* %20
%219 = load i32*, i32** %15
%220 = getelementptr inbounds i32, i32* %219, i64 0
%221 = load i32, i32* %220
%222 = load i32*, i32** %15
%223 = getelementptr inbounds i32, i32* %222, i64 1
%224 = load i32, i32* %223
%225 = mul nsw i32 %221, %224
%226 = mul nsw i32 %218, %225
%227 = load i32, i32* %18
%228 = load i32*, i32** %15
%229 = getelementptr inbounds i32, i32* %228, i64 0
%230 = load i32, i32* %229
%231 = mul nsw i32 %227, %230
%232 = add nsw i32 %226, %231
%233 = load i32, i32* %19
%234 = add nsw i32 %232, %233
%235 = sext i32 %234 to i64
%236 = getelementptr inbounds float, float* %217, i64 %235
store float %216, float* %236
br label %237
237: ; preds = %210
%238 = load i32, i32* %19
%239 = add nsw i32 %238, 1
store i32 %239, i32* %19
br label %44
240: ; preds = %44
br label %241
241: ; preds = %240
%242 = load i32, i32* %18
%243 = add nsw i32 %242, 1
store i32 %243, i32* %18
br label %37
244: ; preds = %37
br label %245
245: ; preds = %244
%246 = load i32, i32* %20
%247 = add nsw i32 %246, 1
store i32 %247, i32* %20
br label %30
248: ; preds = %30
store i32 0, i32* %21
br label %249
249: ; preds = %545, %248
%250 = load i32, i32* %21
%251 = load i32*, i32** %15
%252 = getelementptr inbounds i32, i32* %251, i64 3
%253 = load i32, i32* %252
%254 = icmp slt i32 %250, %253
br i1 %254, label %255, label %548
255: ; preds = %249
store i32 0, i32* %20
br label %256
256: ; preds = %541, %255
%257 = load i32, i32* %20
%258 = load i32*, i32** %15
%259 = getelementptr inbounds i32, i32* %258, i64 2
%260 = load i32, i32* %259
%261 = icmp slt i32 %257, %260
br i1 %261, label %262, label %544
262: ; preds = %256
store i32 0, i32* %18
br label %263
263: ; preds = %537, %262
%264 = load i32, i32* %18
%265 = load i32*, i32** %15
%266 = getelementptr inbounds i32, i32* %265, i64 1
%267 = load i32, i32* %266
%268 = icmp slt i32 %264, %267
br i1 %268, label %269, label %540
269: ; preds = %263
store i32 0, i32* %19
br label %270
270: ; preds = %533, %269
%271 = load i32, i32* %19
%272 = load i32*, i32** %15
%273 = getelementptr inbounds i32, i32* %272, i64 0
%274 = load i32, i32* %273
%275 = icmp slt i32 %271, %274
br i1 %275, label %276, label %536
276: ; preds = %270
store float 0.000000e+00, float* %28
store i32 0, i32* %17
%277 = load i32, i32* %16
%278 = sub nsw i32 0, %277
store i32 %278, i32* %22
br label %279
279: ; preds = %489, %276
%280 = load i32, i32* %22
%281 = load i32, i32* %16
%282 = icmp sle i32 %280, %281
br i1 %282, label %283, label %492
283: ; preds = %279
%284 = load i32, i32* %16
%285 = sub nsw i32 0, %284
store i32 %285, i32* %23
br label %286
286: ; preds = %485, %283
%287 = load i32, i32* %23
%288 = load i32, i32* %16
%289 = icmp sle i32 %287, %288
br i1 %289, label %290, label %488
290: ; preds = %286
%291 = load i32, i32* %16
%292 = sub nsw i32 0, %291
store i32 %292, i32* %24
br label %293
293: ; preds = %481, %290
%294 = load i32, i32* %24
%295 = load i32, i32* %16
%296 = icmp sle i32 %294, %295
br i1 %296, label %297, label %484
297: ; preds = %293
%298 = load i32, i32* %18
%299 = load i32, i32* %22
%300 = add nsw i32 %298, %299
store i32 %300, i32* %25
%301 = load i32, i32* %19
%302 = load i32, i32* %23
%303 = add nsw i32 %301, %302
store i32 %303, i32* %26
%304 = load i32, i32* %20
%305 = load i32, i32* %24
%306 = add nsw i32 %304, %305
store i32 %306, i32* %27
%307 = load i32, i32* %25
%308 = icmp sge i32 %307, 0
br i1 %308, label %309, label %480
309: ; preds = %297
%310 = load i32, i32* %26
%311 = icmp sge i32 %310, 0
br i1 %311, label %312, label %480
312: ; preds = %309
%313 = load i32, i32* %27
%314 = icmp sgt i32 %313, 0
br i1 %314, label %315, label %480
315: ; preds = %312
%316 = load i32, i32* %25
%317 = load i32*, i32** %15
%318 = getelementptr inbounds i32, i32* %317, i64 1
%319 = load i32, i32* %318
%320 = icmp slt i32 %316, %319
br i1 %320, label %321, label %480
321: ; preds = %315
%322 = load i32, i32* %26
%323 = load i32*, i32** %15
%324 = getelementptr inbounds i32, i32* %323, i64 0
%325 = load i32, i32* %324
%326 = icmp slt i32 %322, %325
br i1 %326, label %327, label %480
327: ; preds = %321
%328 = load i32, i32* %27
%329 = load i32*, i32** %15
%330 = getelementptr inbounds i32, i32* %329, i64 2
%331 = load i32, i32* %330
%332 = icmp slt i32 %328, %331
br i1 %332, label %333, label %480
333: ; preds = %327
%334 = load float, float* %28
%335 = load float*, float** %9
%336 = load i32, i32* %21
%337 = load i32*, i32** %15
%338 = getelementptr inbounds i32, i32* %337, i64 0
%339 = load i32, i32* %338
%340 = load i32*, i32** %15
%341 = getelementptr inbounds i32, i32* %340, i64 1
%342 = load i32, i32* %341
%343 = mul nsw i32 %339, %342
%344 = load i32*, i32** %15
%345 = getelementptr inbounds i32, i32* %344, i64 2
%346 = load i32, i32* %345
%347 = mul nsw i32 %343, %346
%348 = mul nsw i32 %336, %347
%349 = load i32, i32* %27
%350 = load i32*, i32** %15
%351 = getelementptr inbounds i32, i32* %350, i64 0
%352 = load i32, i32* %351
%353 = load i32*, i32** %15
%354 = getelementptr inbounds i32, i32* %353, i64 1
%355 = load i32, i32* %354
%356 = mul nsw i32 %352, %355
%357 = mul nsw i32 %349, %356
%358 = add nsw i32 %348, %357
%359 = load i32, i32* %25
%360 = load i32*, i32** %15
%361 = getelementptr inbounds i32, i32* %360, i64 0
%362 = load i32, i32* %361
%363 = mul nsw i32 %359, %362
%364 = add nsw i32 %358, %363
%365 = load i32, i32* %26
%366 = add nsw i32 %364, %365
%367 = sext i32 %366 to i64
%368 = getelementptr inbounds float, float* %335, i64 %367
%369 = load float, float* %368
%370 = load float*, float** %11
%371 = load i32, i32* %21
%372 = load i32*, i32** %15
%373 = getelementptr inbounds i32, i32* %372, i64 0
%374 = load i32, i32* %373
%375 = load i32*, i32** %15
%376 = getelementptr inbounds i32, i32* %375, i64 1
%377 = load i32, i32* %376
%378 = mul nsw i32 %374, %377
%379 = load i32*, i32** %15
%380 = getelementptr inbounds i32, i32* %379, i64 2
%381 = load i32, i32* %380
%382 = mul nsw i32 %378, %381
%383 = mul nsw i32 %371, %382
%384 = load i32, i32* %20
%385 = load i32*, i32** %15
%386 = getelementptr inbounds i32, i32* %385, i64 0
%387 = load i32, i32* %386
%388 = load i32*, i32** %15
%389 = getelementptr inbounds i32, i32* %388, i64 1
%390 = load i32, i32* %389
%391 = mul nsw i32 %387, %390
%392 = mul nsw i32 %384, %391
%393 = add nsw i32 %383, %392
%394 = load i32, i32* %18
%395 = load i32*, i32** %15
%396 = getelementptr inbounds i32, i32* %395, i64 0
%397 = load i32, i32* %396
%398 = mul nsw i32 %394, %397
%399 = add nsw i32 %393, %398
%400 = load i32, i32* %19
%401 = add nsw i32 %399, %400
%402 = sext i32 %401 to i64
%403 = getelementptr inbounds float, float* %370, i64 %402
%404 = load float, float* %403
%405 = fsub float %369, %404
%406 = load float*, float** %9
%407 = load i32, i32* %21
%408 = load i32*, i32** %15
%409 = getelementptr inbounds i32, i32* %408, i64 0
%410 = load i32, i32* %409
%411 = load i32*, i32** %15
%412 = getelementptr inbounds i32, i32* %411, i64 1
%413 = load i32, i32* %412
%414 = mul nsw i32 %410, %413
%415 = load i32*, i32** %15
%416 = getelementptr inbounds i32, i32* %415, i64 2
%417 = load i32, i32* %416
%418 = mul nsw i32 %414, %417
%419 = mul nsw i32 %407, %418
%420 = load i32, i32* %27
%421 = load i32*, i32** %15
%422 = getelementptr inbounds i32, i32* %421, i64 0
%423 = load i32, i32* %422
%424 = load i32*, i32** %15
%425 = getelementptr inbounds i32, i32* %424, i64 1
%426 = load i32, i32* %425
%427 = mul nsw i32 %423, %426
%428 = mul nsw i32 %420, %427
%429 = add nsw i32 %419, %428
%430 = load i32, i32* %25
%431 = load i32*, i32** %15
%432 = getelementptr inbounds i32, i32* %431, i64 0
%433 = load i32, i32* %432
%434 = mul nsw i32 %430, %433
%435 = add nsw i32 %429, %434
%436 = load i32, i32* %26
%437 = add nsw i32 %435, %436
%438 = sext i32 %437 to i64
%439 = getelementptr inbounds float, float* %406, i64 %438
%440 = load float, float* %439
%441 = load float*, float** %11
%442 = load i32, i32* %21
%443 = load i32*, i32** %15
%444 = getelementptr inbounds i32, i32* %443, i64 0
%445 = load i32, i32* %444
%446 = load i32*, i32** %15
%447 = getelementptr inbounds i32, i32* %446, i64 1
%448 = load i32, i32* %447
%449 = mul nsw i32 %445, %448
%450 = load i32*, i32** %15
%451 = getelementptr inbounds i32, i32* %450, i64 2
%452 = load i32, i32* %451
%453 = mul nsw i32 %449, %452
%454 = mul nsw i32 %442, %453
%455 = load i32, i32* %20
%456 = load i32*, i32** %15
%457 = getelementptr inbounds i32, i32* %456, i64 0
%458 = load i32, i32* %457
%459 = load i32*, i32** %15
%460 = getelementptr inbounds i32, i32* %459, i64 1
%461 = load i32, i32* %460
%462 = mul nsw i32 %458, %461
%463 = mul nsw i32 %455, %462
%464 = add nsw i32 %454, %463
%465 = load i32, i32* %18
%466 = load i32*, i32** %15
%467 = getelementptr inbounds i32, i32* %466, i64 0
%468 = load i32, i32* %467
%469 = mul nsw i32 %465, %468
%470 = add nsw i32 %464, %469
%471 = load i32, i32* %19
%472 = add nsw i32 %470, %471
%473 = sext i32 %472 to i64
%474 = getelementptr inbounds float, float* %441, i64 %473
%475 = load float, float* %474
%476 = fsub float %440, %475
%477 = call float @llvm.fmuladd.f32(float %405, float %476, float %334)
store float %477, float* %28
%478 = load i32, i32* %17
%479 = add nsw i32 %478, 1
store i32 %479, i32* %17
br label %480
480: ; preds = %333, %327, %321, %315, %312, %309, %297
br label %481
481: ; preds = %480
%482 = load i32, i32* %24
%483 = add nsw i32 %482, 1
store i32 %483, i32* %24
br label %293
484: ; preds = %293
br label %485
485: ; preds = %484
%486 = load i32, i32* %23
%487 = add nsw i32 %486, 1
store i32 %487, i32* %23
br label %286
488: ; preds = %286
br label %489
489: ; preds = %488
%490 = load i32, i32* %22
%491 = add nsw i32 %490, 1
store i32 %491, i32* %22
br label %279
492: ; preds = %279
%493 = load float, float* %28
%494 = load i32, i32* %17
%495 = sub nsw i32 %494, 1
%496 = sitofp i32 %495 to float
%497 = fdiv float %493, %496
store float %497, float* %28
%498 = load float, float* %28
%499 = load float*, float** %13
%500 = load i32, i32* %21
%501 = load i32*, i32** %15
%502 = getelementptr inbounds i32, i32* %501, i64 0
%503 = load i32, i32* %502
%504 = load i32*, i32** %15
%505 = getelementptr inbounds i32, i32* %504, i64 1
%506 = load i32, i32* %505
%507 = mul nsw i32 %503, %506
%508 = load i32*, i32** %15
%509 = getelementptr inbounds i32, i32* %508, i64 2
%510 = load i32, i32* %509
%511 = mul nsw i32 %507, %510
%512 = mul nsw i32 %500, %511
%513 = load i32, i32* %20
%514 = load i32*, i32** %15
%515 = getelementptr inbounds i32, i32* %514, i64 0
%516 = load i32, i32* %515
%517 = load i32*, i32** %15
%518 = getelementptr inbounds i32, i32* %517, i64 1
%519 = load i32, i32* %518
%520 = mul nsw i32 %516, %519
%521 = mul nsw i32 %513, %520
%522 = add nsw i32 %512, %521
%523 = load i32, i32* %18
%524 = load i32*, i32** %15
%525 = getelementptr inbounds i32, i32* %524, i64 0
%526 = load i32, i32* %525
%527 = mul nsw i32 %523, %526
%528 = add nsw i32 %522, %527
%529 = load i32, i32* %19
%530 = add nsw i32 %528, %529
%531 = sext i32 %530 to i64
%532 = getelementptr inbounds float, float* %499, i64 %531
store float %498, float* %532
br label %533
533: ; preds = %492
%534 = load i32, i32* %19
%535 = add nsw i32 %534, 1
store i32 %535, i32* %19
br label %270
536: ; preds = %270
br label %537
537: ; preds = %536
%538 = load i32, i32* %18
%539 = add nsw i32 %538, 1
store i32 %539, i32* %18
br label %263
540: ; preds = %263
br label %541
541: ; preds = %540
%542 = load i32, i32* %20
%543 = add nsw i32 %542, 1
store i32 %543, i32* %20
br label %256
544: ; preds = %256
br label %545
545: ; preds = %544
%546 = load i32, i32* %21
%547 = add nsw i32 %546, 1
store i32 %547, i32* %21
br label %249
548: ; preds = %249
ret void
}

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