#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 sqrt(int);
void ComputeHsIext(int *Iext, int K, double *h0, double *h1) {
  int i, M;
  double iC, iL, iR;
  M = Iext[K + 1];
  for (i = 1; i < K; i++) {
    iL = Iext[i - 1];
    iC = Iext[i];
    iR = (double)Iext[i + 1];
    h0[i - 1] = ((M - iC) * iC / M) * (iR - iL) / ((iR - iC) * ((iC - iL)));
    h1[i - 1] = -sqrt((M - iC) * iC * (M - iR) * iR) / (M * (iR - iC));
  }
  iL = Iext[i - 1];
  iC = Iext[i];
  iR = (double)Iext[i + 1];
  h0[i - 1] = ((M - iC) * iC / M) * (iR - iL) / ((iR - iC) * ((iC - iL)));
}


x86 assembly (clang, -O0):
<ComputeHsIext>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x28(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    124d <ComputeHsIext+0x13d>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0x28(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
subsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x28(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
subsd  -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sqrt@plt>
mov    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm2
subsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1141 <ComputeHsIext+0x31>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0x28(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
subsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ComputeHsIext>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rcx,%r15
mov    %rdx,%r14
mov    %rdi,%rbx
movslq %esi,%rcx
cvtsi2sdl 0x4(%rdi,%rcx,4),%xmm4
mov    $0x1,%eax
cmp    $0x2,%ecx
jl     11f2 <ComputeHsIext+0xe2>
mov    %esi,%r12d
mov    %esi,%r13d
add    $0xffffffffffffffff,%r13
xor    %ebp,%ebp
movsd  %xmm4,0x8(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl (%rbx,%rbp,4),%xmm0
xorps  %xmm1,%xmm1
cvtsi2sdl 0x4(%rbx,%rbp,4),%xmm1
xorps  %xmm2,%xmm2
cvtsi2sdl 0x8(%rbx,%rbp,4),%xmm2
movapd %xmm4,%xmm3
subsd  %xmm1,%xmm3
mulsd  %xmm1,%xmm3
movsd  0x8(%rsp),%xmm4
subsd  %xmm2,%xmm4
mulsd  %xmm3,%xmm4
divsd  0x8(%rsp),%xmm3
movapd %xmm2,%xmm5
mulsd  %xmm2,%xmm4
subsd  %xmm0,%xmm2
mulsd  %xmm3,%xmm2
subsd  %xmm1,%xmm5
movsd  %xmm5,0x10(%rsp)
subsd  %xmm0,%xmm1
mulsd  %xmm5,%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,(%r14,%rbp,8)
cvttsd2si %xmm4,%edi
call   1030 <sqrt@plt>
movsd  0x8(%rsp),%xmm4
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  %xmm4,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r15,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%r13
jne    1150 <ComputeHsIext+0x40>
mov    %r12d,%eax
cltq
xorps  %xmm0,%xmm0
cvtsi2sdl -0x4(%rbx,%rax,4),%xmm0
mov    %eax,%ecx
xorps  %xmm1,%xmm1
cvtsi2sdl (%rbx,%rcx,4),%xmm1
xorps  %xmm2,%xmm2
cvtsi2sdl 0x4(%rbx,%rcx,4),%xmm2
movapd %xmm4,%xmm3
subsd  %xmm1,%xmm3
mulsd  %xmm1,%xmm3
divsd  %xmm4,%xmm3
movapd %xmm2,%xmm4
subsd  %xmm0,%xmm4
mulsd  %xmm3,%xmm4
subsd  %xmm1,%xmm2
subsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm4
movsd  %xmm4,-0x8(%r14,%rax,8)
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @ComputeHsIext(i32* %0, i32 %1, double* %2, double* %3) {
%5 = alloca i32*
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca double
%13 = alloca double
store i32* %0, i32** %5
store i32 %1, i32* %6
store double* %2, double** %7
store double* %3, double** %8
%14 = load i32*, i32** %5
%15 = load i32, i32* %6
%16 = add nsw i32 %15, 1
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds i32, i32* %14, i64 %17
%19 = load i32, i32* %18
store i32 %19, i32* %10
store i32 1, i32* %9
br label %20
20: ; preds = %100, %4
%21 = load i32, i32* %9
%22 = load i32, i32* %6
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %103
24: ; preds = %20
%25 = load i32*, i32** %5
%26 = load i32, i32* %9
%27 = sub nsw i32 %26, 1
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds i32, i32* %25, i64 %28
%30 = load i32, i32* %29
%31 = sitofp i32 %30 to double
store double %31, double* %12
%32 = load i32*, i32** %5
%33 = load i32, i32* %9
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds i32, i32* %32, i64 %34
%36 = load i32, i32* %35
%37 = sitofp i32 %36 to double
store double %37, double* %11
%38 = load i32*, i32** %5
%39 = load i32, i32* %9
%40 = add nsw i32 %39, 1
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds i32, i32* %38, i64 %41
%43 = load i32, i32* %42
%44 = sitofp i32 %43 to double
store double %44, double* %13
%45 = load i32, i32* %10
%46 = sitofp i32 %45 to double
%47 = load double, double* %11
%48 = fsub double %46, %47
%49 = load double, double* %11
%50 = fmul double %48, %49
%51 = load i32, i32* %10
%52 = sitofp i32 %51 to double
%53 = fdiv double %50, %52
%54 = load double, double* %13
%55 = load double, double* %12
%56 = fsub double %54, %55
%57 = fmul double %53, %56
%58 = load double, double* %13
%59 = load double, double* %11
%60 = fsub double %58, %59
%61 = load double, double* %11
%62 = load double, double* %12
%63 = fsub double %61, %62
%64 = fmul double %60, %63
%65 = fdiv double %57, %64
%66 = load double*, double** %7
%67 = load i32, i32* %9
%68 = sub nsw i32 %67, 1
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %66, i64 %69
store double %65, double* %70
%71 = load i32, i32* %10
%72 = sitofp i32 %71 to double
%73 = load double, double* %11
%74 = fsub double %72, %73
%75 = load double, double* %11
%76 = fmul double %74, %75
%77 = load i32, i32* %10
%78 = sitofp i32 %77 to double
%79 = load double, double* %13
%80 = fsub double %78, %79
%81 = fmul double %76, %80
%82 = load double, double* %13
%83 = fmul double %81, %82
%84 = fptosi double %83 to i32
%85 = call i32 @sqrt(i32 %84)
%86 = sub nsw i32 0, %85
%87 = sitofp i32 %86 to double
%88 = load i32, i32* %10
%89 = sitofp i32 %88 to double
%90 = load double, double* %13
%91 = load double, double* %11
%92 = fsub double %90, %91
%93 = fmul double %89, %92
%94 = fdiv double %87, %93
%95 = load double*, double** %8
%96 = load i32, i32* %9
%97 = sub nsw i32 %96, 1
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %95, i64 %98
store double %94, double* %99
br label %100
100: ; preds = %24
%101 = load i32, i32* %9
%102 = add nsw i32 %101, 1
store i32 %102, i32* %9
br label %20
103: ; preds = %20
%104 = load i32*, i32** %5
%105 = load i32, i32* %9
%106 = sub nsw i32 %105, 1
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds i32, i32* %104, i64 %107
%109 = load i32, i32* %108
%110 = sitofp i32 %109 to double
store double %110, double* %12
%111 = load i32*, i32** %5
%112 = load i32, i32* %9
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds i32, i32* %111, i64 %113
%115 = load i32, i32* %114
%116 = sitofp i32 %115 to double
store double %116, double* %11
%117 = load i32*, i32** %5
%118 = load i32, i32* %9
%119 = add nsw i32 %118, 1
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds i32, i32* %117, i64 %120
%122 = load i32, i32* %121
%123 = sitofp i32 %122 to double
store double %123, double* %13
%124 = load i32, i32* %10
%125 = sitofp i32 %124 to double
%126 = load double, double* %11
%127 = fsub double %125, %126
%128 = load double, double* %11
%129 = fmul double %127, %128
%130 = load i32, i32* %10
%131 = sitofp i32 %130 to double
%132 = fdiv double %129, %131
%133 = load double, double* %13
%134 = load double, double* %12
%135 = fsub double %133, %134
%136 = fmul double %132, %135
%137 = load double, double* %13
%138 = load double, double* %11
%139 = fsub double %137, %138
%140 = load double, double* %11
%141 = load double, double* %12
%142 = fsub double %140, %141
%143 = fmul double %139, %142
%144 = fdiv double %136, %143
%145 = load double*, double** %7
%146 = load i32, i32* %9
%147 = sub nsw i32 %146, 1
%148 = sext i32 %147 to i64
%149 = getelementptr inbounds double, double* %145, i64 %148
store double %144, double* %149
ret void
}

LLVM IR (-O2):
define void @ComputeHsIext(i32* nocapture %0, i32 %1, double* nocapture writeonly %2, double* nocapture writeonly %3) {
%5 = add nsw i32 %1, 1
%6 = sext i32 %5 to i64
%7 = getelementptr inbounds i32, i32* %0, i64 %6
%8 = load i32, i32* %7
%9 = icmp sgt i32 %1, 1
%10 = sitofp i32 %8 to double
br i1 %9, label %11, label %47
11: ; preds = %4
%12 = zext i32 %1 to i64
br label %13
13: ; preds = %11, %13
%14 = phi i64 [ 1, %11 ], [ %22, %13 ]
%15 = add nsw i64 %14, -1
%16 = getelementptr inbounds i32, i32* %0, i64 %15
%17 = load i32, i32* %16
%18 = sitofp i32 %17 to double
%19 = getelementptr inbounds i32, i32* %0, i64 %14
%20 = load i32, i32* %19
%21 = sitofp i32 %20 to double
%22 = add nuw nsw i64 %14, 1
%23 = getelementptr inbounds i32, i32* %0, i64 %22
%24 = load i32, i32* %23
%25 = sitofp i32 %24 to double
%26 = fsub double %10, %21
%27 = fmul double %26, %21
%28 = fdiv double %27, %10
%29 = fsub double %25, %18
%30 = fmul double %29, %28
%31 = fsub double %25, %21
%32 = fsub double %21, %18
%33 = fmul double %32, %31
%34 = fdiv double %30, %33
%35 = getelementptr inbounds double, double* %2, i64 %15
store double %34, double* %35
%36 = fsub double %10, %25
%37 = fmul double %27, %36
%38 = fmul double %37, %25
%39 = fptosi double %38 to i32
%40 = tail call i32 @sqrt(i32 %39)
%41 = sub nsw i32 0, %40
%42 = sitofp i32 %41 to double
%43 = fmul double %31, %10
%44 = fdiv double %42, %43
%45 = getelementptr inbounds double, double* %3, i64 %15
store double %44, double* %45
%46 = icmp eq i64 %22, %12
br i1 %46, label %47, label %13
47: ; preds = %13, %4
%48 = phi i32 [ 1, %4 ], [ %1, %13 ]
%49 = add nsw i32 %48, -1
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds i32, i32* %0, i64 %50
%52 = load i32, i32* %51
%53 = sitofp i32 %52 to double
%54 = zext i32 %48 to i64
%55 = getelementptr inbounds i32, i32* %0, i64 %54
%56 = load i32, i32* %55
%57 = sitofp i32 %56 to double
%58 = add nuw nsw i32 %48, 1
%59 = zext i32 %58 to i64
%60 = getelementptr inbounds i32, i32* %0, i64 %59
%61 = load i32, i32* %60
%62 = sitofp i32 %61 to double
%63 = fsub double %10, %57
%64 = fmul double %63, %57
%65 = fdiv double %64, %10
%66 = fsub double %62, %53
%67 = fmul double %66, %65
%68 = fsub double %62, %57
%69 = fsub double %57, %53
%70 = fmul double %69, %68
%71 = fdiv double %67, %70
%72 = getelementptr inbounds double, double* %2, i64 %50
store double %71, double* %72
ret void
}