#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 DBL_EPSILON;
double fabs(double);
double sqrt(double);
void dsolveSymmetric22(double a[2][2], double v[2][2], double d[2]) {
  double a00 = a[0][0];
  double a01 = a[0][1], a11 = a[1][1];
  double v00 = 1.0, v01 = 0.0;
  double v10 = 0.0, v11 = 1.0;
  double tiny;
  double c, r, s, t, u, vpr, vqr;
  double dt, vt[2];
  if (a01 != 0.0) {
    tiny = 0.1 * sqrt(DBL_EPSILON);
    u = a11 - a00;
    if (fabs(a01) < tiny * fabs(u)) {
      t = a01 / u;
    } else {
      r = 0.5 * u / a01;
      t = (r >= 0.0) ? 1.0 / (r + sqrt(1.0 + r * r))
                     : 1.0 / (r - sqrt(1.0 + r * r));
    }
    c = 1.0 / sqrt(1.0 + t * t);
    s = t * c;
    u = s / (1.0 + c);
    r = t * a01;
    a00 -= r;
    a11 += r;
    vpr = v00;
    vqr = v10;
    v00 = vpr - s * (vqr + vpr * u);
    v10 = vqr + s * (vpr - vqr * u);
    vpr = v01;
    vqr = v11;
    v01 = vpr - s * (vqr + vpr * u);
    v11 = vqr + s * (vpr - vqr * u);
  }
  d[0] = a00;
  d[1] = a11;
  v[0][0] = v00;
  v[0][1] = v01;
  v[1][0] = v10;
  v[1][1] = v11;
  if (d[0] < d[1]) {
    dt = d[1];
    d[1] = d[0];
    d[0] = dt;
    vt[0] = v[1][0];
    vt[1] = v[1][1];
    v[1][0] = v[0][0];
    v[1][1] = v[0][1];
    v[0][0] = vt[0];
    v[0][1] = vt[1];
  }
}


