#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 *dvector(int, int);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, ...);
int /*<<< orphan*/ free_dvector(double *, int, int);
int /*<<< orphan*/ ldl_dcmp(double **, int, double *, double *, double *, int,
                            int /*<<< orphan*/, int *);
double sqrt(double);
int /*<<< orphan*/ stderr;
int /*<<< orphan*/ stdout;
int sturm(double **K, double **M, int n, int m, double shift, double ws,
          int verbose) {
  double ws_shift, *d;
  int ok = 0, i, j, modes;
  d = dvector(1, n);
  modes = (int)((float)(0.5 * m) > (float)(m - 8.0) ? (int)(m / 2.0) : m - 8);
  ws_shift = ws + shift;
  for (i = 1; i <= n; i++)
    for (j = i; j <= n; j++)
      K[i][j] -= ws_shift * M[i][j];
  ldl_dcmp(K, n, d, d, d, 1, 0, &ok);
  if (verbose)
    fprintf(stdout, "  There are %d modes below %f Hz.", -ok,
            sqrt(ws) /
                (2.0 * 3.14159265358979323846264338327950288419716939937510));
  if (-ok > modes) {
    fprintf(stderr, " ... %d modes were not found.\n", -ok - modes);
    fprintf(stderr, " Try increasing the number of modes in \n");
    fprintf(stderr, " order to get the missing modes below %f Hz.\n",
            sqrt(ws) /
                (2.0 * 3.14159265358979323846264338327950288419716939937510));
  } else if (verbose)
    fprintf(stdout, "  All %d modes were found.\n", modes);
  for (i = 1; i <= n; i++)
    for (j = i; j <= n; j++)
      K[i][j] += ws_shift * M[i][j];
  free_dvector(d, 1, n);
  return ok;
}


