double GaussTransform(const double *A, const double *B, int m, int n, int dim,
                      double scale, double *grad) {
  int i, j, d;
  int id, jd;
  double dist_ij, cross_term = 0;
  double cost_ij;
  for (i = 0; i < m * dim; ++i)
    grad[i] = 0;
  for (i = 0; i < m; ++i) {
    for (j = 0; j < n; ++j) {
      dist_ij = 0;
      for (d = 0; d < dim; ++d) {
        id = i * dim + d;
        jd = j * dim + d;
        dist_ij = dist_ij + ((A[id] - B[jd]) * (A[id] - B[jd]));
      }
      cost_ij = exp(-dist_ij / ((scale) * (scale)));
      for (d = 0; d < dim; ++d) {
        id = i * dim + d;
        jd = j * dim + d;
        grad[id] += -cost_ij * 2 * (A[id] - B[jd]);
      }
      cross_term += cost_ij;
    }
  }
  for (i = 0; i < m * dim; ++i) {
    grad[i] /= (m * n * ((scale) * (scale)));
  }
  return cross_term / (m * n);
}


x86 assembly (clang, -O0):
<GaussTransform>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    1172 <GaussTransform+0x62>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1142 <GaussTransform+0x32>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    131c <GaussTransform+0x20c>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1309 <GaussTransform+0x1f9>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1221 <GaussTransform+0x111>
mov    -0x34(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x44(%rbp)
movsd  -0x50(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11a7 <GaussTransform+0x97>
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1030 <exp@plt>
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12ec <GaussTransform+0x1dc>
mov    -0x34(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x44(%rbp)
movsd  -0x60(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xd5a(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    125c <GaussTransform+0x14c>
movsd  -0x60(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    118c <GaussTransform+0x7c>
jmp    130e <GaussTransform+0x1fe>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1179 <GaussTransform+0x69>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    1372 <GaussTransform+0x262>
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  -0x28(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
mulsd  %xmm0,%xmm1
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1323 <GaussTransform+0x213>
movsd  -0x58(%rbp),%xmm0
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<GaussTransform>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
mov    %r9,%r12
movapd %xmm0,%xmm4
mov    %ecx,0x8(%rsp)
mov    %edx,%ebp
mov    %rsi,0x20(%rsp)
mov    %rdi,0x18(%rsp)
mov    %r8d,%ebx
imul   %edx,%ebx
mov    %ebx,%r14d
test   %ebx,%ebx
movapd %xmm0,0x50(%rsp)
mov    %r8d,0xc(%rsp)
jle    117d <GaussTransform+0x5d>
lea    0x0(,%r14,8),%rdx
mov    %r12,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    0xc(%rsp),%r8d
movapd 0x50(%rsp),%xmm4
test   %ebp,%ebp
mov    %r12,0x40(%rsp)
mov    %r14,0x88(%rsp)
mov    %ebx,0x4c(%rsp)
mov    %ebp,0x48(%rsp)
jle    1629 <GaussTransform+0x509>
movapd %xmm4,%xmm0
mulsd  %xmm4,%xmm0
movslq %r8d,%rax
mov    %ebp,%ecx
mov    %rcx,0xa0(%rsp)
mov    0x8(%rsp),%ecx
mov    %rcx,0xf0(%rsp)
movsd  0xe3d(%rip),%xmm1
movsd  %xmm0,0xe8(%rsp)
divsd  %xmm0,%xmm1
movsd  %xmm1,0xe0(%rsp)
mov    %eax,%r10d
mov    %r10d,%ecx
and    $0xfffffffe,%ecx
mov    %rcx,0xd0(%rsp)
mov    %r10d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,0xa8(%rsp)
mov    %r10,%rcx
neg    %rcx
mov    %rcx,0xd8(%rsp)
mov    0x18(%rsp),%r13
lea    0x8(%r13),%r14
mov    %rax,0x70(%rsp)
lea    0x0(,%rax,8),%rax
mov    %rax,0x80(%rsp)
mov    0x20(%rsp),%rax
lea    0x8(%rax),%rcx
mov    %rcx,0x98(%rsp)
lea    0x10(%r12),%r15
lea    0x10(%r13),%rdi
add    $0x10,%rax
mov    %rax,0x90(%rsp)
lea    0x0(,%r10,8),%rax
and    $0xffffffffffffffe0,%rax
mov    %rax,0xb8(%rsp)
xorpd  %xmm5,%xmm5
xor    %ecx,%ecx
mov    %r12,%rbp
mov    %r10,0x28(%rsp)
jmp    12a1 <GaussTransform+0x181>
xchg   %ax,%ax
mov    0xb0(%rsp),%rcx
add    $0x1,%rcx
mov    0x80(%rsp),%rax
add    %rax,%r14
add    %rax,%r15
add    %rax,%rdi
add    %rax,%r13
add    %rax,%rbp
cmp    0xa0(%rsp),%rcx
je     162d <GaussTransform+0x50d>
cmpl   $0x0,0x8(%rsp)
mov    %rcx,0xb0(%rsp)
jle    1270 <GaussTransform+0x150>
imul   0x70(%rsp),%rcx
lea    (%rcx,%r10,1),%rax
mov    0x40(%rsp),%rdx
lea    (%rdx,%rcx,8),%rbx
lea    (%rdx,%rax,8),%r9
mov    0x18(%rsp),%rdx
mov    %rcx,0x68(%rsp)
lea    (%rdx,%rcx,8),%rcx
lea    (%rdx,%rax,8),%rax
mov    %rbx,0xc8(%rsp)
cmp    %rax,%rbx
setb   %al
mov    %r9,0xc0(%rsp)
cmp    %r9,%rcx
setb   %cl
and    %al,%cl
mov    %cl,0x7(%rsp)
mov    0x90(%rsp),%rbx
mov    0x98(%rsp),%r12
xor    %esi,%esi
mov    %rdi,0x30(%rsp)
jmp    1372 <GaussTransform+0x252>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  0xe0(%rsp),%xmm0
call   1040 <exp@plt>
mov    0x38(%rsp),%rsi
mov    0x30(%rsp),%rdi
mov    0x28(%rsp),%r10
mov    0xc(%rsp),%r8d
movsd  0x10(%rsp),%xmm5
movapd 0x50(%rsp),%xmm4
addsd  %xmm0,%xmm5
add    $0x1,%rsi
mov    0x80(%rsp),%rax
add    %rax,%r12
add    %rax,%rbx
cmp    0xf0(%rsp),%rsi
je     1270 <GaussTransform+0x150>
test   %r8d,%r8d
movsd  %xmm5,0x10(%rsp)
mov    %rsi,0x38(%rsp)
jle    1320 <GaussTransform+0x200>
cmp    $0x1,%r8d
jne    1390 <GaussTransform+0x270>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    13d7 <GaussTransform+0x2b7>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
mov    0xd0(%rsp),%rcx
xchg   %ax,%ax
movsd  -0x8(%r14,%rax,8),%xmm2
movsd  (%r14,%rax,8),%xmm0
subsd  -0x8(%r12,%rax,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
subsd  (%r12,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
add    $0x2,%rax
movapd %xmm0,%xmm1
cmp    %rax,%rcx
jne    13a0 <GaussTransform+0x280>
imul   0x70(%rsp),%rsi
test   $0x1,%r10b
je     140a <GaussTransform+0x2ea>
mov    0x68(%rsp),%rcx
add    %rax,%rcx
add    %rsi,%rax
mov    0x18(%rsp),%rdx
movsd  (%rdx,%rcx,8),%xmm1
mov    0x20(%rsp),%rcx
subsd  (%rcx,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
mov    %rsi,0x78(%rsp)
xorpd  0xbf9(%rip),%xmm0
divsd  0xe8(%rsp),%xmm0
call   1040 <exp@plt>
mov    0xc(%rsp),%r8d
test   %r8d,%r8d
jle    1470 <GaussTransform+0x350>
movapd %xmm0,%xmm1
mulsd  0xbcd(%rip),%xmm1
cmp    $0x4,%r8d
jae    1490 <GaussTransform+0x370>
xor    %ecx,%ecx
movapd 0x50(%rsp),%xmm4
movsd  0x10(%rsp),%xmm5
mov    0x28(%rsp),%r10
mov    0x30(%rsp),%rdi
mov    0x38(%rsp),%rsi
mov    0x78(%rsp),%r9
jmp    1580 <GaussTransform+0x460>
nopl   0x0(%rax,%rax,1)
movapd 0x50(%rsp),%xmm4
movsd  0x10(%rsp),%xmm5
mov    0x28(%rsp),%r10
mov    0x30(%rsp),%rdi
mov    0x38(%rsp),%rsi
jmp    134e <GaussTransform+0x22e>
mov    0x28(%rsp),%r10
mov    0x78(%rsp),%r9
lea    (%r9,%r10,1),%rax
mov    0x20(%rsp),%rdx
lea    (%rdx,%r9,8),%rcx
lea    (%rdx,%rax,8),%rax
cmp    %rax,0xc8(%rsp)
setb   %al
cmp    0xc0(%rsp),%rcx
setb   %cl
and    %al,%cl
or     0x7(%rsp),%cl
mov    0x30(%rsp),%rdi
mov    0xb8(%rsp),%rcx
mov    0x38(%rsp),%rsi
je     14ee <GaussTransform+0x3ce>
xor    %ecx,%ecx
movapd 0x50(%rsp),%xmm4
movsd  0x10(%rsp),%xmm5
jmp    1580 <GaussTransform+0x460>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movupd -0x10(%rdi,%rax,1),%xmm3
movupd (%rdi,%rax,1),%xmm4
movupd -0x10(%rbx,%rax,1),%xmm5
subpd  %xmm5,%xmm3
movupd (%rbx,%rax,1),%xmm5
subpd  %xmm5,%xmm4
movupd -0x10(%r15,%rax,1),%xmm5
movupd (%r15,%rax,1),%xmm6
mulpd  %xmm2,%xmm3
addpd  %xmm5,%xmm3
mulpd  %xmm2,%xmm4
addpd  %xmm6,%xmm4
movupd %xmm3,-0x10(%r15,%rax,1)
movupd %xmm4,(%r15,%rax,1)
add    $0x20,%rax
cmp    %rax,%rcx
jne    1500 <GaussTransform+0x3e0>
mov    0xa8(%rsp),%rax
mov    %rax,%rcx
cmp    %r10,%rax
movapd 0x50(%rsp),%xmm4
movsd  0x10(%rsp),%xmm5
je     134e <GaussTransform+0x22e>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rcx,%rax
test   $0x1,%r10b
je     15c3 <GaussTransform+0x4a3>
mov    0x68(%rsp),%rax
add    %rcx,%rax
add    %rcx,%r9
mov    0x18(%rsp),%rdx
movsd  (%rdx,%rax,8),%xmm2
mov    0x20(%rsp),%rdx
subsd  (%rdx,%r9,8),%xmm2
mulsd  %xmm1,%xmm2
mov    0x40(%rsp),%rdx
addsd  (%rdx,%rax,8),%xmm2
movsd  %xmm2,(%rdx,%rax,8)
mov    %rcx,%rax
or     $0x1,%rax
not    %rcx
cmp    0xd8(%rsp),%rcx
je     134e <GaussTransform+0x22e>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  0x0(%r13,%rax,8),%xmm2
subsd  -0x8(%r12,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x0(%rbp,%rax,8),%xmm2
movsd  %xmm2,0x0(%rbp,%rax,8)
movsd  0x8(%r13,%rax,8),%xmm2
subsd  (%r12,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x8(%rbp,%rax,8),%xmm2
movsd  %xmm2,0x8(%rbp,%rax,8)
add    $0x2,%rax
cmp    %rax,%r10
jne    15e0 <GaussTransform+0x4c0>
jmp    134e <GaussTransform+0x22e>
xorpd  %xmm5,%xmm5
mov    0x8(%rsp),%eax
imul   0x48(%rsp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    0x4c(%rsp),%eax
test   %eax,%eax
mov    0x40(%rsp),%rdi
mov    0x88(%rsp),%rbp
jle    16f7 <GaussTransform+0x5d7>
mulsd  %xmm4,%xmm4
mulsd  %xmm0,%xmm4
cmp    $0x1,%eax
jne    1667 <GaussTransform+0x547>
xor    %eax,%eax
jmp    16e0 <GaussTransform+0x5c0>
mov    %ebp,%eax
and    $0xfffffffe,%eax
movapd %xmm4,%xmm1
unpcklpd %xmm4,%xmm1
lea    -0x2(%rax),%rdx
mov    %rdx,%rcx
shr    %rcx
add    $0x1,%rcx
test   %rdx,%rdx
je     1711 <GaussTransform+0x5f1>
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%rdi,%rdx,8),%xmm2
movupd 0x10(%rdi,%rdx,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdi,%rdx,8)
divpd  %xmm1,%xmm3
movupd %xmm3,0x10(%rdi,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    16a0 <GaussTransform+0x580>
test   $0x1,%cl
je     16db <GaussTransform+0x5bb>
movupd (%rdi,%rdx,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdi,%rdx,8)
cmp    %rbp,%rax
je     16f7 <GaussTransform+0x5d7>
movsd  (%rdi,%rax,8),%xmm1
divsd  %xmm4,%xmm1
movsd  %xmm1,(%rdi,%rax,8)
add    $0x1,%rax
cmp    %rax,%rbp
jne    16e0 <GaussTransform+0x5c0>
divsd  %xmm0,%xmm5
movapd %xmm5,%xmm0
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%cl
jne    16cd <GaussTransform+0x5ad>
jmp    16db <GaussTransform+0x5bb>

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

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