// File: /workspace/github_cpp_makefile_star2_O1/Cell-veto[P]postlhc/jellium4.cpp
// Function: Jellium4<2u, Monodisperse2D>::setup_probe_distr(AbstractStorage const*)
void setup_probe_distr (const AbstractStorage *stor_)
    {
        // FIXME introduce CellStorage
        auto stor = dynamic_cast <const Storage <ENCODING> *> (stor_);
        if (!stor)
            std::cerr << "jellium4: incompatible Storage passed in" << ABORT;

        vector_t cwid = stor->cell_widths ();
        cell_volume = fproduct (cwid.begin (), cwid.end ());

        GridIndices <DIM> cell_grid = this->get_cell_grid (*stor);

        std::cerr << "jellium4_init cell_volume " << cwid << " = " << cell_volume << "\n";
        std::cerr << "jellium4_init cell_count(DIM) " << stor->cell_count (DIM) << "\n";

        for (unsigned direction = 0; direction != DIM; ++direction)
        {
            prober_t *pr = &prober[direction];
            pr->clear ();

            auto it = cell_grid.begin (), it_end = cell_grid.end ();
            for (; it != it_end; ++it)
            {
                vector_t cell_vector = elementwise_product (cwid, *it);
                int discret = (norm (cell_vector) < 2 * sr_lr_split) ? 50 : 6;
                GridIndices <DIM> discret_grid = GridIndices <DIM> (-discret, discret+1);

                double max_event_rate_found = 0.;

                auto iit = discret_grid.begin (), iit_end = discret_grid.end ();
                for (; iit != iit_end; ++iit)
                {
                    vector_t dr = elementwise_product (cwid, *iit) / discret;
                    double er = lr_event_rate (cell_vector + dr, direction);
                    update_max (&max_event_rate_found, er);
                }

                // report maximum to prober
                pr->add (max_event_rate_found * SAFETY_CUSHION, cell_vector);
            }

            pr->finish ();
        }
    }

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x108, %rsp            # imm = 0x108
movq %rdi, 0x28(%rsp)
testq %rsi, %rsi
je 0x41c767
leaq 0x22de62(%rip), %rax    # 0x64a5b0
leaq 0x22deab(%rip), %rdx    # 0x64a600
movq %rsi, %rdi
movq %rax, %rsi
xorl %ecx, %ecx
callq 0x4440f0
movq %rax, %r14
jmp 0x41c76a
xorl %r14d, %r14d
testq %r14, %r14
je 0x41cbc1
movq %r14, %rdi
callq 0x4133c6
movsd %xmm0, 0x68(%rsp)
movsd %xmm1, 0x70(%rsp)
movsd 0x1ad889(%rip), %xmm0   # 0x5ca018
xorl %eax, %eax
mulsd 0x68(%rsp,%rax), %xmm0
addq $0x8, %rax
cmpq $0x10, %rax
jne 0x41c791
movq 0x28(%rsp), %rbx
movsd %xmm0, 0x38(%rbx)
movq %r14, %rdi
callq 0x41cbe4
movq %rax, %r13
movq %rdx, 0x98(%rsp)
movq $0x656aa0, %r15         # imm = 0x656AA0
leaq 0x1b112d(%rip), %rsi    # 0x5cd8f9
movl $0x1a, %edx
movq %r15, %rdi
callq 0x4b14d0
leaq 0x68(%rsp), %rsi
movq %r15, %rdi
callq 0x40e5ab
movq %rax, %r12
leaq 0x1aec5e(%rip), %rsi    # 0x5cb44e
movl $0x3, %edx
movq %rax, %rdi
callq 0x4b14d0
movsd 0x38(%rbx), %xmm0
movq %r12, %rdi
callq 0x4b2730
leaq 0x1bf3d9(%rip), %r12    # 0x5dbbea
movl $0x1, %edx
movq %rax, %rdi
movq %r12, %rsi
callq 0x4b14d0
leaq 0x1b10ec(%rip), %rsi    # 0x5cd914
movl $0x1e, %edx
movq %r15, %rdi
callq 0x4b14d0
movq (%r14), %rax
movq %r14, %rdi
movl $0x2, %esi
callq *0x48(%rax)
movq %r15, %rdi
movq %rax, %rsi
callq 0x4b1dd0
movl $0x1, %edx
movq %rax, %rdi
movq %r12, %rsi
callq 0x4b14d0
leaq 0x48(%rbx), %rax
movq %rax, 0x90(%rsp)
movq %r13, 0xa0(%rsp)
shrq $0x20, %r13
movq %r13, 0x88(%rsp)
xorl %r12d, %r12d
leaq 0x58(%rsp), %rbp
imulq $0xe0, %r12, %rcx
movq 0x90(%rsp), %rdx
leaq (%rdx,%rcx), %rbx
movq (%rdx,%rcx), %rax
cmpq %rax, 0x8(%rdx,%rcx)
je 0x41c8a8
movq %rax, 0x8(%rbx)
movq 0x20(%rbx), %rax
cmpq %rax, 0x28(%rbx)
je 0x41c8b6
movq %rax, 0x28(%rbx)
movq 0xa0(%rsp), %rax
movq %rax, 0x38(%rsp)
movq 0x98(%rsp), %rcx
movq %rcx, 0x40(%rsp)
movq %rax, 0x30(%rsp)
movq %rax, 0xd0(%rsp)
movq %rcx, 0xd8(%rsp)
movl %ecx, 0xc8(%rsp)
movq 0x88(%rsp), %rcx
movl %ecx, 0xcc(%rsp)
cmpq 0xc8(%rsp), %rax
je 0x41cb9a
movl 0x30(%rsp), %eax
movl %eax, 0xc(%rsp)
movl 0x34(%rsp), %r14d
movl 0x44(%rsp), %eax
movl %eax, 0x4c(%rsp)
movl 0x3c(%rsp), %eax
movl %eax, 0x48(%rsp)
movq %rbx, 0xa8(%rsp)
xorl %eax, %eax
xorps %xmm0, %xmm0
cvtsi2sdl 0x30(%rsp,%rax,4), %xmm0
mulsd 0x68(%rsp,%rax,8), %xmm0
movsd %xmm0, 0x10(%rsp,%rax,8)
incq %rax
cmpq $0x2, %rax
jne 0x41c930
movsd 0x10(%rsp), %xmm1
movsd 0x18(%rsp), %xmm0
movsd %xmm1, 0xb8(%rsp)
movsd %xmm0, 0xc0(%rsp)
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
ucomisd 0x1ae668(%rip), %xmm0   # 0x5cafe8
jb 0x41c988
sqrtsd %xmm0, %xmm0
jmp 0x41c98d
callq 0x4ee410
movq 0x28(%rsp), %rax
movsd 0x18(%rax), %xmm1
addsd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
movl $0xfffffffa, %ecx       # imm = 0xFFFFFFFA
movl $0xffffffce, %eax       # imm = 0xFFFFFFCE
cmoval %eax, %ecx
movl $0x6, %eax
movl $0x32, %edx
cmoval %edx, %eax
movl %eax, 0x54(%rsp)
xorl %edx, %edx
movl %ecx, 0x78(%rsp,%rdx)
addq $0x4, %rdx
cmpq $0x8, %rdx
jne 0x41c9bf
orl $0x1, %eax
movl $0x8, %ecx
movl %eax, 0x78(%rsp,%rcx)
addq $0x4, %rcx
cmpq $0x10, %rcx
jne 0x41c9d5
movq 0x78(%rsp), %rax
movq 0x7c(%rsp), %rcx
movq %rax, 0x18(%rsp)
movq 0x80(%rsp), %rdx
movq %rdx, 0x20(%rsp)
movq %rax, 0x10(%rsp)
movq %rax, 0xf8(%rsp)
movq %rdx, 0x100(%rsp)
rolq $0x20, %rcx
movq %rcx, 0xf0(%rsp)
cmpq %rcx, %rax
je 0x41cb4b
movl %r14d, 0x50(%rsp)
movl 0x10(%rsp), %r13d
movl 0x14(%rsp), %r15d
movl 0x24(%rsp), %r14d
movl 0x1c(%rsp), %ebx
xorpd %xmm0, %xmm0
movsd %xmm0, 0xb0(%rsp)
xorl %eax, %eax
xorps %xmm0, %xmm0
cvtsi2sdl 0x10(%rsp,%rax,4), %xmm0
mulsd 0x68(%rsp,%rax,8), %xmm0
movsd %xmm0, 0x58(%rsp,%rax,8)
incq %rax
cmpq $0x2, %rax
jne 0x41ca50
movups 0x58(%rsp), %xmm0
movups %xmm0, 0x58(%rsp)
movq %rbp, %rdi
leaq 0x54(%rsp), %rsi
callq 0x41cc2f
movsd %xmm0, 0xe0(%rsp)
movsd %xmm1, 0xe8(%rsp)
xorl %eax, %eax
movsd 0xb8(%rsp,%rax,8), %xmm0
addsd 0xe0(%rsp,%rax,8), %xmm0
movsd %xmm0, 0x58(%rsp,%rax,8)
incq %rax
cmpq $0x2, %rax
jne 0x41ca99
movupd 0x58(%rsp), %xmm0
movupd %xmm0, 0x58(%rsp)
movq 0x28(%rsp), %rdi
movq %rbp, %rsi
movl %r12d, %edx
movl $0x1, %ecx
callq 0x41cc94
incl %r15d
movl %r15d, 0x14(%rsp)
cmpl %r14d, %r15d
jne 0x41caf7
movl %ebx, 0x14(%rsp)
incl %r13d
movl %r13d, 0x10(%rsp)
movl %ebx, %r15d
movsd 0xb0(%rsp), %xmm3
movapd %xmm3, %xmm1
cmpunordsd %xmm3, %xmm1
movapd %xmm1, %xmm2
andpd %xmm0, %xmm2
maxsd %xmm3, %xmm0
andnpd %xmm0, %xmm1
orpd %xmm2, %xmm1
movapd %xmm1, %xmm0
movq 0x10(%rsp), %rax
cmpq 0xf0(%rsp), %rax
jne 0x41ca45
mulsd 0x1b0d1c(%rip), %xmm0   # 0x5cd858
movq 0xa8(%rsp), %rbx
movl 0x50(%rsp), %r14d
jmp 0x41cb4f
xorpd %xmm0, %xmm0
movq %rbx, %rdi
leaq 0xb8(%rsp), %rsi
callq 0x41cf16
incl %r14d
movl %r14d, 0x34(%rsp)
cmpl 0x4c(%rsp), %r14d
jne 0x41cb87
movl 0x48(%rsp), %ecx
movl %ecx, 0x34(%rsp)
movl 0xc(%rsp), %eax
incl %eax
movl %eax, 0xc(%rsp)
movl %eax, 0x30(%rsp)
movl %ecx, %r14d
movq 0x30(%rsp), %rax
cmpq 0xc8(%rsp), %rax
jne 0x41c92e
movq %rbx, %rdi
callq 0x41cf84
incq %r12
cmpq $0x2, %r12
jne 0x41c886
addq $0x108, %rsp            # imm = 0x108
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq $0x656aa0, %rdi         # imm = 0x656AA0
leaq 0x1b0d01(%rip), %rsi    # 0x5cd8d0
callq 0x4b1a00
leaq 0x238a7e(%rip), %rsi    # 0x655659
movq %rax, %rdi
callq 0x440cc2
nop