x86 assembly (clang, -O0):
<sturm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
call   1050 <dvector@plt>
mov    %rax,-0x40(%rbp)
cvtsi2sdl -0x18(%rbp),%xmm1
movsd  0xe6f(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
movsd  0xe52(%rip),%xmm2
subsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    11dc <sturm+0x8c>
cvtsi2sdl -0x18(%rbp),%xmm0
movsd  0xe44(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x54(%rbp)
jmp    11e5 <sturm+0x95>
mov    -0x18(%rbp),%eax
sub    $0x8,%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x50(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1293 <sturm+0x143>
mov    -0x48(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1280 <sturm+0x130>
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1213 <sturm+0xc3>
jmp    1285 <sturm+0x135>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1201 <sturm+0xb1>
mov    -0x8(%rbp),%rdi
mov    -0x14(%rbp),%esi
mov    -0x40(%rbp),%rdx
mov    -0x40(%rbp),%rcx
mov    -0x40(%rbp),%r8
mov    $0x1,%r9d
xor    %eax,%eax
lea    -0x44(%rbp),%rax
movl   $0x0,(%rsp)
mov    %rax,0x8(%rsp)
call   1040 <ldl_dcmp@plt>
cmpl   $0x0,-0x2c(%rbp)
je     130b <sturm+0x1bb>
mov    0x2d04(%rip),%rax
mov    (%rax),%eax
mov    %eax,-0x5c(%rbp)
xor    %eax,%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0x58(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1070 <sqrt@plt>
mov    -0x5c(%rbp),%edi
mov    -0x58(%rbp),%edx
movsd  0xd1f(%rip),%xmm1
divsd  %xmm1,%xmm0
lea    0xd1c(%rip),%rsi
mov    $0x1,%al
call   1030 <fprintf@plt>
xor    %eax,%eax
sub    -0x44(%rbp),%eax
cmp    -0x50(%rbp),%eax
jle    1387 <sturm+0x237>
mov    0x2cc0(%rip),%rax
mov    (%rax),%edi
xor    %edx,%edx
sub    -0x44(%rbp),%edx
sub    -0x50(%rbp),%edx
lea    0xd11(%rip),%rsi
mov    $0x0,%al
call   1030 <fprintf@plt>
mov    0x2ca1(%rip),%rax
mov    (%rax),%edi
lea    0xd19(%rip),%rsi
mov    $0x0,%al
call   1030 <fprintf@plt>
mov    0x2c8a(%rip),%rax
mov    (%rax),%eax
mov    %eax,-0x60(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1070 <sqrt@plt>
mov    -0x60(%rbp),%edi
movsd  0xca8(%rip),%xmm1
divsd  %xmm1,%xmm0
lea    0xd0f(%rip),%rsi
mov    $0x1,%al
call   1030 <fprintf@plt>
jmp    13b0 <sturm+0x260>
cmpl   $0x0,-0x2c(%rbp)
je     13ab <sturm+0x25b>
mov    0x2c40(%rip),%rax
mov    (%rax),%edi
mov    -0x50(%rbp),%edx
lea    0xd14(%rip),%rsi
mov    $0x0,%al
call   1030 <fprintf@plt>
jmp    13b0 <sturm+0x260>
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1432 <sturm+0x2e2>
mov    -0x48(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     141f <sturm+0x2cf>
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    13c9 <sturm+0x279>
jmp    1424 <sturm+0x2d4>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    13b7 <sturm+0x267>
mov    -0x40(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
call   1060 <free_dvector@plt>
mov    -0x44(%rbp),%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sturm>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r8d,0x2c(%rsp)
movsd  %xmm1,0x18(%rsp)
movaps %xmm0,0x40(%rsp)
mov    %ecx,%ebx
mov    %edx,%ebp
mov    %rsi,0x38(%rsp)
mov    %rdi,0x20(%rsp)
movl   $0x0,0xc(%rsp)
mov    $0x1,%edi
mov    %edx,%esi
call   1050 <dvector@plt>
movapd 0x40(%rsp),%xmm7
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  0xe5b(%rip),%xmm1
mulsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm2
addsd  0xe53(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
add    $0xfffffff8,%ebx
ucomiss %xmm0,%xmm2
cvttsd2si %xmm1,%ecx
cmovbe %ebx,%ecx
mov    %ecx,0x10(%rsp)
mov    %rax,0x30(%rsp)
addsd  0x18(%rsp),%xmm7
mov    %ebp,%r14d
mov    %ebp,0x14(%rsp)
test   %ebp,%ebp
jle    13f7 <sturm+0x2a7>
movapd 0xe34(%rip),%xmm0
xorpd  %xmm7,%xmm0
lea    0x1(%r14),%rax
pshufd $0x44,%xmm0,%xmm1
mov    %r14,%r15
not    %r15
lea    -0x4(%r14),%r11
mov    $0x1,%r13d
mov    $0x7,%r8d
xor    %ebp,%ebp
jmp    1239 <sturm+0xe9>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r13
add    $0x1,%rbp
add    $0x1,%r8
add    $0xffffffffffffffff,%r11
cmp    %rax,%r13
je     13f7 <sturm+0x2a7>
mov    %rbp,%rdx
not    %rdx
lea    (%rdx,%rax,1),%r9
mov    0x38(%rsp),%rcx
mov    (%rcx,%r13,8),%rcx
mov    0x20(%rsp),%rsi
mov    (%rsi,%r13,8),%rbx
mov    %r13,%rsi
cmp    $0x4,%r9
jb     1390 <sturm+0x240>
lea    (%rbx,%r13,8),%rsi
lea    (%rcx,%rax,8),%rdi
cmp    %rdi,%rsi
jae    1283 <sturm+0x133>
lea    (%rbx,%rax,8),%r10
lea    (%rcx,%r13,8),%rdi
mov    %r13,%rsi
cmp    %r10,%rdi
jb     1390 <sturm+0x240>
mov    %r15,%rdi
add    %rax,%rdx
add    $0xfffffffffffffffc,%rdx
mov    %rdx,%r10
shr    $0x2,%r10
add    $0x1,%r10
cmp    $0x4,%rdx
jae    12a5 <sturm+0x155>
xor    %edx,%edx
jmp    133b <sturm+0x1eb>
mov    %r11,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
and    $0xfffffffffffffffe,%rsi
lea    (%rbx,%r8,8),%r12
lea    (%rcx,%r8,8),%r15
xor    %edx,%edx
xchg   %ax,%ax
movupd -0x30(%r15,%rdx,8),%xmm2
movupd -0x20(%r15,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
movupd -0x30(%r12,%rdx,8),%xmm4
addpd  %xmm2,%xmm4
movupd -0x20(%r12,%rdx,8),%xmm2
movupd -0x10(%r12,%rdx,8),%xmm5
movupd (%r12,%rdx,8),%xmm6
mulpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm4,-0x30(%r12,%rdx,8)
movupd %xmm3,-0x20(%r12,%rdx,8)
movupd -0x10(%r15,%rdx,8),%xmm2
movupd (%r15,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
addpd  %xmm5,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm6,%xmm3
movupd %xmm2,-0x10(%r12,%rdx,8)
movupd %xmm3,(%r12,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rsi
jne    12c0 <sturm+0x170>
mov    %r9,%rsi
and    $0xfffffffffffffffc,%rsi
test   $0x1,%r10b
je     137c <sturm+0x22c>
add    %r13,%rdx
movupd (%rcx,%rdx,8),%xmm2
movupd 0x10(%rcx,%rdx,8),%xmm3
movupd (%rbx,%rdx,8),%xmm4
movupd 0x10(%rbx,%rdx,8),%xmm5
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,(%rbx,%rdx,8)
movupd %xmm3,0x10(%rbx,%rdx,8)
cmp    %rsi,%r9
mov    %rdi,%r15
je     1220 <sturm+0xd0>
add    %r13,%rsi
nopl   0x0(%rax,%rax,1)
mov    %eax,%edi
sub    %esi,%edi
mov    %rsi,%rdx
test   $0x1,%dil
je     13b4 <sturm+0x264>
movsd  (%rcx,%rsi,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  (%rbx,%rsi,8),%xmm2
movsd  %xmm2,(%rbx,%rsi,8)
lea    0x1(%rsi),%rdx
not    %rsi
cmp    %r15,%rsi
je     1220 <sturm+0xd0>
movsd  (%rcx,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  (%rbx,%rdx,8),%xmm2
movsd  %xmm2,(%rbx,%rdx,8)
movsd  0x8(%rcx,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  0x8(%rbx,%rdx,8),%xmm2
movsd  %xmm2,0x8(%rbx,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rax
jne    13c0 <sturm+0x270>
jmp    1220 <sturm+0xd0>
movapd %xmm7,0x40(%rsp)
lea    0xc(%rsp),%rax
mov    0x20(%rsp),%rdi
mov    0x14(%rsp),%esi
mov    0x30(%rsp),%rcx
mov    %rcx,%rdx
mov    %rcx,%r8
mov    $0x1,%r9d
push   %rax
push   $0x0
call   1040 <ldl_dcmp@plt>
add    $0x10,%rsp
cmpl   $0x0,0x2c(%rsp)
je     1458 <sturm+0x308>
mov    0x2ba2(%rip),%r15
mov    (%r15),%ebp
xor    %r13d,%r13d
xor    %ebx,%ebx
xorpd  %xmm1,%xmm1
sub    0xc(%rsp),%ebx
movsd  0x18(%rsp),%xmm0
ucomisd %xmm1,%xmm0
jb     146e <sturm+0x31e>
sqrtsd %xmm0,%xmm0
jmp    1473 <sturm+0x323>
xor    %r13d,%r13d
sub    0xc(%rsp),%r13d
mov    0x10(%rsp),%edx
cmp    %r13d,%edx
jl     149b <sturm+0x34b>
jmp    1510 <sturm+0x3c0>
call   1070 <sqrt@plt>
divsd  0xb95(%rip),%xmm0
lea    0xbae(%rip),%rsi
mov    %ebp,%edi
mov    %ebx,%edx
mov    $0x1,%al
call   1030 <fprintf@plt>
sub    0xc(%rsp),%r13d
mov    0x10(%rsp),%edx
cmp    %r13d,%edx
jge    14e0 <sturm+0x390>
mov    0x2b3e(%rip),%rbx
mov    (%rbx),%edi
sub    %edx,%r13d
lea    0xba4(%rip),%rsi
mov    %r13d,%edx
xor    %eax,%eax
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xbb0(%rip),%rsi
xor    %eax,%eax
call   1030 <fprintf@plt>
mov    (%rbx),%ebx
xorpd  %xmm1,%xmm1
movsd  0x18(%rsp),%xmm0
ucomisd %xmm1,%xmm0
jb     14f3 <sturm+0x3a3>
sqrtsd %xmm0,%xmm0
jmp    14f8 <sturm+0x3a8>
mov    (%r15),%edi
lea    0xbde(%rip),%rsi
xor    %eax,%eax
call   1030 <fprintf@plt>
jmp    1510 <sturm+0x3c0>
call   1070 <sqrt@plt>
divsd  0xb10(%rip),%xmm0
lea    0xb93(%rip),%rsi
mov    %ebx,%edi
mov    $0x1,%al
call   1030 <fprintf@plt>
mov    0x14(%rsp),%edx
test   %edx,%edx
movapd 0x40(%rsp),%xmm6
jle    1738 <sturm+0x5e8>
lea    0x1(%r14),%r13
movapd %xmm6,%xmm0
unpcklpd %xmm6,%xmm0
mov    %r14,%r8
not    %r8
add    $0xfffffffffffffffc,%r14
mov    $0x1,%r15d
mov    $0x7,%r11d
xor    %r12d,%r12d
jmp    1569 <sturm+0x419>
nopl   0x0(%rax)
add    $0x1,%r15
add    $0x1,%r12
add    $0x1,%r11
add    $0xffffffffffffffff,%r14
cmp    %r13,%r15
je     1738 <sturm+0x5e8>
mov    %r12,%rcx
not    %rcx
lea    (%rcx,%r13,1),%r10
mov    0x38(%rsp),%rax
mov    (%rax,%r15,8),%rbp
mov    0x20(%rsp),%rax
mov    (%rax,%r15,8),%rbx
mov    %r15,%rdi
cmp    $0x4,%r10
jb     16c0 <sturm+0x570>
lea    (%rbx,%r15,8),%rax
lea    0x0(,%r13,8),%rsi
add    %rbp,%rsi
cmp    %rsi,%rax
jae    15c1 <sturm+0x471>
lea    (%rbx,%r13,8),%rsi
lea    0x0(,%r15,8),%rax
add    %rbp,%rax
mov    %r15,%rdi
cmp    %rsi,%rax
jb     16c0 <sturm+0x570>
add    %r13,%rcx
add    $0xfffffffffffffffc,%rcx
mov    %rcx,%r9
shr    $0x2,%r9
add    $0x1,%r9
cmp    $0x4,%rcx
jae    15e0 <sturm+0x490>
xor    %ecx,%ecx
jmp    166f <sturm+0x51f>
mov    %r14,%rax
shr    $0x2,%rax
add    $0x1,%rax
and    $0xfffffffffffffffe,%rax
lea    (%rbx,%r11,8),%rdi
lea    0x0(,%r11,8),%rsi
add    %rbp,%rsi
xor    %ecx,%ecx
movupd -0x30(%rsi,%rcx,8),%xmm1
movupd -0x20(%rsi,%rcx,8),%xmm2
mulpd  %xmm0,%xmm1
movupd -0x30(%rdi,%rcx,8),%xmm3
addpd  %xmm1,%xmm3
movupd -0x20(%rdi,%rcx,8),%xmm1
movupd -0x10(%rdi,%rcx,8),%xmm4
movupd (%rdi,%rcx,8),%xmm5
mulpd  %xmm0,%xmm2
addpd  %xmm1,%xmm2
movupd %xmm3,-0x30(%rdi,%rcx,8)
movupd %xmm2,-0x20(%rdi,%rcx,8)
movupd -0x10(%rsi,%rcx,8),%xmm1
movupd (%rsi,%rcx,8),%xmm2
mulpd  %xmm0,%xmm1
addpd  %xmm4,%xmm1
mulpd  %xmm0,%xmm2
addpd  %xmm5,%xmm2
movupd %xmm1,-0x10(%rdi,%rcx,8)
movupd %xmm2,(%rdi,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rax
jne    1600 <sturm+0x4b0>
mov    %r10,%rdi
and    $0xfffffffffffffffc,%rdi
test   $0x1,%r9b
je     16b1 <sturm+0x561>
add    %r15,%rcx
movupd 0x0(%rbp,%rcx,8),%xmm1
movupd 0x10(%rbp,%rcx,8),%xmm2
movupd (%rbx,%rcx,8),%xmm3
movupd 0x10(%rbx,%rcx,8),%xmm4
mulpd  %xmm0,%xmm1
addpd  %xmm3,%xmm1
mulpd  %xmm0,%xmm2
addpd  %xmm4,%xmm2
movupd %xmm1,(%rbx,%rcx,8)
movupd %xmm2,0x10(%rbx,%rcx,8)
cmp    %rdi,%r10
je     1550 <sturm+0x400>
add    %r15,%rdi
nopl   (%rax)
mov    %r13d,%esi
sub    %edi,%esi
mov    %rdi,%rcx
test   $0x1,%sil
je     16e6 <sturm+0x596>
movsd  0x0(%rbp,%rdi,8),%xmm1
mulsd  %xmm6,%xmm1
addsd  (%rbx,%rdi,8),%xmm1
movsd  %xmm1,(%rbx,%rdi,8)
lea    0x1(%rdi),%rcx
not    %rdi
cmp    %r8,%rdi
je     1550 <sturm+0x400>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%rbp,%rcx,8),%xmm1
mulsd  %xmm6,%xmm1
addsd  (%rbx,%rcx,8),%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
movsd  0x8(%rbp,%rcx,8),%xmm1
mulsd  %xmm6,%xmm1
addsd  0x8(%rbx,%rcx,8),%xmm1
movsd  %xmm1,0x8(%rbx,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r13
jne    1700 <sturm+0x5b0>
jmp    1550 <sturm+0x400>
mov    0x30(%rsp),%rdi
mov    $0x1,%esi
call   1060 <free_dvector@plt>
mov    0xc(%rsp),%eax
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@stdout = global i32 0
@.str = [34 x i8] c" There are %d modes below %f Hz.\00"
@stderr = global i32 0
@.str.1 = [31 x i8] c" ... %d modes were not found.\0A\00"
@.str.2 = [41 x i8] c" Try increasing the number of modes in \0A\00"
@.str.3 = [46 x i8] c" order to get the missing modes below %f Hz.\0A\00"
@.str.4 = [28 x i8] c" All %d modes were found.\0A\00"
define i32 @sturm(double** %0, double** %1, i32 %2, i32 %3, double %4, double %5, i32 %6) {
%8 = alloca double**
%9 = alloca double**
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca i32
%15 = alloca double
%16 = alloca double*
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
store double** %0, double*** %8
store double** %1, double*** %9
store i32 %2, i32* %10
store i32 %3, i32* %11
store double %4, double* %12
store double %5, double* %13
store i32 %6, i32* %14
store i32 0, i32* %17
%21 = load i32, i32* %10
%22 = call double* @dvector(i32 1, i32 %21)
store double* %22, double** %16
%23 = load i32, i32* %11
%24 = sitofp i32 %23 to double
%25 = fmul double 5.000000e-01, %24
%26 = fptrunc double %25 to float
%27 = load i32, i32* %11
%28 = sitofp i32 %27 to double
%29 = fsub double %28, 8.000000e+00
%30 = fptrunc double %29 to float
%31 = fcmp ogt float %26, %30
br i1 %31, label %32, label %37
32: ; preds = %7
%33 = load i32, i32* %11
%34 = sitofp i32 %33 to double
%35 = fdiv double %34, 2.000000e+00
%36 = fptosi double %35 to i32
br label %40
37: ; preds = %7
%38 = load i32, i32* %11
%39 = sub nsw i32 %38, 8
br label %40
40: ; preds = %37, %32
%41 = phi i32 [ %36, %32 ], [ %39, %37 ]
store i32 %41, i32* %20
%42 = load double, double* %13
%43 = load double, double* %12
%44 = fadd double %42, %43
store double %44, double* %15
store i32 1, i32* %18
br label %45
45: ; preds = %81, %40
%46 = load i32, i32* %18
%47 = load i32, i32* %10
%48 = icmp sle i32 %46, %47
br i1 %48, label %49, label %84
49: ; preds = %45
%50 = load i32, i32* %18
store i32 %50, i32* %19
br label %51
51: ; preds = %77, %49
%52 = load i32, i32* %19
%53 = load i32, i32* %10
%54 = icmp sle i32 %52, %53
br i1 %54, label %55, label %80
55: ; preds = %51
%56 = load double, double* %15
%57 = load double**, double*** %9
%58 = load i32, i32* %18
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double*, double** %57, i64 %59
%61 = load double*, double** %60
%62 = load i32, i32* %19
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %61, i64 %63
%65 = load double, double* %64
%66 = load double**, double*** %8
%67 = load i32, i32* %18
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double*, double** %66, i64 %68
%70 = load double*, double** %69
%71 = load i32, i32* %19
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %70, i64 %72
%74 = load double, double* %73
%75 = fneg double %56
%76 = call double @llvm.fmuladd.f64(double %75, double %65, double %74)
store double %76, double* %73
br label %77
77: ; preds = %55
%78 = load i32, i32* %19
%79 = add nsw i32 %78, 1
store i32 %79, i32* %19
br label %51
80: ; preds = %51
br label %81
81: ; preds = %80
%82 = load i32, i32* %18
%83 = add nsw i32 %82, 1
store i32 %83, i32* %18
br label %45
84: ; preds = %45
%85 = load double**, double*** %8
%86 = load i32, i32* %10
%87 = load double*, double** %16
%88 = load double*, double** %16
%89 = load double*, double** %16
%90 = call i32 @ldl_dcmp(double** %85, i32 %86, double* %87, double* %88, double* %89, i32 1, i32 0, i32* %17)
%91 = load i32, i32* %14
%92 = icmp ne i32 %91, 0
br i1 %92, label %93, label %101
93: ; preds = %84
%94 = load i32, i32* @stdout
%95 = load i32, i32* %17
%96 = sub nsw i32 0, %95
%97 = load double, double* %13
%98 = call double @sqrt(double %97)
%99 = fdiv double %98, 0x401921FB54442D18
%100 = call i32 (i32, i8*, ...) @fprintf(i32 %94, i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str, i64 0, i64 0), i32 %96, double %99)
br label %101
101: ; preds = %93, %84
%102 = load i32, i32* %17
%103 = sub nsw i32 0, %102
%104 = load i32, i32* %20
%105 = icmp sgt i32 %103, %104
br i1 %105, label %106, label %120
106: ; preds = %101
%107 = load i32, i32* @stderr
%108 = load i32, i32* %17
%109 = sub nsw i32 0, %108
%110 = load i32, i32* %20
%111 = sub nsw i32 %109, %110
%112 = call i32 (i32, i8*, ...) @fprintf(i32 %107, i8* getelementptr inbounds ([31 x i8], [31 x i8]* @.str.1, i64 0, i64 0), i32 %111)
%113 = load i32, i32* @stderr
%114 = call i32 (i32, i8*, ...) @fprintf(i32 %113, i8* getelementptr inbounds ([41 x i8], [41 x i8]* @.str.2, i64 0, i64 0))
%115 = load i32, i32* @stderr
%116 = load double, double* %13
%117 = call double @sqrt(double %116)
%118 = fdiv double %117, 0x401921FB54442D18
%119 = call i32 (i32, i8*, ...) @fprintf(i32 %115, i8* getelementptr inbounds ([46 x i8], [46 x i8]* @.str.3, i64 0, i64 0), double %118)
br label %128
120: ; preds = %101
%121 = load i32, i32* %14
%122 = icmp ne i32 %121, 0
br i1 %122, label %123, label %127
123: ; preds = %120
%124 = load i32, i32* @stdout
%125 = load i32, i32* %20
%126 = call i32 (i32, i8*, ...) @fprintf(i32 %124, i8* getelementptr inbounds ([28 x i8], [28 x i8]* @.str.4, i64 0, i64 0), i32 %125)
br label %127
127: ; preds = %123, %120
br label %128
128: ; preds = %127, %106
store i32 1, i32* %18
br label %129
129: ; preds = %164, %128
%130 = load i32, i32* %18
%131 = load i32, i32* %10
%132 = icmp sle i32 %130, %131
br i1 %132, label %133, label %167
133: ; preds = %129
%134 = load i32, i32* %18
store i32 %134, i32* %19
br label %135
135: ; preds = %160, %133
%136 = load i32, i32* %19
%137 = load i32, i32* %10
%138 = icmp sle i32 %136, %137
br i1 %138, label %139, label %163
139: ; preds = %135
%140 = load double, double* %15
%141 = load double**, double*** %9
%142 = load i32, i32* %18
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double*, double** %141, i64 %143
%145 = load double*, double** %144
%146 = load i32, i32* %19
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = load double**, double*** %8
%151 = load i32, i32* %18
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double*, double** %150, i64 %152
%154 = load double*, double** %153
%155 = load i32, i32* %19
%156 = sext i32 %155 to i64
%157 = getelementptr inbounds double, double* %154, i64 %156
%158 = load double, double* %157
%159 = call double @llvm.fmuladd.f64(double %140, double %149, double %158)
store double %159, double* %157
br label %160
160: ; preds = %139
%161 = load i32, i32* %19
%162 = add nsw i32 %161, 1
store i32 %162, i32* %19
br label %135
163: ; preds = %135
br label %164
164: ; preds = %163
%165 = load i32, i32* %18
%166 = add nsw i32 %165, 1
store i32 %166, i32* %18
br label %129
167: ; preds = %129
%168 = load double*, double** %16
%169 = load i32, i32* %10
%170 = call i32 @free_dvector(double* %168, i32 1, i32 %169)
%171 = load i32, i32* %17
ret i32 %171
}

LLVM IR (-O2):
@stdout = global i32 0
@.str = [34 x i8] c" There are %d modes below %f Hz.\00"
@stderr = global i32 0
@.str.1 = [31 x i8] c" ... %d modes were not found.\0A\00"
@.str.2 = [41 x i8] c" Try increasing the number of modes in \0A\00"
@.str.3 = [46 x i8] c" order to get the missing modes below %f Hz.\0A\00"
@.str.4 = [28 x i8] c" All %d modes were found.\0A\00"
define i32 @sturm(double** %0, double** nocapture %1, i32 %2, i32 %3, double %4, double %5, i32 %6) {
%8 = alloca i32
%9 = bitcast i32* %8 to i8*
store i32 0, i32* %8
%10 = tail call double* @dvector(i32 1, i32 %2)
%11 = sitofp i32 %3 to double
%12 = fmul double %11, 5.000000e-01
%13 = fptrunc double %12 to float
%14 = fadd double %11, -8.000000e+00
%15 = fptrunc double %14 to float
%16 = fcmp ogt float %13, %15
%17 = fptosi double %12 to i32
%18 = add nsw i32 %3, -8
%19 = select i1 %16, i32 %17, i32 %18
%20 = fadd double %4, %5
%21 = icmp slt i32 %2, 1
br i1 %21, label %160, label %22
22: ; preds = %7
%23 = fneg double %20
%24 = add nuw i32 %2, 1
%25 = zext i32 %24 to i64
%26 = insertelement <2 x double> poison, double %23, i64 0
%27 = shufflevector <2 x double> %26, <2 x double> poison, <2 x i32> zeroinitializer
%28 = insertelement <2 x double> poison, double %23, i64 0
%29 = shufflevector <2 x double> %28, <2 x double> poison, <2 x i32> zeroinitializer
%30 = sub nsw i64 0, %25
br label %31
31: ; preds = %156, %22
%32 = phi i64 [ %159, %156 ], [ 0, %22 ]
%33 = phi i64 [ %157, %156 ], [ 1, %22 ]
%34 = xor i64 %32, -1
%35 = add i64 %34, %25
%36 = add i64 %35, -4
%37 = lshr i64 %36, 2
%38 = add nuw nsw i64 %37, 1
%39 = xor i64 %32, -1
%40 = add i64 %39, %25
%41 = getelementptr inbounds double*, double** %1, i64 %33
%42 = load double*, double** %41
%43 = getelementptr inbounds double*, double** %0, i64 %33
%44 = load double*, double** %43
%45 = icmp ult i64 %40, 4
br i1 %45, label %125, label %46
46: ; preds = %31
%47 = getelementptr double, double* %44, i64 %33
%48 = getelementptr double, double* %44, i64 %25
%49 = getelementptr double, double* %42, i64 %33
%50 = getelementptr double, double* %42, i64 %25
%51 = icmp ult double* %47, %50
%52 = icmp ult double* %49, %48
%53 = and i1 %51, %52
br i1 %53, label %125, label %54
54: ; preds = %46
%55 = and i64 %40, -4
%56 = add i64 %33, %55
%57 = and i64 %38, 1
%58 = icmp ult i64 %36, 4
br i1 %58, label %102, label %59
59: ; preds = %54
%60 = and i64 %38, 9223372036854775806
br label %61
61: ; preds = %61, %59
%62 = phi i64 [ 0, %59 ], [ %99, %61 ]
%63 = phi i64 [ 0, %59 ], [ %100, %61 ]
%64 = add i64 %33, %62
%65 = getelementptr inbounds double, double* %42, i64 %64
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = getelementptr inbounds double, double* %65, i64 2
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = getelementptr inbounds double, double* %44, i64 %64
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = getelementptr inbounds double, double* %71, i64 2
%75 = bitcast double* %74 to <2 x double>*
%76 = load <2 x double>, <2 x double>* %75
%77 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %27, <2 x double> %67, <2 x double> %73)
%78 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %29, <2 x double> %70, <2 x double> %76)
%79 = bitcast double* %71 to <2 x double>*
store <2 x double> %77, <2 x double>* %79
%80 = bitcast double* %74 to <2 x double>*
store <2 x double> %78, <2 x double>* %80
%81 = or i64 %62, 4
%82 = add i64 %33, %81
%83 = getelementptr inbounds double, double* %42, i64 %82
%84 = bitcast double* %83 to <2 x double>*
%85 = load <2 x double>, <2 x double>* %84
%86 = getelementptr inbounds double, double* %83, i64 2
%87 = bitcast double* %86 to <2 x double>*
%88 = load <2 x double>, <2 x double>* %87
%89 = getelementptr inbounds double, double* %44, i64 %82
%90 = bitcast double* %89 to <2 x double>*
%91 = load <2 x double>, <2 x double>* %90
%92 = getelementptr inbounds double, double* %89, i64 2
%93 = bitcast double* %92 to <2 x double>*
%94 = load <2 x double>, <2 x double>* %93
%95 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %27, <2 x double> %85, <2 x double> %91)
%96 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %29, <2 x double> %88, <2 x double> %94)
%97 = bitcast double* %89 to <2 x double>*
store <2 x double> %95, <2 x double>* %97
%98 = bitcast double* %92 to <2 x double>*
store <2 x double> %96, <2 x double>* %98
%99 = add nuw i64 %62, 8
%100 = add i64 %63, 2
%101 = icmp eq i64 %100, %60
br i1 %101, label %102, label %61
102: ; preds = %61, %54
%103 = phi i64 [ 0, %54 ], [ %99, %61 ]
%104 = icmp eq i64 %57, 0
br i1 %104, label %123, label %105
105: ; preds = %102
%106 = add i64 %33, %103
%107 = getelementptr inbounds double, double* %42, i64 %106
%108 = bitcast double* %107 to <2 x double>*
%109 = load <2 x double>, <2 x double>* %108
%110 = getelementptr inbounds double, double* %107, i64 2
%111 = bitcast double* %110 to <2 x double>*
%112 = load <2 x double>, <2 x double>* %111
%113 = getelementptr inbounds double, double* %44, i64 %106
%114 = bitcast double* %113 to <2 x double>*
%115 = load <2 x double>, <2 x double>* %114
%116 = getelementptr inbounds double, double* %113, i64 2
%117 = bitcast double* %116 to <2 x double>*
%118 = load <2 x double>, <2 x double>* %117
%119 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %27, <2 x double> %109, <2 x double> %115)
%120 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %29, <2 x double> %112, <2 x double> %118)
%121 = bitcast double* %113 to <2 x double>*
store <2 x double> %119, <2 x double>* %121
%122 = bitcast double* %116 to <2 x double>*
store <2 x double> %120, <2 x double>* %122
br label %123
123: ; preds = %102, %105
%124 = icmp eq i64 %40, %55
br i1 %124, label %156, label %125
125: ; preds = %46, %31, %123
%126 = phi i64 [ %33, %46 ], [ %33, %31 ], [ %56, %123 ]
%127 = sub i64 %25, %126
%128 = xor i64 %126, -1
%129 = and i64 %127, 1
%130 = icmp eq i64 %129, 0
br i1 %130, label %138, label %131
131: ; preds = %125
%132 = getelementptr inbounds double, double* %42, i64 %126
%133 = load double, double* %132
%134 = getelementptr inbounds double, double* %44, i64 %126
%135 = load double, double* %134
%136 = tail call double @llvm.fmuladd.f64(double %23, double %133, double %135)
store double %136, double* %134
%137 = add nuw nsw i64 %126, 1
br label %138
138: ; preds = %131, %125
%139 = phi i64 [ %126, %125 ], [ %137, %131 ]
%140 = icmp eq i64 %128, %30
br i1 %140, label %156, label %141
141: ; preds = %138, %141
%142 = phi i64 [ %154, %141 ], [ %139, %138 ]
%143 = getelementptr inbounds double, double* %42, i64 %142
%144 = load double, double* %143
%145 = getelementptr inbounds double, double* %44, i64 %142
%146 = load double, double* %145
%147 = tail call double @llvm.fmuladd.f64(double %23, double %144, double %146)
store double %147, double* %145
%148 = add nuw nsw i64 %142, 1
%149 = getelementptr inbounds double, double* %42, i64 %148
%150 = load double, double* %149
%151 = getelementptr inbounds double, double* %44, i64 %148
%152 = load double, double* %151
%153 = tail call double @llvm.fmuladd.f64(double %23, double %150, double %152)
store double %153, double* %151
%154 = add nuw nsw i64 %142, 2
%155 = icmp eq i64 %154, %25
br i1 %155, label %156, label %141
156: ; preds = %138, %141, %123
%157 = add nuw nsw i64 %33, 1
%158 = icmp eq i64 %157, %25
%159 = add i64 %32, 1
br i1 %158, label %160, label %31
160: ; preds = %156, %7
%161 = call i32 @ldl_dcmp(double** %0, i32 %2, double* %10, double* %10, double* %10, i32 1, i32 0, i32* %8)
%162 = icmp eq i32 %6, 0
br i1 %162, label %163, label %167
163: ; preds = %160
%164 = load i32, i32* %8
%165 = sub nsw i32 0, %164
%166 = icmp slt i32 %19, %165
br i1 %166, label %177, label %191
167: ; preds = %160
%168 = load i32, i32* @stdout
%169 = load i32, i32* %8
%170 = sub nsw i32 0, %169
%171 = call double @sqrt(double %5)
%172 = fdiv double %171, 0x401921FB54442D18
%173 = call i32 (i32, i8*, ...) @fprintf(i32 %168, i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str, i64 0, i64 0), i32 %170, double %172)
%174 = load i32, i32* %8
%175 = sub nsw i32 0, %174
%176 = icmp slt i32 %19, %175
br i1 %176, label %177, label %188
177: ; preds = %167, %163
%178 = phi i32 [ %175, %167 ], [ %165, %163 ]
%179 = load i32, i32* @stderr
%180 = sub nsw i32 %178, %19
%181 = call i32 (i32, i8*, ...) @fprintf(i32 %179, i8* getelementptr inbounds ([31 x i8], [31 x i8]* @.str.1, i64 0, i64 0), i32 %180)
%182 = load i32, i32* @stderr
%183 = call i32 (i32, i8*, ...) @fprintf(i32 %182, i8* getelementptr inbounds ([41 x i8], [41 x i8]* @.str.2, i64 0, i64 0))
%184 = load i32, i32* @stderr
%185 = call double @sqrt(double %5)
%186 = fdiv double %185, 0x401921FB54442D18
%187 = call i32 (i32, i8*, ...) @fprintf(i32 %184, i8* getelementptr inbounds ([46 x i8], [46 x i8]* @.str.3, i64 0, i64 0), double %186)
br label %191
188: ; preds = %167
%189 = load i32, i32* @stdout
%190 = call i32 (i32, i8*, ...) @fprintf(i32 %189, i8* getelementptr inbounds ([28 x i8], [28 x i8]* @.str.4, i64 0, i64 0), i32 %19)
br label %191
191: ; preds = %163, %188, %177
br i1 %21, label %329, label %192
192: ; preds = %191
%193 = add nuw i32 %2, 1
%194 = zext i32 %193 to i64
%195 = insertelement <2 x double> poison, double %20, i64 0
%196 = shufflevector <2 x double> %195, <2 x double> poison, <2 x i32> zeroinitializer
%197 = insertelement <2 x double> poison, double %20, i64 0
%198 = shufflevector <2 x double> %197, <2 x double> poison, <2 x i32> zeroinitializer
%199 = sub nsw i64 0, %194
br label %200
200: ; preds = %325, %192
%201 = phi i64 [ %328, %325 ], [ 0, %192 ]
%202 = phi i64 [ %326, %325 ], [ 1, %192 ]
%203 = xor i64 %201, -1
%204 = add i64 %203, %194
%205 = add i64 %204, -4
%206 = lshr i64 %205, 2
%207 = add nuw nsw i64 %206, 1
%208 = xor i64 %201, -1
%209 = add i64 %208, %194
%210 = getelementptr inbounds double*, double** %1, i64 %202
%211 = load double*, double** %210
%212 = getelementptr inbounds double*, double** %0, i64 %202
%213 = load double*, double** %212
%214 = icmp ult i64 %209, 4
br i1 %214, label %294, label %215
215: ; preds = %200
%216 = getelementptr double, double* %213, i64 %202
%217 = getelementptr double, double* %213, i64 %194
%218 = getelementptr double, double* %211, i64 %202
%219 = getelementptr double, double* %211, i64 %194
%220 = icmp ult double* %216, %219
%221 = icmp ult double* %218, %217
%222 = and i1 %220, %221
br i1 %222, label %294, label %223
223: ; preds = %215
%224 = and i64 %209, -4
%225 = add i64 %202, %224
%226 = and i64 %207, 1
%227 = icmp ult i64 %205, 4
br i1 %227, label %271, label %228
228: ; preds = %223
%229 = and i64 %207, 9223372036854775806
br label %230
230: ; preds = %230, %228
%231 = phi i64 [ 0, %228 ], [ %268, %230 ]
%232 = phi i64 [ 0, %228 ], [ %269, %230 ]
%233 = add i64 %202, %231
%234 = getelementptr inbounds double, double* %211, i64 %233
%235 = bitcast double* %234 to <2 x double>*
%236 = load <2 x double>, <2 x double>* %235
%237 = getelementptr inbounds double, double* %234, i64 2
%238 = bitcast double* %237 to <2 x double>*
%239 = load <2 x double>, <2 x double>* %238
%240 = getelementptr inbounds double, double* %213, i64 %233
%241 = bitcast double* %240 to <2 x double>*
%242 = load <2 x double>, <2 x double>* %241
%243 = getelementptr inbounds double, double* %240, i64 2
%244 = bitcast double* %243 to <2 x double>*
%245 = load <2 x double>, <2 x double>* %244
%246 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %196, <2 x double> %236, <2 x double> %242)
%247 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %198, <2 x double> %239, <2 x double> %245)
%248 = bitcast double* %240 to <2 x double>*
store <2 x double> %246, <2 x double>* %248
%249 = bitcast double* %243 to <2 x double>*
store <2 x double> %247, <2 x double>* %249
%250 = or i64 %231, 4
%251 = add i64 %202, %250
%252 = getelementptr inbounds double, double* %211, i64 %251
%253 = bitcast double* %252 to <2 x double>*
%254 = load <2 x double>, <2 x double>* %253
%255 = getelementptr inbounds double, double* %252, i64 2
%256 = bitcast double* %255 to <2 x double>*
%257 = load <2 x double>, <2 x double>* %256
%258 = getelementptr inbounds double, double* %213, i64 %251
%259 = bitcast double* %258 to <2 x double>*
%260 = load <2 x double>, <2 x double>* %259
%261 = getelementptr inbounds double, double* %258, i64 2
%262 = bitcast double* %261 to <2 x double>*
%263 = load <2 x double>, <2 x double>* %262
%264 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %196, <2 x double> %254, <2 x double> %260)
%265 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %198, <2 x double> %257, <2 x double> %263)
%266 = bitcast double* %258 to <2 x double>*
store <2 x double> %264, <2 x double>* %266
%267 = bitcast double* %261 to <2 x double>*
store <2 x double> %265, <2 x double>* %267
%268 = add nuw i64 %231, 8
%269 = add i64 %232, 2
%270 = icmp eq i64 %269, %229
br i1 %270, label %271, label %230
271: ; preds = %230, %223
%272 = phi i64 [ 0, %223 ], [ %268, %230 ]
%273 = icmp eq i64 %226, 0
br i1 %273, label %292, label %274
274: ; preds = %271
%275 = add i64 %202, %272
%276 = getelementptr inbounds double, double* %211, i64 %275
%277 = bitcast double* %276 to <2 x double>*
%278 = load <2 x double>, <2 x double>* %277
%279 = getelementptr inbounds double, double* %276, i64 2
%280 = bitcast double* %279 to <2 x double>*
%281 = load <2 x double>, <2 x double>* %280
%282 = getelementptr inbounds double, double* %213, i64 %275
%283 = bitcast double* %282 to <2 x double>*
%284 = load <2 x double>, <2 x double>* %283
%285 = getelementptr inbounds double, double* %282, i64 2
%286 = bitcast double* %285 to <2 x double>*
%287 = load <2 x double>, <2 x double>* %286
%288 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %196, <2 x double> %278, <2 x double> %284)
%289 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %198, <2 x double> %281, <2 x double> %287)
%290 = bitcast double* %282 to <2 x double>*
store <2 x double> %288, <2 x double>* %290
%291 = bitcast double* %285 to <2 x double>*
store <2 x double> %289, <2 x double>* %291
br label %292
292: ; preds = %271, %274
%293 = icmp eq i64 %209, %224
br i1 %293, label %325, label %294
294: ; preds = %215, %200, %292
%295 = phi i64 [ %202, %215 ], [ %202, %200 ], [ %225, %292 ]
%296 = sub i64 %194, %295
%297 = xor i64 %295, -1
%298 = and i64 %296, 1
%299 = icmp eq i64 %298, 0
br i1 %299, label %307, label %300
300: ; preds = %294
%301 = getelementptr inbounds double, double* %211, i64 %295
%302 = load double, double* %301
%303 = getelementptr inbounds double, double* %213, i64 %295
%304 = load double, double* %303
%305 = call double @llvm.fmuladd.f64(double %20, double %302, double %304)
store double %305, double* %303
%306 = add nuw nsw i64 %295, 1
br label %307
307: ; preds = %300, %294
%308 = phi i64 [ %295, %294 ], [ %306, %300 ]
%309 = icmp eq i64 %297, %199
br i1 %309, label %325, label %310
310: ; preds = %307, %310
%311 = phi i64 [ %323, %310 ], [ %308, %307 ]
%312 = getelementptr inbounds double, double* %211, i64 %311
%313 = load double, double* %312
%314 = getelementptr inbounds double, double* %213, i64 %311
%315 = load double, double* %314
%316 = call double @llvm.fmuladd.f64(double %20, double %313, double %315)
store double %316, double* %314
%317 = add nuw nsw i64 %311, 1
%318 = getelementptr inbounds double, double* %211, i64 %317
%319 = load double, double* %318
%320 = getelementptr inbounds double, double* %213, i64 %317
%321 = load double, double* %320
%322 = call double @llvm.fmuladd.f64(double %20, double %319, double %321)
store double %322, double* %320
%323 = add nuw nsw i64 %311, 2
%324 = icmp eq i64 %323, %194
br i1 %324, label %325, label %310
325: ; preds = %307, %310, %292
%326 = add nuw nsw i64 %202, 1
%327 = icmp eq i64 %326, %194
%328 = add i64 %201, 1
br i1 %327, label %329, label %200
329: ; preds = %325, %191
%330 = call i32 @free_dvector(double* %10, i32 1, i32 %2)
%331 = load i32, i32* %8
ret i32 %331
}