// File: /workspace/github_c_makefile_star2_O1/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: hyperspherical_WKB
int hyperspherical_WKB(int K,int l,double beta,double y, double *Phi){
  double e, w, w2, alpha, alpha2, CscK, ytp, t;
  double S, Q, C, argu, Ai;
  int airy_sign = 1, phisign = 1, dphisign = 1, intbeta;
  double ldbl = l;

  if (K==1){
    //Limit range to [0; pi/2]:
    intbeta = (int)(beta+0.4); //Round to nearest integer (just to be sure)
    ClosedModY(l, intbeta, &y, &phisign, &dphisign);
  }
  e = 1.0/sqrt(ldbl*(ldbl+1.0));
  alpha = beta*e;
  if (K==-1){
    CscK = 1.0/sinh(y);
    ytp = asinh(1.0/alpha);
  }
  else if (K==1){
    CscK = 1.0/sin(y);
    ytp = asin(1.0/alpha);
  }
  else{
    return _FAILURE_;
  }
  w = alpha/CscK;
  w2 = w*w;
  alpha2 = alpha*alpha;
  if (K==-1){
    if (y>ytp){
      S = alpha*log((sqrt(w2-1.0)+sqrt(w2+alpha2))/sqrt(1.0+alpha2))+
        atan(1.0/alpha*sqrt((w2+alpha2)/(w2-1.0)))-0.5*_PI_;
      airy_sign = -1;
    }
    else{
      t = sqrt(1.0-w2)/sqrt(1.0+w2/alpha2);
      S = atanh(t)-alpha*atan(t/alpha);
      airy_sign = 1;
    }
  }
  else if (K==1){
    if (y>ytp){
      t = sqrt(1-w2/alpha2)/sqrt(w2-1.0);
        S = atan(t)+alpha*atan(1.0/(t*alpha))-0.5*_PI_;
        airy_sign = -1;
    }
    else{
      S = atanh(sqrt(1.0-w2)/sqrt(1.0-w2/alpha2))-
        alpha*log((sqrt(alpha2-w2)+sqrt(1.0-w2))/sqrt(alpha2-1.0));
      airy_sign = 1;
    }
  }
  else{
    return _FAILURE_;
  }
  argu = 3.0*S/(2.0*e);
  Q = CscK*CscK-alpha2;
  C = 0.5*sqrt(alpha)/beta;
  Ai = airy_cheb_approx(airy_sign*pow(argu,2.0/3.0));
  *Phi = phisign*2.0*sqrt(_PI_)*C*pow(argu,1.0/6.0)*pow(fabs(Q),-0.25)*Ai*CscK;
  return _SUCCESS_;
}

