#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ daxpy_(int *, double *, double *, int *, double *, int *);
double dnrm2(int, double *, int);
void bestdist(int n, int ncol, double *restrict x, double *restrict a,
              double *restrict b, double *restrict work) {
  int i, j;
  int ia = 0;
  int intone = 1;
  double tmp, best;
  double one = 1.;
  best = 0.;
  for (i = 0; i < n; i++) {
    for (j = 0; j < ncol; j++)
      work[j] = -1. * x[i + n * j];
    daxpy_(&ncol, &one, b, &intone, work, &intone);
    tmp = dnrm2(ncol, work, 1);
    if (tmp > best) {
      best = tmp;
      ia = i;
    }
  }
  for (j = 0; j < ncol; j++)
    a[j] = x[ia + n * j];
}


x86 assembly (clang, -O0):
<bestdist>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x1,-0x38(%rbp)
movsd  0xeac(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1227 <bestdist+0x107>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11c2 <bestdist+0xa2>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x30(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xe66(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    117b <bestdist+0x5b>
mov    -0x20(%rbp),%rdx
mov    -0x28(%rbp),%r8
lea    -0x8(%rbp),%rdi
lea    -0x50(%rbp),%rsi
lea    -0x38(%rbp),%r9
mov    %r9,%rcx
call   1040 <daxpy_@plt>
mov    -0x8(%rbp),%edi
mov    -0x28(%rbp),%rsi
mov    $0x1,%edx
call   1030 <dnrm2@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
ucomisd -0x48(%rbp),%xmm0
jbe    1214 <bestdist+0xf4>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
jmp    1219 <bestdist+0xf9>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1168 <bestdist+0x48>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    126d <bestdist+0x14d>
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x30(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    122e <bestdist+0x10e>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bestdist>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8,0x28(%rsp)
mov    %esi,0x4(%rsp)
movl   $0x1,0xc(%rsp)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x30(%rsp)
xor    %r13d,%r13d
test   %edi,%edi
jle    1280 <bestdist+0x160>
mov    %r9,%rbp
mov    %rcx,0x20(%rsp)
mov    %edi,0x8(%rsp)
mov    %edi,%r15d
mov    %r15,%r12
shl    $0x4,%r12
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x18(%rsp)
movapd 0xe7f(%rip),%xmm1
mov    %rdx,0x10(%rsp)
mov    %rdx,%r14
xor    %ebx,%ebx
test   %esi,%esi
jle    1200 <bestdist+0xe0>
nop
mov    %esi,%eax
cmp    $0x1,%esi
jne    11a0 <bestdist+0x80>
xor    %ecx,%ecx
jmp    11da <bestdist+0xba>
nopl   0x0(%rax,%rax,1)
mov    %eax,%edx
and    $0xfffffffe,%edx
mov    %r14,%rsi
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movsd  (%rsi),%xmm0
xorpd  %xmm1,%xmm0
movlpd %xmm0,0x0(%rbp,%rcx,8)
movsd  (%rsi,%r15,8),%xmm0
xorpd  %xmm1,%xmm0
movlpd %xmm0,0x8(%rbp,%rcx,8)
add    $0x2,%rcx
add    %r12,%rsi
cmp    %rcx,%rdx
jne    11b0 <bestdist+0x90>
test   $0x1,%al
je     1200 <bestdist+0xe0>
mov    %rcx,%rax
imul   %r15,%rax
add    %rbx,%rax
mov    0x10(%rsp),%rdx
movsd  (%rdx,%rax,8),%xmm0
xorpd  0xe06(%rip),%xmm0
movlpd %xmm0,0x0(%rbp,%rcx,8)
lea    0x4(%rsp),%rdi
lea    0x30(%rsp),%rsi
mov    0x28(%rsp),%rdx
lea    0xc(%rsp),%rcx
mov    %rbp,%r8
mov    %rcx,%r9
call   1040 <daxpy_@plt>
mov    0x4(%rsp),%edi
mov    %rbp,%rsi
mov    $0x1,%edx
call   1030 <dnrm2@plt>
movsd  0x18(%rsp),%xmm1
ucomisd %xmm1,%xmm0
cmova  %ebx,%r13d
add    $0x1,%rbx
cmp    %r15,%rbx
je     126b <bestdist+0x14b>
maxsd  %xmm1,%xmm0
mov    0x4(%rsp),%esi
add    $0x8,%r14
movsd  %xmm0,0x18(%rsp)
movapd 0xd9f(%rip),%xmm1
test   %esi,%esi
jg     1190 <bestdist+0x70>
jmp    1200 <bestdist+0xe0>
mov    0x4(%rsp),%esi
movslq %r13d,%r13
mov    0x10(%rsp),%rdx
mov    0x20(%rsp),%rcx
mov    0x8(%rsp),%edi
test   %esi,%esi
jle    1355 <bestdist+0x235>
movslq %edi,%r8
mov    %esi,%r10d
lea    -0x1(%r10),%rax
mov    %r10d,%r9d
and    $0x3,%r9d
cmp    $0x3,%rax
jae    12a3 <bestdist+0x183>
xor    %esi,%esi
jmp    131d <bestdist+0x1fd>
and    $0xfffffffc,%r10d
mov    %r8,%rax
shl    $0x4,%rax
lea    0x0(,%r13,8),%r14
lea    (%rax,%r13,8),%r11
mov    %r8,%rax
shl    $0x5,%rax
lea    (%r8,%r13,1),%rbx
lea    (%r8,%r8,2),%rsi
lea    (%r14,%rsi,8),%rdi
xor    %esi,%esi
mov    %rdx,%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%rbp,%r14,1),%xmm0
movsd  %xmm0,(%rcx,%rsi,8)
movsd  0x0(%rbp,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rsi,8)
movsd  0x0(%rbp,%r11,1),%xmm0
movsd  %xmm0,0x10(%rcx,%rsi,8)
movsd  0x0(%rbp,%rdi,1),%xmm0
movsd  %xmm0,0x18(%rcx,%rsi,8)
add    $0x4,%rsi
add    %rax,%rbp
cmp    %rsi,%r10
jne    12e0 <bestdist+0x1c0>
test   %r9,%r9
je     1355 <bestdist+0x235>
lea    (%rcx,%rsi,8),%rax
imul   %r8,%rsi
add    %rsi,%r13
lea    (%rdx,%r13,8),%rcx
shl    $0x3,%r8
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movsd  (%rcx),%xmm0
movsd  %xmm0,(%rax,%rdx,8)
add    $0x1,%rdx
add    %r8,%rcx
cmp    %rdx,%r9
jne    1340 <bestdist+0x220>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @bestdist(i32 %0, i32 %1, double* noalias %2, double* noalias %3, double* noalias %4, double* noalias %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
%12 = alloca double*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca double
%18 = alloca double
%19 = alloca double
store i32 %0, i32* %7
store i32 %1, i32* %8
store double* %2, double** %9
store double* %3, double** %10
store double* %4, double** %11
store double* %5, double** %12
store i32 0, i32* %15
store i32 1, i32* %16
store double 1.000000e+00, double* %19
store double 0.000000e+00, double* %18
store i32 0, i32* %13
br label %20
20: ; preds = %61, %6
%21 = load i32, i32* %13
%22 = load i32, i32* %7
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %64
24: ; preds = %20
store i32 0, i32* %14
br label %25
25: ; preds = %44, %24
%26 = load i32, i32* %14
%27 = load i32, i32* %8
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %47
29: ; preds = %25
%30 = load double*, double** %9
%31 = load i32, i32* %13
%32 = load i32, i32* %7
%33 = load i32, i32* %14
%34 = mul nsw i32 %32, %33
%35 = add nsw i32 %31, %34
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %30, i64 %36
%38 = load double, double* %37
%39 = fmul double -1.000000e+00, %38
%40 = load double*, double** %12
%41 = load i32, i32* %14
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %40, i64 %42
store double %39, double* %43
br label %44
44: ; preds = %29
%45 = load i32, i32* %14
%46 = add nsw i32 %45, 1
store i32 %46, i32* %14
br label %25
47: ; preds = %25
%48 = load double*, double** %11
%49 = load double*, double** %12
%50 = call i32 @daxpy_(i32* %8, double* %19, double* %48, i32* %16, double* %49, i32* %16)
%51 = load i32, i32* %8
%52 = load double*, double** %12
%53 = call double @dnrm2(i32 %51, double* %52, i32 1)
store double %53, double* %17
%54 = load double, double* %17
%55 = load double, double* %18
%56 = fcmp ogt double %54, %55
br i1 %56, label %57, label %60
57: ; preds = %47
%58 = load double, double* %17
store double %58, double* %18
%59 = load i32, i32* %13
store i32 %59, i32* %15
br label %60
60: ; preds = %57, %47
br label %61
61: ; preds = %60
%62 = load i32, i32* %13
%63 = add nsw i32 %62, 1
store i32 %63, i32* %13
br label %20
64: ; preds = %20
store i32 0, i32* %14
br label %65
65: ; preds = %83, %64
%66 = load i32, i32* %14
%67 = load i32, i32* %8
%68 = icmp slt i32 %66, %67
br i1 %68, label %69, label %86
69: ; preds = %65
%70 = load double*, double** %9
%71 = load i32, i32* %15
%72 = load i32, i32* %7
%73 = load i32, i32* %14
%74 = mul nsw i32 %72, %73
%75 = add nsw i32 %71, %74
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %70, i64 %76
%78 = load double, double* %77
%79 = load double*, double** %10
%80 = load i32, i32* %14
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %79, i64 %81
store double %78, double* %82
br label %83
83: ; preds = %69
%84 = load i32, i32* %14
%85 = add nsw i32 %84, 1
store i32 %85, i32* %14
br label %65
86: ; preds = %65
ret void
}

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