int ijkl2intindex(int i, int j, int k, int l) {
  int tmp, ij, kl;
  if (i < j)
    return ijkl2intindex(j, i, k, l);
  if (k < l)
    return ijkl2intindex(i, j, l, k);
  ij = i * (i + 1) / 2 + j;
  kl = k * (k + 1) / 2 + l;
  if (ij < kl) {
    tmp = ij;
    ij = kl;
    kl = tmp;
  }
  return ij * (ij + 1) / 2 + kl;
}


x86 assembly (clang, -O0):
<ijkl2intindex>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1149 <ijkl2intindex+0x39>
mov    -0xc(%rbp),%edi
mov    -0x8(%rbp),%esi
mov    -0x10(%rbp),%edx
mov    -0x14(%rbp),%ecx
call   1030 <ijkl2intindex@plt>
mov    %eax,-0x4(%rbp)
jmp    11da <ijkl2intindex+0xca>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    116e <ijkl2intindex+0x5e>
mov    -0x8(%rbp),%edi
mov    -0xc(%rbp),%esi
mov    -0x14(%rbp),%edx
mov    -0x10(%rbp),%ecx
call   1030 <ijkl2intindex@plt>
mov    %eax,-0x4(%rbp)
jmp    11da <ijkl2intindex+0xca>
mov    -0x8(%rbp),%eax
mov    -0x8(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0xc(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x14(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11c0 <ijkl2intindex+0xb0>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x20(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ijkl2intindex>:
mov    %edi,%r8d
mov    %edx,%eax
mov    %ecx,%edi
nopw   0x0(%rax,%rax,1)
cmp    %esi,%r8d
jl     1130 <ijkl2intindex+0x30>
mov    %edi,%r9d
mov    %eax,%edi
mov    %r9d,%eax
cmp    %r9d,%edi
jl     1110 <ijkl2intindex+0x10>
jmp    1137 <ijkl2intindex+0x37>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %esi,%edi
mov    %r8d,%esi
jmp    1100 <ijkl2intindex>
lea    0x1(%r8),%eax
imul   %r8d,%eax
mov    %eax,%ecx
shr    $0x1f,%ecx
add    %eax,%ecx
sar    %ecx
add    %esi,%ecx
lea    0x1(%rdi),%eax
imul   %edi,%eax
mov    %eax,%edx
shr    $0x1f,%edx
add    %eax,%edx
sar    %edx
add    %r9d,%edx
cmp    %edx,%ecx
mov    %ecx,%eax
cmovl  %edx,%eax
cmovl  %ecx,%edx
lea    0x1(%rax),%ecx
imul   %eax,%ecx
mov    %ecx,%eax
shr    $0x1f,%eax
add    %ecx,%eax
sar    %eax
add    %edx,%eax
ret

LLVM IR (-O0):
define i32 @ijkl2intindex(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%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 i32 %3, i32* %9
%13 = load i32, i32* %6
%14 = load i32, i32* %7
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %22
16: ; preds = %4
%17 = load i32, i32* %7
%18 = load i32, i32* %6
%19 = load i32, i32* %8
%20 = load i32, i32* %9
%21 = call i32 @ijkl2intindex(i32 %17, i32 %18, i32 %19, i32 %20)
store i32 %21, i32* %5
br label %62
22: ; preds = %4
%23 = load i32, i32* %8
%24 = load i32, i32* %9
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %32
26: ; preds = %22
%27 = load i32, i32* %6
%28 = load i32, i32* %7
%29 = load i32, i32* %9
%30 = load i32, i32* %8
%31 = call i32 @ijkl2intindex(i32 %27, i32 %28, i32 %29, i32 %30)
store i32 %31, i32* %5
br label %62
32: ; preds = %22
%33 = load i32, i32* %6
%34 = load i32, i32* %6
%35 = add nsw i32 %34, 1
%36 = mul nsw i32 %33, %35
%37 = sdiv i32 %36, 2
%38 = load i32, i32* %7
%39 = add nsw i32 %37, %38
store i32 %39, i32* %11
%40 = load i32, i32* %8
%41 = load i32, i32* %8
%42 = add nsw i32 %41, 1
%43 = mul nsw i32 %40, %42
%44 = sdiv i32 %43, 2
%45 = load i32, i32* %9
%46 = add nsw i32 %44, %45
store i32 %46, i32* %12
%47 = load i32, i32* %11
%48 = load i32, i32* %12
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %54
50: ; preds = %32
%51 = load i32, i32* %11
store i32 %51, i32* %10
%52 = load i32, i32* %12
store i32 %52, i32* %11
%53 = load i32, i32* %10
store i32 %53, i32* %12
br label %54
54: ; preds = %50, %32
%55 = load i32, i32* %11
%56 = load i32, i32* %11
%57 = add nsw i32 %56, 1
%58 = mul nsw i32 %55, %57
%59 = sdiv i32 %58, 2
%60 = load i32, i32* %12
%61 = add nsw i32 %59, %60
store i32 %61, i32* %5
br label %62
62: ; preds = %54, %26, %16
%63 = load i32, i32* %5
ret i32 %63
}

LLVM IR (-O2):
define i32 @ijkl2intindex(i32 %0, i32 %1, i32 %2, i32 %3) {
br label %5
5: ; preds = %9, %4
%6 = phi i32 [ %0, %4 ], [ %7, %9 ]
%7 = phi i32 [ %1, %4 ], [ %6, %9 ]
%8 = icmp slt i32 %6, %7
br label %9
9: ; preds = %5, %12
%10 = phi i32 [ %11, %12 ], [ %2, %5 ]
%11 = phi i32 [ %10, %12 ], [ %3, %5 ]
br i1 %8, label %5, label %12
12: ; preds = %9
%13 = icmp slt i32 %10, %11
br i1 %13, label %9, label %14
14: ; preds = %12
%15 = add nsw i32 %6, 1
%16 = mul nsw i32 %15, %6
%17 = sdiv i32 %16, 2
%18 = add nsw i32 %17, %7
%19 = add nsw i32 %10, 1
%20 = mul nsw i32 %19, %10
%21 = sdiv i32 %20, 2
%22 = add nsw i32 %21, %11
%23 = icmp slt i32 %18, %22
%24 = select i1 %23, i32 %22, i32 %18
%25 = select i1 %23, i32 %18, i32 %22
%26 = add nsw i32 %24, 1
%27 = mul nsw i32 %26, %24
%28 = sdiv i32 %27, 2
%29 = add nsw i32 %28, %25
ret i32 %29
}