#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 */
int /*<<< orphan*/ epipole(double *, double const *);
int sqrt(double);
int exFDso(const double *u, const double *F, double *p, double *w, int len,
           double th, int *inl_buff, int **inls) {
  double rx, ry, rwc, ryc, rxc, r;
  int i, po = 0, ne = 0;
  double ec[3];
  double sx, sy, sgn;
  epipole(ec, F);
  for (i = 1; i <= len; i++) {
    rxc = (*F) * (*(u + 3)) + (*(F + 3)) * (*(u + 4)) + (*(F + 6)) * (*(u + 5));
    ryc = (*(F + 1)) * (*(u + 3)) + (*(F + 4)) * (*(u + 4)) +
          (*(F + 7)) * (*(u + 5));
    rwc = (*(F + 2)) * (*(u + 3)) + (*(F + 5)) * (*(u + 4)) +
          (*(F + 8)) * (*(u + 5));
    r = ((*(u)) * rxc + (*(u + 1)) * ryc + (*(u + 2)) * rwc);
    rx = (*F) * (*(u)) + (*(F + 1)) * (*(u + 1)) + (*(F + 2)) * (*(u + 2));
    ry =
        (*(F + 3)) * (*(u)) + (*(F + 4)) * (*(u + 1)) + (*(F + 5)) * (*(u + 2));
    *w = rxc * rxc + ryc * ryc + rx * rx + ry * ry;
    *p = r * r / *w;
    if (*p < th) {
      sx = ec[1] * u[2] - ec[2] * u[1];
      sy = ec[2] * u[0] - ec[0] * u[2];
      sgn = sx * rxc + sy * ryc;
      if (sgn > 0) {
        inl_buff[po] = i;
        po++;
      } else {
        ne++;
        inl_buff[len - ne] = i;
      }
      *w = 1 / sqrt(*w);
    }
    p++;
    w++;
    u += 6;
  }
  if (po >= ne) {
    *inls = inl_buff;
    return po;
  }
  *inls = inl_buff + len - ne;
  return ne;
}


