#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 cubic(double, double, double);
double sqrt(double);
int quarticroots(double a, double b, double c, double d, double *rts) {
  double y, p, q, r, A, B, C, m, n1, n2, asq, dis, rtdis;
  double A_4, B_m;
  int roots = 0;
  asq = a * a;
  A = b - asq * 3.0 / 8.0;
  B = c + a * (asq / 8.0 - b / 2.0);
  C = d + asq * (b / 16.0 - asq * 3.0 / 256.0) - a * c / 4.0;
  p = 2.0 * A;
  q = A * A - 4.0 * C;
  r = -B * B;
  y = cubic(p, q, r);
  if (y <= 0.0)
    return 0;
  m = sqrt(y);
  A_4 = 0.25 * A;
  B_m = B / m;
  n1 = 0.5 * (y + A + B_m);
  n2 = 0.5 * (y + A - B_m);
  dis = y - 4.0 * n1;
  if (dis >= 0.0) {
    rtdis = sqrt(dis);
    rts[roots++] = 0.5 * (-m - rtdis) - A_4;
    rts[roots++] = 0.5 * (-m + rtdis) - A_4;
  }
  dis = y - 4.0 * n2;
  if (dis >= 0.0) {
    rtdis = sqrt(dis);
    rts[roots++] = 0.5 * (-m - rtdis) - A_4;
    rts[roots++] = 0.5 * (-m + rtdis) - A_4;
  }
  return roots;
}


