// File: ekg[P]sgd2/src/sgd2.hpp
// Function: sgd2::schedule_convergent(std::vector<sgd2::term, std::allocator<sgd2::term>> const&, unsigned long, double, unsigned long)
std::vector<double> schedule_convergent(const std::vector<term> &terms, uint64_t t_max, double eps, uint64_t t_maxmax)
{
    double w_min = terms[0].w, w_max = terms[0].w;
    for (uint64_t i=1; i<terms.size(); i++)
    {
        double w = terms[i].w;
        if (w < w_min) w_min = w;
        if (w > w_max) w_max = w;
    }
    double eta_max = 1.0 / w_min;
    double eta_min = eps / w_max;

    double lambda = log(eta_max/eta_min) / ((double)t_max-1);

    // initialize step sizes
    std::vector<double> etas;
    etas.reserve(t_maxmax);
    double eta_switch = 1.0 / w_max;
    for (uint64_t t=0; t<t_maxmax; t++)
    {
        double eta = eta_max * exp(-lambda * t);
        if (eta < eta_switch)
            break;

        etas.push_back(eta);
    }
    uint64_t tau = etas.size();
    for (uint64_t t=tau; t<t_maxmax; t++)
    {
        double eta = eta_switch / (1 + lambda*((double)t-tau));
        etas.push_back(eta);
    }

    return etas;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xd0, %rsp
movq %rdi, -0xb8(%rbp)
movq %rdi, %rax
movq %rax, -0xb0(%rbp)
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq %rdx, -0x18(%rbp)
movsd %xmm0, -0x20(%rbp)
movq %rcx, -0x28(%rbp)
movq -0x10(%rbp), %rdi
xorl %eax, %eax
movl %eax, %esi
callq 0xab00
movsd 0x18(%rax), %xmm0
movsd %xmm0, -0x30(%rbp)
movq -0x10(%rbp), %rdi
xorl %eax, %eax
movl %eax, %esi
callq 0xab00
movsd 0x18(%rax), %xmm0
movsd %xmm0, -0x38(%rbp)
movq $0x1, -0x40(%rbp)
movq -0x40(%rbp), %rax
movq %rax, -0xc0(%rbp)
movq -0x10(%rbp), %rdi
callq 0x9f50
movq %rax, %rcx
movq -0xc0(%rbp), %rax
cmpq %rcx, %rax
jae 0x3e51
movq -0x10(%rbp), %rdi
movq -0x40(%rbp), %rsi
callq 0xab00
movsd 0x18(%rax), %xmm0
movsd %xmm0, -0x48(%rbp)
movsd -0x48(%rbp), %xmm1
movsd -0x30(%rbp), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x3e2b
movsd -0x48(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
movsd -0x48(%rbp), %xmm0
ucomisd -0x38(%rbp), %xmm0
jbe 0x3e41
movsd -0x48(%rbp), %xmm0
movsd %xmm0, -0x38(%rbp)
jmp 0x3e43
movq -0x40(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x40(%rbp)
jmp 0x3dd7
movsd -0x30(%rbp), %xmm1
movsd 0x1a1aa(%rip), %xmm0    # 0x1e008
movsd %xmm0, -0xc8(%rbp)
divsd %xmm1, %xmm0
movsd %xmm0, -0x50(%rbp)
movsd -0x20(%rbp), %xmm0
movsd -0x38(%rbp), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x58(%rbp)
movsd -0x50(%rbp), %xmm0
movsd -0x58(%rbp), %xmm1
divsd %xmm1, %xmm0
callq 0x2130
movsd -0xc8(%rbp), %xmm2
movq -0xb8(%rbp), %rdi
movaps 0x1a195(%rip), %xmm3    # 0x1e040
movq -0x18(%rbp), %xmm1
punpckldq %xmm3, %xmm1    # xmm1 = xmm1[0],xmm3[0],xmm1[1],xmm3[1]
movapd 0x1a194(%rip), %xmm3    # 0x1e050
subpd %xmm3, %xmm1
movaps %xmm1, %xmm3
unpckhpd %xmm1, %xmm1            # xmm1 = xmm1[1,1]
addsd %xmm3, %xmm1
subsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x60(%rbp)
movb $0x0, -0x61(%rbp)
callq 0xab20
movq -0xb8(%rbp), %rdi
movq -0x28(%rbp), %rsi
callq 0xab40
jmp 0x3ef3
movsd 0x1a10d(%rip), %xmm0    # 0x1e008
divsd -0x38(%rbp), %xmm0
movsd %xmm0, -0x80(%rbp)
movq $0x0, -0x88(%rbp)
movq -0x88(%rbp), %rax
cmpq -0x28(%rbp), %rax
jae 0x3fe8
movsd -0x50(%rbp), %xmm0
movsd %xmm0, -0xd0(%rbp)
movaps 0x1a12b(%rip), %xmm1    # 0x1e060
movsd -0x60(%rbp), %xmm0
pxor %xmm1, %xmm0
movaps 0x1a0fb(%rip), %xmm2    # 0x1e040
movq -0x88(%rbp), %xmm1
punpckldq %xmm2, %xmm1    # xmm1 = xmm1[0],xmm2[0],xmm1[1],xmm2[1]
movapd 0x1a0f7(%rip), %xmm2    # 0x1e050
subpd %xmm2, %xmm1
movaps %xmm1, %xmm2
unpckhpd %xmm1, %xmm1            # xmm1 = xmm1[1,1]
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
callq 0x22c0
movaps %xmm0, %xmm1
movsd -0xd0(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x90(%rbp)
movsd -0x90(%rbp), %xmm1
movsd -0x80(%rbp), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x3fba
jmp 0x3fe8
movq -0xb8(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x70(%rbp)
movl %eax, -0x74(%rbp)
callq 0xaa80
jmp 0x40dd
movq -0xb8(%rbp), %rdi
leaq -0x90(%rbp), %rsi
callq 0xac70
jmp 0x3fcf
jmp 0x3fd1
movq -0x88(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x88(%rbp)
jmp 0x3f10
movq -0xb8(%rbp), %rdi
callq 0xace0
movq %rax, -0x98(%rbp)
movq -0x98(%rbp), %rax
movq %rax, -0xa0(%rbp)
movq -0xa0(%rbp), %rax
cmpq -0x28(%rbp), %rax
jae 0x40b7
movq -0xb8(%rbp), %rdi
movsd -0x80(%rbp), %xmm0
movsd -0x60(%rbp), %xmm1
movq -0xa0(%rbp), %xmm2
movaps 0x1a006(%rip), %xmm5    # 0x1e040
punpckldq %xmm5, %xmm2    # xmm2 = xmm2[0],xmm5[0],xmm2[1],xmm5[1]
movapd 0x1a00a(%rip), %xmm4    # 0x1e050
subpd %xmm4, %xmm2
movaps %xmm2, %xmm3
unpckhpd %xmm2, %xmm2            # xmm2 = xmm2[1,1]
addsd %xmm3, %xmm2
movq -0x98(%rbp), %xmm3
punpckldq %xmm5, %xmm3    # xmm3 = xmm3[0],xmm5[0],xmm3[1],xmm5[1]
subpd %xmm4, %xmm3
movaps %xmm3, %xmm4
unpckhpd %xmm3, %xmm3            # xmm3 = xmm3[1,1]
addsd %xmm4, %xmm3
subsd %xmm3, %xmm2
mulsd %xmm2, %xmm1
movsd 0x19f88(%rip), %xmm2    # 0x1e008
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0xa8(%rbp)
leaq -0xa8(%rbp), %rsi
callq 0xac70
jmp 0x409e
jmp 0x40a0
movq -0xa0(%rbp), %rax
addq $0x1, %rax
movq %rax, -0xa0(%rbp)
jmp 0x4009
movb $0x1, -0x61(%rbp)
testb $0x1, -0x61(%rbp)
jne 0x40cd
movq -0xb8(%rbp), %rdi
callq 0xaa80
movq -0xb0(%rbp), %rax
addq $0xd0, %rsp
popq %rbp
retq
movq -0x70(%rbp), %rdi
callq 0x22a0
nopw %cs:(%rax,%rax)