#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 */
scalar_t__ pow(float, int);
float sqrt(scalar_t__);
void fill_euclid_dist(float (*prot1_coords)[3], float (*prot2_coords)[3],
                      float **dist, int *prot1_len, int *prot2_len) {
  int iter1, iter2;
  float x1, y1, z1, x2, y2, z2;
  for (iter1 = 0; iter1 < *prot1_len; iter1++) {
    x1 = prot1_coords[iter1][0];
    y1 = prot1_coords[iter1][1];
    z1 = prot1_coords[iter1][2];
    for (iter2 = 0; iter2 < *prot2_len; iter2++) {
      x2 = prot2_coords[iter2][0];
      y2 = prot2_coords[iter2][1];
      z2 = prot2_coords[iter2][2];
      dist[iter1][iter2] =
          sqrt((pow(x2 - x1, 2) + pow(y2 - y1, 2) + pow(z2 - z1, 2)));
    }
  }
}


x86 assembly (clang, -O0):
<fill_euclid_dist>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
jge    128c <fill_euclid_dist+0x16c>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movss  0x8(%rax),%xmm0
movss  %xmm0,-0x3c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x28(%rbp),%rcx
cmp    (%rcx),%eax
jge    1279 <fill_euclid_dist+0x159>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movss  0x8(%rax),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x40(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
mov    %rax,-0x58(%rbp)
movss  -0x44(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
mov    %rax,%rcx
mov    -0x58(%rbp),%rax
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
movss  -0x48(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
mov    -0x50(%rbp),%rdi
add    %rax,%rdi
call   1040 <sqrt@plt>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11a3 <fill_euclid_dist+0x83>
jmp    127e <fill_euclid_dist+0x15e>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1143 <fill_euclid_dist+0x23>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fill_euclid_dist>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rsi,0x20(%rsp)
mov    %rdi,0x28(%rsp)
mov    %rcx,0x30(%rsp)
mov    (%rcx),%ecx
test   %ecx,%ecx
jle    1247 <fill_euclid_dist+0x127>
mov    %r8,%r15
mov    %rdx,%rbx
mov    (%r8),%eax
addq   $0x8,0x20(%rsp)
xor    %r14d,%r14d
jmp    1170 <fill_euclid_dist+0x50>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r14
movslq %ecx,%rdx
cmp    %rdx,%r14
jge    1247 <fill_euclid_dist+0x127>
test   %eax,%eax
jle    1160 <fill_euclid_dist+0x40>
lea    (%r14,%r14,2),%rax
mov    0x28(%rsp),%rcx
movss  (%rcx,%rax,4),%xmm0
movss  %xmm0,0x14(%rsp)
movss  0x4(%rcx,%rax,4),%xmm0
movss  %xmm0,0x10(%rsp)
movss  0x8(%rcx,%rax,4),%xmm0
movss  %xmm0,0xc(%rsp)
mov    0x20(%rsp),%r12
xor    %r13d,%r13d
nopl   0x0(%rax,%rax,1)
movss  -0x8(%r12),%xmm0
subss  0x14(%rsp),%xmm0
movss  -0x4(%r12),%xmm1
subss  0x10(%rsp),%xmm1
movss  %xmm1,0x18(%rsp)
movss  (%r12),%xmm1
subss  0xc(%rsp),%xmm1
movss  %xmm1,0x1c(%rsp)
mov    $0x2,%edi
call   1030 <pow@plt>
mov    %rax,%rbp
movss  0x18(%rsp),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
add    %rax,%rbp
movss  0x1c(%rsp),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
add    %rax,%rbp
mov    %rbp,%rdi
call   1040 <sqrt@plt>
mov    (%rbx,%r14,8),%rax
movss  %xmm0,(%rax,%r13,4)
add    $0x1,%r13
movslq (%r15),%rax
add    $0xc,%r12
cmp    %rax,%r13
jl     11b0 <fill_euclid_dist+0x90>
mov    0x30(%rsp),%rcx
mov    (%rcx),%ecx
jmp    1160 <fill_euclid_dist+0x40>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @fill_euclid_dist([3 x float]* %0, [3 x float]* %1, float** %2, i32* %3, i32* %4) {
%6 = alloca [3 x float]*
%7 = alloca [3 x float]*
%8 = alloca float**
%9 = alloca i32*
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca float
store [3 x float]* %0, [3 x float]** %6
store [3 x float]* %1, [3 x float]** %7
store float** %2, float*** %8
store i32* %3, i32** %9
store i32* %4, i32** %10
store i32 0, i32* %11
br label %19
19: ; preds = %94, %5
%20 = load i32, i32* %11
%21 = load i32*, i32** %9
%22 = load i32, i32* %21
%23 = icmp slt i32 %20, %22
br i1 %23, label %24, label %97
24: ; preds = %19
%25 = load [3 x float]*, [3 x float]** %6
%26 = load i32, i32* %11
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds [3 x float], [3 x float]* %25, i64 %27
%29 = getelementptr inbounds [3 x float], [3 x float]* %28, i64 0, i64 0
%30 = load float, float* %29
store float %30, float* %13
%31 = load [3 x float]*, [3 x float]** %6
%32 = load i32, i32* %11
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds [3 x float], [3 x float]* %31, i64 %33
%35 = getelementptr inbounds [3 x float], [3 x float]* %34, i64 0, i64 1
%36 = load float, float* %35
store float %36, float* %14
%37 = load [3 x float]*, [3 x float]** %6
%38 = load i32, i32* %11
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds [3 x float], [3 x float]* %37, i64 %39
%41 = getelementptr inbounds [3 x float], [3 x float]* %40, i64 0, i64 2
%42 = load float, float* %41
store float %42, float* %15
store i32 0, i32* %12
br label %43
43: ; preds = %90, %24
%44 = load i32, i32* %12
%45 = load i32*, i32** %10
%46 = load i32, i32* %45
%47 = icmp slt i32 %44, %46
br i1 %47, label %48, label %93
48: ; preds = %43
%49 = load [3 x float]*, [3 x float]** %7
%50 = load i32, i32* %12
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds [3 x float], [3 x float]* %49, i64 %51
%53 = getelementptr inbounds [3 x float], [3 x float]* %52, i64 0, i64 0
%54 = load float, float* %53
store float %54, float* %16
%55 = load [3 x float]*, [3 x float]** %7
%56 = load i32, i32* %12
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds [3 x float], [3 x float]* %55, i64 %57
%59 = getelementptr inbounds [3 x float], [3 x float]* %58, i64 0, i64 1
%60 = load float, float* %59
store float %60, float* %17
%61 = load [3 x float]*, [3 x float]** %7
%62 = load i32, i32* %12
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds [3 x float], [3 x float]* %61, i64 %63
%65 = getelementptr inbounds [3 x float], [3 x float]* %64, i64 0, i64 2
%66 = load float, float* %65
store float %66, float* %18
%67 = load float, float* %16
%68 = load float, float* %13
%69 = fsub float %67, %68
%70 = call i64 @pow(float %69, i32 2)
%71 = load float, float* %17
%72 = load float, float* %14
%73 = fsub float %71, %72
%74 = call i64 @pow(float %73, i32 2)
%75 = add nsw i64 %70, %74
%76 = load float, float* %18
%77 = load float, float* %15
%78 = fsub float %76, %77
%79 = call i64 @pow(float %78, i32 2)
%80 = add nsw i64 %75, %79
%81 = call float @sqrt(i64 %80)
%82 = load float**, float*** %8
%83 = load i32, i32* %11
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float*, float** %82, i64 %84
%86 = load float*, float** %85
%87 = load i32, i32* %12
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds float, float* %86, i64 %88
store float %81, float* %89
br label %90
90: ; preds = %48
%91 = load i32, i32* %12
%92 = add nsw i32 %91, 1
store i32 %92, i32* %12
br label %43
93: ; preds = %43
br label %94
94: ; preds = %93
%95 = load i32, i32* %11
%96 = add nsw i32 %95, 1
store i32 %96, i32* %11
br label %19
97: ; preds = %19
ret void
}

LLVM IR (-O2):
define void @fill_euclid_dist([3 x float]* nocapture %0, [3 x float]* nocapture %1, float** nocapture %2, i32* nocapture %3, i32* nocapture %4) {
%6 = load i32, i32* %3
%7 = icmp sgt i32 %6, 0
br i1 %7, label %8, label %54
8: ; preds = %5
%9 = load i32, i32* %4
br label %10
10: ; preds = %8, %48
%11 = phi i32 [ %6, %8 ], [ %49, %48 ]
%12 = phi i32 [ %9, %8 ], [ %50, %48 ]
%13 = phi i64 [ 0, %8 ], [ %51, %48 ]
%14 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %13, i64 0
%15 = load float, float* %14
%16 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %13, i64 1
%17 = load float, float* %16
%18 = getelementptr inbounds [3 x float], [3 x float]* %0, i64 %13, i64 2
%19 = load float, float* %18
%20 = icmp sgt i32 %12, 0
br i1 %20, label %21, label %48
21: ; preds = %10
%22 = getelementptr inbounds float*, float** %2, i64 %13
br label %23
23: ; preds = %21, %23
%24 = phi i64 [ 0, %21 ], [ %42, %23 ]
%25 = getelementptr inbounds [3 x float], [3 x float]* %1, i64 %24, i64 0
%26 = load float, float* %25
%27 = getelementptr inbounds [3 x float], [3 x float]* %1, i64 %24, i64 1
%28 = load float, float* %27
%29 = getelementptr inbounds [3 x float], [3 x float]* %1, i64 %24, i64 2
%30 = load float, float* %29
%31 = fsub float %26, %15
%32 = tail call i64 @pow(float %31, i32 2)
%33 = fsub float %28, %17
%34 = tail call i64 @pow(float %33, i32 2)
%35 = add nsw i64 %34, %32
%36 = fsub float %30, %19
%37 = tail call i64 @pow(float %36, i32 2)
%38 = add nsw i64 %35, %37
%39 = tail call float @sqrt(i64 %38)
%40 = load float*, float** %22
%41 = getelementptr inbounds float, float* %40, i64 %24
store float %39, float* %41
%42 = add nuw nsw i64 %24, 1
%43 = load i32, i32* %4
%44 = sext i32 %43 to i64
%45 = icmp slt i64 %42, %44
br i1 %45, label %23, label %46
46: ; preds = %23
%47 = load i32, i32* %3
br label %48
48: ; preds = %46, %10
%49 = phi i32 [ %47, %46 ], [ %11, %10 ]
%50 = phi i32 [ %43, %46 ], [ %12, %10 ]
%51 = add nuw nsw i64 %13, 1
%52 = sext i32 %49 to i64
%53 = icmp slt i64 %51, %52
br i1 %53, label %10, label %54
54: ; preds = %48, %5
ret void
}