void transpose_so3(double *arrayIn, double *arrayOut, int m, int n) {
  int k, d, mud;
  double *inPtr, *outPtr;
  inPtr = arrayIn;
  outPtr = arrayOut;
  if (n >= m) {
    for (d = n - m + 1; d < n; d++) {
      mud = n - d;
      for (k = 0; k < mud; k++) {
        outPtr[m * (k + d) + k] = inPtr[n * k + k + d];
      }
    }
    for (d = 0; d < n - m + 1; d++) {
      mud = m;
      for (k = 0; k < mud; k++) {
        outPtr[m * (k + d) + k] = inPtr[n * k + k + d];
      }
    }
    for (d = -m + 1; d < 0; d++) {
      mud = m + d;
      for (k = 0; k < mud; k++) {
        outPtr[m * k + k - d] = inPtr[n * (k - d) + k];
      }
    }
  } else {
    for (d = 0; d < n; d++) {
      mud = n - d;
      for (k = 0; k < mud; k++) {
        outPtr[m * (k + d) + k] = inPtr[n * k + k + d];
      }
    }
    for (d = n - m; d < 0; d++) {
      mud = n;
      for (k = 0; k < mud; k++) {
        outPtr[m * (k) + k - d] = inPtr[n * (k - d) + k];
      }
    }
    for (d = -m + 1; d < n - m; d++) {
      mud = m + d;
      for (k = 0; k < mud; k++) {
        outPtr[m * k + k - d] = inPtr[n * (k - d) + k];
      }
    }
  }
}


