// File: /workspace/github_cpp_makefile_star2_O3/Cell-veto[P]postlhc/./adap_prober.hpp
// Function: AdaptiveProber::calib_add(double, double, double)
void calib_add (double r, double rate, double rtol = 0.)
    {
        if (rate == 0.) return;

        if (! (r >= 0.))
            std::cerr << "calib_add: negative distance" << ABORT;

        if (!std::isfinite (rate))
            std::cerr << "calib_add: nonfinite lifting rate" << ABORT;

        // range of bins affected (including imax)
        unsigned i =    bin_from_radius (r-rtol);
        unsigned imax = bin_from_radius (r+rtol);
        assert (imax < N_);

        // zero bin (const probe rate)
        if (i == 0)
        {
            binrate_[i] = fmax (binrate_[i], rate);
            ++i;
        }

        // outer bins (correct for appropriate exponent)
        for (; i <= imax; ++i)
        {
            double ir = fmin (binrad_[i], r+rtol); // FIXME??
            double irate = rate / pow (ir, -(1.+binexp_[i]));
            binrate_[i] = fmax (binrate_[i], irate);
        }
    }

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
subq $0x20, %rsp
movapd %xmm0, %xmm4
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
jne 0x41067f
jnp 0x410816
ucomisd %xmm0, %xmm4
jb 0x410820
movq %xmm1, %rax
btrq $0x3f, %rax
movabsq $0x7ff0000000000000, %rcx # imm = 0x7FF0000000000000
cmpq %rcx, %rax
jge 0x410830
movapd %xmm4, %xmm0
subsd %xmm2, %xmm0
leaq 0x200(%rdi), %rcx
movl $0x40, %esi
movq %rcx, %rax
movq %rsi, %rdx
shrq %rdx
ucomisd (%rax,%rdx,8), %xmm0
jbe 0x4106d8
leaq (%rax,%rdx,8), %rax
addq $0x8, %rax
notq %rdx
addq %rsi, %rdx
movq %rdx, %rsi
testq %rdx, %rdx
jg 0x4106c0
subq %rcx, %rax
shrq $0x3, %rax
addsd %xmm2, %xmm4
movl $0x40, %esi
movq %rcx, %rbx
movq %rsi, %rdx
shrq %rdx
ucomisd (%rbx,%rdx,8), %xmm4
jbe 0x41070e
leaq (%rbx,%rdx,8), %rbx
addq $0x8, %rbx
notq %rdx
addq %rsi, %rdx
movq %rdx, %rsi
testq %rdx, %rdx
jg 0x4106f6
subq %rcx, %rbx
shrq $0x3, %rbx
cmpl $0x40, %ebx
jae 0x410852
movl %eax, %ecx
testl %eax, %eax
jne 0x41075e
movsd 0x600(%rdi), %xmm0
movapd %xmm1, %xmm3
maxsd %xmm0, %xmm3
cmpunordsd %xmm0, %xmm0
movapd %xmm0, %xmm2
andpd %xmm1, %xmm2
andnpd %xmm3, %xmm0
orpd %xmm2, %xmm0
movlpd %xmm0, 0x600(%rdi)
movl $0x1, %ecx
cmpl %ebx, %ecx
ja 0x410816
movl %eax, %eax
cmpq $0x1, %rax
adcq $0x0, %rax
leaq (%rdi,%rax,8), %r14
addq $0x600, %r14            # imm = 0x600
subl %ecx, %ebx
incl %ebx
xorl %r15d, %r15d
movaps %xmm1, 0x10(%rsp)
movapd %xmm4, (%rsp)
movsd -0x400(%r14,%r15,8), %xmm0
movapd %xmm4, %xmm1
minsd %xmm0, %xmm1
cmpunordsd %xmm0, %xmm0
movapd %xmm0, %xmm2
andnpd %xmm1, %xmm2
andpd %xmm4, %xmm0
orpd %xmm2, %xmm0
movsd -0x200(%r14,%r15,8), %xmm1
addsd 0x1b8853(%rip), %xmm1   # 0x5c9018
xorpd 0x1b9873(%rip), %xmm1   # 0x5ca040
callq 0x4ecf40
movapd 0x10(%rsp), %xmm4
divsd %xmm0, %xmm4
movsd (%r14,%r15,8), %xmm0
movapd %xmm0, %xmm2
cmpunordsd %xmm0, %xmm2
movapd %xmm2, %xmm3
andpd %xmm4, %xmm3
maxsd %xmm0, %xmm4
andnpd %xmm4, %xmm2
movapd (%rsp), %xmm4
orpd %xmm3, %xmm2
movlpd %xmm2, (%r14,%r15,8)
incq %r15
cmpl %r15d, %ebx
jne 0x41078c
addq $0x20, %rsp
popq %rbx
popq %r14
popq %r15
retq
movq $0x655aa0, %rdi         # imm = 0x655AA0
leaq 0x1b9e46(%rip), %rsi    # 0x5ca674
jmp 0x41083e
movq $0x655aa0, %rdi         # imm = 0x655AA0
leaq 0x1b9e53(%rip), %rsi    # 0x5ca691
callq 0x4b0650
leaq 0x243db7(%rip), %rsi    # 0x654601
movq %rax, %rdi
callq 0x43f8ce
leaq 0x1b9e5a(%rip), %rdi    # 0x5ca6b3
leaq 0x1b9d5e(%rip), %rsi    # 0x5ca5be
leaq 0x1b9e56(%rip), %rcx    # 0x5ca6bd
movl $0x45, %edx
callq 0x50b050
nop