#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 fabs(double const);
double sec_root(double *, double *, double const *, int, int const);
void tdeig(double *lambda, double *y, double *d, const double *v, const int n) {
  int i;
  double v1norm = 0, min = v[0], max = v[0];
  for (i = 1; i <= n; ++i) {
    double vi = fabs(v[i]), a = d[i] - vi, b = d[i] + vi;
    v1norm += vi;
    if (a < min)
      min = a;
    if (b > max)
      max = b;
  }
  d[0] = v[0] - v1norm < min ? v[0] - v1norm : min;
  d[n + 1] = v[0] + v1norm > max ? v[0] + v1norm : max;
  for (i = 0; i <= n; ++i)
    lambda[i] = sec_root(&y[i], d, v, i, n);
}


x86 assembly (clang, -O0):
<tdeig>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     1206 <tdeig+0xf6>
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe8b(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x50(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11d9 <tdeig+0xc9>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x58(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    11f3 <tdeig+0xe3>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    11f8 <tdeig+0xe8>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1155 <tdeig+0x45>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1239 <tdeig+0x129>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    1243 <tdeig+0x133>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
addsd  -0x30(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    127f <tdeig+0x16f>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    1289 <tdeig+0x179>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     12f1 <tdeig+0x1e1>
mov    -0x10(%rbp),%rdi
movslq -0x28(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x28(%rbp),%ecx
mov    -0x24(%rbp),%r8d
call   1030 <sec_root@plt>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12a7 <tdeig+0x197>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tdeig>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r8d,%r15d
mov    %rcx,%r12
mov    %rdx,%r11
mov    %rsi,%rbx
mov    %rdi,%r13
movsd  (%rcx),%xmm8
lea    0x1(%r15),%eax
mov    %eax,%r14d
test   %r8d,%r8d
jle    115c <tdeig+0x4c>
lea    -0x1(%r14),%rcx
cmp    $0x2,%eax
jne    116f <tdeig+0x5f>
xorpd  %xmm1,%xmm1
mov    $0x1,%edx
movapd %xmm8,%xmm3
movapd %xmm8,%xmm2
jmp    120e <tdeig+0xfe>
xorpd  %xmm1,%xmm1
movapd %xmm8,%xmm2
movapd %xmm8,%xmm3
jmp    1247 <tdeig+0x137>
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
xorpd  %xmm5,%xmm5
mov    $0x1,%edx
movapd 0xe75(%rip),%xmm9
movapd %xmm8,%xmm3
movapd %xmm8,%xmm6
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r12,%rdx,8),%xmm7
movsd  0x8(%r12,%rdx,8),%xmm1
andpd  %xmm9,%xmm7
movsd  (%r11,%rdx,8),%xmm0
movsd  0x8(%r11,%rdx,8),%xmm2
movapd %xmm0,%xmm4
subsd  %xmm7,%xmm4
minsd  %xmm3,%xmm4
addsd  %xmm7,%xmm0
maxsd  %xmm6,%xmm0
addsd  %xmm5,%xmm7
andpd  %xmm9,%xmm1
movapd %xmm2,%xmm3
subsd  %xmm1,%xmm3
minsd  %xmm4,%xmm3
addsd  %xmm1,%xmm2
maxsd  %xmm0,%xmm2
addsd  %xmm7,%xmm1
lea    (%rsi,%rdx,1),%rdi
add    $0x2,%rdi
add    $0x2,%rdx
movapd %xmm1,%xmm5
movapd %xmm2,%xmm6
cmp    $0x1,%rdi
jne    11a0 <tdeig+0x90>
test   $0x1,%cl
je     1247 <tdeig+0x137>
movsd  (%r12,%rdx,8),%xmm0
andpd  0xddf(%rip),%xmm0
movsd  (%r11,%rdx,8),%xmm4
movapd %xmm4,%xmm5
subsd  %xmm0,%xmm5
minsd  %xmm3,%xmm5
addsd  %xmm0,%xmm4
maxsd  %xmm2,%xmm4
addsd  %xmm0,%xmm1
movapd %xmm4,%xmm2
movapd %xmm5,%xmm3
subsd  %xmm1,%xmm8
minsd  %xmm3,%xmm8
movsd  %xmm8,(%r11)
addsd  (%r12),%xmm1
maxsd  %xmm2,%xmm1
cltq
movsd  %xmm1,(%r11,%rax,8)
test   %r15d,%r15d
js     12a5 <tdeig+0x195>
xor    %ebp,%ebp
nop
mov    %rbx,(%rsp)
mov    %rbx,%rdi
mov    %r11,%rsi
mov    %r12,%rdx
mov    %ebp,%ecx
mov    %r15d,%r8d
mov    %r11,%rbx
call   1030 <sec_root@plt>
mov    %rbx,%r11
mov    (%rsp),%rbx
movsd  %xmm0,0x0(%r13,%rbp,8)
add    $0x1,%rbp
add    $0x8,%rbx
cmp    %rbp,%r14
jne    1270 <tdeig+0x160>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @tdeig(double* nocapture writeonly %0, double* %1, double* %2, double* %3, i32 %4) {
%6 = load double, double* %3
%7 = icmp slt i32 %4, 1
%8 = add i32 %4, 1
br i1 %7, label %72, label %9
9: ; preds = %5
%10 = zext i32 %8 to i64
%11 = add nsw i64 %10, -1
%12 = and i64 %11, 1
%13 = icmp eq i32 %8, 2
br i1 %13, label %50, label %14
14: ; preds = %9
%15 = and i64 %11, -2
br label %16
16: ; preds = %16, %14
%17 = phi i64 [ 1, %14 ], [ %47, %16 ]
%18 = phi double [ 0.000000e+00, %14 ], [ %42, %16 ]
%19 = phi double [ %6, %14 ], [ %44, %16 ]
%20 = phi double [ %6, %14 ], [ %46, %16 ]
%21 = phi i64 [ 0, %14 ], [ %48, %16 ]
%22 = getelementptr inbounds double, double* %3, i64 %17
%23 = load double, double* %22
%24 = tail call double @llvm.fabs.f64(double %23)
%25 = getelementptr inbounds double, double* %2, i64 %17
%26 = load double, double* %25
%27 = fsub double %26, %24
%28 = fadd double %26, %24
%29 = fadd double %18, %24
%30 = fcmp olt double %27, %19
%31 = select i1 %30, double %27, double %19
%32 = fcmp ogt double %28, %20
%33 = select i1 %32, double %28, double %20
%34 = add nuw nsw i64 %17, 1
%35 = getelementptr inbounds double, double* %3, i64 %34
%36 = load double, double* %35
%37 = tail call double @llvm.fabs.f64(double %36)
%38 = getelementptr inbounds double, double* %2, i64 %34
%39 = load double, double* %38
%40 = fsub double %39, %37
%41 = fadd double %39, %37
%42 = fadd double %29, %37
%43 = fcmp olt double %40, %31
%44 = select i1 %43, double %40, double %31
%45 = fcmp ogt double %41, %33
%46 = select i1 %45, double %41, double %33
%47 = add nuw nsw i64 %17, 2
%48 = add i64 %21, 2
%49 = icmp eq i64 %48, %15
br i1 %49, label %50, label %16
50: ; preds = %16, %9
%51 = phi double [ undef, %9 ], [ %42, %16 ]
%52 = phi double [ undef, %9 ], [ %44, %16 ]
%53 = phi double [ undef, %9 ], [ %46, %16 ]
%54 = phi i64 [ 1, %9 ], [ %47, %16 ]
%55 = phi double [ 0.000000e+00, %9 ], [ %42, %16 ]
%56 = phi double [ %6, %9 ], [ %44, %16 ]
%57 = phi double [ %6, %9 ], [ %46, %16 ]
%58 = icmp eq i64 %12, 0
br i1 %58, label %72, label %59
59: ; preds = %50
%60 = getelementptr inbounds double, double* %3, i64 %54
%61 = load double, double* %60
%62 = tail call double @llvm.fabs.f64(double %61)
%63 = getelementptr inbounds double, double* %2, i64 %54
%64 = load double, double* %63
%65 = fsub double %64, %62
%66 = fadd double %64, %62
%67 = fadd double %55, %62
%68 = fcmp olt double %65, %56
%69 = select i1 %68, double %65, double %56
%70 = fcmp ogt double %66, %57
%71 = select i1 %70, double %66, double %57
br label %72
72: ; preds = %59, %50, %5
%73 = phi double [ %6, %5 ], [ %53, %50 ], [ %71, %59 ]
%74 = phi double [ %6, %5 ], [ %52, %50 ], [ %69, %59 ]
%75 = phi double [ 0.000000e+00, %5 ], [ %51, %50 ], [ %67, %59 ]
%76 = fsub double %6, %75
%77 = fcmp olt double %76, %74
%78 = select i1 %77, double %76, double %74
store double %78, double* %2
%79 = load double, double* %3
%80 = fadd double %75, %79
%81 = fcmp ogt double %80, %73
%82 = select i1 %81, double %80, double %73
%83 = sext i32 %8 to i64
%84 = getelementptr inbounds double, double* %2, i64 %83
store double %82, double* %84
%85 = icmp slt i32 %4, 0
br i1 %85, label %96, label %86
86: ; preds = %72
%87 = zext i32 %8 to i64
br label %88
88: ; preds = %86, %88
%89 = phi i64 [ 0, %86 ], [ %94, %88 ]
%90 = getelementptr inbounds double, double* %1, i64 %89
%91 = trunc i64 %89 to i32
%92 = tail call double @sec_root(double* %90, double* %2, double* %3, i32 %91, i32 %4)
%93 = getelementptr inbounds double, double* %0, i64 %89
store double %92, double* %93
%94 = add nuw nsw i64 %89, 1
%95 = icmp eq i64 %94, %87
br i1 %95, label %96, label %88
96: ; preds = %88, %72
ret void
}