// File: /workspace/github_cpp_makefile_star2_O2/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 $0xe8, %rsp
testq %rsi, %rsi
je 0x418503
movq %rdi, %rbx
leaq 0x22a4c6(%rip), %rax    # 0x642678
leaq 0x22a50f(%rip), %rdx    # 0x6426c8
movq %rsi, %rdi
movq %rax, %rsi
xorl %ecx, %ecx
callq 0x437a70
testq %rax, %rax
je 0x418503
movq %rax, %r14
movq %rax, %rdi
callq 0x4100fc
leaq 0x98(%rsp), %rax
movsd %xmm0, -0x10(%rax)
movsd %xmm1, -0x8(%rax)
leaq 0x28(%rsp), %rsi
movq %rax, (%rsi)
leaq 0x88(%rsp), %r12
movq %r12, %rdi
callq 0x418525
movsd %xmm0, 0x38(%rbx)
movq %r14, %rdi
callq 0x418540
leaq 0x78(%rsp), %rcx
movq %rax, (%rcx)
movq %rdx, 0x8(%rcx)
movq $0x64eaa0, %r15         # imm = 0x64EAA0
leaq 0x1a997e(%rip), %rsi    # 0x5c1ba9
movq %r15, %rdi
callq 0x4a4910
movq %rax, %rdi
movq %r12, %rsi
callq 0x40b917
leaq 0x1a7336(%rip), %rsi    # 0x5bf57b
movq %rax, %rdi
callq 0x4a4910
movsd 0x38(%rbx), %xmm0
movq %rax, %rdi
callq 0x4a5640
leaq 0x1b81a9(%rip), %r12    # 0x5d040a
movq %rax, %rdi
movq %r12, %rsi
callq 0x4a4910
leaq 0x1a9951(%rip), %rsi    # 0x5c1bc4
movq %r15, %rdi
callq 0x4a4910
movq %rax, %r15
movq (%r14), %rax
pushq $0x2
popq %rsi
movq %r14, %rdi
callq *0x48(%rax)
movq %r15, %rdi
movq %rax, %rsi
callq 0x4a4ce0
movq %rax, %rdi
movq %r12, %rsi
callq 0x4a4910
leaq 0x48(%rbx), %rax
movq %rax, 0x70(%rsp)
xorl %r13d, %r13d
leaq 0x58(%rsp), %r12
leaq 0x98(%rsp), %rbp
leaq 0x48(%rsp), %r14
leaq 0xd8(%rsp), %r15
cmpq $0x2, %r13
je 0x4184f1
imulq $0xe0, %r13, %rdi
addq 0x70(%rsp), %rdi
movq %rdi, 0x40(%rsp)
callq 0x41858e
leaq 0x28(%rsp), %rdi
leaq 0x78(%rsp), %rsi
callq 0x4185aa
movq 0x78(%rsp), %rax
movq 0x80(%rsp), %rcx
movq %rax, 0xb8(%rsp)
movq %rcx, 0xc0(%rsp)
movl %ecx, 0xb0(%rsp)
shrq $0x20, %rax
movl %eax, 0xb4(%rsp)
leaq 0x28(%rsp), %rdi
leaq 0xb0(%rsp), %rsi
callq 0x418903
testb %al, %al
je 0x4184df
leaq 0x88(%rsp), %rdi
leaq 0x28(%rsp), %rsi
callq 0x4185c8
movsd %xmm0, 0xd8(%rsp)
movsd %xmm1, 0xe0(%rsp)
mulsd %xmm0, %xmm0
mulsd %xmm1, %xmm1
addsd %xmm0, %xmm1
sqrtsd %xmm1, %xmm0
movsd 0x18(%rbx), %xmm1
addsd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
pushq $-0x6
popq %rax
pushq $-0x32
popq %rcx
cmoval %ecx, %eax
pushq $0x6
popq %rcx
pushq $0x32
popq %rdx
cmoval %edx, %ecx
movl %ecx, 0x14(%rsp)
incl %ecx
movl %eax, 0x18(%rsp)
movl %eax, 0x1c(%rsp)
movl %ecx, 0x20(%rsp)
movl %ecx, 0x24(%rsp)
movq %r12, %rdi
leaq 0x18(%rsp), %rsi
callq 0x4185aa
movq 0x18(%rsp), %rax
movq 0x20(%rsp), %rcx
movq %rax, 0xa0(%rsp)
movq %rcx, 0xa8(%rsp)
movl %ecx, 0x98(%rsp)
shrq $0x20, %rax
movl %eax, 0x9c(%rsp)
xorpd %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
movq %r12, %rdi
movq %rbp, %rsi
callq 0x418903
testb %al, %al
je 0x41849f
leaq 0x88(%rsp), %rdi
movq %r12, %rsi
callq 0x4185c8
movsd %xmm0, 0x48(%rsp)
movsd %xmm1, 0x50(%rsp)
movq %r14, %rdi
leaq 0x14(%rsp), %rsi
callq 0x4185f2
movsd %xmm0, 0xc8(%rsp)
movsd %xmm1, 0xd0(%rsp)
movq %r15, %rdi
leaq 0xc8(%rsp), %rsi
callq 0x40ba98
movsd %xmm0, 0x48(%rsp)
movsd %xmm1, 0x50(%rsp)
movq %rbx, %rdi
movq %r14, %rsi
movl %r13d, %edx
pushq $0x1
popq %rcx
callq 0x418656
movapd %xmm0, %xmm1
movsd 0x8(%rsp), %xmm0
callq 0x4e1c90
movsd %xmm0, 0x8(%rsp)
movl 0x5c(%rsp), %eax
incl %eax
movl %eax, 0x5c(%rsp)
cmpl 0x6c(%rsp), %eax
jne 0x4183ea
movl 0x64(%rsp), %eax
movl %eax, 0x5c(%rsp)
incl 0x58(%rsp)
jmp 0x4183ea
movsd 0x8(%rsp), %xmm0
mulsd 0x1a965b(%rip), %xmm0   # 0x5c1b08
movq 0x40(%rsp), %rdi
movq %r15, %rsi
callq 0x418898
movl 0x2c(%rsp), %eax
incl %eax
movl %eax, 0x2c(%rsp)
cmpl 0x3c(%rsp), %eax
jne 0x418324
movl 0x34(%rsp), %eax
movl %eax, 0x2c(%rsp)
incl 0x28(%rsp)
jmp 0x418324
movq 0x40(%rsp), %rdi
callq 0x4188d2
incq %r13
jmp 0x4182c6
addq $0xe8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq $0x64eaa0, %rdi         # imm = 0x64EAA0
leaq 0x1a966f(%rip), %rsi    # 0x5c1b80
callq 0x4a4910
leaq 0x23513c(%rip), %rsi    # 0x64d659
movq %rax, %rdi
callq 0x434f95