// File: /workspace/github_c_makefile_star2_O3/c-zhou[P]yahs/binomlite.c
// Function: do_search
static double do_search(double y, double *z, double p, double n, double pr, double incr)
{
    if (*z >= p) {
        /* search to the left */
        for (;;) {
            double newz;
            if (y == 0 || (newz = pbinom(y - incr, n, pr, 1, 0)) < p)
                return y;
            y = fmax2(0, y - incr);
            *z = newz;
        }
    } else {
        /* search to the right */
        for (;;) {
            y = fmin2(y + incr, n);
            if (y == n || (*z = pbinom(y, n, pr, 1, 0)) >= p)
                return y;
        }
    }
}

// Assembly:
pushq %rbp
pushq %rbx
subq $0x68, %rsp
movsd %xmm4, 0x38(%rsp)
movapd %xmm3, 0x40(%rsp)
movq %rdi, %rbx
movapd %xmm0, %xmm4
movsd (%rdi), %xmm0
movsd %xmm1, 0x30(%rsp)
ucomisd %xmm1, %xmm0
jae 0x1bfb7
movaps %xmm2, %xmm0
movapd %xmm4, 0x10(%rsp)
movaps %xmm2, 0x20(%rsp)
callq 0x2090
movapd 0x20(%rsp), %xmm5
movapd 0x40(%rsp), %xmm3
movsd %xmm0, 0x8(%rsp)
movsd 0x67a6(%rip), %xmm0     # 0x22660
cmplepd %xmm3, %xmm0
xorpd %xmm1, %xmm1
movapd %xmm3, %xmm2
cmplepd %xmm1, %xmm2
orpd %xmm0, %xmm2
movd %xmm2, %ebp
movapd 0x10(%rsp), %xmm4
cmplesd %xmm1, %xmm3
movapd %xmm3, 0x50(%rsp)
addsd 0x38(%rsp), %xmm4
movapd %xmm4, %xmm0
cmplesd %xmm5, %xmm0
andpd %xmm0, %xmm4
andnpd %xmm5, %xmm0
orpd %xmm0, %xmm4
ucomisd %xmm5, %xmm4
jne 0x1bf0c
jnp 0x1c0f0
xorpd %xmm0, %xmm0
ucomisd %xmm4, %xmm0
ja 0x1bfa2
movapd %xmm4, 0x10(%rsp)
addsd 0x7d30(%rip), %xmm4     # 0x23c58
movapd %xmm4, %xmm0
callq 0x2470
movapd %xmm0, %xmm1
ucomisd 0x8(%rsp), %xmm0
setae %al
movsd 0x671a(%rip), %xmm0     # 0x22660
jae 0x1bf5a
movapd 0x50(%rsp), %xmm0
movsd 0x670a(%rip), %xmm2     # 0x22660
andpd %xmm2, %xmm0
orb %bpl, %al
testb $0x1, %al
movapd 0x10(%rsp), %xmm4
movapd 0x40(%rsp), %xmm3
movapd 0x20(%rsp), %xmm5
jne 0x1bfa2
movsd 0x8(%rsp), %xmm2
subsd %xmm1, %xmm2
movsd 0x66db(%rip), %xmm0     # 0x22660
addsd %xmm0, %xmm1
movapd %xmm3, %xmm0
xorl %edi, %edi
xorl %esi, %esi
callq 0x1b899
movapd 0x20(%rsp), %xmm5
movapd 0x10(%rsp), %xmm4
movsd %xmm0, (%rbx)
ucomisd 0x30(%rsp), %xmm0
jb 0x1bee5
jmp 0x1c0f0
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm4
jne 0x1bfc7
jnp 0x1c0f0
movapd %xmm2, %xmm0
movapd %xmm4, 0x10(%rsp)
callq 0x2090
xorpd %xmm2, %xmm2
movapd 0x40(%rsp), %xmm3
movapd 0x10(%rsp), %xmm4
movsd %xmm0, 0x8(%rsp)
movsd 0x666c(%rip), %xmm0     # 0x22660
cmplepd %xmm3, %xmm0
movapd %xmm3, %xmm1
cmplepd %xmm2, %xmm1
orpd %xmm0, %xmm1
movd %xmm1, %ebp
cmplesd %xmm2, %xmm3
movapd %xmm3, 0x50(%rsp)
movapd %xmm4, %xmm5
subsd 0x38(%rsp), %xmm5
xorpd %xmm0, %xmm0
ucomisd %xmm5, %xmm0
ja 0x1c0c3
movapd %xmm4, 0x10(%rsp)
movapd %xmm5, 0x20(%rsp)
movapd %xmm5, %xmm0
addsd 0x7c13(%rip), %xmm0     # 0x23c58
callq 0x2470
movapd %xmm0, %xmm1
ucomisd 0x8(%rsp), %xmm0
setae %al
movsd 0x6601(%rip), %xmm0     # 0x22660
jae 0x1c073
movapd 0x50(%rsp), %xmm0
movsd 0x65f1(%rip), %xmm2     # 0x22660
andpd %xmm2, %xmm0
orb %bpl, %al
testb $0x1, %al
movapd 0x10(%rsp), %xmm4
movapd 0x40(%rsp), %xmm3
xorpd %xmm2, %xmm2
movapd 0x20(%rsp), %xmm5
jne 0x1c0c3
movsd 0x8(%rsp), %xmm2
subsd %xmm1, %xmm2
movsd 0x65be(%rip), %xmm0     # 0x22660
addsd %xmm0, %xmm1
movapd %xmm3, %xmm0
xorl %edi, %edi
xorl %esi, %esi
callq 0x1b899
movapd 0x20(%rsp), %xmm5
xorpd %xmm2, %xmm2
movapd 0x10(%rsp), %xmm4
movsd 0x30(%rsp), %xmm1
ucomisd %xmm0, %xmm1
ja 0x1c0f0
movsd %xmm0, (%rbx)
movapd %xmm5, %xmm4
cmplesd %xmm2, %xmm4
andnpd %xmm5, %xmm4
ucomisd %xmm2, %xmm4
jne 0x1c015
jp 0x1c015
movapd %xmm4, %xmm0
addq $0x68, %rsp
popq %rbx
popq %rbp
retq