#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 interpxy(float, float *, float *, int);
int /*<<< orphan*/ willus_mem_alloc(double **, int, char *);
int /*<<< orphan*/ willus_mem_free(double **, char *);
int arrayf_linearize_xyz(float *x, float *y, float *z, int n, int newn) {
  float *x0, *y0, *z0;
  int i;
  float xx;
  if (!willus_mem_alloc((double **)&x0, sizeof(float) * 3 * newn,
                        "arrayf_linearize_xyz"))
    return (1);
  y0 = &x0[newn];
  z0 = &y0[newn];
  for (i = 0; i < newn; i++) {
    xx = x[0] + (x[n - 1] - x[0]) * (double)i / (newn - 1);
    x0[i] = xx;
    y0[i] = interpxy(xx, x, y, n);
    z0[i] = interpxy(xx, x, z, n);
  }
  for (i = 0; i < newn; i++) {
    x[i] = x0[i];
    y[i] = y0[i];
    z[i] = z0[i];
  }
  willus_mem_free((double **)&x0, "arrayf_linearize_xyz");
  return (0);
}


x86 assembly (clang, -O0):
<arrayf_linearize_xyz>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8d,-0x28(%rbp)
lea    -0x30(%rbp),%rdi
movslq -0x28(%rbp),%rax
imul   $0xc,%rax,%rax
mov    %eax,%esi
lea    0xea0(%rip),%rdx
call   1030 <willus_mem_alloc@plt>
cmp    $0x0,%eax
jne    117a <arrayf_linearize_xyz+0x4a>
movl   $0x1,-0x4(%rbp)
jmp    12eb <arrayf_linearize_xyz+0x1bb>
mov    -0x30(%rbp),%rax
movslq -0x28(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    1265 <arrayf_linearize_xyz+0x135>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
subss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm1
cvtsi2sdl -0x44(%rbp),%xmm2
mulsd  %xmm2,%xmm1
mov    -0x28(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x24(%rbp),%edx
call   1040 <interpxy@plt>
mov    -0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x24(%rbp),%edx
call   1040 <interpxy@plt>
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11a7 <arrayf_linearize_xyz+0x77>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    12d4 <arrayf_linearize_xyz+0x1a4>
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    126c <arrayf_linearize_xyz+0x13c>
lea    -0x30(%rbp),%rdi
lea    0xd21(%rip),%rsi
call   1050 <willus_mem_free@plt>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<arrayf_linearize_xyz>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r8d,%r14d
mov    %ecx,%r13d
mov    %rdx,%r15
mov    %rsi,%rbp
mov    %rdi,%rbx
lea    0x0(,%r14,4),%eax
lea    (%rax,%rax,2),%esi
lea    0xea1(%rip),%rdx
lea    0x8(%rsp),%rdi
call   1030 <willus_mem_alloc@plt>
test   %eax,%eax
je     1288 <arrayf_linearize_xyz+0x158>
test   %r14d,%r14d
jle    133e <arrayf_linearize_xyz+0x20e>
mov    %rbp,0x28(%rsp)
movslq %r14d,%rdx
mov    %r14d,%ecx
mov    0x8(%rsp),%rsi
lea    (%rsi,%rdx,4),%rax
mov    %rax,0x18(%rsp)
mov    %rdx,0x38(%rsp)
lea    (%rax,%rdx,4),%rax
mov    %rax,0x20(%rsp)
lea    -0x1(%r13),%eax
movslq %eax,%rbp
mov    %r14,0x48(%rsp)
lea    -0x1(%r14),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x50(%rsp)
mov    %rcx,0x40(%rsp)
lea    -0x1(%rcx),%r12
xor    %r14d,%r14d
mov    %rsi,0x30(%rsp)
mov    %rsi,%rax
nop
movss  (%rbx),%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
movss  (%rbx,%rbp,4),%xmm2
subss  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
xorps  %xmm2,%xmm2
cvtsi2sd %r14d,%xmm2
mulsd  %xmm0,%xmm2
divsd  0x50(%rsp),%xmm2
addsd  %xmm1,%xmm2
xorps  %xmm0,%xmm0
cvtsd2ss %xmm2,%xmm0
movss  %xmm0,0x14(%rsp)
movss  %xmm0,(%rax,%r14,4)
mov    %rbx,%rdi
mov    0x28(%rsp),%rsi
mov    %r13d,%edx
call   1040 <interpxy@plt>
mov    0x18(%rsp),%rax
movss  %xmm0,(%rax,%r14,4)
movss  0x14(%rsp),%xmm0
mov    %rbx,%rdi
mov    %r15,%rsi
mov    %r13d,%edx
call   1040 <interpxy@plt>
mov    0x20(%rsp),%rax
movss  %xmm0,(%rax,%r14,4)
cmp    %r14,%r12
je     1261 <arrayf_linearize_xyz+0x131>
mov    0x8(%rsp),%rax
add    $0x1,%r14
jmp    11d0 <arrayf_linearize_xyz+0xa0>
cmpl   $0x0,0x48(%rsp)
mov    0x28(%rsp),%rbp
jle    133e <arrayf_linearize_xyz+0x20e>
mov    0x8(%rsp),%rax
mov    0x40(%rsp),%r8
cmp    $0x1,%r8d
jne    1292 <arrayf_linearize_xyz+0x162>
xor    %ecx,%ecx
jmp    130e <arrayf_linearize_xyz+0x1de>
mov    $0x1,%eax
jmp    1351 <arrayf_linearize_xyz+0x221>
mov    %r8d,%edx
and    $0xfffffffe,%edx
mov    0x38(%rsp),%rcx
mov    0x30(%rsp),%rdi
lea    (%rdi,%rcx,4),%rsi
add    $0x4,%rsi
lea    (%rdi,%rcx,8),%rdi
add    $0x4,%rdi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rbx,%rcx,4)
movss  -0x4(%rsi,%rcx,4),%xmm0
movss  %xmm0,0x0(%rbp,%rcx,4)
movss  -0x4(%rdi,%rcx,4),%xmm0
movss  %xmm0,(%r15,%rcx,4)
movss  0x4(%rax,%rcx,4),%xmm0
movss  %xmm0,0x4(%rbx,%rcx,4)
movss  (%rsi,%rcx,4),%xmm0
movss  %xmm0,0x4(%rbp,%rcx,4)
movss  (%rdi,%rcx,4),%xmm0
movss  %xmm0,0x4(%r15,%rcx,4)
add    $0x2,%rcx
cmp    %rcx,%rdx
jne    12c0 <arrayf_linearize_xyz+0x190>
test   $0x1,%r8b
je     133e <arrayf_linearize_xyz+0x20e>
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rbx,%rcx,4)
mov    0x18(%rsp),%rax
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,0x0(%rbp,%rcx,4)
mov    0x20(%rsp),%rax
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%r15,%rcx,4)
lea    0xcbb(%rip),%rsi
lea    0x8(%rsp),%rdi
call   1050 <willus_mem_free@plt>
xor    %eax,%eax
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [21 x i8] c"arrayf_linearize_xyz\00"
define i32 @arrayf_linearize_xyz(float* %0, float* %1, float* %2, i32 %3, i32 %4) {
%6 = alloca i32
%7 = alloca float*
%8 = alloca float*
%9 = alloca float*
%10 = alloca i32
%11 = alloca i32
%12 = alloca float*
%13 = alloca float*
%14 = alloca float*
%15 = alloca i32
%16 = alloca float
store float* %0, float** %7
store float* %1, float** %8
store float* %2, float** %9
store i32 %3, i32* %10
store i32 %4, i32* %11
%17 = bitcast float** %12 to double**
%18 = load i32, i32* %11
%19 = sext i32 %18 to i64
%20 = mul i64 12, %19
%21 = trunc i64 %20 to i32
%22 = call i32 @willus_mem_alloc(double** %17, i32 %21, i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str, i64 0, i64 0))
%23 = icmp ne i32 %22, 0
br i1 %23, label %25, label %24
24: ; preds = %5
store i32 1, i32* %6
br label %128
25: ; preds = %5
%26 = load float*, float** %12
%27 = load i32, i32* %11
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %26, i64 %28
store float* %29, float** %13
%30 = load float*, float** %13
%31 = load i32, i32* %11
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %30, i64 %32
store float* %33, float** %14
store i32 0, i32* %15
br label %34
34: ; preds = %86, %25
%35 = load i32, i32* %15
%36 = load i32, i32* %11
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %89
38: ; preds = %34
%39 = load float*, float** %7
%40 = getelementptr inbounds float, float* %39, i64 0
%41 = load float, float* %40
%42 = fpext float %41 to double
%43 = load float*, float** %7
%44 = load i32, i32* %10
%45 = sub nsw i32 %44, 1
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds float, float* %43, i64 %46
%48 = load float, float* %47
%49 = load float*, float** %7
%50 = getelementptr inbounds float, float* %49, i64 0
%51 = load float, float* %50
%52 = fsub float %48, %51
%53 = fpext float %52 to double
%54 = load i32, i32* %15
%55 = sitofp i32 %54 to double
%56 = fmul double %53, %55
%57 = load i32, i32* %11
%58 = sub nsw i32 %57, 1
%59 = sitofp i32 %58 to double
%60 = fdiv double %56, %59
%61 = fadd double %42, %60
%62 = fptrunc double %61 to float
store float %62, float* %16
%63 = load float, float* %16
%64 = load float*, float** %12
%65 = load i32, i32* %15
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds float, float* %64, i64 %66
store float %63, float* %67
%68 = load float, float* %16
%69 = load float*, float** %7
%70 = load float*, float** %8
%71 = load i32, i32* %10
%72 = call float @interpxy(float %68, float* %69, float* %70, i32 %71)
%73 = load float*, float** %13
%74 = load i32, i32* %15
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds float, float* %73, i64 %75
store float %72, float* %76
%77 = load float, float* %16
%78 = load float*, float** %7
%79 = load float*, float** %9
%80 = load i32, i32* %10
%81 = call float @interpxy(float %77, float* %78, float* %79, i32 %80)
%82 = load float*, float** %14
%83 = load i32, i32* %15
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
store float %81, float* %85
br label %86
86: ; preds = %38
%87 = load i32, i32* %15
%88 = add nsw i32 %87, 1
store i32 %88, i32* %15
br label %34
89: ; preds = %34
store i32 0, i32* %15
br label %90
90: ; preds = %122, %89
%91 = load i32, i32* %15
%92 = load i32, i32* %11
%93 = icmp slt i32 %91, %92
br i1 %93, label %94, label %125
94: ; preds = %90
%95 = load float*, float** %12
%96 = load i32, i32* %15
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds float, float* %95, i64 %97
%99 = load float, float* %98
%100 = load float*, float** %7
%101 = load i32, i32* %15
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds float, float* %100, i64 %102
store float %99, float* %103
%104 = load float*, float** %13
%105 = load i32, i32* %15
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds float, float* %104, i64 %106
%108 = load float, float* %107
%109 = load float*, float** %8
%110 = load i32, i32* %15
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds float, float* %109, i64 %111
store float %108, float* %112
%113 = load float*, float** %14
%114 = load i32, i32* %15
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds float, float* %113, i64 %115
%117 = load float, float* %116
%118 = load float*, float** %9
%119 = load i32, i32* %15
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %118, i64 %120
store float %117, float* %121
br label %122
122: ; preds = %94
%123 = load i32, i32* %15
%124 = add nsw i32 %123, 1
store i32 %124, i32* %15
br label %90
125: ; preds = %90
%126 = bitcast float** %12 to double**
%127 = call i32 @willus_mem_free(double** %126, i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str, i64 0, i64 0))
store i32 0, i32* %6
br label %128
128: ; preds = %125, %24
%129 = load i32, i32* %6
ret i32 %129
}