x86 assembly (clang, -O0):
<quarticroots>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  %xmm3,-0x28(%rbp)
mov    %rdi,-0x30(%rbp)
movl   $0x0,-0xac(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  0xeac(%rip),%xmm1
mulsd  -0x88(%rbp),%xmm1
ff
movsd  0xeac(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x20(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
movsd  0xe85(%rip),%xmm3
divsd  %xmm3,%xmm2
movsd  -0x18(%rbp),%xmm3
movsd  0xe54(%rip),%xmm4
divsd  %xmm4,%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x18(%rbp),%xmm2
movsd  0xe3d(%rip),%xmm3
divsd  %xmm3,%xmm2
movsd  0xe29(%rip),%xmm3
mulsd  -0x88(%rbp),%xmm3
ff
movsd  0xe11(%rip),%xmm4
divsd  %xmm4,%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x10(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
movsd  0xddf(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  0xdd2(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  0xdae(%rip),%xmm1
mulsd  -0x68(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x60(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm2
call   1030 <cubic@plt>
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x38(%rbp),%xmm0
jb     12d4 <quarticroots+0x1b4>
movl   $0x0,-0x4(%rbp)
jmp    15b1 <quarticroots+0x491>
movsd  -0x38(%rbp),%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x70(%rbp)
movsd  0xd55(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x38(%rbp),%xmm1
addsd  -0x58(%rbp),%xmm1
addsd  -0xa8(%rbp),%xmm1
ff
movsd  0xd14(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x38(%rbp),%xmm1
addsd  -0x58(%rbp),%xmm1
subsd  -0xa8(%rbp),%xmm1
ff
movsd  0xcf1(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x78(%rbp),%xmm0
movsd  0xcce(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     147c <quarticroots+0x35c>
movsd  -0x90(%rbp),%xmm0
ff
call   1040 <sqrt@plt>
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xc51(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mov    -0xac(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0xac(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xbe1(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mov    -0xac(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0xac(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm0
movsd  0xba2(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     15a8 <quarticroots+0x488>
movsd  -0x90(%rbp),%xmm0
ff
call   1040 <sqrt@plt>
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xb25(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mov    -0xac(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0xac(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xab5(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mov    -0xac(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0xac(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xac(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<quarticroots>:
push   %rbp
push   %rbx
sub    $0x38,%rsp
mov    %rdi,%rbx
movapd %xmm0,%xmm4
mulsd  %xmm0,%xmm4
movsd  0xec7(%rip),%xmm5
mulsd  %xmm4,%xmm5
movsd  0xec3(%rip),%xmm6
movapd %xmm5,%xmm7
mulsd  %xmm6,%xmm7
movapd %xmm1,%xmm8
subsd  %xmm7,%xmm8
mulsd  %xmm4,%xmm6
movsd  0xead(%rip),%xmm7
mulsd  %xmm1,%xmm7
addsd  %xmm6,%xmm7
mulsd  %xmm0,%xmm7
mulsd  0xea1(%rip),%xmm1
mulsd  0xea1(%rip),%xmm5
addsd  %xmm2,%xmm7
addsd  %xmm5,%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
mulsd  0xe8d(%rip),%xmm0
addsd  %xmm0,%xmm1
movapd %xmm8,%xmm0
addsd  %xmm8,%xmm0
mulsd  0xe7f(%rip),%xmm1
movsd  %xmm8,(%rsp)
movapd %xmm8,%xmm2
mulsd  %xmm8,%xmm2
addsd  %xmm2,%xmm1
movapd 0xe83(%rip),%xmm2
xorpd  %xmm7,%xmm2
movapd %xmm7,0x10(%rsp)
mulsd  %xmm7,%xmm2
call   1030 <cubic@plt>
movapd %xmm0,%xmm5
xor    %ebp,%ebp
xorpd  %xmm0,%xmm0
ucomisd %xmm5,%xmm0
jae    135b <quarticroots+0x23b>
ucomisd %xmm0,%xmm5
jb     1203 <quarticroots+0xe3>
xorps  %xmm4,%xmm4
sqrtsd %xmm5,%xmm4
jmp    121c <quarticroots+0xfc>
movapd %xmm5,%xmm0
movsd  %xmm5,0x8(%rsp)
call   1040 <sqrt@plt>
movsd  0x8(%rsp),%xmm5
movapd %xmm0,%xmm4
movsd  (%rsp),%xmm2
movsd  0xe0f(%rip),%xmm3
mulsd  %xmm2,%xmm3
movapd 0x10(%rsp),%xmm6
divsd  %xmm4,%xmm6
addsd  %xmm5,%xmm2
movapd %xmm2,%xmm0
addsd  %xmm6,%xmm0
movsd  0xdf5(%rip),%xmm1
mulsd  %xmm1,%xmm0
subsd  %xmm6,%xmm2
mulsd  %xmm1,%xmm2
mulsd  0xdd1(%rip),%xmm0
addsd  %xmm5,%xmm0
xor    %ebp,%ebp
xorpd  %xmm6,%xmm6
ucomisd %xmm6,%xmm0
jb     12e3 <quarticroots+0x1c3>
jb     1277 <quarticroots+0x157>
sqrtsd %xmm0,%xmm0
jmp    12ae <quarticroots+0x18e>
movsd  %xmm2,(%rsp)
movapd %xmm3,0x10(%rsp)
movapd %xmm4,0x20(%rsp)
movsd  %xmm5,0x8(%rsp)
call   1040 <sqrt@plt>
xorpd  %xmm6,%xmm6
movsd  0x8(%rsp),%xmm5
movapd 0x20(%rsp),%xmm4
movapd 0x10(%rsp),%xmm3
movsd  (%rsp),%xmm2
movapd 0xd9a(%rip),%xmm1
xorpd  %xmm4,%xmm1
unpcklpd %xmm0,%xmm1
unpcklpd %xmm4,%xmm0
subpd  %xmm0,%xmm1
movapd %xmm3,%xmm0
mulpd  0xd8e(%rip),%xmm1
unpcklpd %xmm3,%xmm0
subpd  %xmm0,%xmm1
movupd %xmm1,(%rbx)
mov    $0x2,%ebp
mulsd  0xd45(%rip),%xmm2
addsd  %xmm5,%xmm2
ucomisd %xmm6,%xmm2
jb     135b <quarticroots+0x23b>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm2
jb     1308 <quarticroots+0x1e8>
xorps  %xmm0,%xmm0
sqrtsd %xmm2,%xmm0
jmp    1329 <quarticroots+0x209>
movapd %xmm2,%xmm0
movapd %xmm3,0x10(%rsp)
movapd %xmm4,0x20(%rsp)
call   1040 <sqrt@plt>
movapd 0x20(%rsp),%xmm4
movapd 0x10(%rsp),%xmm3
movapd 0xd1f(%rip),%xmm1
xorpd  %xmm4,%xmm1
mov    %ebp,%eax
unpcklpd %xmm0,%xmm1
unpcklpd %xmm4,%xmm0
subpd  %xmm0,%xmm1
mulpd  0xd15(%rip),%xmm1
unpcklpd %xmm3,%xmm3
subpd  %xmm3,%xmm1
add    $0x2,%ebp
movupd %xmm1,(%rbx,%rax,8)
mov    %ebp,%eax
add    $0x38,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define i32 @quarticroots(double %0, double %1, double %2, double %3, double* %4) {
%6 = alloca i32
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double*
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca i32
store double %0, double* %7
store double %1, double* %8
store double %2, double* %9
store double %3, double* %10
store double* %4, double** %11
store i32 0, i32* %27
%28 = load double, double* %7
%29 = load double, double* %7
%30 = fmul double %28, %29
store double %30, double* %22
%31 = load double, double* %8
%32 = load double, double* %22
%33 = fmul double %32, 3.000000e+00
%34 = fdiv double %33, 8.000000e+00
%35 = fsub double %31, %34
store double %35, double* %16
%36 = load double, double* %9
%37 = load double, double* %7
%38 = load double, double* %22
%39 = fdiv double %38, 8.000000e+00
%40 = load double, double* %8
%41 = fdiv double %40, 2.000000e+00
%42 = fsub double %39, %41
%43 = call double @llvm.fmuladd.f64(double %37, double %42, double %36)
store double %43, double* %17
%44 = load double, double* %10
%45 = load double, double* %22
%46 = load double, double* %8
%47 = fdiv double %46, 1.600000e+01
%48 = load double, double* %22
%49 = fmul double %48, 3.000000e+00
%50 = fdiv double %49, 2.560000e+02
%51 = fsub double %47, %50
%52 = call double @llvm.fmuladd.f64(double %45, double %51, double %44)
%53 = load double, double* %7
%54 = load double, double* %9
%55 = fmul double %53, %54
%56 = fdiv double %55, 4.000000e+00
%57 = fsub double %52, %56
store double %57, double* %18
%58 = load double, double* %16
%59 = fmul double 2.000000e+00, %58
store double %59, double* %13
%60 = load double, double* %16
%61 = load double, double* %16
%62 = load double, double* %18
%63 = fmul double 4.000000e+00, %62
%64 = fneg double %63
%65 = call double @llvm.fmuladd.f64(double %60, double %61, double %64)
store double %65, double* %14
%66 = load double, double* %17
%67 = fneg double %66
%68 = load double, double* %17
%69 = fmul double %67, %68
store double %69, double* %15
%70 = load double, double* %13
%71 = load double, double* %14
%72 = load double, double* %15
%73 = call double @cubic(double %70, double %71, double %72)
store double %73, double* %12
%74 = load double, double* %12
%75 = fcmp ole double %74, 0.000000e+00
br i1 %75, label %76, label %77
76: ; preds = %5
store i32 0, i32* %6
br label %164
77: ; preds = %5
%78 = load double, double* %12
%79 = call double @sqrt(double %78)
store double %79, double* %19
%80 = load double, double* %16
%81 = fmul double 2.500000e-01, %80
store double %81, double* %25
%82 = load double, double* %17
%83 = load double, double* %19
%84 = fdiv double %82, %83
store double %84, double* %26
%85 = load double, double* %12
%86 = load double, double* %16
%87 = fadd double %85, %86
%88 = load double, double* %26
%89 = fadd double %87, %88
%90 = fmul double 5.000000e-01, %89
store double %90, double* %20
%91 = load double, double* %12
%92 = load double, double* %16
%93 = fadd double %91, %92
%94 = load double, double* %26
%95 = fsub double %93, %94
%96 = fmul double 5.000000e-01, %95
store double %96, double* %21
%97 = load double, double* %12
%98 = load double, double* %20
%99 = call double @llvm.fmuladd.f64(double -4.000000e+00, double %98, double %97)
store double %99, double* %23
%100 = load double, double* %23
%101 = fcmp oge double %100, 0.000000e+00
br i1 %101, label %102, label %129
102: ; preds = %77
%103 = load double, double* %23
%104 = call double @sqrt(double %103)
store double %104, double* %24
%105 = load double, double* %19
%106 = fneg double %105
%107 = load double, double* %24
%108 = fsub double %106, %107
%109 = load double, double* %25
%110 = fneg double %109
%111 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %108, double %110)
%112 = load double*, double** %11
%113 = load i32, i32* %27
%114 = add nsw i32 %113, 1
store i32 %114, i32* %27
%115 = sext i32 %113 to i64
%116 = getelementptr inbounds double, double* %112, i64 %115
store double %111, double* %116
%117 = load double, double* %19
%118 = fneg double %117
%119 = load double, double* %24
%120 = fadd double %118, %119
%121 = load double, double* %25
%122 = fneg double %121
%123 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %120, double %122)
%124 = load double*, double** %11
%125 = load i32, i32* %27
%126 = add nsw i32 %125, 1
store i32 %126, i32* %27
%127 = sext i32 %125 to i64
%128 = getelementptr inbounds double, double* %124, i64 %127
store double %123, double* %128
br label %129
129: ; preds = %102, %77
%130 = load double, double* %12
%131 = load double, double* %21
%132 = call double @llvm.fmuladd.f64(double -4.000000e+00, double %131, double %130)
store double %132, double* %23
%133 = load double, double* %23
%134 = fcmp oge double %133, 0.000000e+00
br i1 %134, label %135, label %162
135: ; preds = %129
%136 = load double, double* %23
%137 = call double @sqrt(double %136)
store double %137, double* %24
%138 = load double, double* %19
%139 = fneg double %138
%140 = load double, double* %24
%141 = fsub double %139, %140
%142 = load double, double* %25
%143 = fneg double %142
%144 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %141, double %143)
%145 = load double*, double** %11
%146 = load i32, i32* %27
%147 = add nsw i32 %146, 1
store i32 %147, i32* %27
%148 = sext i32 %146 to i64
%149 = getelementptr inbounds double, double* %145, i64 %148
store double %144, double* %149
%150 = load double, double* %19
%151 = fneg double %150
%152 = load double, double* %24
%153 = fadd double %151, %152
%154 = load double, double* %25
%155 = fneg double %154
%156 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %153, double %155)
%157 = load double*, double** %11
%158 = load i32, i32* %27
%159 = add nsw i32 %158, 1
store i32 %159, i32* %27
%160 = sext i32 %158 to i64
%161 = getelementptr inbounds double, double* %157, i64 %160
store double %156, double* %161
br label %162
162: ; preds = %135, %129
%163 = load i32, i32* %27
store i32 %163, i32* %6
br label %164
164: ; preds = %162, %76
%165 = load i32, i32* %6
ret i32 %165
}

LLVM IR (-O2):
define i32 @quarticroots(double %0, double %1, double %2, double %3, double* nocapture writeonly %4) {
%6 = fmul double %0, %0
%7 = fmul double %6, 3.000000e+00
%8 = fmul double %7, 1.250000e-01
%9 = fsub double %1, %8
%10 = fmul double %6, 1.250000e-01
%11 = fmul double %1, 5.000000e-01
%12 = fsub double %10, %11
%13 = tail call double @llvm.fmuladd.f64(double %0, double %12, double %2)
%14 = fmul double %1, 6.250000e-02
%15 = fmul double %7, 3.906250e-03
%16 = fsub double %14, %15
%17 = tail call double @llvm.fmuladd.f64(double %6, double %16, double %3)
%18 = fmul double %0, %2
%19 = fmul double %18, 2.500000e-01
%20 = fsub double %17, %19
%21 = fmul double %9, 2.000000e+00
%22 = fmul double %20, -4.000000e+00
%23 = tail call double @llvm.fmuladd.f64(double %9, double %9, double %22)
%24 = fneg double %13
%25 = fmul double %13, %24
%26 = tail call double @cubic(double %21, double %23, double %25)
%27 = fcmp ugt double %26, 0.000000e+00
br i1 %27, label %28, label %72
28: ; preds = %5
%29 = tail call double @sqrt(double %26)
%30 = fmul double %9, 2.500000e-01
%31 = fdiv double %13, %29
%32 = fadd double %9, %26
%33 = fadd double %32, %31
%34 = fmul double %33, 5.000000e-01
%35 = fsub double %32, %31
%36 = fmul double %35, 5.000000e-01
%37 = tail call double @llvm.fmuladd.f64(double %34, double -4.000000e+00, double %26)
%38 = fcmp ult double %37, 0.000000e+00
br i1 %38, label %52, label %39
39: ; preds = %28
%40 = tail call double @sqrt(double %37)
%41 = fneg double %29
%42 = fneg double %30
%43 = insertelement <2 x double> poison, double %41, i64 0
%44 = insertelement <2 x double> %43, double %40, i64 1
%45 = insertelement <2 x double> poison, double %40, i64 0
%46 = insertelement <2 x double> %45, double %29, i64 1
%47 = fsub <2 x double> %44, %46
%48 = insertelement <2 x double> poison, double %42, i64 0
%49 = shufflevector <2 x double> %48, <2 x double> poison, <2 x i32> zeroinitializer
%50 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %47, <2 x double> <double 5.000000e-01, double 5.000000e-01>, <2 x double> %49)
%51 = bitcast double* %4 to <2 x double>*
store <2 x double> %50, <2 x double>* %51
br label %52
52: ; preds = %39, %28
%53 = phi i32 [ 2, %39 ], [ 0, %28 ]
%54 = tail call double @llvm.fmuladd.f64(double %36, double -4.000000e+00, double %26)
%55 = fcmp ult double %54, 0.000000e+00
br i1 %55, label %72, label %56
56: ; preds = %52
%57 = tail call double @sqrt(double %54)
%58 = fneg double %29
%59 = fneg double %30
%60 = zext i32 %53 to i64
%61 = getelementptr inbounds double, double* %4, i64 %60
%62 = insertelement <2 x double> poison, double %58, i64 0
%63 = insertelement <2 x double> %62, double %57, i64 1
%64 = insertelement <2 x double> poison, double %29, i64 1
%65 = insertelement <2 x double> %64, double %57, i64 0
%66 = fsub <2 x double> %63, %65
%67 = insertelement <2 x double> poison, double %59, i64 0
%68 = shufflevector <2 x double> %67, <2 x double> poison, <2 x i32> zeroinitializer
%69 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %66, <2 x double> <double 5.000000e-01, double 5.000000e-01>, <2 x double> %68)
%70 = add nuw nsw i32 %53, 2
%71 = bitcast double* %61 to <2 x double>*
store <2 x double> %69, <2 x double>* %71
br label %72
72: ; preds = %52, %56, %5
%73 = phi i32 [ 0, %5 ], [ %70, %56 ], [ %53, %52 ]
ret i32 %73
}