x86 assembly (clang, -O0):
<transpose_so3>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     12cb <transpose_so3+0x1cb>
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11b7 <transpose_so3+0xb7>
mov    -0x18(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11a4 <transpose_so3+0xa4>
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    -0x20(%rbp),%edx
imul   %edx,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1156 <transpose_so3+0x56>
jmp    11a9 <transpose_so3+0xa9>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    113a <transpose_so3+0x3a>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1240 <transpose_so3+0x140>
mov    -0x14(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    122d <transpose_so3+0x12d>
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    -0x20(%rbp),%edx
imul   %edx,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11df <transpose_so3+0xdf>
jmp    1232 <transpose_so3+0x132>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11be <transpose_so3+0xbe>
xor    %eax,%eax
sub    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jge    12c6 <transpose_so3+0x1c6>
mov    -0x14(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12b3 <transpose_so3+0x1b3>
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    -0x20(%rbp),%edx
imul   %edx,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
sub    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1265 <transpose_so3+0x165>
jmp    12b8 <transpose_so3+0x1b8>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    124b <transpose_so3+0x14b>
jmp    1462 <transpose_so3+0x362>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    134f <transpose_so3+0x24f>
mov    -0x18(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    133c <transpose_so3+0x23c>
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    -0x20(%rbp),%edx
imul   %edx,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    12ee <transpose_so3+0x1ee>
jmp    1341 <transpose_so3+0x241>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    12d2 <transpose_so3+0x1d2>
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jge    13d0 <transpose_so3+0x2d0>
mov    -0x18(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    13bd <transpose_so3+0x2bd>
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    -0x20(%rbp),%edx
imul   %edx,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
sub    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    136f <transpose_so3+0x26f>
jmp    13c2 <transpose_so3+0x2c2>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1358 <transpose_so3+0x258>
xor    %eax,%eax
sub    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
cmp    %ecx,%eax
jge    145d <transpose_so3+0x35d>
mov    -0x14(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    144a <transpose_so3+0x34a>
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    -0x20(%rbp),%edx
imul   %edx,%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x1c(%rbp),%ecx
sub    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    13fc <transpose_so3+0x2fc>
jmp    144f <transpose_so3+0x34f>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    13db <transpose_so3+0x2db>
jmp    1462 <transpose_so3+0x362>
pop    %rbp
ret

x86 assembly (clang, -O2):
<transpose_so3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %ecx,%ebp
sub    %edx,%ebp
mov    %rcx,-0x40(%rsp)
mov    %rdx,-0x38(%rsp)
jge    1234 <transpose_so3+0x134>
test   %ecx,%ecx
jle    1385 <transpose_so3+0x285>
lea    (%rdx,%rdx,1),%r13d
lea    0x1(%rcx),%r11d
lea    (%rcx,%rcx,1),%r15d
add    $0x2,%r15d
xor    %r8d,%r8d
mov    %ecx,%r10d
mov    %edx,%r14d
xor    %r12d,%r12d
jmp    1167 <transpose_so3+0x67>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r12d
add    %edx,%r8d
add    %edx,%r14d
add    $0xffffffff,%r10d
cmp    %ecx,%r12d
je     1385 <transpose_so3+0x285>
mov    %ecx,%eax
sub    %r12d,%eax
jle    1150 <transpose_so3+0x50>
cmp    $0x1,%eax
mov    %r8,-0x28(%rsp)
mov    %r14d,-0x30(%rsp)
mov    %rax,-0x20(%rsp)
jne    1190 <transpose_so3+0x90>
xor    %ebx,%ebx
jmp    11f0 <transpose_so3+0xf0>
cs nopw 0x0(%rax,%rax,1)
mov    %r10d,%r9d
and    $0xfffffffe,%r9d
mov    %r14d,%ecx
mov    %r12d,%edx
mov    %r8d,%eax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %ecx,%rcx
lea    (%rcx,%rbx,1),%r14
cltq
movslq %edx,%rdx
movsd  (%rdi,%rdx,8),%xmm0
lea    (%rbx,%rax,1),%r8
movsd  %xmm0,(%rsi,%r8,8)
lea    (%r11,%rdx,1),%ebp
movslq %ebp,%rbp
movsd  (%rdi,%rbp,8),%xmm0
movsd  %xmm0,0x8(%rsi,%r14,8)
add    $0x2,%rbx
add    %r13d,%eax
add    %r15d,%edx
add    %r13d,%ecx
cmp    %rbx,%r9
jne    11b0 <transpose_so3+0xb0>
testb  $0x1,-0x20(%rsp)
mov    -0x40(%rsp),%rcx
mov    -0x38(%rsp),%rdx
mov    -0x28(%rsp),%r8
mov    -0x30(%rsp),%r14d
je     1150 <transpose_so3+0x50>
lea    (%r12,%rbx,1),%eax
mov    %ebx,%ebp
imul   %ecx,%ebp
add    %eax,%ebp
movslq %ebp,%rbp
movsd  (%rdi,%rbp,8),%xmm0
imul   %edx,%eax
cltq
add    %rbx,%rax
movsd  %xmm0,(%rsi,%rax,8)
jmp    1150 <transpose_so3+0x50>
lea    0x1(%rbp),%r9d
cmp    %ecx,%r9d
mov    %rbp,-0x30(%rsp)
jge    14c6 <transpose_so3+0x3c6>
lea    0x2(%rbp),%r14d
imul   %edx,%r14d
lea    (%rdx,%rdx,1),%r11d
lea    (%rcx,%rcx,1),%eax
add    $0x2,%eax
mov    %eax,%r15d
sub    %edx,%r15d
mov    %edx,%r10d
imul   %r9d,%r10d
lea    -0x1(%rdx),%r13d
xor    %r12d,%r12d
jmp    128f <transpose_so3+0x18f>
xchg   %ax,%ax
add    $0x1,%r9d
add    $0x1,%r12d
add    %edx,%r14d
add    $0x1,%r15d
add    %edx,%r10d
add    $0xffffffff,%r13d
cmp    %ecx,%r9d
je     14c6 <transpose_so3+0x3c6>
cmp    %ecx,%r9d
jge    1270 <transpose_so3+0x170>
mov    %r12,-0x10(%rsp)
not    %r12d
add    %edx,%r12d
mov    %r12,-0x8(%rsp)
cmp    $0x1,%r12d
mov    %r15d,-0x28(%rsp)
mov    %r10d,-0x20(%rsp)
mov    %r13,-0x18(%rsp)
jne    12c0 <transpose_so3+0x1c0>
mov    %r14d,%ebp
xor    %ebx,%ebx
jmp    1320 <transpose_so3+0x220>
and    $0xfffffffe,%r13d
mov    %r10d,%r8d
mov    %r9d,%r10d
mov    %r15d,%ecx
mov    %r14d,%ebp
mov    %r14d,%r12d
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movslq %r12d,%r12
lea    (%r12,%rbx,1),%r14
movslq %r8d,%r8
movslq %r10d,%r10
movsd  (%rdi,%r10,8),%xmm0
lea    (%rbx,%r8,1),%r15
movsd  %xmm0,(%rsi,%r15,8)
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rsi,%r14,8)
add    $0x2,%rbx
add    %r11d,%r12d
add    %eax,%ecx
add    %eax,%r10d
add    %r11d,%r8d
cmp    %rbx,%r13
jne    12e0 <transpose_so3+0x1e0>
testb  $0x1,-0x8(%rsp)
mov    -0x40(%rsp),%rcx
mov    %ebp,%r14d
mov    -0x30(%rsp),%rbp
mov    -0x28(%rsp),%r15d
mov    -0x20(%rsp),%r10d
mov    -0x18(%rsp),%r13
mov    -0x10(%rsp),%r12
je     1270 <transpose_so3+0x170>
lea    (%r9,%rbx,1),%r8d
mov    %ebx,%ecx
imul   -0x40(%rsp),%ecx
add    %r8d,%ecx
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm0
imul   %edx,%r8d
movslq %r8d,%rbp
mov    -0x40(%rsp),%rcx
add    %rbx,%rbp
movsd  %xmm0,(%rsi,%rbp,8)
mov    -0x30(%rsp),%rbp
mov    -0x38(%rsp),%rdx
jmp    1270 <transpose_so3+0x170>
mov    %ecx,%eax
sub    %edx,%eax
mov    %eax,-0x18(%rsp)
js     15d5 <transpose_so3+0x4d5>
mov    $0x1,%r15d
mov    $0x1,%r9d
sub    -0x38(%rsp),%r9d
cmp    -0x18(%rsp),%r9d
jge    15ca <transpose_so3+0x4ca>
mov    -0x38(%rsp),%rcx
lea    (%rcx,%rcx,1),%r8d
lea    (%rcx,%rcx,1),%r12d
add    $0x2,%r12d
mov    -0x40(%rsp),%rax
mov    %eax,%r10d
imul   %ecx,%r10d
lea    (%rax,%rax,1),%r13d
add    $0xffffffff,%ecx
mov    %eax,%r14d
imul   %ecx,%r14d
jmp    13fe <transpose_so3+0x2fe>
nopl   0x0(%rax)
add    $0x1,%r9d
add    $0x1,%r15d
add    $0xffffffff,%r8d
sub    %ebp,%r10d
add    $0xffffffff,%ecx
sub    %ebp,%r14d
cmp    %ebp,%r15d
je     15ca <transpose_so3+0x4ca>
mov    %r15d,%r15d
mov    %r9d,%eax
add    -0x38(%rsp),%eax
mov    -0x40(%rsp),%rbp
jle    13e0 <transpose_so3+0x2e0>
cmp    $0x1,%r15d
mov    %r8d,-0x28(%rsp)
mov    %r10d,-0x30(%rsp)
mov    %ecx,-0x20(%rsp)
jne    1430 <transpose_so3+0x330>
xor    %eax,%eax
jmp    1481 <transpose_so3+0x381>
nopw   0x0(%rax,%rax,1)
mov    %r15d,%ebx
and    $0xfffffffe,%ebx
mov    %r14d,%edx
mov    %r10d,%ebp
xor    %eax,%eax
xchg   %ax,%ax
movslq %ebp,%rbp
lea    (%rax,%rbp,1),%r10
movslq %edx,%rdx
lea    (%rax,%rdx,1),%r11
movsd  (%rdi,%r11,8),%xmm0
movslq %ecx,%rcx
movsd  %xmm0,(%rsi,%rcx,8)
movsd  0x8(%rdi,%r10,8),%xmm0
movslq %r8d,%r8
movsd  %xmm0,(%rsi,%r8,8)
add    $0x2,%rax
add    %r12d,%r8d
add    %r13d,%ebp
add    %r12d,%ecx
add    %r13d,%edx
cmp    %rax,%rbx
jne    1440 <transpose_so3+0x340>
test   $0x1,%r15b
mov    -0x40(%rsp),%rbp
mov    -0x28(%rsp),%r8d
mov    -0x30(%rsp),%r10d
mov    -0x20(%rsp),%ecx
je     13e0 <transpose_so3+0x2e0>
mov    %eax,%ebx
sub    %r9d,%ebx
mov    %ebx,%edx
imul   %ebp,%edx
movslq %edx,%rdx
add    %rax,%rdx
movsd  (%rdi,%rdx,8),%xmm0
imul   -0x38(%rsp),%eax
add    %ebx,%eax
cltq
movsd  %xmm0,(%rsi,%rax,8)
jmp    13e0 <transpose_so3+0x2e0>
mov    %edx,%eax
mov    %rax,-0x28(%rsp)
test   %ebp,%ebp
js     15bb <transpose_so3+0x4bb>
mov    -0x28(%rsp),%r10
and    $0xfffffffe,%r10d
lea    (%rdx,%rdx,1),%r9d
mov    -0x40(%rsp),%rax
lea    0x1(%rax),%r12d
lea    (%rax,%rax,1),%r8d
add    $0x2,%r8d
xor    %r11d,%r11d
mov    %edx,%r14d
xor    %r15d,%r15d
jmp    1516 <transpose_so3+0x416>
xchg   %ax,%ax
lea    0x1(%r15),%eax
add    %edx,%r11d
add    %edx,%r14d
cmp    %ebp,%r15d
mov    %eax,%r15d
je     15bb <transpose_so3+0x4bb>
test   %edx,%edx
jle    1500 <transpose_so3+0x400>
cmp    $0x1,%edx
jne    1530 <transpose_so3+0x430>
xor    %ebx,%ebx
jmp    157f <transpose_so3+0x47f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r14d,%ecx
mov    %r15d,%edx
mov    %r11d,%eax
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movslq %ecx,%rcx
lea    (%rcx,%rbx,1),%r13
cltq
movslq %edx,%rdx
movsd  (%rdi,%rdx,8),%xmm0
lea    (%rbx,%rax,1),%rbp
movsd  %xmm0,(%rsi,%rbp,8)
lea    (%r12,%rdx,1),%ebp
movslq %ebp,%rbp
movsd  (%rdi,%rbp,8),%xmm0
movsd  %xmm0,0x8(%rsi,%r13,8)
add    $0x2,%rbx
add    %r9d,%eax
add    %r8d,%edx
add    %r9d,%ecx
cmp    %rbx,%r10
jne    1540 <transpose_so3+0x440>
testb  $0x1,-0x28(%rsp)
mov    -0x38(%rsp),%rdx
mov    -0x30(%rsp),%rbp
je     1500 <transpose_so3+0x400>
lea    (%r15,%rbx,1),%eax
mov    %ebx,%ecx
imul   -0x40(%rsp),%ecx
add    %eax,%ecx
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm0
imul   %edx,%eax
cltq
add    %rbx,%rax
movsd  %xmm0,(%rsi,%rax,8)
jmp    1500 <transpose_so3+0x400>
mov    $0x1,%r10d
sub    %edx,%r10d
js     16e4 <transpose_so3+0x5e4>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %ecx,%r13d
mov    %r13,-0x20(%rsp)
and    $0xfffffffe,%r13d
mov    -0x38(%rsp),%rax
lea    (%rax,%rax,1),%r9d
sub    %ecx,%r9d
add    $0x1,%r9d
lea    (%rax,%rax,1),%ebp
add    $0x2,%ebp
mov    %eax,%r10d
sub    %ecx,%r10d
lea    0x1(%r10),%r14d
imul   %ecx,%r14d
lea    (%rcx,%rcx,1),%ebx
mov    %ecx,%r15d
imul   %r10d,%r15d
mov    -0x18(%rsp),%r12d
jmp    1637 <transpose_so3+0x537>
cs nopw 0x0(%rax,%rax,1)
add    $0xffffffff,%r9d
sub    %ecx,%r14d
add    $0xffffffff,%r10d
sub    %ecx,%r15d
inc    %r12d
je     1393 <transpose_so3+0x293>
test   %ecx,%ecx
jle    1620 <transpose_so3+0x520>
cmp    $0x1,%ecx
mov    %r9d,-0x28(%rsp)
mov    %r10,-0x30(%rsp)
jne    1650 <transpose_so3+0x550>
xor    %eax,%eax
jmp    169f <transpose_so3+0x59f>
xchg   %ax,%ax
mov    %r15d,%edx
mov    %r10d,%ecx
mov    %r14d,%r8d
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movslq %r8d,%r8
lea    (%r8,%rax,1),%r10
movslq %edx,%rdx
lea    (%rax,%rdx,1),%r11
movsd  (%rdi,%r11,8),%xmm0
movslq %ecx,%rcx
movsd  %xmm0,(%rsi,%rcx,8)
movsd  0x8(%rdi,%r10,8),%xmm0
movslq %r9d,%r9
movsd  %xmm0,(%rsi,%r9,8)
add    $0x2,%rax
add    %ebp,%r9d
add    %ebx,%r8d
add    %ebp,%ecx
add    %ebx,%edx
cmp    %rax,%r13
jne    1660 <transpose_so3+0x560>
testb  $0x1,-0x20(%rsp)
mov    -0x40(%rsp),%rcx
mov    -0x28(%rsp),%r9d
mov    -0x30(%rsp),%r10
je     1620 <transpose_so3+0x520>
mov    %eax,%r8d
sub    %r12d,%r8d
mov    %r8d,%edx
imul   %ecx,%edx
movslq %edx,%rdx
add    %rax,%rdx
movsd  (%rdi,%rdx,8),%xmm0
imul   -0x38(%rsp),%eax
add    %r8d,%eax
cltq
movsd  %xmm0,(%rsi,%rax,8)
jmp    1620 <transpose_so3+0x520>
mov    -0x38(%rsp),%rdx
lea    (%rdx,%rdx,1),%r11d
lea    (%rdx,%rdx,1),%eax
add    $0x2,%eax
mov    -0x40(%rsp),%rcx
mov    %ecx,%r9d
imul   %edx,%r9d
lea    (%rcx,%rcx,1),%ebp
lea    -0x1(%rdx),%r14d
mov    %ecx,%r15d
imul   %r14d,%r15d
mov    $0x1,%r8d
xor    %r12d,%r12d
jmp    174a <transpose_so3+0x64a>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r10d
add    $0x1,%r8
add    $0x1,%r12
add    $0xffffffff,%r11d
mov    -0x40(%rsp),%rcx
sub    %ecx,%r9d
add    $0xffffffff,%r14d
sub    %ecx,%r15d
cmp    -0x28(%rsp),%r8
je     15ca <transpose_so3+0x4ca>
mov    %r10d,%ecx
add    -0x38(%rsp),%ecx
jle    1720 <transpose_so3+0x620>
test   %r12,%r12
mov    %r11d,-0x30(%rsp)
mov    %r9d,-0x20(%rsp)
mov    %r14d,-0x18(%rsp)
mov    %r15d,-0x10(%rsp)
je     17d0 <transpose_so3+0x6d0>
mov    %r8,%r13
and    $0xfffffffffffffffe,%r13
mov    %r15d,%edx
mov    %r14d,%ecx
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movslq %r9d,%r9
lea    (%r9,%rbx,1),%r14
movslq %edx,%rdx
lea    (%rbx,%rdx,1),%r15
movsd  (%rdi,%r15,8),%xmm0
movslq %ecx,%rcx
movsd  %xmm0,(%rsi,%rcx,8)
movsd  0x8(%rdi,%r14,8),%xmm0
movslq %r11d,%r11
movsd  %xmm0,(%rsi,%r11,8)
add    $0x2,%rbx
add    %eax,%r11d
add    %ebp,%r9d
add    %eax,%ecx
add    %ebp,%edx
cmp    %rbx,%r13
jne    1780 <transpose_so3+0x680>
jmp    17d2 <transpose_so3+0x6d2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xor    %ebx,%ebx
test   $0x1,%r8b
mov    -0x30(%rsp),%r11d
mov    -0x20(%rsp),%r9d
mov    -0x18(%rsp),%r14d
mov    -0x10(%rsp),%r15d
je     1720 <transpose_so3+0x620>
mov    %ebx,%ecx
sub    %r10d,%ecx
mov    %ecx,%edx
imul   -0x40(%rsp),%edx
movslq %edx,%rdx
add    %rbx,%rdx
movsd  (%rdi,%rdx,8),%xmm0
imul   -0x38(%rsp),%ebx
add    %ecx,%ebx
movslq %ebx,%rcx
movsd  %xmm0,(%rsi,%rcx,8)
jmp    1720 <transpose_so3+0x620>

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

LLVM IR (-O2):
define void @transpose_so3(double* nocapture %0, double* nocapture writeonly %1, i32 %2, i32 %3) {
%5 = icmp slt i32 %3, %2
br i1 %5, label %6, label %8
6: ; preds = %4
%7 = icmp sgt i32 %3, 0
br i1 %7, label %204, label %264
8: ; preds = %4
%9 = sub nsw i32 %3, %2
%10 = add i32 %9, 1
%11 = icmp slt i32 %10, %3
br i1 %11, label %41, label %32
12: ; preds = %53, %48
%13 = phi i64 [ 0, %48 ], [ %81, %53 ]
%14 = icmp eq i64 %49, 0
br i1 %14, label %28, label %15
15: ; preds = %12
%16 = trunc i64 %13 to i32
%17 = add nsw i32 %43, %16
%18 = trunc i64 %13 to i32
%19 = mul i32 %18, %3
%20 = add i32 %17, %19
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %0, i64 %21
%23 = load double, double* %22
%24 = mul nsw i32 %17, %2
%25 = sext i32 %24 to i64
%26 = add nsw i64 %13, %25
%27 = getelementptr inbounds double, double* %1, i64 %26
store double %23, double* %27
br label %28
28: ; preds = %15, %12, %41
%29 = add nsw i32 %43, 1
%30 = icmp eq i32 %29, %3
%31 = add i32 %42, 1
br i1 %30, label %32, label %41
32: ; preds = %28, %8
%33 = icmp slt i32 %9, 0
br i1 %33, label %136, label %34
34: ; preds = %32
%35 = icmp sgt i32 %2, 0
%36 = zext i32 %2 to i64
%37 = and i64 %36, 1
%38 = icmp eq i32 %2, 1
%39 = and i64 %36, 4294967294
%40 = icmp eq i64 %37, 0
br label %84
41: ; preds = %8, %28
%42 = phi i32 [ %31, %28 ], [ 0, %8 ]
%43 = phi i32 [ %29, %28 ], [ %10, %8 ]
%44 = xor i32 %42, -1
%45 = add i32 %44, %2
%46 = zext i32 %45 to i64
%47 = icmp slt i32 %43, %3
br i1 %47, label %48, label %28
48: ; preds = %41
%49 = and i64 %46, 1
%50 = icmp eq i32 %45, 1
br i1 %50, label %12, label %51
51: ; preds = %48
%52 = and i64 %46, 4294967294
br label %53
53: ; preds = %53, %51
%54 = phi i64 [ 0, %51 ], [ %81, %53 ]
%55 = phi i64 [ 0, %51 ], [ %82, %53 ]
%56 = trunc i64 %54 to i32
%57 = add nsw i32 %43, %56
%58 = trunc i64 %54 to i32
%59 = mul i32 %58, %3
%60 = add i32 %57, %59
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = load double, double* %62
%64 = mul nsw i32 %57, %2
%65 = sext i32 %64 to i64
%66 = add nsw i64 %54, %65
%67 = getelementptr inbounds double, double* %1, i64 %66
store double %63, double* %67
%68 = or i64 %54, 1
%69 = trunc i64 %68 to i32
%70 = add nsw i32 %43, %69
%71 = trunc i64 %68 to i32
%72 = mul i32 %71, %3
%73 = add i32 %70, %72
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %0, i64 %74
%76 = load double, double* %75
%77 = mul nsw i32 %70, %2
%78 = sext i32 %77 to i64
%79 = add nsw i64 %68, %78
%80 = getelementptr inbounds double, double* %1, i64 %79
store double %76, double* %80
%81 = add nuw nsw i64 %54, 2
%82 = add i64 %55, 2
%83 = icmp eq i64 %82, %52
br i1 %83, label %12, label %53
84: ; preds = %34, %133
%85 = phi i32 [ 0, %34 ], [ %134, %133 ]
br i1 %35, label %86, label %133
86: ; preds = %84
br i1 %38, label %118, label %87
87: ; preds = %86, %87
%88 = phi i64 [ %115, %87 ], [ 0, %86 ]
%89 = phi i64 [ %116, %87 ], [ 0, %86 ]
%90 = trunc i64 %88 to i32
%91 = add nuw nsw i32 %85, %90
%92 = trunc i64 %88 to i32
%93 = mul i32 %92, %3
%94 = add i32 %91, %93
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %0, i64 %95
%97 = load double, double* %96
%98 = mul nsw i32 %91, %2
%99 = sext i32 %98 to i64
%100 = add nsw i64 %88, %99
%101 = getelementptr inbounds double, double* %1, i64 %100
store double %97, double* %101
%102 = or i64 %88, 1
%103 = trunc i64 %102 to i32
%104 = add nuw nsw i32 %85, %103
%105 = trunc i64 %102 to i32
%106 = mul i32 %105, %3
%107 = add i32 %104, %106
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds double, double* %0, i64 %108
%110 = load double, double* %109
%111 = mul nsw i32 %104, %2
%112 = sext i32 %111 to i64
%113 = add nsw i64 %102, %112
%114 = getelementptr inbounds double, double* %1, i64 %113
store double %110, double* %114
%115 = add nuw nsw i64 %88, 2
%116 = add i64 %89, 2
%117 = icmp eq i64 %116, %39
br i1 %117, label %118, label %87
118: ; preds = %87, %86
%119 = phi i64 [ 0, %86 ], [ %115, %87 ]
br i1 %40, label %133, label %120
120: ; preds = %118
%121 = trunc i64 %119 to i32
%122 = add nuw nsw i32 %85, %121
%123 = trunc i64 %119 to i32
%124 = mul i32 %123, %3
%125 = add i32 %122, %124
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds double, double* %0, i64 %126
%128 = load double, double* %127
%129 = mul nsw i32 %122, %2
%130 = sext i32 %129 to i64
%131 = add nsw i64 %119, %130
%132 = getelementptr inbounds double, double* %1, i64 %131
store double %128, double* %132
br label %133
133: ; preds = %120, %118, %84
%134 = add nuw i32 %85, 1
%135 = icmp eq i32 %85, %9
br i1 %135, label %136, label %84
136: ; preds = %133, %32
%137 = sub i32 1, %2
%138 = icmp slt i32 %137, 0
br i1 %138, label %139, label %390
139: ; preds = %136
%140 = zext i32 %2 to i64
br label %141
141: ; preds = %139, %199
%142 = phi i64 [ 0, %139 ], [ %203, %199 ]
%143 = phi i64 [ 1, %139 ], [ %201, %199 ]
%144 = phi i32 [ %137, %139 ], [ %200, %199 ]
%145 = add nsw i32 %144, %2
%146 = icmp sgt i32 %145, 0
br i1 %146, label %147, label %199
147: ; preds = %141
%148 = and i64 %143, 1
%149 = icmp eq i64 %142, 0
br i1 %149, label %183, label %150
150: ; preds = %147
%151 = and i64 %143, 9223372036854775806
br label %152
152: ; preds = %152, %150
%153 = phi i64 [ 0, %150 ], [ %180, %152 ]
%154 = phi i64 [ 0, %150 ], [ %181, %152 ]
%155 = trunc i64 %153 to i32
%156 = sub i32 %155, %144
%157 = mul nsw i32 %156, %3
%158 = sext i32 %157 to i64
%159 = add nsw i64 %153, %158
%160 = getelementptr inbounds double, double* %0, i64 %159
%161 = load double, double* %160
%162 = trunc i64 %153 to i32
%163 = mul i32 %162, %2
%164 = add i32 %156, %163
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds double, double* %1, i64 %165
store double %161, double* %166
%167 = or i64 %153, 1
%168 = trunc i64 %167 to i32
%169 = sub i32 %168, %144
%170 = mul nsw i32 %169, %3
%171 = sext i32 %170 to i64
%172 = add nsw i64 %167, %171
%173 = getelementptr inbounds double, double* %0, i64 %172
%174 = load double, double* %173
%175 = trunc i64 %167 to i32
%176 = mul i32 %175, %2
%177 = add i32 %169, %176
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds double, double* %1, i64 %178
store double %174, double* %179
%180 = add nuw nsw i64 %153, 2
%181 = add i64 %154, 2
%182 = icmp eq i64 %181, %151
br i1 %182, label %183, label %152
183: ; preds = %152, %147
%184 = phi i64 [ 0, %147 ], [ %180, %152 ]
%185 = icmp eq i64 %148, 0
br i1 %185, label %199, label %186
186: ; preds = %183
%187 = trunc i64 %184 to i32
%188 = sub i32 %187, %144
%189 = mul nsw i32 %188, %3
%190 = sext i32 %189 to i64
%191 = add nsw i64 %184, %190
%192 = getelementptr inbounds double, double* %0, i64 %191
%193 = load double, double* %192
%194 = trunc i64 %184 to i32
%195 = mul i32 %194, %2
%196 = add i32 %188, %195
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds double, double* %1, i64 %197
store double %193, double* %198
br label %199
199: ; preds = %186, %183, %141
%200 = add nsw i32 %144, 1
%201 = add nuw nsw i64 %143, 1
%202 = icmp eq i64 %201, %140
%203 = add i64 %142, 1
br i1 %202, label %390, label %141
204: ; preds = %6, %261
%205 = phi i32 [ %262, %261 ], [ 0, %6 ]
%206 = sub i32 %3, %205
%207 = zext i32 %206 to i64
%208 = icmp slt i32 %205, %3
br i1 %208, label %209, label %261
209: ; preds = %204
%210 = and i64 %207, 1
%211 = icmp eq i32 %206, 1
br i1 %211, label %245, label %212
212: ; preds = %209
%213 = and i64 %207, 4294967294
br label %214
214: ; preds = %214, %212
%215 = phi i64 [ 0, %212 ], [ %242, %214 ]
%216 = phi i64 [ 0, %212 ], [ %243, %214 ]
%217 = trunc i64 %215 to i32
%218 = add nuw nsw i32 %205, %217
%219 = trunc i64 %215 to i32
%220 = mul i32 %219, %3
%221 = add i32 %218, %220
%222 = sext i32 %221 to i64
%223 = getelementptr inbounds double, double* %0, i64 %222
%224 = load double, double* %223
%225 = mul nsw i32 %218, %2
%226 = sext i32 %225 to i64
%227 = add nsw i64 %215, %226
%228 = getelementptr inbounds double, double* %1, i64 %227
store double %224, double* %228
%229 = or i64 %215, 1
%230 = trunc i64 %229 to i32
%231 = add nuw nsw i32 %205, %230
%232 = trunc i64 %229 to i32
%233 = mul i32 %232, %3
%234 = add i32 %231, %233
%235 = sext i32 %234 to i64
%236 = getelementptr inbounds double, double* %0, i64 %235
%237 = load double, double* %236
%238 = mul nsw i32 %231, %2
%239 = sext i32 %238 to i64
%240 = add nsw i64 %229, %239
%241 = getelementptr inbounds double, double* %1, i64 %240
store double %237, double* %241
%242 = add nuw nsw i64 %215, 2
%243 = add i64 %216, 2
%244 = icmp eq i64 %243, %213
br i1 %244, label %245, label %214
245: ; preds = %214, %209
%246 = phi i64 [ 0, %209 ], [ %242, %214 ]
%247 = icmp eq i64 %210, 0
br i1 %247, label %261, label %248
248: ; preds = %245
%249 = trunc i64 %246 to i32
%250 = add nuw nsw i32 %205, %249
%251 = trunc i64 %246 to i32
%252 = mul i32 %251, %3
%253 = add i32 %250, %252
%254 = sext i32 %253 to i64
%255 = getelementptr inbounds double, double* %0, i64 %254
%256 = load double, double* %255
%257 = mul nsw i32 %250, %2
%258 = sext i32 %257 to i64
%259 = add nsw i64 %246, %258
%260 = getelementptr inbounds double, double* %1, i64 %259
store double %256, double* %260
br label %261
261: ; preds = %248, %245, %204
%262 = add nuw nsw i32 %205, 1
%263 = icmp eq i32 %262, %3
br i1 %263, label %264, label %204
264: ; preds = %261, %6
%265 = sub nsw i32 %3, %2
%266 = icmp slt i32 %265, 0
br i1 %266, label %267, label %325
267: ; preds = %264
%268 = zext i32 %3 to i64
%269 = and i64 %268, 1
%270 = icmp eq i32 %3, 1
%271 = and i64 %268, 4294967294
%272 = icmp eq i64 %269, 0
br label %273
273: ; preds = %267, %322
%274 = phi i32 [ %265, %267 ], [ %323, %322 ]
br i1 %7, label %275, label %322
275: ; preds = %273
br i1 %270, label %307, label %276
276: ; preds = %275, %276
%277 = phi i64 [ %304, %276 ], [ 0, %275 ]
%278 = phi i64 [ %305, %276 ], [ 0, %275 ]
%279 = trunc i64 %277 to i32
%280 = sub i32 %279, %274
%281 = mul nsw i32 %280, %3
%282 = sext i32 %281 to i64
%283 = add nsw i64 %277, %282
%284 = getelementptr inbounds double, double* %0, i64 %283
%285 = load double, double* %284
%286 = trunc i64 %277 to i32
%287 = mul i32 %286, %2
%288 = add i32 %280, %287
%289 = sext i32 %288 to i64
%290 = getelementptr inbounds double, double* %1, i64 %289
store double %285, double* %290
%291 = or i64 %277, 1
%292 = trunc i64 %291 to i32
%293 = sub i32 %292, %274
%294 = mul nsw i32 %293, %3
%295 = sext i32 %294 to i64
%296 = add nsw i64 %291, %295
%297 = getelementptr inbounds double, double* %0, i64 %296
%298 = load double, double* %297
%299 = trunc i64 %291 to i32
%300 = mul i32 %299, %2
%301 = add i32 %293, %300
%302 = sext i32 %301 to i64
%303 = getelementptr inbounds double, double* %1, i64 %302
store double %298, double* %303
%304 = add nuw nsw i64 %277, 2
%305 = add i64 %278, 2
%306 = icmp eq i64 %305, %271
br i1 %306, label %307, label %276
307: ; preds = %276, %275
%308 = phi i64 [ 0, %275 ], [ %304, %276 ]
br i1 %272, label %322, label %309
309: ; preds = %307
%310 = trunc i64 %308 to i32
%311 = sub i32 %310, %274
%312 = mul nsw i32 %311, %3
%313 = sext i32 %312 to i64
%314 = add nsw i64 %308, %313
%315 = getelementptr inbounds double, double* %0, i64 %314
%316 = load double, double* %315
%317 = trunc i64 %308 to i32
%318 = mul i32 %317, %2
%319 = add i32 %311, %318
%320 = sext i32 %319 to i64
%321 = getelementptr inbounds double, double* %1, i64 %320
store double %316, double* %321
br label %322
322: ; preds = %309, %307, %273
%323 = add i32 %274, 1
%324 = icmp eq i32 %323, 0
br i1 %324, label %325, label %273
325: ; preds = %322, %264
%326 = sub i32 1, %2
%327 = icmp slt i32 %326, %265
br i1 %327, label %328, label %390
328: ; preds = %325, %386
%329 = phi i32 [ %388, %386 ], [ 1, %325 ]
%330 = phi i32 [ %387, %386 ], [ %326, %325 ]
%331 = zext i32 %329 to i64
%332 = add nsw i32 %330, %2
%333 = icmp sgt i32 %332, 0
br i1 %333, label %334, label %386
334: ; preds = %328
%335 = and i64 %331, 1
%336 = icmp eq i32 %329, 1
br i1 %336, label %370, label %337
337: ; preds = %334
%338 = and i64 %331, 4294967294
br label %339
339: ; preds = %339, %337
%340 = phi i64 [ 0, %337 ], [ %367, %339 ]
%341 = phi i64 [ 0, %337 ], [ %368, %339 ]
%342 = trunc i64 %340 to i32
%343 = sub i32 %342, %330
%344 = mul nsw i32 %343, %3
%345 = sext i32 %344 to i64
%346 = add nsw i64 %340, %345
%347 = getelementptr inbounds double, double* %0, i64 %346
%348 = load double, double* %347
%349 = trunc i64 %340 to i32
%350 = mul i32 %349, %2
%351 = add i32 %343, %350
%352 = sext i32 %351 to i64
%353 = getelementptr inbounds double, double* %1, i64 %352
store double %348, double* %353
%354 = or i64 %340, 1
%355 = trunc i64 %354 to i32
%356 = sub i32 %355, %330
%357 = mul nsw i32 %356, %3
%358 = sext i32 %357 to i64
%359 = add nsw i64 %354, %358
%360 = getelementptr inbounds double, double* %0, i64 %359
%361 = load double, double* %360
%362 = trunc i64 %354 to i32
%363 = mul i32 %362, %2
%364 = add i32 %356, %363
%365 = sext i32 %364 to i64
%366 = getelementptr inbounds double, double* %1, i64 %365
store double %361, double* %366
%367 = add nuw nsw i64 %340, 2
%368 = add i64 %341, 2
%369 = icmp eq i64 %368, %338
br i1 %369, label %370, label %339
370: ; preds = %339, %334
%371 = phi i64 [ 0, %334 ], [ %367, %339 ]
%372 = icmp eq i64 %335, 0
br i1 %372, label %386, label %373
373: ; preds = %370
%374 = trunc i64 %371 to i32
%375 = sub i32 %374, %330
%376 = mul nsw i32 %375, %3
%377 = sext i32 %376 to i64
%378 = add nsw i64 %371, %377
%379 = getelementptr inbounds double, double* %0, i64 %378
%380 = load double, double* %379
%381 = trunc i64 %371 to i32
%382 = mul i32 %381, %2
%383 = add i32 %375, %382
%384 = sext i32 %383 to i64
%385 = getelementptr inbounds double, double* %1, i64 %384
store double %380, double* %385
br label %386
386: ; preds = %373, %370, %328
%387 = add nsw i32 %330, 1
%388 = add i32 %329, 1
%389 = icmp eq i32 %388, %3
br i1 %389, label %390, label %328
390: ; preds = %199, %386, %136, %325
ret void
}