// File: /workspace/github_c_makefile_star2_O0/yui0[P]img2vec/./potracelib.h
// Function: tangent
static double tangent(dpoint_t p0, dpoint_t p1, dpoint_t p2, dpoint_t p3, dpoint_t q0, dpoint_t q1) {
  double A, B, C;   /* (1-t)^2 A + 2(1-t)t B + t^2 C = 0 */
  double a, b, c;   /* a t^2 + b t + c = 0 */
  double d, s, r1, r2;
  A = cprod(p0, p1, q0, q1);
  B = cprod(p1, p2, q0, q1);
  C = cprod(p2, p3, q0, q1);
  a = A - 2*B + C;
  b = -2*A + 2*B;
  c = A;
  
  d = b*b - 4*a*c;
  if (a==0 || d<0) {
    return -1.0;
  }
  s = sqrt(d);
  r1 = (-b + s) / (2 * a);
  r2 = (-b - s) / (2 * a);
  if (r1 >= 0 && r1 <= 1) {
    return r1;
  } else if (r2 >= 0 && r2 <= 1) {
    return r2;
  } else {
    return -1.0;
  }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xb0, %rsp
leaq 0x20(%rbp), %rax
movq %rax, -0xa0(%rbp)
leaq 0x10(%rbp), %rcx
movq %rcx, -0xa8(%rbp)
movsd %xmm0, -0x18(%rbp)
movsd %xmm1, -0x10(%rbp)
movsd %xmm2, -0x28(%rbp)
movsd %xmm3, -0x20(%rbp)
movsd %xmm4, -0x38(%rbp)
movsd %xmm5, -0x30(%rbp)
movsd %xmm6, -0x48(%rbp)
movsd %xmm7, -0x40(%rbp)
movsd -0x18(%rbp), %xmm0
movsd -0x10(%rbp), %xmm1
movsd -0x28(%rbp), %xmm2
movsd -0x20(%rbp), %xmm3
movsd (%rcx), %xmm4
movsd 0x8(%rcx), %xmm5
movsd (%rax), %xmm6
movsd 0x8(%rax), %xmm7
callq 0x2e7f0
movq -0xa8(%rbp), %rcx
movq -0xa0(%rbp), %rax
movsd %xmm0, -0x50(%rbp)
movsd -0x28(%rbp), %xmm0
movsd -0x20(%rbp), %xmm1
movsd -0x38(%rbp), %xmm2
movsd -0x30(%rbp), %xmm3
movsd (%rcx), %xmm4
movsd 0x8(%rcx), %xmm5
movsd (%rax), %xmm6
movsd 0x8(%rax), %xmm7
callq 0x2e7f0
movq -0xa8(%rbp), %rcx
movq -0xa0(%rbp), %rax
movsd %xmm0, -0x58(%rbp)
movsd -0x38(%rbp), %xmm0
movsd -0x30(%rbp), %xmm1
movsd -0x48(%rbp), %xmm2
movsd -0x40(%rbp), %xmm3
movsd (%rcx), %xmm4
movsd 0x8(%rcx), %xmm5
movsd (%rax), %xmm6
movsd 0x8(%rax), %xmm7
callq 0x2e7f0
movsd %xmm0, -0x60(%rbp)
movsd -0x50(%rbp), %xmm0
movsd -0x58(%rbp), %xmm1
addsd %xmm1, %xmm1
subsd %xmm1, %xmm0
addsd -0x60(%rbp), %xmm0
movsd %xmm0, -0x68(%rbp)
movsd -0x50(%rbp), %xmm1
movsd 0x360a(%rip), %xmm0     # 0x32048
mulsd -0x58(%rbp), %xmm0
addsd %xmm1, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, -0x70(%rbp)
movsd -0x50(%rbp), %xmm0
movsd %xmm0, -0x78(%rbp)
movsd -0x70(%rbp), %xmm0
movsd -0x70(%rbp), %xmm2
movsd 0x35e4(%rip), %xmm1     # 0x32050
mulsd -0x68(%rbp), %xmm1
mulsd -0x78(%rbp), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x80(%rbp)
movsd -0x68(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jne 0x2eaac
jp 0x2eaac
jmp 0x2eab6
xorps %xmm0, %xmm0
ucomisd -0x80(%rbp), %xmm0
jbe 0x2eac8
movsd 0x3562(%rip), %xmm0     # 0x32020
movsd %xmm0, -0x8(%rbp)
jmp 0x2ebc5
movsd -0x80(%rbp), %xmm0
callq 0x1220
movsd %xmm0, -0x88(%rbp)
movsd -0x70(%rbp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
addsd -0x88(%rbp), %xmm0
movsd 0x3542(%rip), %xmm1     # 0x32048
mulsd -0x68(%rbp), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x90(%rbp)
movsd -0x70(%rbp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
subsd -0x88(%rbp), %xmm0
movsd 0x3505(%rip), %xmm1     # 0x32048
mulsd -0x68(%rbp), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x90(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x2eb86
movsd 0x34ab(%rip), %xmm0     # 0x32018
ucomisd -0x90(%rbp), %xmm0
jb 0x2eb86
movsd -0x90(%rbp), %xmm0
movsd %xmm0, -0x8(%rbp)
jmp 0x2ebc5
movsd -0x98(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x2ebb8
movsd 0x3479(%rip), %xmm0     # 0x32018
ucomisd -0x98(%rbp), %xmm0
jb 0x2ebb8
movsd -0x98(%rbp), %xmm0
movsd %xmm0, -0x8(%rbp)
jmp 0x2ebc5
movsd 0x3460(%rip), %xmm0     # 0x32020
movsd %xmm0, -0x8(%rbp)
movsd -0x8(%rbp), %xmm0
addq $0xb0, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)