// File: Crawler995[P]Cpp-Ray-Tracing/libs/sphere.cpp
// Function: Sphere::intersect(Ray3)
IntersectResult Sphere::intersect(Ray3 ray) {
    Vector3 v = ray.get_origin().subtract(center);
    double a = ray.get_dir().square_length();
    double b = 2 * v.dot(ray.get_dir());
    double c = v.square_length() - radius * radius;

    double delta = b * b - 4 * a * c;

    IntersectResult res(0, 0, Vector3(0, 0, 0), Vector3(0, 0, 0));
    if(delta >= 0 && b <= 0) {
        res.geometry = this;
        res.distance = (-b - sqrt(delta)) / (2 * a);
        res.position = ray.get_point(res.distance);
        res.normal = res.position.subtract(center).normalize();
    }

    return res;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x190, %rsp            # imm = 0x190
movq %rsi, %rax
movq %rdi, -0x138(%rbp)
movq %rdi, -0x150(%rbp)
leaq 0x10(%rbp), %rsi
movq %rsi, -0x148(%rbp)
movq %rax, -0x8(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x140(%rbp)
leaq -0x38(%rbp), %rdi
callq 0x3440
movq -0x140(%rbp), %rax
movq 0x10(%rax), %rcx
movq %rcx, -0x50(%rbp)
movq 0x18(%rax), %rcx
movq %rcx, -0x48(%rbp)
movq 0x20(%rax), %rax
movq %rax, -0x40(%rbp)
leaq -0x20(%rbp), %rdi
leaq -0x38(%rbp), %rsi
leaq -0x50(%rbp), %rax
movq (%rax), %rcx
movq %rcx, (%rsp)
movq 0x8(%rax), %rcx
movq %rcx, 0x8(%rsp)
movq 0x10(%rax), %rax
movq %rax, 0x10(%rsp)
callq 0x2f80
movq -0x148(%rbp), %rsi
leaq -0x70(%rbp), %rdi
callq 0x3470
leaq -0x70(%rbp), %rdi
callq 0x2d50
movq -0x148(%rbp), %rsi
movsd %xmm0, -0x58(%rbp)
leaq -0x90(%rbp), %rdi
callq 0x3470
leaq -0x20(%rbp), %rdi
leaq -0x90(%rbp), %rax
movq (%rax), %rcx
movq %rcx, (%rsp)
movq 0x8(%rax), %rcx
movq %rcx, 0x8(%rsp)
movq 0x10(%rax), %rax
movq %rax, 0x10(%rsp)
callq 0x30d0
movaps %xmm0, %xmm1
movsd 0x58b3(%rip), %xmm0     # 0x90e0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
leaq -0x20(%rbp), %rdi
callq 0x2d50
movq -0x140(%rbp), %rax
movaps %xmm0, %xmm1
movsd 0x28(%rax), %xmm0
movsd 0x28(%rax), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x78(%rbp), %xmm0
movsd -0x78(%rbp), %xmm2
movsd 0x584c(%rip), %xmm1     # 0x90d8
mulsd -0x58(%rbp), %xmm1
mulsd -0x98(%rbp), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
leaq -0xb8(%rbp), %rdi
xorps %xmm2, %xmm2
movaps %xmm2, %xmm0
movaps %xmm2, %xmm1
callq 0x2c50
leaq -0xd0(%rbp), %rdi
xorps %xmm2, %xmm2
movaps %xmm2, %xmm0
movaps %xmm2, %xmm1
callq 0x2c50
movq -0x138(%rbp), %rdi
xorl %eax, %eax
movl %eax, %esi
xorps %xmm0, %xmm0
leaq -0xb8(%rbp), %rcx
leaq -0xd0(%rbp), %rax
movq (%rcx), %rdx
movq %rdx, (%rsp)
movq 0x8(%rcx), %rdx
movq %rdx, 0x8(%rsp)
movq 0x10(%rcx), %rcx
movq %rcx, 0x10(%rsp)
movq (%rax), %rcx
movq %rcx, 0x18(%rsp)
movq 0x8(%rax), %rcx
movq %rcx, 0x20(%rsp)
movq 0x10(%rax), %rax
movq %rax, 0x28(%rsp)
callq 0x3bc0
movsd -0xa0(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x3aa2
xorps %xmm0, %xmm0
ucomisd -0x78(%rbp), %xmm0
jb 0x3aa2
movq -0x138(%rbp), %rax
movq -0x140(%rbp), %rcx
movq %rcx, (%rax)
movsd -0x78(%rbp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
movsd %xmm0, -0x158(%rbp)
movsd -0xa0(%rbp), %xmm0
callq 0x20d0
movq -0x148(%rbp), %rsi
movq -0x138(%rbp), %rax
movaps %xmm0, %xmm1
movsd -0x158(%rbp), %xmm0
subsd %xmm1, %xmm0
movsd 0x5718(%rip), %xmm1     # 0x90e0
mulsd -0x58(%rbp), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rax)
movsd 0x8(%rax), %xmm0
leaq -0xe8(%rbp), %rdi
callq 0x33d0
movq -0x140(%rbp), %rax
movq -0x138(%rbp), %rsi
movq -0xe8(%rbp), %rcx
movq %rcx, 0x10(%rsi)
movq -0xe0(%rbp), %rcx
movq %rcx, 0x18(%rsi)
movq -0xd8(%rbp), %rcx
movq %rcx, 0x20(%rsi)
addq $0x10, %rsi
movq 0x10(%rax), %rcx
movq %rcx, -0x130(%rbp)
movq 0x18(%rax), %rcx
movq %rcx, -0x128(%rbp)
movq 0x20(%rax), %rax
movq %rax, -0x120(%rbp)
leaq -0x118(%rbp), %rdi
leaq -0x130(%rbp), %rax
movq (%rax), %rcx
movq %rcx, (%rsp)
movq 0x8(%rax), %rcx
movq %rcx, 0x8(%rsp)
movq 0x10(%rax), %rax
movq %rax, 0x10(%rsp)
callq 0x2f80
leaq -0x100(%rbp), %rdi
leaq -0x118(%rbp), %rsi
callq 0x2d90
movq -0x138(%rbp), %rax
movq -0x100(%rbp), %rcx
movq %rcx, 0x28(%rax)
movq -0xf8(%rbp), %rcx
movq %rcx, 0x30(%rax)
movq -0xf0(%rbp), %rcx
movq %rcx, 0x38(%rax)
movq -0x150(%rbp), %rax
addq $0x190, %rsp            # imm = 0x190
popq %rbp
retq
nopw %cs:(%rax,%rax)