// File: Cell-veto[P]postlhc/./ipl.hpp
// Function: PowerLawInteraction::PowerLawInteraction()
PowerLawInteraction ()
    {
        // we set beta = 1.  scale density accordingly!
        gamma = 1.;
        // sr_lr_split = cutoff for a truncated IPL.  otherwise, it's the
        // handshake distance between probing events and direct enumeration.
        sr_lr_split = 5.8;

        // built-in self-test, will leave exponent at 32 at the end
        double rsq_list[] = { 1e-6, 1e-2, 1., 1e2, 1e6 };
        double ex_list[] = { 1024., 6., 0., .5, /*-.5,*/ 1., 32. };
        for (double rsq : rsq_list)
        for (double ex : ex_list)
        {
            set_parameter ("exponent", ex);
            double delta = repulsive_lift (rsq, 1e-2);
            if (! (delta <= 0.))
            {
                std::cerr << "ipl: self-test error\n" << ABORT;
            }
        }
    }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xf0, %rsp
movq %rdi, -0x8(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0xe8(%rbp)
movsd 0x1e3726(%rip), %xmm0   # 0x5fb018
movsd %xmm0, (%rax)
movsd 0x1e49ea(%rip), %xmm0   # 0x5fc2e8
movsd %xmm0, 0x18(%rax)
leaq -0x30(%rbp), %rdi
xorl %esi, %esi
movl $0x28, %edx
callq 0x401080
movsd 0x1e496d(%rip), %xmm0   # 0x5fc288
movsd %xmm0, -0x30(%rbp)
movsd 0x1e36e8(%rip), %xmm0   # 0x5fb010
movsd %xmm0, -0x28(%rbp)
movsd 0x1e36e3(%rip), %xmm0   # 0x5fb018
movsd %xmm0, -0x20(%rbp)
movsd 0x1e36c6(%rip), %xmm0   # 0x5fb008
movsd %xmm0, -0x18(%rbp)
movsd 0x1e4949(%rip), %xmm0   # 0x5fc298
movsd %xmm0, -0x10(%rbp)
leaq -0x60(%rbp), %rdi
xorl %esi, %esi
movl $0x30, %edx
callq 0x401080
movsd 0x1e4974(%rip), %xmm0   # 0x5fc2e0
movsd %xmm0, -0x60(%rbp)
movsd 0x1e495f(%rip), %xmm0   # 0x5fc2d8
movsd %xmm0, -0x58(%rbp)
movsd 0x1e491a(%rip), %xmm0   # 0x5fc2a0
movsd %xmm0, -0x48(%rbp)
movsd 0x1e3685(%rip), %xmm0   # 0x5fb018
movsd %xmm0, -0x40(%rbp)
movsd 0x1e4930(%rip), %xmm0   # 0x5fc2d0
movsd %xmm0, -0x38(%rbp)
leaq -0x30(%rbp), %rax
movq %rax, -0x68(%rbp)
movq -0x68(%rbp), %rax
movq %rax, -0x70(%rbp)
movq -0x68(%rbp), %rax
addq $0x28, %rax
movq %rax, -0x78(%rbp)
movq -0x70(%rbp), %rax
cmpq -0x78(%rbp), %rax
je 0x417b4e
movq -0x70(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0x80(%rbp)
leaq -0x60(%rbp), %rax
movq %rax, -0x88(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0x90(%rbp)
movq -0x88(%rbp), %rax
addq $0x30, %rax
movq %rax, -0x98(%rbp)
movq -0x90(%rbp), %rax
cmpq -0x98(%rbp), %rax
je 0x417b3b
movq -0x90(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0xa0(%rbp)
leaq -0xc1(%rbp), %rdi
movq %rdi, -0xf0(%rbp)
callq 0x4771f0
movq -0xf0(%rbp), %rdx
leaq 0x1e49b8(%rip), %rsi    # 0x5fc407
leaq -0xc0(%rbp), %rdi
callq 0x407bd0
jmp 0x417a5d
movq -0xe8(%rbp), %rdi
movsd -0xa0(%rbp), %xmm0
leaq -0xc0(%rbp), %rsi
callq 0x417b70
jmp 0x417a7a
leaq -0xc0(%rbp), %rdi
callq 0x4f1e40
leaq -0xc1(%rbp), %rdi
callq 0x477210
movq -0xe8(%rbp), %rdi
movsd -0x80(%rbp), %xmm0
movsd 0x1e356a(%rip), %xmm1   # 0x5fb010
callq 0x417db0
movsd %xmm0, -0xe0(%rbp)
xorps %xmm0, %xmm0
ucomisd -0xe0(%rbp), %xmm0
jae 0x417b22
movq $0x698ac0, %rdi         # imm = 0x698AC0
leaq 0x1e4942(%rip), %rsi    # 0x5fc410
callq 0x4e1b60
movq %rax, %rdi
leaq 0x27fb44(%rip), %rsi    # 0x697621
callq 0x46d700
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0xd0(%rbp)
movl %eax, -0xd4(%rbp)
jmp 0x417b14
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0xd0(%rbp)
movl %eax, -0xd4(%rbp)
leaq -0xc0(%rbp), %rdi
callq 0x4f1e40
leaq -0xc1(%rbp), %rdi
callq 0x477210
jmp 0x417b57
jmp 0x417b24
movq -0x90(%rbp), %rax
addq $0x8, %rax
movq %rax, -0x90(%rbp)
jmp 0x417a07
jmp 0x417b3d
movq -0x70(%rbp), %rax
addq $0x8, %rax
movq %rax, -0x70(%rbp)
jmp 0x4179c1
addq $0xf0, %rsp
popq %rbp
retq
movq -0xd0(%rbp), %rdi
callq 0x5283c0
nopw %cs:(%rax,%rax)
nopl (%rax)