#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__ floorf(float);
void ambe_resample_mag(float *mag_dst, int L_dst, float *mag_src, int L_src) {
  float avg, step, pos;
  int i;
  avg = 0.0f;
  step = (float)L_src / (float)L_dst;
  pos = step;
  for (i = 0; i < L_dst; i++) {
    int posi = (int)floorf(pos);
    if (posi == 0) {
      mag_dst[i] = mag_src[0];
    } else if (posi >= L_src) {
      mag_dst[i] = mag_src[L_src - 1];
    } else {
      float alpha = pos - posi;
      mag_dst[i] = mag_src[posi - 1] * (1.0f - alpha) + mag_src[posi] * alpha;
    }
    avg += mag_dst[i];
    pos += step;
  }
  avg /= L_dst;
  for (i = 0; i < L_dst; i++)
    mag_dst[i] -= avg;
}


x86 assembly (clang, -O0):
<ambe_resample_mag>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
cvtsi2ssl -0x1c(%rbp),%xmm0
cvtsi2ssl -0xc(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1251 <ambe_resample_mag+0x141>
movss  -0x28(%rbp),%xmm0
call   1030 <floorf@plt>
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x30(%rbp)
jne    118f <ambe_resample_mag+0x7f>
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    121d <ambe_resample_mag+0x10d>
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     11bf <ambe_resample_mag+0xaf>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1218 <ambe_resample_mag+0x108>
movss  -0x28(%rbp),%xmm0
cvtsi2ssl -0x30(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  0xe14(%rip),%xmm2
subss  -0x34(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x34(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    121d <ambe_resample_mag+0x10d>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x24(%rbp),%xmm0
addss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1152 <ambe_resample_mag+0x42>
cvtsi2ssl -0xc(%rbp),%xmm1
movss  -0x20(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12a0 <ambe_resample_mag+0x190>
movss  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    126b <ambe_resample_mag+0x15b>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ambe_resample_mag>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
test   %esi,%esi
jle    12c7 <ambe_resample_mag+0x1b7>
mov    %ecx,%r15d
mov    %rdx,%rbp
mov    %rdi,%rbx
cvtsi2ss %ecx,%xmm2
cvtsi2ss %esi,%xmm0
movss  %xmm0,0x4(%rsp)
divss  %xmm0,%xmm2
lea    -0x1(%r15),%eax
movslq %eax,%r13
mov    %esi,0x8(%rsp)
mov    %esi,%r12d
xorps  %xmm5,%xmm5
xor    %r14d,%r14d
movss  %xmm2,0xc(%rsp)
jmp    11c6 <ambe_resample_mag+0xb6>
nopl   (%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  (%rsp),%xmm2
movaps %xmm2,%xmm1
subss  %xmm0,%xmm1
shl    $0x20,%rax
movabs $0xffffffff00000000,%rcx
ff ff ff
add    %rax,%rcx
movss  0xe74(%rip),%xmm0
subss  %xmm1,%xmm0
sar    $0x1e,%rax
mulss  0x0(%rbp,%rax,1),%xmm1
sar    $0x1e,%rcx
mulss  0x0(%rbp,%rcx,1),%xmm0
addss  %xmm1,%xmm0
movaps 0x10(%rsp),%xmm5
movss  %xmm0,(%rbx,%r14,4)
addss  %xmm0,%xmm5
addss  0xc(%rsp),%xmm2
add    $0x1,%r14
cmp    %r14,%r12
je     1201 <ambe_resample_mag+0xf1>
movaps %xmm5,0x10(%rsp)
movss  %xmm2,(%rsp)
movaps %xmm2,%xmm0
call   1030 <floorf@plt>
test   %eax,%eax
je     11f0 <ambe_resample_mag+0xe0>
cmp    %r15d,%eax
jl     1160 <ambe_resample_mag+0x50>
movss  0x0(%rbp,%r13,4),%xmm0
jmp    11f5 <ambe_resample_mag+0xe5>
xchg   %ax,%ax
movss  0x0(%rbp),%xmm0
movaps 0x10(%rsp),%xmm5
movss  (%rsp),%xmm2
jmp    11ad <ambe_resample_mag+0x9d>
mov    0x8(%rsp),%eax
test   %eax,%eax
jle    12c7 <ambe_resample_mag+0x1b7>
divss  0x4(%rsp),%xmm5
cmp    $0x8,%eax
jae    121f <ambe_resample_mag+0x10f>
xor    %eax,%eax
jmp    12b0 <ambe_resample_mag+0x1a0>
mov    %r12d,%eax
and    $0xfffffff8,%eax
movaps %xmm5,%xmm0
shufps $0x0,%xmm5,%xmm0
lea    -0x8(%rax),%rcx
mov    %rcx,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
test   %rcx,%rcx
je     12d6 <ambe_resample_mag+0x1c6>
mov    %rdx,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ecx,%ecx
nopl   (%rax)
movups (%rbx,%rcx,4),%xmm1
movups 0x10(%rbx,%rcx,4),%xmm2
movups 0x20(%rbx,%rcx,4),%xmm3
movups 0x30(%rbx,%rcx,4),%xmm4
subps  %xmm0,%xmm1
subps  %xmm0,%xmm2
movups %xmm1,(%rbx,%rcx,4)
movups %xmm2,0x10(%rbx,%rcx,4)
subps  %xmm0,%xmm3
subps  %xmm0,%xmm4
movups %xmm3,0x20(%rbx,%rcx,4)
movups %xmm4,0x30(%rbx,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1250 <ambe_resample_mag+0x140>
test   $0x1,%dl
je     12a9 <ambe_resample_mag+0x199>
movups (%rbx,%rcx,4),%xmm1
movups 0x10(%rbx,%rcx,4),%xmm2
subps  %xmm0,%xmm1
subps  %xmm0,%xmm2
movups %xmm1,(%rbx,%rcx,4)
movups %xmm2,0x10(%rbx,%rcx,4)
cmp    %r12,%rax
je     12c7 <ambe_resample_mag+0x1b7>
xchg   %ax,%ax
movss  (%rbx,%rax,4),%xmm0
subss  %xmm5,%xmm0
movss  %xmm0,(%rbx,%rax,4)
add    $0x1,%rax
cmp    %rax,%r12
jne    12b0 <ambe_resample_mag+0x1a0>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ecx,%ecx
test   $0x1,%dl
jne    1291 <ambe_resample_mag+0x181>
jmp    12a9 <ambe_resample_mag+0x199>

LLVM IR (-O0):
define void @ambe_resample_mag(float* %0, i32 %1, float* %2, i32 %3) {
%5 = alloca float*
%6 = alloca i32
%7 = alloca float*
%8 = alloca i32
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca i32
%13 = alloca i32
%14 = alloca float
store float* %0, float** %5
store i32 %1, i32* %6
store float* %2, float** %7
store i32 %3, i32* %8
store float 0.000000e+00, float* %9
%15 = load i32, i32* %8
%16 = sitofp i32 %15 to float
%17 = load i32, i32* %6
%18 = sitofp i32 %17 to float
%19 = fdiv float %16, %18
store float %19, float* %10
%20 = load float, float* %10
store float %20, float* %11
store i32 0, i32* %12
br label %21
21: ; preds = %91, %4
%22 = load i32, i32* %12
%23 = load i32, i32* %6
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %94
25: ; preds = %21
%26 = load float, float* %11
%27 = call i64 @floorf(float %26)
%28 = trunc i64 %27 to i32
store i32 %28, i32* %13
%29 = load i32, i32* %13
%30 = icmp eq i32 %29, 0
br i1 %30, label %31, label %39
31: ; preds = %25
%32 = load float*, float** %7
%33 = getelementptr inbounds float, float* %32, i64 0
%34 = load float, float* %33
%35 = load float*, float** %5
%36 = load i32, i32* %12
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds float, float* %35, i64 %37
store float %34, float* %38
br label %80
39: ; preds = %25
%40 = load i32, i32* %13
%41 = load i32, i32* %8
%42 = icmp sge i32 %40, %41
br i1 %42, label %43, label %54
43: ; preds = %39
%44 = load float*, float** %7
%45 = load i32, i32* %8
%46 = sub nsw i32 %45, 1
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds float, float* %44, i64 %47
%49 = load float, float* %48
%50 = load float*, float** %5
%51 = load i32, i32* %12
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds float, float* %50, i64 %52
store float %49, float* %53
br label %79
54: ; preds = %39
%55 = load float, float* %11
%56 = load i32, i32* %13
%57 = sitofp i32 %56 to float
%58 = fsub float %55, %57
store float %58, float* %14
%59 = load float*, float** %7
%60 = load i32, i32* %13
%61 = sub nsw i32 %60, 1
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %59, i64 %62
%64 = load float, float* %63
%65 = load float, float* %14
%66 = fsub float 1.000000e+00, %65
%67 = load float*, float** %7
%68 = load i32, i32* %13
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds float, float* %67, i64 %69
%71 = load float, float* %70
%72 = load float, float* %14
%73 = fmul float %71, %72
%74 = call float @llvm.fmuladd.f32(float %64, float %66, float %73)
%75 = load float*, float** %5
%76 = load i32, i32* %12
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float, float* %75, i64 %77
store float %74, float* %78
br label %79
79: ; preds = %54, %43
br label %80
80: ; preds = %79, %31
%81 = load float*, float** %5
%82 = load i32, i32* %12
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds float, float* %81, i64 %83
%85 = load float, float* %84
%86 = load float, float* %9
%87 = fadd float %86, %85
store float %87, float* %9
%88 = load float, float* %10
%89 = load float, float* %11
%90 = fadd float %89, %88
store float %90, float* %11
br label %91
91: ; preds = %80
%92 = load i32, i32* %12
%93 = add nsw i32 %92, 1
store i32 %93, i32* %12
br label %21
94: ; preds = %21
%95 = load i32, i32* %6
%96 = sitofp i32 %95 to float
%97 = load float, float* %9
%98 = fdiv float %97, %96
store float %98, float* %9
store i32 0, i32* %12
br label %99
99: ; preds = %111, %94
%100 = load i32, i32* %12
%101 = load i32, i32* %6
%102 = icmp slt i32 %100, %101
br i1 %102, label %103, label %114
103: ; preds = %99
%104 = load float, float* %9
%105 = load float*, float** %5
%106 = load i32, i32* %12
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds float, float* %105, i64 %107
%109 = load float, float* %108
%110 = fsub float %109, %104
store float %110, float* %108
br label %111
111: ; preds = %103
%112 = load i32, i32* %12
%113 = add nsw i32 %112, 1
store i32 %113, i32* %12
br label %99
114: ; preds = %99
ret void
}

LLVM IR (-O2):
define void @ambe_resample_mag(float* nocapture %0, i32 %1, float* nocapture %2, i32 %3) {
%5 = sitofp i32 %3 to float
%6 = sitofp i32 %1 to float
%7 = fdiv float %5, %6
%8 = icmp sgt i32 %1, 0
br i1 %8, label %9, label %118
9: ; preds = %4
%10 = add nsw i32 %3, -1
%11 = sext i32 %10 to i64
%12 = getelementptr inbounds float, float* %2, i64 %11
%13 = zext i32 %1 to i64
br label %14
14: ; preds = %9, %41
%15 = phi i64 [ 0, %9 ], [ %46, %41 ]
%16 = phi float [ 0.000000e+00, %9 ], [ %44, %41 ]
%17 = phi float [ %7, %9 ], [ %45, %41 ]
%18 = tail call i64 @floorf(float %17)
%19 = trunc i64 %18 to i32
%20 = icmp eq i32 %19, 0
br i1 %20, label %21, label %23
21: ; preds = %14
%22 = load float, float* %2
br label %41
23: ; preds = %14
%24 = icmp slt i32 %19, %3
br i1 %24, label %27, label %25
25: ; preds = %23
%26 = load float, float* %12
br label %41
27: ; preds = %23
%28 = sitofp i32 %19 to float
%29 = fsub float %17, %28
%30 = shl i64 %18, 32
%31 = add i64 %30, -4294967296
%32 = ashr exact i64 %31, 32
%33 = getelementptr inbounds float, float* %2, i64 %32
%34 = load float, float* %33
%35 = fsub float 1.000000e+00, %29
%36 = ashr exact i64 %30, 32
%37 = getelementptr inbounds float, float* %2, i64 %36
%38 = load float, float* %37
%39 = fmul float %29, %38
%40 = tail call float @llvm.fmuladd.f32(float %34, float %35, float %39)
br label %41
41: ; preds = %25, %27, %21
%42 = phi float [ %26, %25 ], [ %40, %27 ], [ %22, %21 ]
%43 = getelementptr inbounds float, float* %0, i64 %15
store float %42, float* %43
%44 = fadd float %16, %42
%45 = fadd float %7, %17
%46 = add nuw nsw i64 %15, 1
%47 = icmp eq i64 %46, %13
br i1 %47, label %48, label %14
48: ; preds = %41
%49 = fdiv float %44, %6
br i1 %8, label %50, label %118
50: ; preds = %48
%51 = zext i32 %1 to i64
%52 = icmp ult i32 %1, 8
br i1 %52, label %109, label %53
53: ; preds = %50
%54 = and i64 %13, 4294967288
%55 = insertelement <4 x float> poison, float %49, i64 0
%56 = shufflevector <4 x float> %55, <4 x float> poison, <4 x i32> zeroinitializer
%57 = insertelement <4 x float> poison, float %49, i64 0
%58 = shufflevector <4 x float> %57, <4 x float> poison, <4 x i32> zeroinitializer
%59 = add nsw i64 %54, -8
%60 = lshr exact i64 %59, 3
%61 = add nuw nsw i64 %60, 1
%62 = and i64 %61, 1
%63 = icmp eq i64 %59, 0
br i1 %63, label %93, label %64
64: ; preds = %53
%65 = and i64 %61, 4611686018427387902
br label %66
66: ; preds = %66, %64
%67 = phi i64 [ 0, %64 ], [ %90, %66 ]
%68 = phi i64 [ 0, %64 ], [ %91, %66 ]
%69 = getelementptr inbounds float, float* %0, i64 %67
%70 = bitcast float* %69 to <4 x float>*
%71 = load <4 x float>, <4 x float>* %70
%72 = getelementptr inbounds float, float* %69, i64 4
%73 = bitcast float* %72 to <4 x float>*
%74 = load <4 x float>, <4 x float>* %73
%75 = fsub <4 x float> %71, %56
%76 = fsub <4 x float> %74, %58
%77 = bitcast float* %69 to <4 x float>*
store <4 x float> %75, <4 x float>* %77
%78 = bitcast float* %72 to <4 x float>*
store <4 x float> %76, <4 x float>* %78
%79 = or i64 %67, 8
%80 = getelementptr inbounds float, float* %0, i64 %79
%81 = bitcast float* %80 to <4 x float>*
%82 = load <4 x float>, <4 x float>* %81
%83 = getelementptr inbounds float, float* %80, i64 4
%84 = bitcast float* %83 to <4 x float>*
%85 = load <4 x float>, <4 x float>* %84
%86 = fsub <4 x float> %82, %56
%87 = fsub <4 x float> %85, %58
%88 = bitcast float* %80 to <4 x float>*
store <4 x float> %86, <4 x float>* %88
%89 = bitcast float* %83 to <4 x float>*
store <4 x float> %87, <4 x float>* %89
%90 = add nuw i64 %67, 16
%91 = add i64 %68, 2
%92 = icmp eq i64 %91, %65
br i1 %92, label %93, label %66
93: ; preds = %66, %53
%94 = phi i64 [ 0, %53 ], [ %90, %66 ]
%95 = icmp eq i64 %62, 0
br i1 %95, label %107, label %96
96: ; preds = %93
%97 = getelementptr inbounds float, float* %0, i64 %94
%98 = bitcast float* %97 to <4 x float>*
%99 = load <4 x float>, <4 x float>* %98
%100 = getelementptr inbounds float, float* %97, i64 4
%101 = bitcast float* %100 to <4 x float>*
%102 = load <4 x float>, <4 x float>* %101
%103 = fsub <4 x float> %99, %56
%104 = fsub <4 x float> %102, %58
%105 = bitcast float* %97 to <4 x float>*
store <4 x float> %103, <4 x float>* %105
%106 = bitcast float* %100 to <4 x float>*
store <4 x float> %104, <4 x float>* %106
br label %107
107: ; preds = %93, %96
%108 = icmp eq i64 %54, %13
br i1 %108, label %118, label %109
109: ; preds = %50, %107
%110 = phi i64 [ 0, %50 ], [ %54, %107 ]
br label %111
111: ; preds = %109, %111
%112 = phi i64 [ %116, %111 ], [ %110, %109 ]
%113 = getelementptr inbounds float, float* %0, i64 %112
%114 = load float, float* %113
%115 = fsub float %114, %49
store float %115, float* %113
%116 = add nuw nsw i64 %112, 1
%117 = icmp eq i64 %116, %51
br i1 %117, label %118, label %111
118: ; preds = %111, %107, %4, %48
ret void
}