// File: mouda[P]LTE_random_access-simulate/alglibinternal.cpp
// Function: alglib_impl::taskgenint1dcheb2(double, double, long, alglib_impl::ae_vector*, alglib_impl::ae_vector*, alglib_impl::ae_state*)
void taskgenint1dcheb2(double a,
     double b,
     ae_int_t n,
     /* Real    */ ae_vector* x,
     /* Real    */ ae_vector* y,
     ae_state *_state)
{
    ae_int_t i;

    ae_vector_clear(x);
    ae_vector_clear(y);

    ae_assert(n>=1, "TaskGenInterpolation1DCheb2: N<1!", _state);
    ae_vector_set_length(x, n, _state);
    ae_vector_set_length(y, n, _state);
    if( n>1 )
    {
        for(i=0; i<=n-1; i++)
        {
            x->ptr.p_double[i] = 0.5*(b+a)+0.5*(b-a)*ae_cos(ae_pi*i/(n-1), _state);
            if( i==0 )
            {
                y->ptr.p_double[i] = 2*ae_randomreal(_state)-1;
            }
            else
            {
                y->ptr.p_double[i] = y->ptr.p_double[i-1]+(2*ae_randomreal(_state)-1)*(x->ptr.p_double[i]-x->ptr.p_double[i-1]);
            }
        }
    }
    else
    {
        x->ptr.p_double[0] = 0.5*(a+b);
        y->ptr.p_double[0] = 2*ae_randomreal(_state)-1;
    }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x50, %rsp
movsd %xmm0, -0x8(%rbp)
movsd %xmm1, -0x10(%rbp)
movq %rdi, -0x18(%rbp)
movq %rsi, -0x20(%rbp)
movq %rdx, -0x28(%rbp)
movq %rcx, -0x30(%rbp)
movq -0x20(%rbp), %rdi
callq 0x3b6f0
movq -0x28(%rbp), %rdi
callq 0x3b6f0
cmpq $0x1, -0x18(%rbp)
setge %al
movq -0x30(%rbp), %rdx
movzbl %al, %edi
andl $0x1, %edi
leaq 0x39ac8(%rip), %rsi     # 0x579c5
callq 0x3b240
movq -0x20(%rbp), %rdi
movq -0x18(%rbp), %rsi
movq -0x30(%rbp), %rdx
callq 0x3b620
movq -0x28(%rbp), %rdi
movq -0x18(%rbp), %rsi
movq -0x30(%rbp), %rdx
callq 0x3b620
cmpq $0x1, -0x18(%rbp)
jle 0x1e093
movq $0x0, -0x38(%rbp)
movq -0x38(%rbp), %rax
movq -0x18(%rbp), %rcx
subq $0x1, %rcx
cmpq %rcx, %rax
jg 0x1e091
movsd -0x10(%rbp), %xmm0
addsd -0x8(%rbp), %xmm0
movsd %xmm0, -0x40(%rbp)
movsd -0x10(%rbp), %xmm1
subsd -0x8(%rbp), %xmm1
movsd 0x374e3(%rip), %xmm0    # 0x55450
mulsd %xmm1, %xmm0
movsd %xmm0, -0x48(%rbp)
cvtsi2sdq -0x38(%rbp), %xmm1
movsd 0x374d4(%rip), %xmm0    # 0x55458
mulsd %xmm1, %xmm0
movq -0x18(%rbp), %rax
subq $0x1, %rax
cvtsi2sd %rax, %xmm1
divsd %xmm1, %xmm0
movq -0x30(%rbp), %rdi
callq 0x3cbc0
movsd -0x48(%rbp), %xmm1
movaps %xmm0, %xmm2
movsd -0x40(%rbp), %xmm0
mulsd %xmm2, %xmm1
movsd 0x37495(%rip), %xmm2    # 0x55450
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq -0x20(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x38(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
cmpq $0x0, -0x38(%rbp)
jne 0x1e007
movq -0x30(%rbp), %rdi
callq 0x3caf0
addsd %xmm0, %xmm0
movsd 0x371e0(%rip), %xmm1    # 0x551d0
addsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x38(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0x1e07e
movq -0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x38(%rbp), %rcx
subq $0x1, %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, -0x50(%rbp)
movq -0x30(%rbp), %rdi
callq 0x3caf0
movsd -0x50(%rbp), %xmm1
addsd %xmm0, %xmm0
movsd 0x37195(%rip), %xmm2    # 0x551d0
addsd %xmm2, %xmm0
movq -0x20(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x38(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm2
movq -0x20(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x38(%rbp), %rcx
subq $0x1, %rcx
subsd (%rax,%rcx,8), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x38(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0x1e080
movq -0x38(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x38(%rbp)
jmp 0x1df37
jmp 0x1e0da
movsd -0x8(%rbp), %xmm1
addsd -0x10(%rbp), %xmm1
movsd 0x373ab(%rip), %xmm0    # 0x55450
mulsd %xmm1, %xmm0
movq -0x20(%rbp), %rax
movq 0x28(%rax), %rax
movsd %xmm0, (%rax)
movq -0x30(%rbp), %rdi
callq 0x3caf0
addsd %xmm0, %xmm0
movsd 0x37106(%rip), %xmm1    # 0x551d0
addsd %xmm1, %xmm0
movq -0x28(%rbp), %rax
movq 0x28(%rax), %rax
movsd %xmm0, (%rax)
addq $0x50, %rsp
popq %rbp
retq