int pcapack_symmetrize(const int triang, const int m, const int n, double *x) {
  int i, j;
  const int k = m < n ? m : n;
  if (m == 0 || n == 0)
    return 0;
  if (triang != 1 && triang != 2)
    return -1;
  if (__builtin_expect((triang == 1), 1)) {
    for (j = 0; j < k; j++) {
      for (i = 0; i < j; i++)
        x[j + m * i] = x[i + m * j];
    }
  } else if (__builtin_expect((triang == 2), 0)) {
    for (j = 0; j < k; j++) {
      for (i = j + 1; i < k; i++)
        x[j + m * i] = x[i + m * j];
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<pcapack_symmetrize>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1128 <pcapack_symmetrize+0x28>
mov    -0xc(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    112e <pcapack_symmetrize+0x2e>
mov    -0x10(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0xc(%rbp)
je     1148 <pcapack_symmetrize+0x48>
cmpl   $0x0,-0x10(%rbp)
jne    1154 <pcapack_symmetrize+0x54>
movl   $0x0,-0x4(%rbp)
jmp    12aa <pcapack_symmetrize+0x1aa>
cmpl   $0x1,-0x8(%rbp)
je     1174 <pcapack_symmetrize+0x74>
cmpl   $0x2,-0x8(%rbp)
je     1174 <pcapack_symmetrize+0x74>
movl   $0xffffffff,-0x4(%rbp)
jmp    12aa <pcapack_symmetrize+0x1aa>
cmpl   $0x1,-0x8(%rbp)
sete   %al
and    $0x1,%al
movzbl %al,%eax
cltq
cmp    $0x0,%rax
je     1208 <pcapack_symmetrize+0x108>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1203 <pcapack_symmetrize+0x103>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11f0 <pcapack_symmetrize+0xf0>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11a6 <pcapack_symmetrize+0xa6>
jmp    11f5 <pcapack_symmetrize+0xf5>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1193 <pcapack_symmetrize+0x93>
jmp    12a3 <pcapack_symmetrize+0x1a3>
cmpl   $0x2,-0x8(%rbp)
sete   %al
and    $0x1,%al
movzbl %al,%eax
cltq
cmp    $0x0,%rax
je     129e <pcapack_symmetrize+0x19e>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1299 <pcapack_symmetrize+0x199>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1286 <pcapack_symmetrize+0x186>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    123c <pcapack_symmetrize+0x13c>
jmp    128b <pcapack_symmetrize+0x18b>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1227 <pcapack_symmetrize+0x127>
jmp    129e <pcapack_symmetrize+0x19e>
jmp    12a3 <pcapack_symmetrize+0x1a3>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<pcapack_symmetrize>:
push   %r14
push   %rbx
cmp    %edx,%esi
mov    %edx,%r9d
cmovl  %esi,%r9d
test   %esi,%esi
je     11e8 <pcapack_symmetrize+0xe8>
test   %edx,%edx
je     11e8 <pcapack_symmetrize+0xe8>
lea    -0x3(%rdi),%edx
mov    $0xffffffff,%eax
cmp    $0xfffffffe,%edx
jb     11ea <pcapack_symmetrize+0xea>
cmp    $0x1,%edi
jne    11e3 <pcapack_symmetrize+0xe3>
mov    $0x0,%eax
test   %r9d,%r9d
jle    11ea <pcapack_symmetrize+0xea>
movslq %esi,%r8
mov    %r9d,%r9d
lea    0x8(%rcx),%rsi
lea    0x0(,%r8,8),%rdi
mov    %r8,%r11
shl    $0x4,%r11
xor    %r14d,%r14d
mov    %rcx,%r10
jmp    1180 <pcapack_symmetrize+0x80>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r14
add    %rdi,%rsi
add    $0x8,%r10
cmp    %r9,%r14
je     11e8 <pcapack_symmetrize+0xe8>
test   %r14,%r14
je     1170 <pcapack_symmetrize+0x70>
cmp    $0x1,%r14
jne    1190 <pcapack_symmetrize+0x90>
xor    %eax,%eax
jmp    11c0 <pcapack_symmetrize+0xc0>
nop
mov    %r14,%rbx
and    $0xfffffffffffffffe,%rbx
mov    %r10,%rdx
xor    %eax,%eax
nopl   0x0(%rax)
movsd  -0x8(%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx)
movsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rdi,1)
add    $0x2,%rax
add    %r11,%rdx
cmp    %rax,%rbx
jne    11a0 <pcapack_symmetrize+0xa0>
test   $0x1,%r14b
je     1170 <pcapack_symmetrize+0x70>
mov    %r14,%rdx
imul   %r8,%rdx
add    %rax,%rdx
movsd  (%rcx,%rdx,8),%xmm0
imul   %r8,%rax
add    %r14,%rax
movsd  %xmm0,(%rcx,%rax,8)
jmp    1170 <pcapack_symmetrize+0x70>
cmp    $0x2,%edi
je     11ee <pcapack_symmetrize+0xee>
xor    %eax,%eax
pop    %rbx
pop    %r14
ret
mov    $0x0,%eax
test   %r9d,%r9d
jle    11ea <pcapack_symmetrize+0xea>
movslq %esi,%rdx
mov    %r9d,%r8d
lea    (%rcx,%rdx,8),%rbx
lea    0x8(,%rdx,8),%r9
shl    $0x3,%rdx
xor    %r10d,%r10d
mov    %r8,%rsi
jmp    1227 <pcapack_symmetrize+0x127>
add    %r9,%rbx
add    %r9,%rcx
add    $0xffffffffffffffff,%rsi
xor    %eax,%eax
cmp    %r8,%r10
je     11ea <pcapack_symmetrize+0xea>
add    $0x1,%r10
cmp    %r8,%r10
jae    1216 <pcapack_symmetrize+0x116>
mov    $0x1,%eax
mov    %rbx,%rdi
movsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,(%rdi)
add    %rdx,%rdi
add    $0x1,%rax
cmp    %rax,%rsi
jne    1238 <pcapack_symmetrize+0x138>
jmp    1216 <pcapack_symmetrize+0x116>

LLVM IR (-O0):
define i32 @pcapack_symmetrize(i32 %0, i32 %1, i32 %2, double* %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double* %3, double** %9
%13 = load i32, i32* %7
%14 = load i32, i32* %8
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %18
16: ; preds = %4
%17 = load i32, i32* %7
br label %20
18: ; preds = %4
%19 = load i32, i32* %8
br label %20
20: ; preds = %18, %16
%21 = phi i32 [ %17, %16 ], [ %19, %18 ]
store i32 %21, i32* %12
%22 = load i32, i32* %7
%23 = icmp eq i32 %22, 0
br i1 %23, label %27, label %24
24: ; preds = %20
%25 = load i32, i32* %8
%26 = icmp eq i32 %25, 0
br i1 %26, label %27, label %28
27: ; preds = %24, %20
store i32 0, i32* %5
br label %123
28: ; preds = %24
%29 = load i32, i32* %6
%30 = icmp ne i32 %29, 1
br i1 %30, label %31, label %35
31: ; preds = %28
%32 = load i32, i32* %6
%33 = icmp ne i32 %32, 2
br i1 %33, label %34, label %35
34: ; preds = %31
store i32 -1, i32* %5
br label %123
35: ; preds = %31, %28
%36 = load i32, i32* %6
%37 = icmp eq i32 %36, 1
%38 = zext i1 %37 to i32
%39 = sext i32 %38 to i64
%40 = icmp ne i64 %39, 0
br i1 %40, label %41, label %77
41: ; preds = %35
store i32 0, i32* %11
br label %42
42: ; preds = %73, %41
%43 = load i32, i32* %11
%44 = load i32, i32* %12
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %76
46: ; preds = %42
store i32 0, i32* %10
br label %47
47: ; preds = %69, %46
%48 = load i32, i32* %10
%49 = load i32, i32* %11
%50 = icmp slt i32 %48, %49
br i1 %50, label %51, label %72
51: ; preds = %47
%52 = load double*, double** %9
%53 = load i32, i32* %10
%54 = load i32, i32* %7
%55 = load i32, i32* %11
%56 = mul nsw i32 %54, %55
%57 = add nsw i32 %53, %56
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %52, i64 %58
%60 = load double, double* %59
%61 = load double*, double** %9
%62 = load i32, i32* %11
%63 = load i32, i32* %7
%64 = load i32, i32* %10
%65 = mul nsw i32 %63, %64
%66 = add nsw i32 %62, %65
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %61, i64 %67
store double %60, double* %68
br label %69
69: ; preds = %51
%70 = load i32, i32* %10
%71 = add nsw i32 %70, 1
store i32 %71, i32* %10
br label %47
72: ; preds = %47
br label %73
73: ; preds = %72
%74 = load i32, i32* %11
%75 = add nsw i32 %74, 1
store i32 %75, i32* %11
br label %42
76: ; preds = %42
br label %122
77: ; preds = %35
%78 = load i32, i32* %6
%79 = icmp eq i32 %78, 2
%80 = zext i1 %79 to i32
%81 = sext i32 %80 to i64
%82 = icmp ne i64 %81, 0
br i1 %82, label %83, label %121
83: ; preds = %77
store i32 0, i32* %11
br label %84
84: ; preds = %117, %83
%85 = load i32, i32* %11
%86 = load i32, i32* %12
%87 = icmp slt i32 %85, %86
br i1 %87, label %88, label %120
88: ; preds = %84
%89 = load i32, i32* %11
%90 = add nsw i32 %89, 1
store i32 %90, i32* %10
br label %91
91: ; preds = %113, %88
%92 = load i32, i32* %10
%93 = load i32, i32* %12
%94 = icmp slt i32 %92, %93
br i1 %94, label %95, label %116
95: ; preds = %91
%96 = load double*, double** %9
%97 = load i32, i32* %10
%98 = load i32, i32* %7
%99 = load i32, i32* %11
%100 = mul nsw i32 %98, %99
%101 = add nsw i32 %97, %100
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %96, i64 %102
%104 = load double, double* %103
%105 = load double*, double** %9
%106 = load i32, i32* %11
%107 = load i32, i32* %7
%108 = load i32, i32* %10
%109 = mul nsw i32 %107, %108
%110 = add nsw i32 %106, %109
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds double, double* %105, i64 %111
store double %104, double* %112
br label %113
113: ; preds = %95
%114 = load i32, i32* %10
%115 = add nsw i32 %114, 1
store i32 %115, i32* %10
br label %91
116: ; preds = %91
br label %117
117: ; preds = %116
%118 = load i32, i32* %11
%119 = add nsw i32 %118, 1
store i32 %119, i32* %11
br label %84
120: ; preds = %84
br label %121
121: ; preds = %120, %77
br label %122
122: ; preds = %121, %76
store i32 0, i32* %5
br label %123
123: ; preds = %122, %34, %27
%124 = load i32, i32* %5
ret i32 %124
}

LLVM IR (-O2):
define i32 @pcapack_symmetrize(i32 %0, i32 %1, i32 %2, double* nocapture %3) {
%5 = icmp slt i32 %1, %2
%6 = select i1 %5, i32 %1, i32 %2
%7 = icmp eq i32 %1, 0
%8 = icmp eq i32 %2, 0
%9 = or i1 %7, %8
br i1 %9, label %86, label %10
10: ; preds = %4
%11 = add i32 %0, -3
%12 = icmp ult i32 %11, -2
br i1 %12, label %86, label %13
13: ; preds = %10
switch i32 %0, label %86 [
i32 1, label %20
i32 2, label %14
]
14: ; preds = %13
%15 = icmp sgt i32 %6, 0
br i1 %15, label %16, label %86
16: ; preds = %14
%17 = sext i32 %1 to i64
%18 = zext i32 %6 to i64
%19 = zext i32 %6 to i64
br label %69
20: ; preds = %13
%21 = icmp sgt i32 %6, 0
br i1 %21, label %22, label %86
22: ; preds = %20
%23 = sext i32 %1 to i64
%24 = zext i32 %6 to i64
br label %25
25: ; preds = %22, %63
%26 = phi i64 [ 0, %22 ], [ %64, %63 ]
%27 = icmp eq i64 %26, 0
br i1 %27, label %63, label %28
28: ; preds = %25
%29 = mul nsw i64 %26, %23
%30 = and i64 %26, 1
%31 = icmp eq i64 %26, 1
br i1 %31, label %53, label %32
32: ; preds = %28
%33 = and i64 %26, 9223372036854775806
br label %34
34: ; preds = %34, %32
%35 = phi i64 [ 0, %32 ], [ %50, %34 ]
%36 = phi i64 [ 0, %32 ], [ %51, %34 ]
%37 = add nsw i64 %35, %29
%38 = getelementptr inbounds double, double* %3, i64 %37
%39 = load double, double* %38
%40 = mul nsw i64 %35, %23
%41 = add nsw i64 %40, %26
%42 = getelementptr inbounds double, double* %3, i64 %41
store double %39, double* %42
%43 = or i64 %35, 1
%44 = add nsw i64 %43, %29
%45 = getelementptr inbounds double, double* %3, i64 %44
%46 = load double, double* %45
%47 = mul nsw i64 %43, %23
%48 = add nsw i64 %47, %26
%49 = getelementptr inbounds double, double* %3, i64 %48
store double %46, double* %49
%50 = add nuw nsw i64 %35, 2
%51 = add i64 %36, 2
%52 = icmp eq i64 %51, %33
br i1 %52, label %53, label %34
53: ; preds = %34, %28
%54 = phi i64 [ 0, %28 ], [ %50, %34 ]
%55 = icmp eq i64 %30, 0
br i1 %55, label %63, label %56
56: ; preds = %53
%57 = add nsw i64 %54, %29
%58 = getelementptr inbounds double, double* %3, i64 %57
%59 = load double, double* %58
%60 = mul nsw i64 %54, %23
%61 = add nsw i64 %60, %26
%62 = getelementptr inbounds double, double* %3, i64 %61
store double %59, double* %62
br label %63
63: ; preds = %56, %53, %25
%64 = add nuw nsw i64 %26, 1
%65 = icmp eq i64 %64, %24
br i1 %65, label %86, label %25
66: ; preds = %76, %69
%67 = add nuw nsw i64 %71, 1
%68 = icmp eq i64 %72, %19
br i1 %68, label %86, label %69
69: ; preds = %16, %66
%70 = phi i64 [ 0, %16 ], [ %72, %66 ]
%71 = phi i64 [ 1, %16 ], [ %67, %66 ]
%72 = add nuw nsw i64 %70, 1
%73 = icmp ult i64 %72, %18
br i1 %73, label %74, label %66
74: ; preds = %69
%75 = mul nsw i64 %70, %17
br label %76
76: ; preds = %74, %76
%77 = phi i64 [ %71, %74 ], [ %84, %76 ]
%78 = add nsw i64 %77, %75
%79 = getelementptr inbounds double, double* %3, i64 %78
%80 = load double, double* %79
%81 = mul nsw i64 %77, %17
%82 = add nsw i64 %81, %70
%83 = getelementptr inbounds double, double* %3, i64 %82
store double %80, double* %83
%84 = add nuw nsw i64 %77, 1
%85 = icmp eq i64 %84, %19
br i1 %85, label %66, label %76
86: ; preds = %66, %63, %14, %20, %13, %10, %4
%87 = phi i32 [ 0, %4 ], [ -1, %10 ], [ 0, %13 ], [ 0, %20 ], [ 0, %14 ], [ 0, %63 ], [ 0, %66 ]
ret i32 %87
}