#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 */
float FLT_MAX;
float edgedf(float, float, float);
float sqrtf(float);
float distaa3(const float *img, const float *gximg, const float *gyimg, int w,
              int c, int xc, int yc, int xi, int yi) {
  float di, df, dx, dy, gx, gy, a;
  int closest;
  closest = c - xc - yc * w;
  a = img[closest];
  gx = gximg[closest];
  gy = gyimg[closest];
  if (a > 1)
    a = 1;
  if (a < 0)
    a = 0;
  if (a == 0)
    return FLT_MAX;
  dx = (float)xi;
  dy = (float)yi;
  di = sqrtf(dx * dx + dy * dy);
  if (di == 0) {
    df = edgedf(gx, gy, a);
  } else {
    df = edgedf(dx, dy, a);
  }
  return di + df;
}


x86 assembly (clang, -O0):
<distaa3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8d,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    0x10(%rbp),%ecx
imul   -0x24(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x48(%rbp),%xmm0
movss  0xe63(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    11b3 <distaa3+0x93>
movss  0xe52(%rip),%xmm0
movss  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
ucomiss -0x48(%rbp),%xmm0
jbe    11c8 <distaa3+0xa8>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    11f4 <distaa3+0xd4>
jp     11f4 <distaa3+0xd4>
mov    0x2e02(%rip),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    128b <distaa3+0x16b>
cvtsi2ssl 0x18(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
cvtsi2ssl 0x20(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x38(%rbp),%xmm0
movss  -0x38(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  -0x3c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
call   1040 <sqrtf@plt>
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1263 <distaa3+0x143>
jp     1263 <distaa3+0x143>
movss  -0x40(%rbp),%xmm0
movss  -0x44(%rbp),%xmm1
movss  -0x48(%rbp),%xmm2
call   1030 <edgedf@plt>
movss  %xmm0,-0x34(%rbp)
jmp    127c <distaa3+0x15c>
movss  -0x38(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm1
movss  -0x48(%rbp),%xmm2
call   1030 <edgedf@plt>
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm0
addss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<distaa3>:
sub    $0x18,%rsp
imul   0x20(%rsp),%ecx
add    %r9d,%ecx
sub    %ecx,%r8d
movslq %r8d,%rax
movss  0xec6(%rip),%xmm0
minss  (%rdi,%rax,4),%xmm0
xorps  %xmm5,%xmm5
xorps  %xmm2,%xmm2
maxss  %xmm0,%xmm2
ucomiss %xmm5,%xmm2
jne    1160 <distaa3+0x40>
jp     1160 <distaa3+0x40>
mov    0x2e91(%rip),%rax
movss  (%rax),%xmm1
jmp    11fe <distaa3+0xde>
movss  (%rsi,%rax,4),%xmm4
movss  (%rdx,%rax,4),%xmm1
cvtsi2ssl 0x28(%rsp),%xmm6
cvtsi2ssl 0x30(%rsp),%xmm7
movaps %xmm7,%xmm3
mulss  %xmm7,%xmm3
movaps %xmm6,%xmm0
mulss  %xmm6,%xmm0
addss  %xmm3,%xmm0
ucomiss %xmm5,%xmm0
jb     1193 <distaa3+0x73>
sqrtss %xmm0,%xmm0
jmp    11d4 <distaa3+0xb4>
movss  %xmm2,0x4(%rsp)
movss  %xmm4,0x14(%rsp)
movss  %xmm1,0x10(%rsp)
movss  %xmm6,0xc(%rsp)
movss  %xmm7,0x8(%rsp)
call   1040 <sqrtf@plt>
movss  0x8(%rsp),%xmm7
movss  0xc(%rsp),%xmm6
movss  0x10(%rsp),%xmm1
movss  0x14(%rsp),%xmm4
movss  0x4(%rsp),%xmm2
xorps  %xmm3,%xmm3
movss  %xmm0,0x4(%rsp)
ucomiss %xmm3,%xmm0
jne    11e9 <distaa3+0xc9>
jp     11e9 <distaa3+0xc9>
movaps %xmm4,%xmm0
jmp    11ef <distaa3+0xcf>
movaps %xmm6,%xmm0
movaps %xmm7,%xmm1
call   1030 <edgedf@plt>
movss  0x4(%rsp),%xmm1
addss  %xmm0,%xmm1
movaps %xmm1,%xmm0
add    $0x18,%rsp
ret

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

LLVM IR (-O2):
@FLT_MAX = global float 0.000000e+00
define float @distaa3(float* nocapture %0, float* nocapture %1, float* nocapture %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, i32 %8) {
%10 = mul nsw i32 %6, %3
%11 = add i32 %10, %5
%12 = sub i32 %4, %11
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds float, float* %0, i64 %13
%15 = load float, float* %14
%16 = getelementptr inbounds float, float* %1, i64 %13
%17 = load float, float* %16
%18 = getelementptr inbounds float, float* %2, i64 %13
%19 = load float, float* %18
%20 = fcmp ogt float %15, 1.000000e+00
%21 = select i1 %20, float 1.000000e+00, float %15
%22 = fcmp olt float %21, 0.000000e+00
%23 = select i1 %22, float 0.000000e+00, float %21
%24 = fcmp oeq float %23, 0.000000e+00
br i1 %24, label %25, label %27
25: ; preds = %9
%26 = load float, float* @FLT_MAX
br label %41
27: ; preds = %9
%28 = sitofp i32 %7 to float
%29 = sitofp i32 %8 to float
%30 = fmul float %29, %29
%31 = tail call float @llvm.fmuladd.f32(float %28, float %28, float %30)
%32 = tail call float @sqrtf(float %31)
%33 = fcmp oeq float %32, 0.000000e+00
br i1 %33, label %34, label %36
34: ; preds = %27
%35 = tail call float @edgedf(float %17, float %19, float %23)
br label %38
36: ; preds = %27
%37 = tail call float @edgedf(float %28, float %29, float %23)
br label %38
38: ; preds = %36, %34
%39 = phi float [ %35, %34 ], [ %37, %36 ]
%40 = fadd float %32, %39
br label %41
41: ; preds = %38, %25
%42 = phi float [ %26, %25 ], [ %40, %38 ]
ret float %42
}