// File: /workspace/github_c_makefile_star2_O0/PoulinV[P]AxiCLASS/build/../external/HMcode/hmcode.c
// Function: hmcode_nowiggle_init
int hmcode_nowiggle_init(
                         struct fourier *pfo,
                         double **lnpk_l,
                         double **ddlnpk_l,
                         int index_pk,
                         struct hmcode_workspace *phw
                         ) {

  int i;
  double loganorm = log(1/(2.*_PI_*_PI_));
  double * pk_nw;

  /* Temporary array to store P_nowiggle (will not be used here) */

  class_alloc(pk_nw, pfo->nk_wiggle*sizeof(double), pfo->error_message);

  /* Get P_wiggle = P_l - P_nowiggle */

  class_call(hmcode_nowiggle(pfo,
                             lnpk_l[index_pk],
                             ddlnpk_l[index_pk],
                             pfo->nk_wiggle,
                             phw->lnk_wiggle,
                             pk_nw,
                             phw->pk_wiggle),
             pfo->error_message,
             pfo->error_message);

  /* Transform P_wiggle into a dimensionless linear spectrum {cal
     P}_wiggle = k^3/(2pi^2) P_l, where {cal P}_l would scale like
     k^(4+n_s) for k<k_eq. Store the result in phw->pk_wiggle. */

  for(i=0;i<pfo->nk_wiggle;++i){
    phw->pk_wiggle[i] = phw->pk_wiggle[i]*exp(3.*phw->lnk_wiggle[i]+loganorm);
  }

  /* Spline phw->pk_wiggle w.r.t. ln(k), to be able to
     interpolate later on. */

  class_call(array_spline_table_columns(phw->lnk_wiggle,
                                        pfo->nk_wiggle,
                                        phw->pk_wiggle,
                                        1,
                                        phw->ddpk_wiggle,
                                        _SPLINE_EST_DERIV_,
                                        pfo->error_message),
             pfo->error_message,
             pfo->error_message);

  free(pk_nw);

  return _SUCCESS_;
}

// Assembly:
subq $0x1868, %rsp           # imm = 0x1868
movq %rdi, 0x1858(%rsp)
movq %rsi, 0x1850(%rsp)
movq %rdx, 0x1848(%rsp)
movl %ecx, 0x1844(%rsp)
movq %r8, 0x1838(%rsp)
movsd 0x54b32(%rip), %xmm0    # 0x210368
callq 0x82f0
movsd %xmm0, 0x1828(%rsp)
movq 0x1858(%rsp), %rax
movslq 0x2c(%rax), %rdi
shlq $0x3, %rdi
callq 0x8460
movq %rax, 0x1820(%rsp)
cmpq $0x0, 0x1820(%rsp)
jne 0x1bb8e8
movq 0x1858(%rsp), %rax
movslq 0x2c(%rax), %rax
shlq $0x3, %rax
movl %eax, 0x181c(%rsp)
leaq 0x1010(%rsp), %rdi
movl 0x181c(%rsp), %ecx
leaq 0x1676b(%rip), %rsi     # 0x1d2004
leaq 0x567a3(%rip), %rdx     # 0x212043
movb $0x0, %al
callq 0x5b9c0
movq 0x1858(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x1010(%rsp), %r8
leaq 0x1676f(%rip), %rsi     # 0x1d2034
leaq 0x56655(%rip), %rdx     # 0x211f21
movl $0x770, %ecx            # imm = 0x770
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x1864(%rsp)
jmp 0x1bbb43
movq 0x1858(%rsp), %rdi
movq 0x1850(%rsp), %rax
movslq 0x1844(%rsp), %rcx
movq (%rax,%rcx,8), %rsi
movq 0x1848(%rsp), %rax
movslq 0x1844(%rsp), %rcx
movq (%rax,%rcx,8), %rdx
movq 0x1858(%rsp), %rax
movl 0x2c(%rax), %ecx
movq 0x1838(%rsp), %rax
movq 0x58(%rax), %r8
movq 0x1820(%rsp), %r9
movq 0x1838(%rsp), %rax
movq 0x60(%rax), %rax
movq %rax, (%rsp)
callq 0x1c02a0
cmpl $0x1, %eax
jne 0x1bb9be
leaq 0x810(%rsp), %rdi
movq 0x1858(%rsp), %rcx
addq $0x192, %rcx            # imm = 0x192
leaq 0x16a69(%rip), %rsi     # 0x1d23d8
leaq 0x565c0(%rip), %rdx     # 0x211f36
movb $0x0, %al
callq 0x5b9c0
movq 0x1858(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x810(%rsp), %r8
leaq 0x16699(%rip), %rsi     # 0x1d2034
leaq 0x5657f(%rip), %rdx     # 0x211f21
movl $0x77c, %ecx            # imm = 0x77C
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x1864(%rsp)
jmp 0x1bbb43
movl $0x0, 0x1834(%rsp)
movl 0x1834(%rsp), %eax
movq 0x1858(%rsp), %rcx
cmpl 0x2c(%rcx), %eax
jge 0x1bba73
movq 0x1838(%rsp), %rax
movq 0x60(%rax), %rax
movslq 0x1834(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, 0x8(%rsp)
movq 0x1838(%rsp), %rax
movq 0x58(%rax), %rax
movslq 0x1834(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x1828(%rsp), %xmm1
movsd 0x17c6e(%rip), %xmm2    # 0x1d3698
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x8670
movaps %xmm0, %xmm1
movsd 0x8(%rsp), %xmm0
mulsd %xmm1, %xmm0
movq 0x1838(%rsp), %rax
movq 0x60(%rax), %rax
movslq 0x1834(%rsp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl 0x1834(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x1834(%rsp)
jmp 0x1bb9c9
movq 0x1838(%rsp), %rax
movq 0x58(%rax), %rdi
movq 0x1858(%rsp), %rax
movl 0x2c(%rax), %esi
movq 0x1838(%rsp), %rax
movq 0x60(%rax), %rdx
movq 0x1838(%rsp), %rax
movq 0x68(%rax), %r8
movq 0x1858(%rsp), %rax
addq $0x192, %rax            # imm = 0x192
movl $0x1, %r9d
movl %r9d, %ecx
movq %rax, (%rsp)
callq 0x1f4b0
cmpl $0x1, %eax
jne 0x1bbb2b
leaq 0x10(%rsp), %rdi
movq 0x1858(%rsp), %rcx
addq $0x192, %rcx            # imm = 0x192
leaq 0x168f6(%rip), %rsi     # 0x1d23d8
leaq 0x564c0(%rip), %rdx     # 0x211fa9
movb $0x0, %al
callq 0x5b9c0
movq 0x1858(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x10(%rsp), %r8
leaq 0x16529(%rip), %rsi     # 0x1d2034
leaq 0x5640f(%rip), %rdx     # 0x211f21
movl $0x791, %ecx            # imm = 0x791
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x1864(%rsp)
jmp 0x1bbb43
movq 0x1820(%rsp), %rdi
callq 0x8420
movl $0x0, 0x1864(%rsp)
movl 0x1864(%rsp), %eax
addq $0x1868, %rsp           # imm = 0x1868
retq
nopw %cs:(%rax,%rax)