void switch_row_and_colMAT(double *A, int n, int i, int j) {
  double tmp;
  int k;
  for (k = 0; k < j; k++) {
    tmp = A[(((j) * ((j) + 1)) / 2 + (k))];
    A[(((j) * ((j) + 1)) / 2 + (k))] = A[(((i) * ((i) + 1)) / 2 + (k))];
    A[(((i) * ((i) + 1)) / 2 + (k))] = tmp;
  }
  for (k = j + 1; k < i; k++) {
    tmp = A[(((k) * ((k) + 1)) / 2 + (j))];
    A[(((k) * ((k) + 1)) / 2 + (j))] = A[(((i) * ((i) + 1)) / 2 + (k))];
    A[(((i) * ((i) + 1)) / 2 + (k))] = tmp;
  }
  for (k = i + 1; k < n; k++) {
    tmp = A[(((k) * ((k) + 1)) / 2 + (j))];
    A[(((k) * ((k) + 1)) / 2 + (j))] = A[(((k) * ((k) + 1)) / 2 + (i))];
    A[(((k) * ((k) + 1)) / 2 + (i))] = tmp;
  }
  tmp = A[(((j) * ((j) + 1)) / 2 + (j))];
  A[(((j) * ((j) + 1)) / 2 + (j))] = A[(((i) * ((i) + 1)) / 2 + (i))];
  A[(((i) * ((i) + 1)) / 2 + (i))] = tmp;
}


