#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 /*<<< orphan*/ memcpy(float *, float *, long);
int ffr4fr4(float *input, long ntodo, double scale, double zero, float *output,
            int *status) {
  long ii;
  if (scale == 1. && zero == 0.) {
    memcpy(output, input, ntodo * sizeof(float));
  } else {
    for (ii = 0; ii < ntodo; ii++)
      output[ii] = (input[ii] - zero) / scale;
  }
  return (*status);
}


x86 assembly (clang, -O0):
<ffr4fr4>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xec1(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1181 <ffr4fr4+0x71>
jp     1181 <ffr4fr4+0x71>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1181 <ffr4fr4+0x71>
jp     1181 <ffr4fr4+0x71>
mov    -0x28(%rbp),%rdi
mov    -0x8(%rbp),%rsi
mov    -0x10(%rbp),%rdx
shl    $0x2,%rdx
call   1030 <memcpy@plt>
jmp    11d9 <ffr4fr4+0xc9>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11d4 <ffr4fr4+0xc4>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  -0x20(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1189 <ffr4fr4+0x79>
jmp    11d9 <ffr4fr4+0xc9>
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ffr4fr4>:
push   %r14
push   %rbx
push   %rax
mov    %rcx,%r14
mov    %rsi,%rax
mov    %rdi,%rsi
ucomisd 0xedb(%rip),%xmm0
jne    1149 <ffr4fr4+0x39>
jp     1149 <ffr4fr4+0x39>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    1149 <ffr4fr4+0x39>
jp     1149 <ffr4fr4+0x39>
shl    $0x2,%rax
mov    %rdx,%rdi
mov    %rax,%rdx
call   1030 <memcpy@plt>
jmp    12c9 <ffr4fr4+0x1b9>
test   %rax,%rax
jle    12c9 <ffr4fr4+0x1b9>
cmp    $0x4,%rax
jb     1172 <ffr4fr4+0x62>
lea    (%rsi,%rax,4),%rcx
cmp    %rdx,%rcx
jbe    11f4 <ffr4fr4+0xe4>
lea    (%rdx,%rax,4),%rcx
cmp    %rsi,%rcx
jbe    11f4 <ffr4fr4+0xe4>
xor    %ecx,%ecx
mov    %rcx,%rdi
not    %rdi
test   $0x1,%al
je     119c <ffr4fr4+0x8c>
movss  (%rsi,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
subsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%rdx,%rcx,4)
or     $0x1,%rcx
add    %rax,%rdi
je     12c9 <ffr4fr4+0x1b9>
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rsi,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
subsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%rdx,%rcx,4)
movss  0x4(%rsi,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
subsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,0x4(%rdx,%rcx,4)
add    $0x2,%rcx
cmp    %rcx,%rax
jne    11b0 <ffr4fr4+0xa0>
jmp    12c9 <ffr4fr4+0x1b9>
mov    %rax,%rcx
and    $0xfffffffffffffffc,%rcx
movapd %xmm1,%xmm3
unpcklpd %xmm1,%xmm3
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
lea    -0x4(%rcx),%rdi
mov    %rdi,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rdi,%rdi
je     12d4 <ffr4fr4+0x1c4>
mov    %r8,%rbx
and    $0xfffffffffffffffe,%rbx
xor    %edi,%edi
nopl   0x0(%rax)
cvtps2pd (%rsi,%rdi,4),%xmm4
cvtps2pd 0x8(%rsi,%rdi,4),%xmm5
subpd  %xmm3,%xmm5
subpd  %xmm3,%xmm4
divpd  %xmm2,%xmm4
divpd  %xmm2,%xmm5
cvtpd2ps %xmm5,%xmm5
cvtpd2ps %xmm4,%xmm4
unpcklpd %xmm5,%xmm4
movupd %xmm4,(%rdx,%rdi,4)
cvtps2pd 0x10(%rsi,%rdi,4),%xmm4
cvtps2pd 0x18(%rsi,%rdi,4),%xmm5
subpd  %xmm3,%xmm5
subpd  %xmm3,%xmm4
divpd  %xmm2,%xmm4
divpd  %xmm2,%xmm5
cvtpd2ps %xmm5,%xmm5
cvtpd2ps %xmm4,%xmm4
unpcklpd %xmm5,%xmm4
movupd %xmm4,0x10(%rdx,%rdi,4)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rbx
jne    1230 <ffr4fr4+0x120>
test   $0x1,%r8b
je     12c0 <ffr4fr4+0x1b0>
cvtps2pd (%rsi,%rdi,4),%xmm4
cvtps2pd 0x8(%rsi,%rdi,4),%xmm5
subpd  %xmm3,%xmm5
subpd  %xmm3,%xmm4
divpd  %xmm2,%xmm4
divpd  %xmm2,%xmm5
cvtpd2ps %xmm5,%xmm2
cvtpd2ps %xmm4,%xmm3
unpcklpd %xmm2,%xmm3
movupd %xmm3,(%rdx,%rdi,4)
cmp    %rax,%rcx
jne    1174 <ffr4fr4+0x64>
mov    (%r14),%eax
add    $0x8,%rsp
pop    %rbx
pop    %r14
ret
xor    %edi,%edi
test   $0x1,%r8b
jne    1296 <ffr4fr4+0x186>
jmp    12c0 <ffr4fr4+0x1b0>

LLVM IR (-O0):
define i32 @ffr4fr4(float* %0, i64 %1, double %2, double %3, float* %4, i32* %5) {
%7 = alloca float*
%8 = alloca i64
%9 = alloca double
%10 = alloca double
%11 = alloca float*
%12 = alloca i32*
%13 = alloca i64
store float* %0, float** %7
store i64 %1, i64* %8
store double %2, double* %9
store double %3, double* %10
store float* %4, float** %11
store i32* %5, i32** %12
%14 = load double, double* %9
%15 = fcmp oeq double %14, 1.000000e+00
br i1 %15, label %16, label %25
16: ; preds = %6
%17 = load double, double* %10
%18 = fcmp oeq double %17, 0.000000e+00
br i1 %18, label %19, label %25
19: ; preds = %16
%20 = load float*, float** %11
%21 = load float*, float** %7
%22 = load i64, i64* %8
%23 = mul i64 %22, 4
%24 = call i32 @memcpy(float* %20, float* %21, i64 %23)
br label %48
25: ; preds = %16, %6
store i64 0, i64* %13
br label %26
26: ; preds = %44, %25
%27 = load i64, i64* %13
%28 = load i64, i64* %8
%29 = icmp slt i64 %27, %28
br i1 %29, label %30, label %47
30: ; preds = %26
%31 = load float*, float** %7
%32 = load i64, i64* %13
%33 = getelementptr inbounds float, float* %31, i64 %32
%34 = load float, float* %33
%35 = fpext float %34 to double
%36 = load double, double* %10
%37 = fsub double %35, %36
%38 = load double, double* %9
%39 = fdiv double %37, %38
%40 = fptrunc double %39 to float
%41 = load float*, float** %11
%42 = load i64, i64* %13
%43 = getelementptr inbounds float, float* %41, i64 %42
store float %40, float* %43
br label %44
44: ; preds = %30
%45 = load i64, i64* %13
%46 = add nsw i64 %45, 1
store i64 %46, i64* %13
br label %26
47: ; preds = %26
br label %48
48: ; preds = %47, %19
%49 = load i32*, i32** %12
%50 = load i32, i32* %49
ret i32 %50
}

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