// Assembly:
pushq %rbp
pushq %rbx
subq $0x68, %rsp
movq %rdx, %rbx
movsd %xmm0, 0x40(%rsp)
movl %edi, %ebp
movsd %xmm1, 0x48(%rsp)
movl $0x1, %eax
movl %eax, 0x50(%rsp)
movl %eax, 0x54(%rsp)
xorps %xmm1, %xmm1
cvtsi2sd %esi, %xmm1
cmpl %eax, %edi
jne 0x255cc
movsd 0xcef3a(%rip), %xmm0    # 0xf44d8
addsd 0x40(%rsp), %xmm0
cvttsd2si %xmm0, %eax
leaq 0x48(%rsp), %rdx
leaq 0x50(%rsp), %rcx
leaq 0x54(%rsp), %r8
movl %esi, %edi
movl %eax, %esi
movsd %xmm1, 0x20(%rsp)
callq 0x241c1
movsd 0x20(%rsp), %xmm1
movsd 0xcbc8c(%rip), %xmm0    # 0xf1260
addsd %xmm1, %xmm0
mulsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x255ec
sqrtsd %xmm0, %xmm0
jmp 0x255f1
callq 0x82b0
movsd 0xcbc67(%rip), %xmm1    # 0xf1260
divsd %xmm0, %xmm1
movsd %xmm1, 0x60(%rsp)
mulsd 0x40(%rsp), %xmm1
cmpl $-0x1, %ebp
je 0x25647
movl $0x1, %eax
cmpl $0x1, %ebp
jne 0x25b9b
movsd 0x48(%rsp), %xmm0
movsd %xmm1, 0x20(%rsp)
callq 0x80a0
movsd %xmm0, (%rsp)
movsd 0xcbc26(%rip), %xmm0    # 0xf1260
divsd 0x20(%rsp), %xmm0
callq 0x85a0
jmp 0x25670
movsd 0x48(%rsp), %xmm0
movsd %xmm1, 0x20(%rsp)
callq 0x8370
movsd %xmm0, (%rsp)
movsd 0xcbbfb(%rip), %xmm0    # 0xf1260
divsd 0x20(%rsp), %xmm0
callq 0x8470
movsd 0xcbbe8(%rip), %xmm4    # 0xf1260
movapd %xmm4, %xmm1
divsd (%rsp), %xmm1
movsd 0x20(%rsp), %xmm3
movapd %xmm3, %xmm5
movsd %xmm1, 0x58(%rsp)
divsd %xmm1, %xmm5
mulsd %xmm5, %xmm5
mulsd %xmm3, %xmm3
cmpl $-0x1, %ebp
movsd %xmm3, 0x10(%rsp)
je 0x256f1
movl $0x1, %eax
cmpl $0x1, %ebp
jne 0x25b9b
movsd 0x48(%rsp), %xmm1
ucomisd %xmm0, %xmm1
jbe 0x2574b
movapd %xmm5, %xmm1
divsd %xmm3, %xmm1
movsd 0xcbb8a(%rip), %xmm0    # 0xf1260
subsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x258f6
sqrtsd %xmm0, %xmm0
jmp 0x2590b
movsd 0x48(%rsp), %xmm1
ucomisd %xmm0, %xmm1
jbe 0x25722
movsd 0xcc0e3(%rip), %xmm0    # 0xf17e8
addsd %xmm5, %xmm0
xorpd %xmm2, %xmm2
ucomisd %xmm2, %xmm0
movsd %xmm0, 0x18(%rsp)
jb 0x2577c
xorps %xmm1, %xmm1
sqrtsd %xmm0, %xmm1
jmp 0x2579b
movsd 0xcbb36(%rip), %xmm1    # 0xf1260
movapd %xmm1, %xmm0
subsd %xmm5, %xmm0
xorpd %xmm2, %xmm2
ucomisd %xmm2, %xmm0
movsd %xmm1, 0x18(%rsp)
jb 0x257b4
sqrtsd %xmm0, %xmm0
jmp 0x257cf
movapd %xmm4, %xmm0
subsd %xmm5, %xmm0
xorpd %xmm6, %xmm6
ucomisd %xmm6, %xmm0
movsd %xmm4, 0x18(%rsp)
movsd %xmm5, 0x8(%rsp)
movsd %xmm0, 0x38(%rsp)
jb 0x25927
sqrtsd %xmm0, %xmm2
jmp 0x25948
movsd %xmm5, 0x8(%rsp)
callq 0x82b0
xorpd %xmm2, %xmm2
movsd 0x8(%rsp), %xmm5
movsd 0x10(%rsp), %xmm3
movapd %xmm0, %xmm1
addsd %xmm3, %xmm5
ucomisd %xmm2, %xmm5
movsd %xmm5, 0x8(%rsp)
jb 0x257f0
xorps %xmm0, %xmm0
sqrtsd %xmm5, %xmm0
jmp 0x25809
movsd %xmm5, 0x8(%rsp)
callq 0x82b0
xorpd %xmm2, %xmm2
movsd 0x8(%rsp), %xmm5
movsd 0x10(%rsp), %xmm3
movsd %xmm0, 0x38(%rsp)
divsd %xmm3, %xmm5
addsd 0xcba7f(%rip), %xmm5    # 0xf1260
ucomisd %xmm2, %xmm5
jb 0x25829
xorps %xmm0, %xmm0
sqrtsd %xmm5, %xmm0
jmp 0x25832
movapd %xmm5, %xmm0
movsd %xmm1, (%rsp)
callq 0x82b0
movsd (%rsp), %xmm1
movsd 0x10(%rsp), %xmm3
addsd %xmm0, %xmm1
movsd 0xcba4b(%rip), %xmm0    # 0xf1260
addsd %xmm3, %xmm0
xorpd %xmm2, %xmm2
ucomisd %xmm2, %xmm0
jb 0x25875
sqrtsd %xmm0, %xmm0
jmp 0x25884
movapd %xmm5, %xmm0
callq 0x82b0
movsd 0x38(%rsp), %xmm1
divsd %xmm0, %xmm1
movsd %xmm1, 0x38(%rsp)
movapd %xmm1, %xmm0
callq 0x8640
movsd %xmm0, (%rsp)
movsd 0x38(%rsp), %xmm0
divsd 0x20(%rsp), %xmm0
callq 0x82e0
movsd (%rsp), %xmm2
movsd 0x20(%rsp), %xmm1
mulsd %xmm1, %xmm0
jmp 0x25a9b
movsd %xmm1, (%rsp)
callq 0x82b0
movsd (%rsp), %xmm1
divsd %xmm0, %xmm1
movapd %xmm1, %xmm0
callq 0x82c0
movsd %xmm0, (%rsp)
movsd 0xcb9c2(%rip), %xmm1    # 0xf1260
divsd 0x20(%rsp), %xmm1
movsd 0x8(%rsp), %xmm0
divsd 0x18(%rsp), %xmm0
ucomisd 0xcb9a0(%rip), %xmm0    # 0xf1258
jb 0x258c0
sqrtsd %xmm0, %xmm0
jmp 0x258d1
movsd %xmm1, 0x18(%rsp)
callq 0x82b0
movsd 0x18(%rsp), %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, %xmm0
callq 0x82e0
movsd 0x20(%rsp), %xmm1
movsd (%rsp), %xmm2
mulsd %xmm1, %xmm2
addsd %xmm0, %xmm2
jmp 0x259b9
movsd %xmm5, 0x8(%rsp)
callq 0x82b0
xorpd %xmm1, %xmm1
movsd 0x8(%rsp), %xmm5
movsd %xmm0, (%rsp)
addsd 0xcbed0(%rip), %xmm5    # 0xf17e8
ucomisd %xmm1, %xmm5
jb 0x25963
xorps %xmm0, %xmm0
sqrtsd %xmm5, %xmm0
jmp 0x2596c
callq 0x82b0
xorpd %xmm6, %xmm6
movsd 0x8(%rsp), %xmm5
movsd 0xcb922(%rip), %xmm4    # 0xf1260
movsd 0x10(%rsp), %xmm3
movapd %xmm0, %xmm2
divsd %xmm3, %xmm5
movapd %xmm4, %xmm0
subsd %xmm5, %xmm0
ucomisd %xmm6, %xmm0
jb 0x259d4
sqrtsd %xmm0, %xmm0
jmp 0x259e3
movapd %xmm5, %xmm0
callq 0x82b0
movsd (%rsp), %xmm1
divsd %xmm0, %xmm1
movsd %xmm1, (%rsp)
movapd %xmm1, %xmm0
callq 0x82e0
movsd %xmm0, 0x18(%rsp)
movsd (%rsp), %xmm1
mulsd 0x20(%rsp), %xmm1
movsd 0xcb8c4(%rip), %xmm0    # 0xf1260
divsd %xmm1, %xmm0
callq 0x82e0
movsd 0x20(%rsp), %xmm1
movapd %xmm0, %xmm2
mulsd %xmm1, %xmm2
addsd 0x18(%rsp), %xmm2
addsd 0xceaff(%rip), %xmm2    # 0xf44c0
movsd 0xcbe1f(%rip), %xmm0    # 0xf17e8
movsd %xmm0, 0x18(%rsp)
jmp 0x25a9f
movsd %xmm2, (%rsp)
callq 0x82b0
movsd (%rsp), %xmm2
divsd %xmm0, %xmm2
movapd %xmm2, %xmm0
callq 0x8640
movapd %xmm0, %xmm1
movsd 0x10(%rsp), %xmm0
subsd 0x8(%rsp), %xmm0
xorpd %xmm3, %xmm3
ucomisd %xmm3, %xmm0
movsd %xmm1, (%rsp)
jb 0x25a18
xorps %xmm2, %xmm2
sqrtsd %xmm0, %xmm2
jmp 0x25a25
callq 0x82b0
xorpd %xmm3, %xmm3
movapd %xmm0, %xmm2
movsd 0x38(%rsp), %xmm0
ucomisd %xmm3, %xmm0
jb 0x25a37
sqrtsd %xmm0, %xmm0
jmp 0x25a48
movsd %xmm2, 0x8(%rsp)
callq 0x82b0
movsd 0x8(%rsp), %xmm2
movsd 0x10(%rsp), %xmm1
addsd %xmm0, %xmm2
movsd 0xcbd8e(%rip), %xmm0    # 0xf17e8
addsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x25a6e
sqrtsd %xmm0, %xmm0
jmp 0x25a7f
movsd %xmm2, 0x8(%rsp)
callq 0x82b0
movsd 0x8(%rsp), %xmm2
divsd %xmm0, %xmm2
movapd %xmm2, %xmm0
callq 0x82c0
movsd 0x20(%rsp), %xmm1
mulsd %xmm1, %xmm0
movsd (%rsp), %xmm2
subsd %xmm0, %xmm2
mulsd 0xccbd9(%rip), %xmm2    # 0xf2680
movsd 0x60(%rsp), %xmm0
addsd %xmm0, %xmm0
divsd %xmm0, %xmm2
movsd 0x58(%rsp), %xmm0
mulsd %xmm0, %xmm0
subsd 0x10(%rsp), %xmm0
movapd %xmm0, 0x20(%rsp)
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
movsd %xmm2, (%rsp)
jb 0x25ae3
xorps %xmm0, %xmm0
sqrtsd %xmm1, %xmm0
jmp 0x25af1
movapd %xmm1, %xmm0
callq 0x82b0
movsd (%rsp), %xmm2
mulsd 0xcc1d7(%rip), %xmm0    # 0xf1cd0
divsd 0x40(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0xce9d3(%rip), %xmm1    # 0xf44e0
movapd %xmm2, %xmm0
callq 0x8620
mulsd 0x18(%rsp), %xmm0
callq 0x25530
movsd %xmm0, 0x40(%rsp)
xorps %xmm0, %xmm0
cvtsi2sdl 0x50(%rsp), %xmm0
addsd %xmm0, %xmm0
mulsd 0xce9ac(%rip), %xmm0    # 0xf44e8
mulsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0xce978(%rip), %xmm1    # 0xf44c8
movsd (%rsp), %xmm0
callq 0x8620
mulsd 0x10(%rsp), %xmm0
movsd %xmm0, (%rsp)
movapd 0x20(%rsp), %xmm0
andpd 0xcb6cd(%rip), %xmm0    # 0xf1240
movsd 0xcb70d(%rip), %xmm1    # 0xf1288
callq 0x8620
mulsd (%rsp), %xmm0
mulsd 0x40(%rsp), %xmm0
movsd 0x58(%rsp), %xmm1
mulsd %xmm0, %xmm1
movsd %xmm1, (%rbx)
xorl %eax, %eax
addq $0x68, %rsp
popq %rbx
popq %rbp
retq