LLVM IR (-O2):
@.str = [21 x i8] c"arrayf_linearize_xyz\00"
define i32 @arrayf_linearize_xyz(float* %0, float* %1, float* %2, i32 %3, i32 %4) {
%6 = alloca float*
%7 = bitcast float** %6 to i8*
%8 = bitcast float** %6 to double**
%9 = sext i32 %4 to i64
%10 = mul i32 %4, 12
%11 = call i32 @willus_mem_alloc(double** %8, i32 %10, i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str, i64 0, i64 0))
%12 = icmp eq i32 %11, 0
br i1 %12, label %95, label %13
13: ; preds = %5
%14 = load float*, float** %6
%15 = getelementptr inbounds float, float* %14, i64 %9
%16 = getelementptr inbounds float, float* %15, i64 %9
%17 = icmp sgt i32 %4, 0
br i1 %17, label %18, label %93
18: ; preds = %13
%19 = add nsw i32 %3, -1
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds float, float* %0, i64 %20
%22 = add nsw i32 %4, -1
%23 = sitofp i32 %22 to double
%24 = zext i32 %4 to i64
br label %32
25: ; preds = %32
br i1 %17, label %26, label %93
26: ; preds = %25
%27 = load float*, float** %6
%28 = and i64 %24, 1
%29 = icmp eq i32 %4, 1
br i1 %29, label %80, label %30
30: ; preds = %26
%31 = and i64 %24, 4294967294
br label %55
32: ; preds = %53, %18
%33 = phi float* [ %14, %18 ], [ %54, %53 ]
%34 = phi i64 [ 0, %18 ], [ %51, %53 ]
%35 = load float, float* %0
%36 = fpext float %35 to double
%37 = load float, float* %21
%38 = fsub float %37, %35
%39 = fpext float %38 to double
%40 = trunc i64 %34 to i32
%41 = sitofp i32 %40 to double
%42 = fmul double %41, %39
%43 = fdiv double %42, %23
%44 = fadd double %43, %36
%45 = fptrunc double %44 to float
%46 = getelementptr inbounds float, float* %33, i64 %34
store float %45, float* %46
%47 = call float @interpxy(float %45, float* %0, float* %1, i32 %3)
%48 = getelementptr inbounds float, float* %15, i64 %34
store float %47, float* %48
%49 = call float @interpxy(float %45, float* %0, float* %2, i32 %3)
%50 = getelementptr inbounds float, float* %16, i64 %34
store float %49, float* %50
%51 = add nuw nsw i64 %34, 1
%52 = icmp eq i64 %51, %24
br i1 %52, label %25, label %53
53: ; preds = %32
%54 = load float*, float** %6
br label %32
55: ; preds = %55, %30
%56 = phi i64 [ 0, %30 ], [ %77, %55 ]
%57 = phi i64 [ 0, %30 ], [ %78, %55 ]
%58 = getelementptr inbounds float, float* %27, i64 %56
%59 = load float, float* %58
%60 = getelementptr inbounds float, float* %0, i64 %56
store float %59, float* %60
%61 = getelementptr inbounds float, float* %15, i64 %56
%62 = load float, float* %61
%63 = getelementptr inbounds float, float* %1, i64 %56
store float %62, float* %63
%64 = getelementptr inbounds float, float* %16, i64 %56
%65 = load float, float* %64
%66 = getelementptr inbounds float, float* %2, i64 %56
store float %65, float* %66
%67 = or i64 %56, 1
%68 = getelementptr inbounds float, float* %27, i64 %67
%69 = load float, float* %68
%70 = getelementptr inbounds float, float* %0, i64 %67
store float %69, float* %70
%71 = getelementptr inbounds float, float* %15, i64 %67
%72 = load float, float* %71
%73 = getelementptr inbounds float, float* %1, i64 %67
store float %72, float* %73
%74 = getelementptr inbounds float, float* %16, i64 %67
%75 = load float, float* %74
%76 = getelementptr inbounds float, float* %2, i64 %67
store float %75, float* %76
%77 = add nuw nsw i64 %56, 2
%78 = add i64 %57, 2
%79 = icmp eq i64 %78, %31
br i1 %79, label %80, label %55
80: ; preds = %55, %26
%81 = phi i64 [ 0, %26 ], [ %77, %55 ]
%82 = icmp eq i64 %28, 0
br i1 %82, label %93, label %83
83: ; preds = %80
%84 = getelementptr inbounds float, float* %27, i64 %81
%85 = load float, float* %84
%86 = getelementptr inbounds float, float* %0, i64 %81
store float %85, float* %86
%87 = getelementptr inbounds float, float* %15, i64 %81
%88 = load float, float* %87
%89 = getelementptr inbounds float, float* %1, i64 %81
store float %88, float* %89
%90 = getelementptr inbounds float, float* %16, i64 %81
%91 = load float, float* %90
%92 = getelementptr inbounds float, float* %2, i64 %81
store float %91, float* %92
br label %93
93: ; preds = %83, %80, %13, %25
%94 = call i32 @willus_mem_free(double** %8, i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str, i64 0, i64 0))
br label %95
95: ; preds = %5, %93
%96 = phi i32 [ 0, %93 ], [ 1, %5 ]
ret i32 %96
}