x86 assembly (clang, -O0):
<exFDso>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
movsd  %xmm0,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
movl   $0x0,-0x78(%rbp)
movl   $0x0,-0x7c(%rbp)
lea    -0xa0(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   1030 <epipole@plt>
movl   $0x1,-0x74(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1504 <exFDso+0x3e4>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x20(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x20(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x38(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x20(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x40(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x68(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mulsd  -0x60(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  -0x58(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x48(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x70(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
mov    -0x28(%rbp),%rax
divsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    14d2 <exFDso+0x3b2>
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mov    -0x10(%rbp),%rax
mulsd  0x8(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0xa0(%rbp),%xmm1
ff
mov    -0x10(%rbp),%rax
mulsd  0x10(%rax),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0x68(%rbp),%xmm2
movsd  -0xb0(%rbp),%xmm1
ff
mulsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1493 <exFDso+0x373>
mov    -0x74(%rbp),%edx
mov    -0x40(%rbp),%rax
movslq -0x78(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x78(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
jmp    14af <exFDso+0x38f>
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
mov    -0x74(%rbp),%edx
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x7c(%rbp),%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
call   1040 <sqrt@plt>
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
add    $0x30,%rax
mov    %rax,-0x10(%rbp)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    1171 <exFDso+0x51>
mov    -0x78(%rbp),%eax
cmp    -0x7c(%rbp),%eax
jl     1526 <exFDso+0x406>
mov    -0x40(%rbp),%rcx
mov    0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x78(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1552 <exFDso+0x432>
mov    -0x40(%rbp),%rcx
movslq -0x2c(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rcx
movslq -0x7c(%rbp),%rdx
xor    %eax,%eax
sub    %rdx,%rax
shl    $0x2,%rax
add    %rax,%rcx
mov    0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x7c(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<exFDso>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9,0x10(%rsp)
movsd  %xmm0,0x28(%rsp)
mov    %r8d,%r15d
mov    %rcx,%rbp
mov    %rdx,0x20(%rsp)
mov    %rsi,%rbx
mov    %rdi,%r14
lea    0x30(%rsp),%rdi
call   1030 <epipole@plt>
test   %r15d,%r15d
mov    %r15d,0xc(%rsp)
jle    1378 <exFDso+0x258>
mov    %r15d,%ecx
xor    %r12d,%r12d
movl   $0x0,0x8(%rsp)
xor    %r13d,%r13d
mov    %rcx,0x18(%rsp)
jmp    11d9 <exFDso+0xb9>
nopw   0x0(%rax,%rax,1)
mov    0x8(%rsp),%ecx
add    $0x1,%ecx
mov    0xc(%rsp),%eax
mov    %ecx,0x8(%rsp)
sub    %ecx,%eax
cltq
mov    0x10(%rsp),%rcx
mov    %r15d,(%rcx,%rax,4)
movsd  0x0(%rbp,%r12,8),%xmm0
call   1040 <sqrt@plt>
lea    0x1(%rax),%ecx
cmp    $0x3,%ecx
mov    $0x0,%ecx
cmovae %ecx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x0(%rbp,%r12,8)
mov    0x18(%rsp),%rcx
add    $0x30,%r14
mov    %r15,%r12
cmp    %r15d,%ecx
je     1383 <exFDso+0x263>
lea    0x1(%r12),%r15
movsd  0x18(%r14),%xmm0
movsd  0x18(%rbx),%xmm11
movsd  0x20(%r14),%xmm5
movapd %xmm11,%xmm1
mulsd  %xmm5,%xmm1
movsd  (%rbx),%xmm4
movsd  0x8(%rbx),%xmm8
movapd %xmm4,%xmm6
mulsd  %xmm0,%xmm6
addsd  %xmm1,%xmm6
movsd  0x28(%r14),%xmm2
movsd  0x30(%rbx),%xmm9
mulsd  %xmm2,%xmm9
addsd  %xmm6,%xmm9
movsd  0x20(%rbx),%xmm10
movapd %xmm5,%xmm1
mulsd  %xmm10,%xmm1
movapd %xmm8,%xmm6
mulsd  %xmm0,%xmm6
addsd  %xmm1,%xmm6
movsd  0x38(%rbx),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm6,%xmm3
movsd  0x10(%rbx),%xmm7
movsd  0x28(%rbx),%xmm6
mulsd  %xmm6,%xmm5
mulsd  %xmm7,%xmm0
addsd  %xmm5,%xmm0
mulsd  0x40(%rbx),%xmm2
addsd  %xmm0,%xmm2
movsd  (%r14),%xmm0
movsd  0x8(%r14),%xmm5
movapd %xmm3,%xmm1
mulsd  %xmm5,%xmm1
mulsd  %xmm0,%xmm4
mulsd  %xmm0,%xmm11
mulsd  %xmm9,%xmm0
addsd  %xmm1,%xmm0
movsd  0x10(%r14),%xmm1
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm5,%xmm8
addsd  %xmm8,%xmm4
mulsd  %xmm1,%xmm7
addsd  %xmm4,%xmm7
mulsd  %xmm5,%xmm10
addsd  %xmm10,%xmm11
mulsd  %xmm1,%xmm6
addsd  %xmm11,%xmm6
movapd %xmm3,%xmm0
mulsd  %xmm3,%xmm0
movapd %xmm9,%xmm1
mulsd  %xmm9,%xmm1
addsd  %xmm0,%xmm1
mulsd  %xmm7,%xmm7
addsd  %xmm1,%xmm7
mulsd  %xmm6,%xmm6
addsd  %xmm7,%xmm6
movsd  %xmm6,0x0(%rbp,%r12,8)
mulsd  %xmm2,%xmm2
divsd  %xmm6,%xmm2
mov    0x20(%rsp),%rax
movsd  %xmm2,(%rax,%r12,8)
movsd  0x28(%rsp),%xmm0
ucomisd %xmm2,%xmm0
jbe    11c9 <exFDso+0xa9>
movsd  0x10(%r14),%xmm0
movsd  0x40(%rsp),%xmm1
movsd  0x8(%r14),%xmm2
mulsd  %xmm1,%xmm2
movsd  0x38(%rsp),%xmm4
mulsd  %xmm0,%xmm4
subsd  %xmm2,%xmm4
mulsd  0x30(%rsp),%xmm0
mulsd  (%r14),%xmm1
subsd  %xmm0,%xmm1
mulsd  %xmm1,%xmm3
mulsd  %xmm9,%xmm4
addsd  %xmm3,%xmm4
ucomisd 0xca3(%rip),%xmm4
jbe    1180 <exFDso+0x60>
movslq %r13d,%r13
mov    0x10(%rsp),%rax
mov    %r15d,(%rax,%r13,4)
add    $0x1,%r13d
jmp    119c <exFDso+0x7c>
xor    %r13d,%r13d
movl   $0x0,0x8(%rsp)
movslq 0xc(%rsp),%rax
mov    0x10(%rsp),%rsi
lea    (%rsi,%rax,4),%rax
movslq 0x8(%rsp),%rcx
lea    0x0(,%rcx,4),%rdx
sub    %rdx,%rax
cmp    %ecx,%r13d
cmovge %rsi,%rax
cmovl  %ecx,%r13d
mov    0x80(%rsp),%rcx
mov    %rax,(%rcx)
mov    %r13d,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @exFDso(double* %0, double* %1, double* %2, double* %3, i32 %4, double %5, i32* %6, i32** %7) {
%9 = alloca i32
%10 = alloca double*
%11 = alloca double*
%12 = alloca double*
%13 = alloca double*
%14 = alloca i32
%15 = alloca double
%16 = alloca i32*
%17 = alloca i32**
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca [3 x double]
%28 = alloca double
%29 = alloca double
%30 = alloca double
store double* %0, double** %10
store double* %1, double** %11
store double* %2, double** %12
store double* %3, double** %13
store i32 %4, i32* %14
store double %5, double* %15
store i32* %6, i32** %16
store i32** %7, i32*** %17
store i32 0, i32* %25
store i32 0, i32* %26
%31 = getelementptr inbounds [3 x double], [3 x double]* %27, i64 0, i64 0
%32 = load double*, double** %11
%33 = call i32 @epipole(double* %31, double* %32)
store i32 1, i32* %24
br label %34
34: ; preds = %245, %8
%35 = load i32, i32* %24
%36 = load i32, i32* %14
%37 = icmp sle i32 %35, %36
br i1 %37, label %38, label %248
38: ; preds = %34
%39 = load double*, double** %11
%40 = load double, double* %39
%41 = load double*, double** %10
%42 = getelementptr inbounds double, double* %41, i64 3
%43 = load double, double* %42
%44 = load double*, double** %11
%45 = getelementptr inbounds double, double* %44, i64 3
%46 = load double, double* %45
%47 = load double*, double** %10
%48 = getelementptr inbounds double, double* %47, i64 4
%49 = load double, double* %48
%50 = fmul double %46, %49
%51 = call double @llvm.fmuladd.f64(double %40, double %43, double %50)
%52 = load double*, double** %11
%53 = getelementptr inbounds double, double* %52, i64 6
%54 = load double, double* %53
%55 = load double*, double** %10
%56 = getelementptr inbounds double, double* %55, i64 5
%57 = load double, double* %56
%58 = call double @llvm.fmuladd.f64(double %54, double %57, double %51)
store double %58, double* %22
%59 = load double*, double** %11
%60 = getelementptr inbounds double, double* %59, i64 1
%61 = load double, double* %60
%62 = load double*, double** %10
%63 = getelementptr inbounds double, double* %62, i64 3
%64 = load double, double* %63
%65 = load double*, double** %11
%66 = getelementptr inbounds double, double* %65, i64 4
%67 = load double, double* %66
%68 = load double*, double** %10
%69 = getelementptr inbounds double, double* %68, i64 4
%70 = load double, double* %69
%71 = fmul double %67, %70
%72 = call double @llvm.fmuladd.f64(double %61, double %64, double %71)
%73 = load double*, double** %11
%74 = getelementptr inbounds double, double* %73, i64 7
%75 = load double, double* %74
%76 = load double*, double** %10
%77 = getelementptr inbounds double, double* %76, i64 5
%78 = load double, double* %77
%79 = call double @llvm.fmuladd.f64(double %75, double %78, double %72)
store double %79, double* %21
%80 = load double*, double** %11
%81 = getelementptr inbounds double, double* %80, i64 2
%82 = load double, double* %81
%83 = load double*, double** %10
%84 = getelementptr inbounds double, double* %83, i64 3
%85 = load double, double* %84
%86 = load double*, double** %11
%87 = getelementptr inbounds double, double* %86, i64 5
%88 = load double, double* %87
%89 = load double*, double** %10
%90 = getelementptr inbounds double, double* %89, i64 4
%91 = load double, double* %90
%92 = fmul double %88, %91
%93 = call double @llvm.fmuladd.f64(double %82, double %85, double %92)
%94 = load double*, double** %11
%95 = getelementptr inbounds double, double* %94, i64 8
%96 = load double, double* %95
%97 = load double*, double** %10
%98 = getelementptr inbounds double, double* %97, i64 5
%99 = load double, double* %98
%100 = call double @llvm.fmuladd.f64(double %96, double %99, double %93)
store double %100, double* %20
%101 = load double*, double** %10
%102 = load double, double* %101
%103 = load double, double* %22
%104 = load double*, double** %10
%105 = getelementptr inbounds double, double* %104, i64 1
%106 = load double, double* %105
%107 = load double, double* %21
%108 = fmul double %106, %107
%109 = call double @llvm.fmuladd.f64(double %102, double %103, double %108)
%110 = load double*, double** %10
%111 = getelementptr inbounds double, double* %110, i64 2
%112 = load double, double* %111
%113 = load double, double* %20
%114 = call double @llvm.fmuladd.f64(double %112, double %113, double %109)
store double %114, double* %23
%115 = load double*, double** %11
%116 = load double, double* %115
%117 = load double*, double** %10
%118 = load double, double* %117
%119 = load double*, double** %11
%120 = getelementptr inbounds double, double* %119, i64 1
%121 = load double, double* %120
%122 = load double*, double** %10
%123 = getelementptr inbounds double, double* %122, i64 1
%124 = load double, double* %123
%125 = fmul double %121, %124
%126 = call double @llvm.fmuladd.f64(double %116, double %118, double %125)
%127 = load double*, double** %11
%128 = getelementptr inbounds double, double* %127, i64 2
%129 = load double, double* %128
%130 = load double*, double** %10
%131 = getelementptr inbounds double, double* %130, i64 2
%132 = load double, double* %131
%133 = call double @llvm.fmuladd.f64(double %129, double %132, double %126)
store double %133, double* %18
%134 = load double*, double** %11
%135 = getelementptr inbounds double, double* %134, i64 3
%136 = load double, double* %135
%137 = load double*, double** %10
%138 = load double, double* %137
%139 = load double*, double** %11
%140 = getelementptr inbounds double, double* %139, i64 4
%141 = load double, double* %140
%142 = load double*, double** %10
%143 = getelementptr inbounds double, double* %142, i64 1
%144 = load double, double* %143
%145 = fmul double %141, %144
%146 = call double @llvm.fmuladd.f64(double %136, double %138, double %145)
%147 = load double*, double** %11
%148 = getelementptr inbounds double, double* %147, i64 5
%149 = load double, double* %148
%150 = load double*, double** %10
%151 = getelementptr inbounds double, double* %150, i64 2
%152 = load double, double* %151
%153 = call double @llvm.fmuladd.f64(double %149, double %152, double %146)
store double %153, double* %19
%154 = load double, double* %22
%155 = load double, double* %22
%156 = load double, double* %21
%157 = load double, double* %21
%158 = fmul double %156, %157
%159 = call double @llvm.fmuladd.f64(double %154, double %155, double %158)
%160 = load double, double* %18
%161 = load double, double* %18
%162 = call double @llvm.fmuladd.f64(double %160, double %161, double %159)
%163 = load double, double* %19
%164 = load double, double* %19
%165 = call double @llvm.fmuladd.f64(double %163, double %164, double %162)
%166 = load double*, double** %13
store double %165, double* %166
%167 = load double, double* %23
%168 = load double, double* %23
%169 = fmul double %167, %168
%170 = load double*, double** %13
%171 = load double, double* %170
%172 = fdiv double %169, %171
%173 = load double*, double** %12
store double %172, double* %173
%174 = load double*, double** %12
%175 = load double, double* %174
%176 = load double, double* %15
%177 = fcmp olt double %175, %176
br i1 %177, label %178, label %238
178: ; preds = %38
%179 = getelementptr inbounds [3 x double], [3 x double]* %27, i64 0, i64 1
%180 = load double, double* %179
%181 = load double*, double** %10
%182 = getelementptr inbounds double, double* %181, i64 2
%183 = load double, double* %182
%184 = getelementptr inbounds [3 x double], [3 x double]* %27, i64 0, i64 2
%185 = load double, double* %184
%186 = load double*, double** %10
%187 = getelementptr inbounds double, double* %186, i64 1
%188 = load double, double* %187
%189 = fmul double %185, %188
%190 = fneg double %189
%191 = call double @llvm.fmuladd.f64(double %180, double %183, double %190)
store double %191, double* %28
%192 = getelementptr inbounds [3 x double], [3 x double]* %27, i64 0, i64 2
%193 = load double, double* %192
%194 = load double*, double** %10
%195 = getelementptr inbounds double, double* %194, i64 0
%196 = load double, double* %195
%197 = getelementptr inbounds [3 x double], [3 x double]* %27, i64 0, i64 0
%198 = load double, double* %197
%199 = load double*, double** %10
%200 = getelementptr inbounds double, double* %199, i64 2
%201 = load double, double* %200
%202 = fmul double %198, %201
%203 = fneg double %202
%204 = call double @llvm.fmuladd.f64(double %193, double %196, double %203)
store double %204, double* %29
%205 = load double, double* %28
%206 = load double, double* %22
%207 = load double, double* %29
%208 = load double, double* %21
%209 = fmul double %207, %208
%210 = call double @llvm.fmuladd.f64(double %205, double %206, double %209)
store double %210, double* %30
%211 = load double, double* %30
%212 = fcmp ogt double %211, 0.000000e+00
br i1 %212, label %213, label %221
213: ; preds = %178
%214 = load i32, i32* %24
%215 = load i32*, i32** %16
%216 = load i32, i32* %25
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds i32, i32* %215, i64 %217
store i32 %214, i32* %218
%219 = load i32, i32* %25
%220 = add nsw i32 %219, 1
store i32 %220, i32* %25
br label %231
221: ; preds = %178
%222 = load i32, i32* %26
%223 = add nsw i32 %222, 1
store i32 %223, i32* %26
%224 = load i32, i32* %24
%225 = load i32*, i32** %16
%226 = load i32, i32* %14
%227 = load i32, i32* %26
%228 = sub nsw i32 %226, %227
%229 = sext i32 %228 to i64
%230 = getelementptr inbounds i32, i32* %225, i64 %229
store i32 %224, i32* %230
br label %231
231: ; preds = %221, %213
%232 = load double*, double** %13
%233 = load double, double* %232
%234 = call i32 @sqrt(double %233)
%235 = sdiv i32 1, %234
%236 = sitofp i32 %235 to double
%237 = load double*, double** %13
store double %236, double* %237
br label %238
238: ; preds = %231, %38
%239 = load double*, double** %12
%240 = getelementptr inbounds double, double* %239, i32 1
store double* %240, double** %12
%241 = load double*, double** %13
%242 = getelementptr inbounds double, double* %241, i32 1
store double* %242, double** %13
%243 = load double*, double** %10
%244 = getelementptr inbounds double, double* %243, i64 6
store double* %244, double** %10
br label %245
245: ; preds = %238
%246 = load i32, i32* %24
%247 = add nsw i32 %246, 1
store i32 %247, i32* %24
br label %34
248: ; preds = %34
%249 = load i32, i32* %25
%250 = load i32, i32* %26
%251 = icmp sge i32 %249, %250
br i1 %251, label %252, label %256
252: ; preds = %248
%253 = load i32*, i32** %16
%254 = load i32**, i32*** %17
store i32* %253, i32** %254
%255 = load i32, i32* %25
store i32 %255, i32* %9
br label %267
256: ; preds = %248
%257 = load i32*, i32** %16
%258 = load i32, i32* %14
%259 = sext i32 %258 to i64
%260 = getelementptr inbounds i32, i32* %257, i64 %259
%261 = load i32, i32* %26
%262 = sext i32 %261 to i64
%263 = sub i64 0, %262
%264 = getelementptr inbounds i32, i32* %260, i64 %263
%265 = load i32**, i32*** %17
store i32* %264, i32** %265
%266 = load i32, i32* %26
store i32 %266, i32* %9
br label %267
267: ; preds = %256, %252
%268 = load i32, i32* %9
ret i32 %268
}

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