x86 assembly (clang, -O0):
<switch_row_and_colMAT>:
push   %rbp
mov    %rsp,%rbp
sub    $0x28,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11f4 <switch_row_and_colMAT+0xf4>
mov    -0x8(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x48(%rbp),%rax
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x38(%rbp),%rax
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x30(%rbp),%rax
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111c <switch_row_and_colMAT+0x1c>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12d5 <switch_row_and_colMAT+0x1d5>
mov    -0x8(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x68(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x60(%rbp),%rax
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x58(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x50(%rbp),%rax
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11fd <switch_row_and_colMAT+0xfd>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13bc <switch_row_and_colMAT+0x2bc>
mov    -0x8(%rbp),%rax
mov    %rax,-0x88(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x88(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x80(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x80(%rbp),%rax
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x78(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x78(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x70(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x70(%rbp),%rax
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12de <switch_row_and_colMAT+0x1de>
mov    -0x8(%rbp),%rax
mov    %rax,-0xa8(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xa8(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0xa0(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xa0(%rbp),%rax
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x98(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x98(%rbp),%rax
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x90(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x90(%rbp),%rax
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x28,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<switch_row_and_colMAT>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    0x1(%rcx),%r8d
test   %ecx,%ecx
jle    12e6 <switch_row_and_colMAT+0x1e6>
mov    %r8d,%ebx
imul   %ecx,%ebx
mov    %ebx,%eax
shr    $0x1f,%eax
add    %ebx,%eax
sar    %eax
lea    0x1(%rdx),%r9d
imul   %edx,%r9d
mov    %r9d,%ebx
shr    $0x1f,%ebx
add    %r9d,%ebx
sar    %ebx
movslq %eax,%r9
movslq %ebx,%r10
mov    %ecx,%r11d
cmp    $0x4,%ecx
jb     1170 <switch_row_and_colMAT+0x70>
lea    (%rdi,%r9,8),%rax
lea    (%r10,%r11,1),%rbx
lea    (%rdi,%rbx,8),%rbx
cmp    %rbx,%rax
jae    11fe <switch_row_and_colMAT+0xfe>
lea    (%r9,%r11,1),%rax
lea    (%rdi,%rax,8),%rax
lea    (%rdi,%r10,8),%rbx
cmp    %rax,%rbx
jae    11fe <switch_row_and_colMAT+0xfe>
xor    %eax,%eax
mov    %rax,%rbx
not    %rbx
test   $0x1,%r11b
je     11a2 <switch_row_and_colMAT+0xa2>
lea    (%rax,%r9,1),%r14
movsd  (%rdi,%r14,8),%xmm0
lea    (%rax,%r10,1),%r15
movsd  (%rdi,%r15,8),%xmm1
movsd  %xmm1,(%rdi,%r14,8)
movsd  %xmm0,(%rdi,%r15,8)
or     $0x1,%rax
add    %r11,%rbx
je     12e6 <switch_row_and_colMAT+0x1e6>
lea    (%rdi,%r10,8),%r10
add    $0x8,%r10
lea    (%rdi,%r9,8),%rbx
add    $0x8,%rbx
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rbx,%rax,8),%xmm0
movsd  -0x8(%r10,%rax,8),%xmm1
movsd  %xmm1,-0x8(%rbx,%rax,8)
movsd  %xmm0,-0x8(%r10,%rax,8)
movsd  (%rbx,%rax,8),%xmm0
movsd  (%r10,%rax,8),%xmm1
movsd  %xmm1,(%rbx,%rax,8)
movsd  %xmm0,(%r10,%rax,8)
add    $0x2,%rax
cmp    %rax,%r11
jne    11c0 <switch_row_and_colMAT+0xc0>
jmp    12e6 <switch_row_and_colMAT+0x1e6>
mov    %r11d,%eax
and    $0xfffffffc,%eax
lea    -0x4(%rax),%rbx
mov    %rbx,%r14
shr    $0x2,%r14
add    $0x1,%r14
test   %rbx,%rbx
je     13c9 <switch_row_and_colMAT+0x2c9>
mov    %r14,%r15
and    $0xfffffffffffffffe,%r15
lea    (%rdi,%r9,8),%r12
add    $0x30,%r12
lea    (%rdi,%r10,8),%r13
add    $0x30,%r13
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movups -0x30(%r12,%rbx,8),%xmm0
movups -0x20(%r12,%rbx,8),%xmm1
movups -0x30(%r13,%rbx,8),%xmm2
movups -0x20(%r13,%rbx,8),%xmm3
movups %xmm2,-0x30(%r12,%rbx,8)
movups %xmm3,-0x20(%r12,%rbx,8)
movups %xmm0,-0x30(%r13,%rbx,8)
movups %xmm1,-0x20(%r13,%rbx,8)
movups -0x10(%r12,%rbx,8),%xmm0
movups (%r12,%rbx,8),%xmm1
movups -0x10(%r13,%rbx,8),%xmm2
movups 0x0(%r13,%rbx,8),%xmm3
movups %xmm2,-0x10(%r12,%rbx,8)
movups %xmm3,(%r12,%rbx,8)
movups %xmm0,-0x10(%r13,%rbx,8)
movups %xmm1,0x0(%r13,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%r15
jne    1240 <switch_row_and_colMAT+0x140>
test   $0x1,%r14b
je     12dd <switch_row_and_colMAT+0x1dd>
lea    (%rbx,%r9,1),%r14
movups (%rdi,%r14,8),%xmm0
movups 0x10(%rdi,%r14,8),%xmm1
add    %r10,%rbx
movups (%rdi,%rbx,8),%xmm2
movups 0x10(%rdi,%rbx,8),%xmm3
movups %xmm2,(%rdi,%r14,8)
movups %xmm3,0x10(%rdi,%r14,8)
movups %xmm0,(%rdi,%rbx,8)
movups %xmm1,0x10(%rdi,%rbx,8)
cmp    %r11,%rax
jne    1172 <switch_row_and_colMAT+0x72>
lea    0x1(%rdx),%r9d
cmp    %edx,%r8d
jge    1347 <switch_row_and_colMAT+0x247>
mov    %r9d,%ebx
imul   %edx,%ebx
mov    %ebx,%eax
shr    $0x1f,%eax
add    %ebx,%eax
sar    %eax
movslq %r8d,%rbx
cltq
lea    (%rdi,%rax,8),%r10
mov    %ecx,%r14d
not    %r14d
add    %edx,%r14d
lea    0x1(%rbx),%r11
mov    %r11d,%eax
imul   %ebx,%eax
mov    %eax,%ebp
shr    $0x1f,%ebp
add    %eax,%ebp
sar    %ebp
add    %ecx,%ebp
movslq %ebp,%rax
movsd  (%rdi,%rax,8),%xmm0
movsd  (%r10,%rbx,8),%xmm1
movsd  %xmm1,(%rdi,%rax,8)
movsd  %xmm0,(%r10,%rbx,8)
mov    %r11,%rbx
add    $0xffffffff,%r14d
jne    1310 <switch_row_and_colMAT+0x210>
cmp    %esi,%r9d
jge    1383 <switch_row_and_colMAT+0x283>
mov    %r9d,%eax
nop
lea    0x1(%rax),%ebp
imul   %ebp,%eax
mov    %eax,%ebx
shr    $0x1f,%ebx
add    %eax,%ebx
sar    %ebx
lea    (%rbx,%rcx,1),%eax
cltq
movsd  (%rdi,%rax,8),%xmm0
add    %edx,%ebx
movslq %ebx,%rbx
movsd  (%rdi,%rbx,8),%xmm1
movsd  %xmm1,(%rdi,%rax,8)
movsd  %xmm0,(%rdi,%rbx,8)
mov    %ebp,%eax
cmp    %ebp,%esi
jne    1350 <switch_row_and_colMAT+0x250>
imul   %ecx,%r8d
mov    %r8d,%eax
shr    $0x1f,%eax
add    %r8d,%eax
sar    %eax
add    %ecx,%eax
cltq
movsd  (%rdi,%rax,8),%xmm0
imul   %edx,%r9d
mov    %r9d,%ecx
shr    $0x1f,%ecx
add    %r9d,%ecx
sar    %ecx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rdi,%rcx,8),%xmm1
movsd  %xmm1,(%rdi,%rax,8)
movsd  %xmm0,(%rdi,%rcx,8)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ebx,%ebx
test   $0x1,%r14b
jne    12ae <switch_row_and_colMAT+0x1ae>
jmp    12dd <switch_row_and_colMAT+0x1dd>

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

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