#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 */
double R_pow_di(double, int);
double log(double);
int sqrt(double);
double cpoly_scale(int n, double *pot, double eps, double BIG, double small,
                   double base) {
  int i, ell;
  double x, high, sc, lo, min_, max_;
  high = sqrt(BIG);
  lo = small / eps;
  max_ = 0.;
  min_ = BIG;
  for (i = 0; i < n; i++) {
    x = pot[i];
    if (x > max_)
      max_ = x;
    if (x != 0. && x < min_)
      min_ = x;
  }
  if (min_ < lo || max_ > high) {
    x = lo / min_;
    if (x <= 1.)
      sc = 1. / (sqrt(max_) * sqrt(min_));
    else {
      sc = x;
      if (BIG / sc > max_)
        sc = 1.0;
    }
    ell = (int)(log(sc) / log(base) + 0.5);
    return R_pow_di(base, ell);
  } else
    return 1.0;
}


x86 assembly (clang, -O0):
<cpoly_scale>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movsd  %xmm2,-0x30(%rbp)
movsd  %xmm3,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1050 <sqrt@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1217 <cpoly_scale+0xe7>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
jbe    11c9 <cpoly_scale+0x99>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11e6 <cpoly_scale+0xb6>
jp     11e6 <cpoly_scale+0xb6>
jmp    1204 <cpoly_scale+0xd4>
movsd  -0x48(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1204 <cpoly_scale+0xd4>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    1209 <cpoly_scale+0xd9>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1191 <cpoly_scale+0x61>
movsd  -0x68(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm0
ucomisd %xmm1,%xmm0
ja     123b <cpoly_scale+0x10b>
movsd  -0x70(%rbp),%xmm0
ucomisd -0x50(%rbp),%xmm0
jbe    1316 <cpoly_scale+0x1e6>
movsd  -0x60(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xdae(%rip),%xmm0
ucomisd -0x48(%rbp),%xmm0
jb     1296 <cpoly_scale+0x166>
movsd  -0x70(%rbp),%xmm0
call   1050 <sqrt@plt>
mov    %eax,-0x74(%rbp)
movsd  -0x68(%rbp),%xmm0
call   1050 <sqrt@plt>
mov    %eax,%ecx
mov    -0x74(%rbp),%eax
imul   %ecx,%eax
cvtsi2sd %eax,%xmm1
movsd  0xd78(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    12c7 <cpoly_scale+0x197>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x28(%rbp),%xmm0
divsd  -0x58(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
jbe    12c2 <cpoly_scale+0x192>
movsd  0xd43(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    12c7 <cpoly_scale+0x197>
movsd  -0x58(%rbp),%xmm0
call   1040 <log@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1040 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  0xd14(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x40(%rbp),%edi
call   1030 <R_pow_di@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    1323 <cpoly_scale+0x1f3>
movsd  0xce2(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cpoly_scale>:
push   %rbp
push   %rbx
sub    $0x28,%rsp
movsd  %xmm3,0x18(%rsp)
movsd  %xmm2,0x20(%rsp)
movsd  %xmm0,(%rsp)
mov    %rsi,%rbx
mov    %edi,%ebp
movsd  %xmm1,0x10(%rsp)
movapd %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  0x20(%rsp),%xmm5
divsd  (%rsp),%xmm5
test   %ebp,%ebp
jle    1186 <cpoly_scale+0x56>
movsd  0x10(%rsp),%xmm6
mov    %ebp,%ecx
cmp    $0x1,%ebp
jne    1199 <cpoly_scale+0x69>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
movapd %xmm6,%xmm7
jmp    120a <cpoly_scale+0xda>
xorpd  %xmm0,%xmm0
movsd  0x10(%rsp),%xmm6
movapd %xmm6,%xmm7
jmp    123d <cpoly_scale+0x10d>
mov    %ecx,%esi
and    $0xfffffffe,%esi
xorpd  %xmm1,%xmm1
xor    %edx,%edx
xorpd  %xmm0,%xmm0
movapd %xmm6,%xmm2
nopl   0x0(%rax)
movsd  (%rbx,%rdx,8),%xmm3
movsd  0x8(%rbx,%rdx,8),%xmm7
movapd %xmm3,%xmm4
maxsd  %xmm0,%xmm4
movapd %xmm3,%xmm0
minsd  %xmm2,%xmm0
cmpneqsd %xmm1,%xmm3
andpd  %xmm3,%xmm0
andnpd %xmm2,%xmm3
orpd   %xmm0,%xmm3
movapd %xmm7,%xmm0
maxsd  %xmm4,%xmm0
movapd %xmm7,%xmm2
cmpneqsd %xmm1,%xmm2
minsd  %xmm3,%xmm7
andpd  %xmm2,%xmm7
andnpd %xmm3,%xmm2
orpd   %xmm2,%xmm7
add    $0x2,%rdx
movapd %xmm7,%xmm2
cmp    %rdx,%rsi
jne    11b0 <cpoly_scale+0x80>
test   $0x1,%cl
je     123d <cpoly_scale+0x10d>
movsd  (%rbx,%rdx,8),%xmm1
movapd %xmm1,%xmm2
xorpd  %xmm3,%xmm3
cmpneqsd %xmm1,%xmm3
maxsd  %xmm0,%xmm1
minsd  %xmm7,%xmm2
andpd  %xmm3,%xmm2
andnpd %xmm7,%xmm3
orpd   %xmm2,%xmm3
movapd %xmm3,%xmm7
movapd %xmm1,%xmm0
ucomisd %xmm7,%xmm5
ja     125f <cpoly_scale+0x12f>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
ucomisd %xmm1,%xmm0
ja     125f <cpoly_scale+0x12f>
movsd  0xda8(%rip),%xmm0
add    $0x28,%rsp
pop    %rbx
pop    %rbp
ret
divsd  %xmm7,%xmm5
movsd  0xd95(%rip),%xmm1
ucomisd %xmm5,%xmm1
jae    1281 <cpoly_scale+0x151>
divsd  %xmm5,%xmm6
ucomisd %xmm0,%xmm6
jbe    12ac <cpoly_scale+0x17c>
movapd %xmm1,%xmm5
jmp    12ac <cpoly_scale+0x17c>
movapd %xmm7,(%rsp)
call   1050 <sqrt@plt>
mov    %eax,%ebx
movaps (%rsp),%xmm0
call   1050 <sqrt@plt>
imul   %ebx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0xd58(%rip),%xmm5
divsd  %xmm0,%xmm5
movapd %xmm5,%xmm0
call   1040 <log@plt>
movsd  %xmm0,(%rsp)
movsd  0x18(%rsp),%xmm0
call   1040 <log@plt>
movsd  (%rsp),%xmm1
divsd  %xmm0,%xmm1
addsd  0xd32(%rip),%xmm1
cvttsd2si %xmm1,%edi
movsd  0x18(%rsp),%xmm0
add    $0x28,%rsp
pop    %rbx
pop    %rbp
jmp    1030 <R_pow_di@plt>

LLVM IR (-O0):
define double @cpoly_scale(i32 %0, double* %1, double %2, double %3, double %4, double %5) {
%7 = alloca double
%8 = alloca i32
%9 = alloca double*
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
store i32 %0, i32* %8
store double* %1, double** %9
store double %2, double* %10
store double %3, double* %11
store double %4, double* %12
store double %5, double* %13
%22 = load double, double* %11
%23 = call i32 @sqrt(double %22)
%24 = sitofp i32 %23 to double
store double %24, double* %17
%25 = load double, double* %12
%26 = load double, double* %10
%27 = fdiv double %25, %26
store double %27, double* %19
store double 0.000000e+00, double* %21
%28 = load double, double* %11
store double %28, double* %20
store i32 0, i32* %14
br label %29
29: ; preds = %54, %6
%30 = load i32, i32* %14
%31 = load i32, i32* %8
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %57
33: ; preds = %29
%34 = load double*, double** %9
%35 = load i32, i32* %14
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
store double %38, double* %16
%39 = load double, double* %16
%40 = load double, double* %21
%41 = fcmp ogt double %39, %40
br i1 %41, label %42, label %44
42: ; preds = %33
%43 = load double, double* %16
store double %43, double* %21
br label %44
44: ; preds = %42, %33
%45 = load double, double* %16
%46 = fcmp une double %45, 0.000000e+00
br i1 %46, label %47, label %53
47: ; preds = %44
%48 = load double, double* %16
%49 = load double, double* %20
%50 = fcmp olt double %48, %49
br i1 %50, label %51, label %53
51: ; preds = %47
%52 = load double, double* %16
store double %52, double* %20
br label %53
53: ; preds = %51, %47, %44
br label %54
54: ; preds = %53
%55 = load i32, i32* %14
%56 = add nsw i32 %55, 1
store i32 %56, i32* %14
br label %29
57: ; preds = %29
%58 = load double, double* %20
%59 = load double, double* %19
%60 = fcmp olt double %58, %59
br i1 %60, label %65, label %61
61: ; preds = %57
%62 = load double, double* %21
%63 = load double, double* %17
%64 = fcmp ogt double %62, %63
br i1 %64, label %65, label %99
65: ; preds = %61, %57
%66 = load double, double* %19
%67 = load double, double* %20
%68 = fdiv double %66, %67
store double %68, double* %16
%69 = load double, double* %16
%70 = fcmp ole double %69, 1.000000e+00
br i1 %70, label %71, label %79
71: ; preds = %65
%72 = load double, double* %21
%73 = call i32 @sqrt(double %72)
%74 = load double, double* %20
%75 = call i32 @sqrt(double %74)
%76 = mul nsw i32 %73, %75
%77 = sitofp i32 %76 to double
%78 = fdiv double 1.000000e+00, %77
store double %78, double* %18
br label %88
79: ; preds = %65
%80 = load double, double* %16
store double %80, double* %18
%81 = load double, double* %11
%82 = load double, double* %18
%83 = fdiv double %81, %82
%84 = load double, double* %21
%85 = fcmp ogt double %83, %84
br i1 %85, label %86, label %87
86: ; preds = %79
store double 1.000000e+00, double* %18
br label %87
87: ; preds = %86, %79
br label %88
88: ; preds = %87, %71
%89 = load double, double* %18
%90 = call double @log(double %89)
%91 = load double, double* %13
%92 = call double @log(double %91)
%93 = fdiv double %90, %92
%94 = fadd double %93, 5.000000e-01
%95 = fptosi double %94 to i32
store i32 %95, i32* %15
%96 = load double, double* %13
%97 = load i32, i32* %15
%98 = call double @R_pow_di(double %96, i32 %97)
store double %98, double* %7
br label %100
99: ; preds = %61
store double 1.000000e+00, double* %7
br label %100
100: ; preds = %99, %88
%101 = load double, double* %7
ret double %101
}

LLVM IR (-O2):
define double @cpoly_scale(i32 %0, double* nocapture %1, double %2, double %3, double %4, double %5) {
%7 = tail call i32 @sqrt(double %3)
%8 = sitofp i32 %7 to double
%9 = fdiv double %4, %2
%10 = icmp sgt i32 %0, 0
br i1 %10, label %11, label %58
11: ; preds = %6
%12 = zext i32 %0 to i64
%13 = and i64 %12, 1
%14 = icmp eq i32 %0, 1
br i1 %14, label %42, label %15
15: ; preds = %11
%16 = and i64 %12, 4294967294
br label %17
17: ; preds = %17, %15
%18 = phi i64 [ 0, %15 ], [ %39, %17 ]
%19 = phi double [ 0.000000e+00, %15 ], [ %34, %17 ]
%20 = phi double [ %3, %15 ], [ %38, %17 ]
%21 = phi i64 [ 0, %15 ], [ %40, %17 ]
%22 = getelementptr inbounds double, double* %1, i64 %18
%23 = load double, double* %22
%24 = fcmp ogt double %23, %19
%25 = select i1 %24, double %23, double %19
%26 = fcmp une double %23, 0.000000e+00
%27 = fcmp olt double %23, %20
%28 = select i1 %26, i1 %27, i1 false
%29 = select i1 %28, double %23, double %20
%30 = or i64 %18, 1
%31 = getelementptr inbounds double, double* %1, i64 %30
%32 = load double, double* %31
%33 = fcmp ogt double %32, %25
%34 = select i1 %33, double %32, double %25
%35 = fcmp une double %32, 0.000000e+00
%36 = fcmp olt double %32, %29
%37 = select i1 %35, i1 %36, i1 false
%38 = select i1 %37, double %32, double %29
%39 = add nuw nsw i64 %18, 2
%40 = add i64 %21, 2
%41 = icmp eq i64 %40, %16
br i1 %41, label %42, label %17
42: ; preds = %17, %11
%43 = phi double [ undef, %11 ], [ %34, %17 ]
%44 = phi double [ undef, %11 ], [ %38, %17 ]
%45 = phi i64 [ 0, %11 ], [ %39, %17 ]
%46 = phi double [ 0.000000e+00, %11 ], [ %34, %17 ]
%47 = phi double [ %3, %11 ], [ %38, %17 ]
%48 = icmp eq i64 %13, 0
br i1 %48, label %58, label %49
49: ; preds = %42
%50 = getelementptr inbounds double, double* %1, i64 %45
%51 = load double, double* %50
%52 = fcmp ogt double %51, %46
%53 = select i1 %52, double %51, double %46
%54 = fcmp une double %51, 0.000000e+00
%55 = fcmp olt double %51, %47
%56 = select i1 %54, i1 %55, i1 false
%57 = select i1 %56, double %51, double %47
br label %58
58: ; preds = %49, %42, %6
%59 = phi double [ %3, %6 ], [ %44, %42 ], [ %57, %49 ]
%60 = phi double [ 0.000000e+00, %6 ], [ %43, %42 ], [ %53, %49 ]
%61 = fcmp olt double %59, %9
%62 = fcmp ogt double %60, %8
%63 = select i1 %61, i1 true, i1 %62
br i1 %63, label %64, label %85
64: ; preds = %58
%65 = fdiv double %9, %59
%66 = fcmp ugt double %65, 1.000000e+00
br i1 %66, label %73, label %67
67: ; preds = %64
%68 = tail call i32 @sqrt(double %60)
%69 = tail call i32 @sqrt(double %59)
%70 = mul nsw i32 %69, %68
%71 = sitofp i32 %70 to double
%72 = fdiv double 1.000000e+00, %71
br label %77
73: ; preds = %64
%74 = fdiv double %3, %65
%75 = fcmp ogt double %74, %60
br i1 %75, label %76, label %77
76: ; preds = %73
br label %77
77: ; preds = %73, %76, %67
%78 = phi double [ %72, %67 ], [ 1.000000e+00, %76 ], [ %65, %73 ]
%79 = tail call double @log(double %78)
%80 = tail call double @log(double %5)
%81 = fdiv double %79, %80
%82 = fadd double %81, 5.000000e-01
%83 = fptosi double %82 to i32
%84 = tail call double @R_pow_di(double %5, i32 %83)
br label %85
85: ; preds = %58, %77
%86 = phi double [ %84, %77 ], [ 1.000000e+00, %58 ]
ret double %86
}