#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 */
double sqrt(double);
double calc_dmin(double **p, size_t ndim, size_t i) {
  size_t j, k;
  double d, di_min = 0;
  for (j = 0; j < ndim; j++) {
    d = 0;
    for (k = 0; k < ndim; k++) {
      if (k != j)
        d += p[i][k] * p[i][k];
      else
        d += (p[i][k] - 1) * (p[i][k] - 1);
    }
    d = sqrt(d);
    if (j == 0 || d < di_min)
      di_min = d;
  }
  return di_min;
}


x86 assembly (clang, -O0):
<calc_dmin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movq   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    1267 <calc_dmin+0x157>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    1219 <calc_dmin+0x109>
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
je     11af <calc_dmin+0x9f>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    1203 <calc_dmin+0xf3>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe34(%rip),%xmm1
subsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  0xe13(%rip),%xmm1
subsd  %xmm1,%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    1208 <calc_dmin+0xf8>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1152 <calc_dmin+0x42>
movsd  -0x30(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x30(%rbp)
cmpq   $0x0,-0x20(%rbp)
je     1247 <calc_dmin+0x137>
movsd  -0x30(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1251 <calc_dmin+0x141>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    1256 <calc_dmin+0x146>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1134 <calc_dmin+0x24>
movsd  -0x38(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_dmin>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
test   %rsi,%rsi
je     1266 <calc_dmin+0x156>
mov    %rdx,%r14
mov    %rsi,%r12
mov    %rdi,%r15
mov    %rsi,%r13
and    $0xfffffffffffffffe,%r13
xorpd  %xmm2,%xmm2
xor    %ebx,%ebx
movsd  0xebc(%rip),%xmm3
xorpd  %xmm4,%xmm4
jmp    1161 <calc_dmin+0x51>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rbx
movapd %xmm0,%xmm4
cmp    %r12,%rbx
je     126a <calc_dmin+0x15a>
cmp    $0x1,%r12
jne    1190 <calc_dmin+0x80>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
test   $0x1,%r12b
jne    11fe <calc_dmin+0xee>
ucomisd %xmm2,%xmm0
jb     1237 <calc_dmin+0x127>
sqrtsd %xmm0,%xmm0
test   %rbx,%rbx
je     1150 <calc_dmin+0x40>
jmp    125d <calc_dmin+0x14d>
nop
mov    (%r15,%r14,8),%rcx
xorpd  %xmm0,%xmm0
xor    %edx,%edx
jmp    11ba <calc_dmin+0xaa>
nopl   0x0(%rax)
movsd  (%rcx,%rbx,8),%xmm0
addsd  %xmm3,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %r13,%rdx
je     11f0 <calc_dmin+0xe0>
mov    %rdx,%rax
cmp    %rdx,%rbx
jne    11d0 <calc_dmin+0xc0>
movsd  (%rcx,%rbx,8),%xmm1
addsd  %xmm3,%xmm1
jmp    11d5 <calc_dmin+0xc5>
nopl   (%rax)
movsd  (%rcx,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
lea    0x1(%rax),%rdx
cmp    %rbx,%rdx
je     11a0 <calc_dmin+0x90>
movsd  0x8(%rcx,%rax,8),%xmm0
jmp    11a9 <calc_dmin+0x99>
xchg   %ax,%ax
add    $0x2,%rax
test   $0x1,%r12b
je     1177 <calc_dmin+0x67>
mov    (%r15,%r14,8),%rcx
cmp    %rbx,%rax
jne    1220 <calc_dmin+0x110>
movsd  (%rcx,%rbx,8),%xmm1
addsd  %xmm3,%xmm1
jmp    1225 <calc_dmin+0x115>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rcx,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
ucomisd %xmm2,%xmm0
jae    1181 <calc_dmin+0x71>
movsd  %xmm4,0x8(%rsp)
call   1030 <sqrt@plt>
movsd  0x8(%rsp),%xmm4
movsd  0xdb0(%rip),%xmm3
xorpd  %xmm2,%xmm2
test   %rbx,%rbx
je     1150 <calc_dmin+0x40>
minsd  %xmm4,%xmm0
jmp    1150 <calc_dmin+0x40>
xorpd  %xmm0,%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define double @calc_dmin(double** %0, i64 %1, i64 %2) {
%4 = alloca double**
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca double
%10 = alloca double
store double** %0, double*** %4
store i64 %1, i64* %5
store i64 %2, i64* %6
store double 0.000000e+00, double* %10
store i64 0, i64* %7
br label %11
11: ; preds = %76, %3
%12 = load i64, i64* %7
%13 = load i64, i64* %5
%14 = icmp ult i64 %12, %13
br i1 %14, label %15, label %79
15: ; preds = %11
store double 0.000000e+00, double* %9
store i64 0, i64* %8
br label %16
16: ; preds = %61, %15
%17 = load i64, i64* %8
%18 = load i64, i64* %5
%19 = icmp ult i64 %17, %18
br i1 %19, label %20, label %64
20: ; preds = %16
%21 = load i64, i64* %8
%22 = load i64, i64* %7
%23 = icmp ne i64 %21, %22
br i1 %23, label %24, label %41
24: ; preds = %20
%25 = load double**, double*** %4
%26 = load i64, i64* %6
%27 = getelementptr inbounds double*, double** %25, i64 %26
%28 = load double*, double** %27
%29 = load i64, i64* %8
%30 = getelementptr inbounds double, double* %28, i64 %29
%31 = load double, double* %30
%32 = load double**, double*** %4
%33 = load i64, i64* %6
%34 = getelementptr inbounds double*, double** %32, i64 %33
%35 = load double*, double** %34
%36 = load i64, i64* %8
%37 = getelementptr inbounds double, double* %35, i64 %36
%38 = load double, double* %37
%39 = load double, double* %9
%40 = call double @llvm.fmuladd.f64(double %31, double %38, double %39)
store double %40, double* %9
br label %60
41: ; preds = %20
%42 = load double**, double*** %4
%43 = load i64, i64* %6
%44 = getelementptr inbounds double*, double** %42, i64 %43
%45 = load double*, double** %44
%46 = load i64, i64* %8
%47 = getelementptr inbounds double, double* %45, i64 %46
%48 = load double, double* %47
%49 = fsub double %48, 1.000000e+00
%50 = load double**, double*** %4
%51 = load i64, i64* %6
%52 = getelementptr inbounds double*, double** %50, i64 %51
%53 = load double*, double** %52
%54 = load i64, i64* %8
%55 = getelementptr inbounds double, double* %53, i64 %54
%56 = load double, double* %55
%57 = fsub double %56, 1.000000e+00
%58 = load double, double* %9
%59 = call double @llvm.fmuladd.f64(double %49, double %57, double %58)
store double %59, double* %9
br label %60
60: ; preds = %41, %24
br label %61
61: ; preds = %60
%62 = load i64, i64* %8
%63 = add i64 %62, 1
store i64 %63, i64* %8
br label %16
64: ; preds = %16
%65 = load double, double* %9
%66 = call double @sqrt(double %65)
store double %66, double* %9
%67 = load i64, i64* %7
%68 = icmp eq i64 %67, 0
br i1 %68, label %73, label %69
69: ; preds = %64
%70 = load double, double* %9
%71 = load double, double* %10
%72 = fcmp olt double %70, %71
br i1 %72, label %73, label %75
73: ; preds = %69, %64
%74 = load double, double* %9
store double %74, double* %10
br label %75
75: ; preds = %73, %69
br label %76
76: ; preds = %75
%77 = load i64, i64* %7
%78 = add i64 %77, 1
store i64 %78, i64* %7
br label %11
79: ; preds = %11
%80 = load double, double* %10
ret double %80
}

LLVM IR (-O2):
define double @calc_dmin(double** nocapture %0, i64 %1, i64 %2) {
%4 = icmp eq i64 %1, 0
br i1 %4, label %74, label %5
5: ; preds = %3
%6 = getelementptr inbounds double*, double** %0, i64 %2
%7 = and i64 %1, 1
%8 = icmp eq i64 %1, 1
%9 = and i64 %1, -2
%10 = icmp eq i64 %7, 0
br label %11
11: ; preds = %5, %65
%12 = phi double [ 0.000000e+00, %5 ], [ %71, %65 ]
%13 = phi i64 [ 0, %5 ], [ %72, %65 ]
br i1 %8, label %49, label %14
14: ; preds = %11
%15 = load double*, double** %6
%16 = getelementptr inbounds double, double* %15, i64 %13
%17 = load double*, double** %6
%18 = getelementptr inbounds double, double* %17, i64 %13
br label %19
19: ; preds = %44, %14
%20 = phi double [ 0.000000e+00, %14 ], [ %45, %44 ]
%21 = phi i64 [ 0, %14 ], [ %46, %44 ]
%22 = phi i64 [ 0, %14 ], [ %47, %44 ]
%23 = icmp eq i64 %21, %13
br i1 %23, label %28, label %24
24: ; preds = %19
%25 = getelementptr inbounds double, double* %15, i64 %21
%26 = load double, double* %25
%27 = tail call double @llvm.fmuladd.f64(double %26, double %26, double %20)
br label %32
28: ; preds = %19
%29 = load double, double* %16
%30 = fadd double %29, -1.000000e+00
%31 = tail call double @llvm.fmuladd.f64(double %30, double %30, double %20)
br label %32
32: ; preds = %24, %28
%33 = phi double [ %27, %24 ], [ %31, %28 ]
%34 = or i64 %21, 1
%35 = icmp eq i64 %34, %13
br i1 %35, label %40, label %36
36: ; preds = %32
%37 = getelementptr inbounds double, double* %17, i64 %34
%38 = load double, double* %37
%39 = tail call double @llvm.fmuladd.f64(double %38, double %38, double %33)
br label %44
40: ; preds = %32
%41 = load double, double* %18
%42 = fadd double %41, -1.000000e+00
%43 = tail call double @llvm.fmuladd.f64(double %42, double %42, double %33)
br label %44
44: ; preds = %40, %36
%45 = phi double [ %39, %36 ], [ %43, %40 ]
%46 = add nuw i64 %21, 2
%47 = add i64 %22, 2
%48 = icmp eq i64 %47, %9
br i1 %48, label %49, label %19
49: ; preds = %44, %11
%50 = phi double [ undef, %11 ], [ %45, %44 ]
%51 = phi double [ 0.000000e+00, %11 ], [ %45, %44 ]
%52 = phi i64 [ 0, %11 ], [ %46, %44 ]
br i1 %10, label %65, label %53
53: ; preds = %49
%54 = icmp eq i64 %52, %13
%55 = load double*, double** %6
br i1 %54, label %60, label %56
56: ; preds = %53
%57 = getelementptr inbounds double, double* %55, i64 %52
%58 = load double, double* %57
%59 = tail call double @llvm.fmuladd.f64(double %58, double %58, double %51)
br label %65
60: ; preds = %53
%61 = getelementptr inbounds double, double* %55, i64 %13
%62 = load double, double* %61
%63 = fadd double %62, -1.000000e+00
%64 = tail call double @llvm.fmuladd.f64(double %63, double %63, double %51)
br label %65
65: ; preds = %56, %60, %49
%66 = phi double [ %50, %49 ], [ %59, %56 ], [ %64, %60 ]
%67 = tail call double @sqrt(double %66)
%68 = icmp eq i64 %13, 0
%69 = fcmp olt double %67, %12
%70 = select i1 %68, i1 true, i1 %69
%71 = select i1 %70, double %67, double %12
%72 = add nuw i64 %13, 1
%73 = icmp eq i64 %72, %1
br i1 %73, label %74, label %11
74: ; preds = %65, %3
%75 = phi double [ 0.000000e+00, %3 ], [ %71, %65 ]
ret double %75
}