x86 assembly (clang, -O0):
<dsolveSymmetric22>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xea8(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xe8b(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1197 <dsolveSymmetric22+0x87>
jp     1197 <dsolveSymmetric22+0x87>
jmp    14de <dsolveSymmetric22+0x3ce>
mov    0x2e4a(%rip),%rax
movsd  (%rax),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xe56(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x28(%rbp),%xmm1
movaps 0xe4a(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movaps 0xe35(%rip),%xmm3
pand   %xmm3,%xmm2
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    1211 <dsolveSymmetric22+0x101>
movsd  -0x28(%rbp),%xmm0
divsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    12e4 <dsolveSymmetric22+0x1d4>
movsd  0xdf7(%rip),%xmm0
mulsd  -0x80(%rbp),%xmm0
divsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     128b <dsolveSymmetric22+0x17b>
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xda3(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  -0xc0(%rbp),%xmm1
ff
addsd  %xmm0,%xmm1
movsd  0xd86(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
jmp    12d7 <dsolveSymmetric22+0x1c7>
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd52(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  -0xc8(%rbp),%xmm1
ff
subsd  %xmm0,%xmm1
movsd  0xd35(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd06(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xcf2(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  0xccd(%rip),%xmm1
addsd  -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x68(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x70(%rbp),%xmm0
movsd  -0x90(%rbp),%xmm3
ff
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x70(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm3
ff
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x70(%rbp),%xmm0
movsd  -0x90(%rbp),%xmm3
ff
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x70(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm3
ff
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x10(%rax)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x18(%rax)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    15e3 <dsolveSymmetric22+0x4d3>
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x10(%rax)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x18(%rax)
movsd  -0xb0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0xa8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dsolveSymmetric22>:
push   %r14
push   %rbx
sub    $0x38,%rsp
mov    %rdx,%rbx
mov    %rsi,%r14
movsd  (%rdi),%xmm4
movsd  0x8(%rdi),%xmm3
movsd  0x18(%rdi),%xmm5
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm3
jne    1137 <dsolveSymmetric22+0x27>
jnp    115f <dsolveSymmetric22+0x4f>
mov    0x2eaa(%rip),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
movsd  %xmm4,0x8(%rsp)
movsd  %xmm5,(%rsp)
movapd %xmm3,0x20(%rsp)
jb     1178 <dsolveSymmetric22+0x68>
sqrtsd %xmm0,%xmm0
jmp    118e <dsolveSymmetric22+0x7e>
xorpd  %xmm3,%xmm3
movhpd 0xe95(%rip),%xmm3
movsd  0xe8d(%rip),%xmm2
jmp    12ff <dsolveSymmetric22+0x1ef>
call   1030 <sqrt@plt>
movapd 0x20(%rsp),%xmm3
movsd  (%rsp),%xmm5
movsd  0x8(%rsp),%xmm4
mulsd  0xe72(%rip),%xmm0
movapd %xmm5,%xmm7
subsd  %xmm4,%xmm7
movapd 0xe7a(%rip),%xmm1
movapd %xmm3,%xmm2
andpd  %xmm1,%xmm2
andpd  %xmm7,%xmm1
mulsd  %xmm0,%xmm1
ucomisd %xmm2,%xmm1
jbe    11c6 <dsolveSymmetric22+0xb6>
movapd %xmm3,%xmm6
divsd  %xmm7,%xmm6
jmp    1245 <dsolveSymmetric22+0x135>
mulsd  0xe42(%rip),%xmm7
divsd  %xmm3,%xmm7
movapd %xmm7,%xmm0
mulsd  %xmm7,%xmm0
addsd  0xe1e(%rip),%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
jb     11f2 <dsolveSymmetric22+0xe2>
sqrtsd %xmm0,%xmm0
jmp    1218 <dsolveSymmetric22+0x108>
movapd %xmm7,0x10(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm2,%xmm2
movapd 0x10(%rsp),%xmm7
movapd 0x20(%rsp),%xmm3
movsd  (%rsp),%xmm5
movsd  0x8(%rsp),%xmm4
movapd 0xe10(%rip),%xmm1
xorpd  %xmm0,%xmm1
cmpnlesd %xmm7,%xmm2
andpd  %xmm2,%xmm1
andnpd %xmm0,%xmm2
orpd   %xmm1,%xmm2
addsd  %xmm7,%xmm2
movsd  0xdbf(%rip),%xmm6
divsd  %xmm2,%xmm6
movapd %xmm6,%xmm0
mulsd  %xmm6,%xmm0
addsd  0xdab(%rip),%xmm0
xorpd  %xmm7,%xmm7
ucomisd %xmm7,%xmm0
jb     1265 <dsolveSymmetric22+0x155>
sqrtsd %xmm0,%xmm0
jmp    128b <dsolveSymmetric22+0x17b>
movapd %xmm6,0x10(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm7,%xmm7
movapd 0x10(%rsp),%xmm6
movapd 0x20(%rsp),%xmm3
movsd  (%rsp),%xmm5
movsd  0x8(%rsp),%xmm4
movsd  0xd6d(%rip),%xmm1
movapd %xmm1,%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm6,%xmm3
mulsd  %xmm2,%xmm6
addsd  %xmm1,%xmm2
movapd %xmm6,%xmm0
divsd  %xmm2,%xmm0
subsd  %xmm3,%xmm4
addsd  %xmm3,%xmm5
movapd %xmm0,%xmm3
addsd  %xmm7,%xmm3
mulsd  %xmm0,%xmm7
addsd  %xmm1,%xmm7
unpcklpd %xmm7,%xmm3
unpcklpd %xmm6,%xmm6
mulpd  %xmm6,%xmm3
movsd  0xd25(%rip),%xmm2
unpcklpd %xmm0,%xmm0
mulpd  0xd59(%rip),%xmm0
addpd  %xmm2,%xmm0
subpd  %xmm3,%xmm2
mulpd  %xmm6,%xmm0
xorpd  %xmm3,%xmm3
unpcklpd %xmm1,%xmm3
addpd  %xmm0,%xmm3
movsd  %xmm4,(%rbx)
movsd  %xmm5,0x8(%rbx)
movupd %xmm2,(%r14)
movupd %xmm3,0x10(%r14)
movsd  (%rbx),%xmm0
movsd  0x8(%rbx),%xmm1
ucomisd %xmm0,%xmm1
jbe    1341 <dsolveSymmetric22+0x231>
movsd  %xmm0,0x8(%rbx)
movsd  %xmm1,(%rbx)
movupd (%r14),%xmm0
movupd 0x10(%r14),%xmm1
movupd %xmm0,0x10(%r14)
movupd %xmm1,(%r14)
add    $0x38,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
@DBL_EPSILON = global double 0.000000e+00
define void @dsolveSymmetric22([2 x double]* %0, [2 x double]* %1, double* %2) {
%4 = alloca [2 x double]*
%5 = alloca [2 x double]*
%6 = alloca double*
%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 [2 x double]
store [2 x double]* %0, [2 x double]** %4
store [2 x double]* %1, [2 x double]** %5
store double* %2, double** %6
%24 = load [2 x double]*, [2 x double]** %4
%25 = getelementptr inbounds [2 x double], [2 x double]* %24, i64 0
%26 = getelementptr inbounds [2 x double], [2 x double]* %25, i64 0, i64 0
%27 = load double, double* %26
store double %27, double* %7
%28 = load [2 x double]*, [2 x double]** %4
%29 = getelementptr inbounds [2 x double], [2 x double]* %28, i64 0
%30 = getelementptr inbounds [2 x double], [2 x double]* %29, i64 0, i64 1
%31 = load double, double* %30
store double %31, double* %8
%32 = load [2 x double]*, [2 x double]** %4
%33 = getelementptr inbounds [2 x double], [2 x double]* %32, i64 1
%34 = getelementptr inbounds [2 x double], [2 x double]* %33, i64 0, i64 1
%35 = load double, double* %34
store double %35, double* %9
store double 1.000000e+00, double* %10
store double 0.000000e+00, double* %11
store double 0.000000e+00, double* %12
store double 1.000000e+00, double* %13
%36 = load double, double* %8
%37 = fcmp une double %36, 0.000000e+00
br i1 %37, label %38, label %139
38: ; preds = %3
%39 = load double, double* @DBL_EPSILON
%40 = call double @sqrt(double %39)
%41 = fmul double 1.000000e-01, %40
store double %41, double* %14
%42 = load double, double* %9
%43 = load double, double* %7
%44 = fsub double %42, %43
store double %44, double* %19
%45 = load double, double* %8
%46 = call double @llvm.fabs.f64(double %45)
%47 = load double, double* %14
%48 = load double, double* %19
%49 = call double @llvm.fabs.f64(double %48)
%50 = fmul double %47, %49
%51 = fcmp olt double %46, %50
br i1 %51, label %52, label %56
52: ; preds = %38
%53 = load double, double* %8
%54 = load double, double* %19
%55 = fdiv double %53, %54
store double %55, double* %18
br label %81
56: ; preds = %38
%57 = load double, double* %19
%58 = fmul double 5.000000e-01, %57
%59 = load double, double* %8
%60 = fdiv double %58, %59
store double %60, double* %16
%61 = load double, double* %16
%62 = fcmp oge double %61, 0.000000e+00
br i1 %62, label %63, label %71
63: ; preds = %56
%64 = load double, double* %16
%65 = load double, double* %16
%66 = load double, double* %16
%67 = call double @llvm.fmuladd.f64(double %65, double %66, double 1.000000e+00)
%68 = call double @sqrt(double %67)
%69 = fadd double %64, %68
%70 = fdiv double 1.000000e+00, %69
br label %79
71: ; preds = %56
%72 = load double, double* %16
%73 = load double, double* %16
%74 = load double, double* %16
%75 = call double @llvm.fmuladd.f64(double %73, double %74, double 1.000000e+00)
%76 = call double @sqrt(double %75)
%77 = fsub double %72, %76
%78 = fdiv double 1.000000e+00, %77
br label %79
79: ; preds = %71, %63
%80 = phi double [ %70, %63 ], [ %78, %71 ]
store double %80, double* %18
br label %81
81: ; preds = %79, %52
%82 = load double, double* %18
%83 = load double, double* %18
%84 = call double @llvm.fmuladd.f64(double %82, double %83, double 1.000000e+00)
%85 = call double @sqrt(double %84)
%86 = fdiv double 1.000000e+00, %85
store double %86, double* %15
%87 = load double, double* %18
%88 = load double, double* %15
%89 = fmul double %87, %88
store double %89, double* %17
%90 = load double, double* %17
%91 = load double, double* %15
%92 = fadd double 1.000000e+00, %91
%93 = fdiv double %90, %92
store double %93, double* %19
%94 = load double, double* %18
%95 = load double, double* %8
%96 = fmul double %94, %95
store double %96, double* %16
%97 = load double, double* %16
%98 = load double, double* %7
%99 = fsub double %98, %97
store double %99, double* %7
%100 = load double, double* %16
%101 = load double, double* %9
%102 = fadd double %101, %100
store double %102, double* %9
%103 = load double, double* %10
store double %103, double* %20
%104 = load double, double* %12
store double %104, double* %21
%105 = load double, double* %20
%106 = load double, double* %17
%107 = load double, double* %21
%108 = load double, double* %20
%109 = load double, double* %19
%110 = call double @llvm.fmuladd.f64(double %108, double %109, double %107)
%111 = fneg double %106
%112 = call double @llvm.fmuladd.f64(double %111, double %110, double %105)
store double %112, double* %10
%113 = load double, double* %21
%114 = load double, double* %17
%115 = load double, double* %20
%116 = load double, double* %21
%117 = load double, double* %19
%118 = fneg double %116
%119 = call double @llvm.fmuladd.f64(double %118, double %117, double %115)
%120 = call double @llvm.fmuladd.f64(double %114, double %119, double %113)
store double %120, double* %12
%121 = load double, double* %11
store double %121, double* %20
%122 = load double, double* %13
store double %122, double* %21
%123 = load double, double* %20
%124 = load double, double* %17
%125 = load double, double* %21
%126 = load double, double* %20
%127 = load double, double* %19
%128 = call double @llvm.fmuladd.f64(double %126, double %127, double %125)
%129 = fneg double %124
%130 = call double @llvm.fmuladd.f64(double %129, double %128, double %123)
store double %130, double* %11
%131 = load double, double* %21
%132 = load double, double* %17
%133 = load double, double* %20
%134 = load double, double* %21
%135 = load double, double* %19
%136 = fneg double %134
%137 = call double @llvm.fmuladd.f64(double %136, double %135, double %133)
%138 = call double @llvm.fmuladd.f64(double %132, double %137, double %131)
store double %138, double* %13
br label %139
139: ; preds = %81, %3
%140 = load double, double* %7
%141 = load double*, double** %6
%142 = getelementptr inbounds double, double* %141, i64 0
store double %140, double* %142
%143 = load double, double* %9
%144 = load double*, double** %6
%145 = getelementptr inbounds double, double* %144, i64 1
store double %143, double* %145
%146 = load double, double* %10
%147 = load [2 x double]*, [2 x double]** %5
%148 = getelementptr inbounds [2 x double], [2 x double]* %147, i64 0
%149 = getelementptr inbounds [2 x double], [2 x double]* %148, i64 0, i64 0
store double %146, double* %149
%150 = load double, double* %11
%151 = load [2 x double]*, [2 x double]** %5
%152 = getelementptr inbounds [2 x double], [2 x double]* %151, i64 0
%153 = getelementptr inbounds [2 x double], [2 x double]* %152, i64 0, i64 1
store double %150, double* %153
%154 = load double, double* %12
%155 = load [2 x double]*, [2 x double]** %5
%156 = getelementptr inbounds [2 x double], [2 x double]* %155, i64 1
%157 = getelementptr inbounds [2 x double], [2 x double]* %156, i64 0, i64 0
store double %154, double* %157
%158 = load double, double* %13
%159 = load [2 x double]*, [2 x double]** %5
%160 = getelementptr inbounds [2 x double], [2 x double]* %159, i64 1
%161 = getelementptr inbounds [2 x double], [2 x double]* %160, i64 0, i64 1
store double %158, double* %161
%162 = load double*, double** %6
%163 = getelementptr inbounds double, double* %162, i64 0
%164 = load double, double* %163
%165 = load double*, double** %6
%166 = getelementptr inbounds double, double* %165, i64 1
%167 = load double, double* %166
%168 = fcmp olt double %164, %167
br i1 %168, label %169, label %215
169: ; preds = %139
%170 = load double*, double** %6
%171 = getelementptr inbounds double, double* %170, i64 1
%172 = load double, double* %171
store double %172, double* %22
%173 = load double*, double** %6
%174 = getelementptr inbounds double, double* %173, i64 0
%175 = load double, double* %174
%176 = load double*, double** %6
%177 = getelementptr inbounds double, double* %176, i64 1
store double %175, double* %177
%178 = load double, double* %22
%179 = load double*, double** %6
%180 = getelementptr inbounds double, double* %179, i64 0
store double %178, double* %180
%181 = load [2 x double]*, [2 x double]** %5
%182 = getelementptr inbounds [2 x double], [2 x double]* %181, i64 1
%183 = getelementptr inbounds [2 x double], [2 x double]* %182, i64 0, i64 0
%184 = load double, double* %183
%185 = getelementptr inbounds [2 x double], [2 x double]* %23, i64 0, i64 0
store double %184, double* %185
%186 = load [2 x double]*, [2 x double]** %5
%187 = getelementptr inbounds [2 x double], [2 x double]* %186, i64 1
%188 = getelementptr inbounds [2 x double], [2 x double]* %187, i64 0, i64 1
%189 = load double, double* %188
%190 = getelementptr inbounds [2 x double], [2 x double]* %23, i64 0, i64 1
store double %189, double* %190
%191 = load [2 x double]*, [2 x double]** %5
%192 = getelementptr inbounds [2 x double], [2 x double]* %191, i64 0
%193 = getelementptr inbounds [2 x double], [2 x double]* %192, i64 0, i64 0
%194 = load double, double* %193
%195 = load [2 x double]*, [2 x double]** %5
%196 = getelementptr inbounds [2 x double], [2 x double]* %195, i64 1
%197 = getelementptr inbounds [2 x double], [2 x double]* %196, i64 0, i64 0
store double %194, double* %197
%198 = load [2 x double]*, [2 x double]** %5
%199 = getelementptr inbounds [2 x double], [2 x double]* %198, i64 0
%200 = getelementptr inbounds [2 x double], [2 x double]* %199, i64 0, i64 1
%201 = load double, double* %200
%202 = load [2 x double]*, [2 x double]** %5
%203 = getelementptr inbounds [2 x double], [2 x double]* %202, i64 1
%204 = getelementptr inbounds [2 x double], [2 x double]* %203, i64 0, i64 1
store double %201, double* %204
%205 = getelementptr inbounds [2 x double], [2 x double]* %23, i64 0, i64 0
%206 = load double, double* %205
%207 = load [2 x double]*, [2 x double]** %5
%208 = getelementptr inbounds [2 x double], [2 x double]* %207, i64 0
%209 = getelementptr inbounds [2 x double], [2 x double]* %208, i64 0, i64 0
store double %206, double* %209
%210 = getelementptr inbounds [2 x double], [2 x double]* %23, i64 0, i64 1
%211 = load double, double* %210
%212 = load [2 x double]*, [2 x double]** %5
%213 = getelementptr inbounds [2 x double], [2 x double]* %212, i64 0
%214 = getelementptr inbounds [2 x double], [2 x double]* %213, i64 0, i64 1
store double %211, double* %214
br label %215
215: ; preds = %169, %139
ret void
}

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