float revcst(float x[], float y[], int iorder[], int ncity, int n[]) {
  float xx[5], yy[5], de;
  int j, ii;
  n[3] = 1 + ((n[1] + ncity - 2) % ncity);
  n[4] = 1 + (n[2] % ncity);
  for (j = 1; j <= 4; j++) {
    ii = iorder[n[j]];
    xx[j] = x[ii];
    yy[j] = y[ii];
  }
  de = -sqrt(((xx[3]) - (xx[1])) * ((xx[3]) - (xx[1])) +
             ((yy[3]) - (yy[1])) * ((yy[3]) - (yy[1])));
  de -= sqrt(((xx[4]) - (xx[2])) * ((xx[4]) - (xx[2])) +
             ((yy[4]) - (yy[2])) * ((yy[4]) - (yy[2])));
  de += sqrt(((xx[4]) - (xx[1])) * ((xx[4]) - (xx[1])) +
             ((yy[4]) - (yy[1])) * ((yy[4]) - (yy[1])));
  de += sqrt(((xx[3]) - (xx[2])) * ((xx[3]) - (xx[2])) +
             ((yy[3]) - (yy[2])) * ((yy[3]) - (yy[2])));
  return de;
}


x86 assembly (clang, -O0):
<revcst>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x28(%rbp),%rax
mov    0x4(%rax),%eax
add    -0x1c(%rbp),%eax
sub    $0x2,%eax
cltd
idivl  -0x1c(%rbp)
add    $0x1,%edx
mov    -0x28(%rbp),%rax
mov    %edx,0xc(%rax)
mov    -0x28(%rbp),%rax
mov    0x8(%rax),%eax
cltd
idivl  -0x1c(%rbp)
add    $0x1,%edx
mov    -0x28(%rbp),%rax
mov    %edx,0x10(%rax)
movl   $0x1,-0x68(%rbp)
cmpl   $0x4,-0x68(%rbp)
jg     11be <revcst+0xae>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x68(%rbp),%rax
movss  %xmm0,-0x40(%rbp,%rax,4)
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x68(%rbp),%rax
movss  %xmm0,-0x60(%rbp,%rax,4)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    1162 <revcst+0x52>
movss  -0x34(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  -0x34(%rbp),%xmm2
subss  -0x3c(%rbp),%xmm2
movss  -0x54(%rbp),%xmm1
subss  -0x5c(%rbp),%xmm1
movss  -0x54(%rbp),%xmm3
subss  -0x5c(%rbp),%xmm3
mulss  %xmm3,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  -0x30(%rbp),%xmm2
subss  -0x38(%rbp),%xmm2
movss  -0x50(%rbp),%xmm1
subss  -0x58(%rbp),%xmm1
movss  -0x50(%rbp),%xmm3
subss  -0x58(%rbp),%xmm3
mulss  %xmm3,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  -0x30(%rbp),%xmm2
subss  -0x3c(%rbp),%xmm2
movss  -0x50(%rbp),%xmm1
subss  -0x5c(%rbp),%xmm1
movss  -0x50(%rbp),%xmm3
subss  -0x5c(%rbp),%xmm3
mulss  %xmm3,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x34(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  -0x34(%rbp),%xmm2
subss  -0x38(%rbp),%xmm2
movss  -0x54(%rbp),%xmm1
subss  -0x58(%rbp),%xmm1
movss  -0x54(%rbp),%xmm3
subss  -0x58(%rbp),%xmm3
mulss  %xmm3,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x64(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<revcst>:
push   %rbx
sub    $0x30,%rsp
mov    %rdx,%r9
movslq 0x4(%r8),%rbx
lea    (%rcx,%rbx,1),%eax
add    $0xfffffffe,%eax
cltd
idiv   %ecx
mov    %edx,%r10d
lea    0x1(%r10),%eax
mov    %eax,0xc(%r8)
movslq 0x8(%r8),%r11
mov    %r11d,%eax
cltd
idiv   %ecx
lea    0x1(%rdx),%eax
mov    %eax,0x10(%r8)
movslq (%r9,%rbx,4),%rax
movss  (%rdi,%rax,4),%xmm3
movss  (%rsi,%rax,4),%xmm2
movslq (%r9,%r11,4),%rax
movss  (%rdi,%rax,4),%xmm0
movss  %xmm0,0x10(%rsp)
movss  (%rsi,%rax,4),%xmm0
movss  %xmm0,(%rsp)
movslq %r10d,%rax
movslq 0x4(%r9,%rax,4),%rax
movss  (%rdi,%rax,4),%xmm0
movss  (%rsi,%rax,4),%xmm1
movslq %edx,%rax
movslq 0x4(%r9,%rax,4),%rax
movss  (%rdi,%rax,4),%xmm4
movss  (%rsi,%rax,4),%xmm5
movss  %xmm0,0x1c(%rsp)
movss  %xmm3,0x14(%rsp)
subss  %xmm3,%xmm0
movss  %xmm1,0x18(%rsp)
movss  %xmm2,0x4(%rsp)
subss  %xmm2,%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
movss  %xmm4,0xc(%rsp)
movss  %xmm5,0x8(%rsp)
jb     11e5 <revcst+0xc5>
xorps  %xmm2,%xmm2
sqrtss %xmm0,%xmm2
jmp    11f9 <revcst+0xd9>
call   1030 <sqrtf@plt>
movss  0x8(%rsp),%xmm5
movss  0xc(%rsp),%xmm4
movaps %xmm0,%xmm2
xorps  0xe00(%rip),%xmm2
movaps %xmm4,%xmm0
subss  0x10(%rsp),%xmm0
movaps %xmm5,%xmm1
subss  (%rsp),%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
jb     1230 <revcst+0x110>
sqrtsd %xmm0,%xmm0
jmp    124e <revcst+0x12e>
movaps %xmm2,0x20(%rsp)
call   1040 <sqrt@plt>
xorps  %xmm3,%xmm3
movaps 0x20(%rsp),%xmm2
movss  0x8(%rsp),%xmm5
movss  0xc(%rsp),%xmm4
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
subsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm1
subss  0x14(%rsp),%xmm4
subss  0x4(%rsp),%xmm5
mulss  %xmm5,%xmm5
mulss  %xmm4,%xmm4
addss  %xmm5,%xmm4
xorps  %xmm0,%xmm0
cvtss2sd %xmm4,%xmm0
ucomisd %xmm3,%xmm0
jb     1288 <revcst+0x168>
sqrtsd %xmm0,%xmm0
jmp    1299 <revcst+0x179>
movss  %xmm1,0x4(%rsp)
call   1040 <sqrt@plt>
movss  0x4(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtsd2ss %xmm1,%xmm2
movss  0x1c(%rsp),%xmm0
subss  0x10(%rsp),%xmm0
movss  0x18(%rsp),%xmm1
subss  (%rsp),%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     12de <revcst+0x1be>
sqrtsd %xmm0,%xmm0
jmp    12ed <revcst+0x1cd>
movss  %xmm2,(%rsp)
call   1040 <sqrt@plt>
movss  (%rsp),%xmm2
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
add    $0x30,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define float @revcst(float* %0, float* %1, i32* %2, i32 %3, i32* %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca i32*
%9 = alloca i32
%10 = alloca i32*
%11 = alloca [5 x float]
%12 = alloca [5 x float]
%13 = alloca float
%14 = alloca i32
%15 = alloca i32
store float* %0, float** %6
store float* %1, float** %7
store i32* %2, i32** %8
store i32 %3, i32* %9
store i32* %4, i32** %10
%16 = load i32*, i32** %10
%17 = getelementptr inbounds i32, i32* %16, i64 1
%18 = load i32, i32* %17
%19 = load i32, i32* %9
%20 = add nsw i32 %18, %19
%21 = sub nsw i32 %20, 2
%22 = load i32, i32* %9
%23 = srem i32 %21, %22
%24 = add nsw i32 1, %23
%25 = load i32*, i32** %10
%26 = getelementptr inbounds i32, i32* %25, i64 3
store i32 %24, i32* %26
%27 = load i32*, i32** %10
%28 = getelementptr inbounds i32, i32* %27, i64 2
%29 = load i32, i32* %28
%30 = load i32, i32* %9
%31 = srem i32 %29, %30
%32 = add nsw i32 1, %31
%33 = load i32*, i32** %10
%34 = getelementptr inbounds i32, i32* %33, i64 4
store i32 %32, i32* %34
store i32 1, i32* %14
br label %35
35: ; preds = %64, %5
%36 = load i32, i32* %14
%37 = icmp sle i32 %36, 4
br i1 %37, label %38, label %67
38: ; preds = %35
%39 = load i32*, i32** %8
%40 = load i32*, i32** %10
%41 = load i32, i32* %14
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds i32, i32* %40, i64 %42
%44 = load i32, i32* %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds i32, i32* %39, i64 %45
%47 = load i32, i32* %46
store i32 %47, i32* %15
%48 = load float*, float** %6
%49 = load i32, i32* %15
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds float, float* %48, i64 %50
%52 = load float, float* %51
%53 = load i32, i32* %14
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 %54
store float %52, float* %55
%56 = load float*, float** %7
%57 = load i32, i32* %15
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %56, i64 %58
%60 = load float, float* %59
%61 = load i32, i32* %14
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 %62
store float %60, float* %63
br label %64
64: ; preds = %38
%65 = load i32, i32* %14
%66 = add nsw i32 %65, 1
store i32 %66, i32* %14
br label %35
67: ; preds = %35
%68 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 3
%69 = load float, float* %68
%70 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 1
%71 = load float, float* %70
%72 = fsub float %69, %71
%73 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 3
%74 = load float, float* %73
%75 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 1
%76 = load float, float* %75
%77 = fsub float %74, %76
%78 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 3
%79 = load float, float* %78
%80 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 1
%81 = load float, float* %80
%82 = fsub float %79, %81
%83 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 3
%84 = load float, float* %83
%85 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 1
%86 = load float, float* %85
%87 = fsub float %84, %86
%88 = fmul float %82, %87
%89 = call float @llvm.fmuladd.f32(float %72, float %77, float %88)
%90 = fpext float %89 to double
%91 = call double @sqrt(double %90)
%92 = fneg double %91
%93 = fptrunc double %92 to float
store float %93, float* %13
%94 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 4
%95 = load float, float* %94
%96 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 2
%97 = load float, float* %96
%98 = fsub float %95, %97
%99 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 4
%100 = load float, float* %99
%101 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 2
%102 = load float, float* %101
%103 = fsub float %100, %102
%104 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 4
%105 = load float, float* %104
%106 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 2
%107 = load float, float* %106
%108 = fsub float %105, %107
%109 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 4
%110 = load float, float* %109
%111 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 2
%112 = load float, float* %111
%113 = fsub float %110, %112
%114 = fmul float %108, %113
%115 = call float @llvm.fmuladd.f32(float %98, float %103, float %114)
%116 = fpext float %115 to double
%117 = call double @sqrt(double %116)
%118 = load float, float* %13
%119 = fpext float %118 to double
%120 = fsub double %119, %117
%121 = fptrunc double %120 to float
store float %121, float* %13
%122 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 4
%123 = load float, float* %122
%124 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 1
%125 = load float, float* %124
%126 = fsub float %123, %125
%127 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 4
%128 = load float, float* %127
%129 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 1
%130 = load float, float* %129
%131 = fsub float %128, %130
%132 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 4
%133 = load float, float* %132
%134 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 1
%135 = load float, float* %134
%136 = fsub float %133, %135
%137 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 4
%138 = load float, float* %137
%139 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 1
%140 = load float, float* %139
%141 = fsub float %138, %140
%142 = fmul float %136, %141
%143 = call float @llvm.fmuladd.f32(float %126, float %131, float %142)
%144 = fpext float %143 to double
%145 = call double @sqrt(double %144)
%146 = load float, float* %13
%147 = fpext float %146 to double
%148 = fadd double %147, %145
%149 = fptrunc double %148 to float
store float %149, float* %13
%150 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 3
%151 = load float, float* %150
%152 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 2
%153 = load float, float* %152
%154 = fsub float %151, %153
%155 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 3
%156 = load float, float* %155
%157 = getelementptr inbounds [5 x float], [5 x float]* %11, i64 0, i64 2
%158 = load float, float* %157
%159 = fsub float %156, %158
%160 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 3
%161 = load float, float* %160
%162 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 2
%163 = load float, float* %162
%164 = fsub float %161, %163
%165 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 3
%166 = load float, float* %165
%167 = getelementptr inbounds [5 x float], [5 x float]* %12, i64 0, i64 2
%168 = load float, float* %167
%169 = fsub float %166, %168
%170 = fmul float %164, %169
%171 = call float @llvm.fmuladd.f32(float %154, float %159, float %170)
%172 = fpext float %171 to double
%173 = call double @sqrt(double %172)
%174 = load float, float* %13
%175 = fpext float %174 to double
%176 = fadd double %175, %173
%177 = fptrunc double %176 to float
store float %177, float* %13
%178 = load float, float* %13
ret float %178
}

LLVM IR (-O2):
define float @revcst(float* nocapture %0, float* nocapture %1, i32* nocapture %2, i32 %3, i32* nocapture %4) {
%6 = getelementptr inbounds i32, i32* %4, i64 1
%7 = load i32, i32* %6
%8 = add i32 %3, -2
%9 = add i32 %8, %7
%10 = srem i32 %9, %3
%11 = add nsw i32 %10, 1
%12 = getelementptr inbounds i32, i32* %4, i64 3
store i32 %11, i32* %12
%13 = getelementptr inbounds i32, i32* %4, i64 2
%14 = load i32, i32* %13
%15 = srem i32 %14, %3
%16 = add nsw i32 %15, 1
%17 = getelementptr inbounds i32, i32* %4, i64 4
store i32 %16, i32* %17
%18 = sext i32 %7 to i64
%19 = getelementptr inbounds i32, i32* %2, i64 %18
%20 = load i32, i32* %19
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds float, float* %0, i64 %21
%23 = load float, float* %22
%24 = getelementptr inbounds float, float* %1, i64 %21
%25 = load float, float* %24
%26 = sext i32 %14 to i64
%27 = getelementptr inbounds i32, i32* %2, i64 %26
%28 = load i32, i32* %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %0, i64 %29
%31 = load float, float* %30
%32 = getelementptr inbounds float, float* %1, i64 %29
%33 = load float, float* %32
%34 = sext i32 %11 to i64
%35 = getelementptr inbounds i32, i32* %2, i64 %34
%36 = load i32, i32* %35
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds float, float* %0, i64 %37
%39 = load float, float* %38
%40 = getelementptr inbounds float, float* %1, i64 %37
%41 = load float, float* %40
%42 = sext i32 %16 to i64
%43 = getelementptr inbounds i32, i32* %2, i64 %42
%44 = load i32, i32* %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds float, float* %0, i64 %45
%47 = load float, float* %46
%48 = getelementptr inbounds float, float* %1, i64 %45
%49 = load float, float* %48
%50 = fsub float %39, %23
%51 = fsub float %41, %25
%52 = fmul float %51, %51
%53 = tail call float @llvm.fmuladd.f32(float %50, float %50, float %52)
%54 = tail call float @sqrtf(float %53)
%55 = fneg float %54
%56 = fsub float %47, %31
%57 = fsub float %49, %33
%58 = fmul float %57, %57
%59 = tail call float @llvm.fmuladd.f32(float %56, float %56, float %58)
%60 = fpext float %59 to double
%61 = tail call double @sqrt(double %60)
%62 = fpext float %55 to double
%63 = fsub double %62, %61
%64 = fptrunc double %63 to float
%65 = fsub float %47, %23
%66 = fsub float %49, %25
%67 = fmul float %66, %66
%68 = tail call float @llvm.fmuladd.f32(float %65, float %65, float %67)
%69 = fpext float %68 to double
%70 = tail call double @sqrt(double %69)
%71 = fpext float %64 to double
%72 = fadd double %70, %71
%73 = fptrunc double %72 to float
%74 = fsub float %39, %31
%75 = fsub float %41, %33
%76 = fmul float %75, %75
%77 = tail call float @llvm.fmuladd.f32(float %74, float %74, float %76)
%78 = fpext float %77 to double
%79 = tail call double @sqrt(double %78)
%80 = fpext float %73 to double
%81 = fadd double %79, %80
%82 = fptrunc double %81 to float
ret float %82
}