double find_closest_pair(int n, double **distmatrix, int *ip, int *jp) {
  int i, j;
  double temp;
  double distance = distmatrix[1][0];
  *ip = 1;
  *jp = 0;
  for (i = 1; i < n; i++) {
    for (j = 0; j < i; j++) {
      temp = distmatrix[i][j];
      if (temp < distance) {
        distance = temp;
        *ip = i;
        *jp = j;
      }
    }
  }
  return distance;
}


x86 assembly (clang, -O0):
<find_closest_pair>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movl   $0x1,(%rax)
mov    -0x20(%rbp),%rax
movl   $0x0,(%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11ce <find_closest_pair+0xce>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11bb <find_closest_pair+0xbb>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11a8 <find_closest_pair+0xa8>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x28(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
jmp    11ad <find_closest_pair+0xad>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1152 <find_closest_pair+0x52>
jmp    11c0 <find_closest_pair+0xc0>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    113f <find_closest_pair+0x3f>
movsd  -0x38(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_closest_pair>:
push   %rbx
mov    0x8(%rsi),%rax
movsd  (%rax),%xmm0
movl   $0x1,(%rdx)
movl   $0x0,(%rcx)
cmp    $0x2,%edi
jl     11b8 <find_closest_pair+0xb8>
mov    %edi,%r8d
mov    $0x1,%r10d
xor    %r9d,%r9d
test   %r9,%r9
je     1180 <find_closest_pair+0x80>
nop
mov    %r10,%r11
and    $0xfffffffffffffffe,%r11
xor    %edi,%edi
jmp    1149 <find_closest_pair+0x49>
nopl   0x0(%rax,%rax,1)
add    $0x2,%rdi
cmp    %rdi,%r11
je     1182 <find_closest_pair+0x82>
movsd  (%rax,%rdi,8),%xmm1
ucomisd %xmm1,%xmm0
jbe    115d <find_closest_pair+0x5d>
mov    %r10d,(%rdx)
mov    %edi,(%rcx)
movapd %xmm1,%xmm0
movsd  0x8(%rax,%rdi,8),%xmm1
ucomisd %xmm1,%xmm0
jbe    1140 <find_closest_pair+0x40>
mov    %r10d,(%rdx)
lea    0x1(%rdi),%ebx
mov    %ebx,(%rcx)
movapd %xmm1,%xmm0
jmp    1140 <find_closest_pair+0x40>
nopw   0x0(%rax,%rax,1)
xor    %edi,%edi
test   $0x1,%r10b
je     119c <find_closest_pair+0x9c>
movsd  (%rax,%rdi,8),%xmm1
ucomisd %xmm1,%xmm0
jbe    119c <find_closest_pair+0x9c>
mov    %r10d,(%rdx)
mov    %edi,(%rcx)
movapd %xmm1,%xmm0
add    $0x1,%r10
cmp    %r8,%r10
je     11b8 <find_closest_pair+0xb8>
mov    (%rsi,%r10,8),%rax
add    $0x1,%r9
test   %r9,%r9
jne    1130 <find_closest_pair+0x30>
jmp    1180 <find_closest_pair+0x80>
pop    %rbx
ret

LLVM IR (-O0):
define double @find_closest_pair(i32 %0, double** %1, i32* %2, i32* %3) {
%5 = alloca i32
%6 = alloca double**
%7 = alloca i32*
%8 = alloca i32*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca double
store i32 %0, i32* %5
store double** %1, double*** %6
store i32* %2, i32** %7
store i32* %3, i32** %8
%13 = load double**, double*** %6
%14 = getelementptr inbounds double*, double** %13, i64 1
%15 = load double*, double** %14
%16 = getelementptr inbounds double, double* %15, i64 0
%17 = load double, double* %16
store double %17, double* %12
%18 = load i32*, i32** %7
store i32 1, i32* %18
%19 = load i32*, i32** %8
store i32 0, i32* %19
store i32 1, i32* %9
br label %20
20: ; preds = %53, %4
%21 = load i32, i32* %9
%22 = load i32, i32* %5
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %56
24: ; preds = %20
store i32 0, i32* %10
br label %25
25: ; preds = %49, %24
%26 = load i32, i32* %10
%27 = load i32, i32* %9
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %52
29: ; preds = %25
%30 = load double**, double*** %6
%31 = load i32, i32* %9
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double*, double** %30, i64 %32
%34 = load double*, double** %33
%35 = load i32, i32* %10
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
store double %38, double* %11
%39 = load double, double* %11
%40 = load double, double* %12
%41 = fcmp olt double %39, %40
br i1 %41, label %42, label %48
42: ; preds = %29
%43 = load double, double* %11
store double %43, double* %12
%44 = load i32, i32* %9
%45 = load i32*, i32** %7
store i32 %44, i32* %45
%46 = load i32, i32* %10
%47 = load i32*, i32** %8
store i32 %46, i32* %47
br label %48
48: ; preds = %42, %29
br label %49
49: ; preds = %48
%50 = load i32, i32* %10
%51 = add nsw i32 %50, 1
store i32 %51, i32* %10
br label %25
52: ; preds = %25
br label %53
53: ; preds = %52
%54 = load i32, i32* %9
%55 = add nsw i32 %54, 1
store i32 %55, i32* %9
br label %20
56: ; preds = %20
%57 = load double, double* %12
ret double %57
}

LLVM IR (-O2):
define double @find_closest_pair(i32 %0, double** nocapture %1, i32* nocapture writeonly %2, i32* nocapture writeonly %3) {
%5 = getelementptr inbounds double*, double** %1, i64 1
%6 = load double*, double** %5
%7 = load double, double* %6
store i32 1, i32* %2
store i32 0, i32* %3
%8 = icmp sgt i32 %0, 1
br i1 %8, label %9, label %62
9: ; preds = %4
%10 = zext i32 %0 to i64
br label %11
11: ; preds = %58, %9
%12 = phi i64 [ %61, %58 ], [ 0, %9 ]
%13 = phi double* [ %60, %58 ], [ %6, %9 ]
%14 = phi i64 [ %56, %58 ], [ 1, %9 ]
%15 = phi double [ %55, %58 ], [ %7, %9 ]
%16 = trunc i64 %14 to i32
%17 = and i64 %14, 1
%18 = icmp eq i64 %12, 0
br i1 %18, label %43, label %19
19: ; preds = %11
%20 = and i64 %14, 9223372036854775806
br label %21
21: ; preds = %38, %19
%22 = phi i64 [ 0, %19 ], [ %40, %38 ]
%23 = phi double [ %15, %19 ], [ %39, %38 ]
%24 = phi i64 [ 0, %19 ], [ %41, %38 ]
%25 = getelementptr inbounds double, double* %13, i64 %22
%26 = load double, double* %25
%27 = fcmp olt double %26, %23
br i1 %27, label %28, label %30
28: ; preds = %21
store i32 %16, i32* %2
%29 = trunc i64 %22 to i32
store i32 %29, i32* %3
br label %30
30: ; preds = %21, %28
%31 = phi double [ %26, %28 ], [ %23, %21 ]
%32 = or i64 %22, 1
%33 = getelementptr inbounds double, double* %13, i64 %32
%34 = load double, double* %33
%35 = fcmp olt double %34, %31
br i1 %35, label %36, label %38
36: ; preds = %30
store i32 %16, i32* %2
%37 = trunc i64 %32 to i32
store i32 %37, i32* %3
br label %38
38: ; preds = %36, %30
%39 = phi double [ %34, %36 ], [ %31, %30 ]
%40 = add nuw nsw i64 %22, 2
%41 = add i64 %24, 2
%42 = icmp eq i64 %41, %20
br i1 %42, label %43, label %21
43: ; preds = %38, %11
%44 = phi double [ undef, %11 ], [ %39, %38 ]
%45 = phi i64 [ 0, %11 ], [ %40, %38 ]
%46 = phi double [ %15, %11 ], [ %39, %38 ]
%47 = icmp eq i64 %17, 0
br i1 %47, label %54, label %48
48: ; preds = %43
%49 = getelementptr inbounds double, double* %13, i64 %45
%50 = load double, double* %49
%51 = fcmp olt double %50, %46
br i1 %51, label %52, label %54
52: ; preds = %48
store i32 %16, i32* %2
%53 = trunc i64 %45 to i32
store i32 %53, i32* %3
br label %54
54: ; preds = %48, %52, %43
%55 = phi double [ %44, %43 ], [ %50, %52 ], [ %46, %48 ]
%56 = add nuw nsw i64 %14, 1
%57 = icmp eq i64 %56, %10
br i1 %57, label %62, label %58
58: ; preds = %54
%59 = getelementptr inbounds double*, double** %1, i64 %56
%60 = load double*, double** %59
%61 = add i64 %12, 1
br label %11
62: ; preds = %54, %4
%63 = phi double [ %7, %4 ], [ %55, %54 ]
ret double %63
}