// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Hyperbolas.cpp
// Function: Hyperbolas::Palphabeta(int, int, double, double, double)
double Palphabeta(int alpha, int beta, double a, double LE, double delta)
{
	assert (alpha >= 0 && alpha <= 2); // alpha valid for 0, 1, 2
	assert (beta >= 0 && beta <= 2); // beta valid for 0, 1, 2

	double L2E = km_per_GeV_to_per_eV2 * LE / 2; // in eV^-2 now

	// See eq. 4.0.1 or eq. 3.12 in arXiv:1505.01826
	std::complex<double> A(0, 0);
	A += Talphabetai(alpha, beta, 0, a, delta) * std::complex<double>(cos(Hyperbolas::l1(a) * L2E), -sin(Hyperbolas::l1(a) * L2E));
	A += Talphabetai(alpha, beta, 1, a, delta) * std::complex<double>(cos(Hyperbolas::l2(a) * L2E), -sin(Hyperbolas::l2(a) * L2E));
	A += Talphabetai(alpha, beta, 2, a, delta) * std::complex<double>(cos(Hyperbolas::l3(a) * L2E), -sin(Hyperbolas::l3(a) * L2E));

	return std::norm(A);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x68, %rsp
movsd %xmm0, (%rsp)
cmpl $0x3, %edi
jae 0x9cdf
movl %esi, %ebx
cmpl $0x3, %esi
jae 0x9cfe
movl %edi, %ebp
leaq 0x4a98(%rip), %rax      # 0xe328
mulsd (%rax), %xmm1
mulsd 0x476c(%rip), %xmm1     # 0xe008
movsd %xmm1, 0x20(%rsp)
movl %ebx, %esi
xorl %edx, %edx
movsd (%rsp), %xmm0
movsd %xmm2, 0x38(%rsp)
movaps %xmm2, %xmm1
callq 0x9721
movsd %xmm0, 0x30(%rsp)
movsd %xmm1, 0x28(%rsp)
leaq 0x4a3c(%rip), %rax      # 0xe308
movsd (%rax), %xmm0
addsd %xmm0, %xmm0
movsd %xmm0, 0x18(%rsp)
callq 0x22c0
leaq 0x98aa(%rip), %r13      # 0x13190
mulsd (%r13), %xmm0
leaq 0x9895(%rip), %r12      # 0x13188
mulsd (%r12), %xmm0
movsd (%rsp), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
mulsd %xmm1, %xmm1
movsd %xmm1, 0x10(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x2070
mulsd (%r13), %xmm0
mulsd (%r12), %xmm0
mulsd %xmm0, %xmm0
addsd 0x10(%rsp), %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x9943
sqrtsd %xmm0, %xmm0
jmp 0x9948
callq 0x2120
movsd 0x8(%rsp), %xmm1
subsd %xmm0, %xmm1
mulsd 0x46ae(%rip), %xmm1     # 0xe008
addsd 0x984e(%rip), %xmm1     # 0x131b0
mulsd 0x20(%rsp), %xmm1
movapd %xmm1, %xmm0
callq 0x22c0
movsd %xmm0, 0x8(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x22c0
mulsd (%r13), %xmm0
mulsd (%r12), %xmm0
movsd (%rsp), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rsp)
mulsd %xmm1, %xmm1
movsd %xmm1, 0x40(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x2070
mulsd (%r13), %xmm0
mulsd (%r12), %xmm0
mulsd %xmm0, %xmm0
addsd 0x40(%rsp), %xmm0
ucomisd 0x4988(%rip), %xmm0     # 0xe358
jb 0x99d8
sqrtsd %xmm0, %xmm0
jmp 0x99dd
callq 0x2120
movsd 0x10(%rsp), %xmm1
subsd %xmm0, %xmm1
mulsd 0x4619(%rip), %xmm1     # 0xe008
addsd 0x97b9(%rip), %xmm1     # 0x131b0
mulsd 0x20(%rsp), %xmm1
movapd %xmm1, %xmm0
callq 0x2070
xorpd 0x4612(%rip), %xmm0     # 0xe020
leaq 0x48(%rsp), %r14
movsd 0x8(%rsp), %xmm1
movsd %xmm1, (%r14)
movlpd %xmm0, 0x8(%r14)
leaq 0x58(%rsp), %r15
movsd 0x30(%rsp), %xmm0
movsd %xmm0, (%r15)
movsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x8(%r15)
movq %r15, %rdi
movq %r14, %rsi
callq 0x3994
xorpd %xmm0, %xmm0
movsd (%r15), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x18(%rsp)
movsd 0x8(%r15), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x30(%rsp)
movl %ebp, %edi
movl %ebx, %esi
movl $0x1, %edx
movsd (%rsp), %xmm0
movsd 0x38(%rsp), %xmm1
callq 0x9721
movsd %xmm0, 0x28(%rsp)
movsd %xmm1, 0x8(%rsp)
movsd (%rsp), %xmm0
callq 0x83c3
mulsd 0x20(%rsp), %xmm0
callq 0x22c0
movsd %xmm0, 0x10(%rsp)
movsd (%rsp), %xmm0
callq 0x83c3
mulsd 0x20(%rsp), %xmm0
callq 0x2070
xorpd 0x4555(%rip), %xmm0     # 0xe020
movsd 0x10(%rsp), %xmm1
movsd %xmm1, (%r14)
movlpd %xmm0, 0x8(%r14)
movsd 0x28(%rsp), %xmm0
movsd %xmm0, (%r15)
movsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x8(%r15)
leaq 0x58(%rsp), %r14
leaq 0x48(%rsp), %rsi
movq %r14, %rdi
callq 0x3994
movsd 0x18(%rsp), %xmm0
addsd (%r14), %xmm0
movsd %xmm0, 0x18(%rsp)
movsd 0x30(%rsp), %xmm0
addsd 0x8(%r14), %xmm0
movsd %xmm0, 0x30(%rsp)
movl %ebp, %edi
movl %ebx, %esi
movl $0x2, %edx
movsd (%rsp), %xmm0
movsd 0x38(%rsp), %xmm1
callq 0x9721
movsd %xmm0, 0x38(%rsp)
leaq 0x475a(%rip), %rax      # 0xe2a8
movsd (%rax), %xmm3
leaq 0x473f(%rip), %rax      # 0xe298
movsd (%rax), %xmm0
movsd (%r13), %xmm2
movsd %xmm0, 0x8(%rsp)
mulsd %xmm0, %xmm2
movsd %xmm3, 0x10(%rsp)
addsd %xmm3, %xmm2
mulsd (%r12), %xmm2
movsd %xmm1, 0x28(%rsp)
addsd (%rsp), %xmm2
movsd 0x9628(%rip), %xmm0     # 0x131b8
movapd %xmm2, %xmm1
addsd %xmm0, %xmm1
subsd %xmm0, %xmm2
mulsd %xmm2, %xmm2
movsd 0x9618(%rip), %xmm0     # 0x131c0
mulsd %xmm0, %xmm0
addsd %xmm2, %xmm0
ucomisd 0x47a0(%rip), %xmm0     # 0xe358
jb 0x9bc0
sqrtsd %xmm0, %xmm0
jmp 0x9bd1
movsd %xmm1, 0x40(%rsp)
callq 0x2120
movsd 0x40(%rsp), %xmm1
addsd %xmm0, %xmm1
mulsd 0x442b(%rip), %xmm1     # 0xe008
mulsd 0x20(%rsp), %xmm1
movapd %xmm1, %xmm0
callq 0x22c0
movsd 0x8(%rsp), %xmm1
mulsd (%r13), %xmm1
addsd 0x10(%rsp), %xmm1
mulsd (%r12), %xmm1
movsd %xmm0, 0x8(%rsp)
addsd (%rsp), %xmm1
movsd 0x95a1(%rip), %xmm0     # 0x131b8
movapd %xmm1, %xmm2
addsd %xmm0, %xmm2
subsd %xmm0, %xmm1
mulsd %xmm1, %xmm1
movsd 0x9591(%rip), %xmm0     # 0x131c0
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
ucomisd 0x4719(%rip), %xmm0     # 0xe358
jb 0x9c47
sqrtsd %xmm0, %xmm0
jmp 0x9c56
movsd %xmm2, (%rsp)
callq 0x2120
movsd (%rsp), %xmm2
addsd %xmm0, %xmm2
mulsd 0x43a6(%rip), %xmm2     # 0xe008
movsd 0x20(%rsp), %xmm0
mulsd %xmm2, %xmm0
callq 0x2070
xorpd 0x43a7(%rip), %xmm0     # 0xe020
leaq 0x48(%rsp), %rsi
movsd 0x8(%rsp), %xmm1
movsd %xmm1, (%rsi)
movlpd %xmm0, 0x8(%rsi)
leaq 0x58(%rsp), %rbx
movsd 0x38(%rsp), %xmm0
movsd %xmm0, (%rbx)
movsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x8(%rbx)
movq %rbx, %rdi
callq 0x3994
movsd 0x18(%rsp), %xmm0
addsd (%rbx), %xmm0
movsd 0x30(%rsp), %xmm1
addsd 0x8(%rbx), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
addq $0x68, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x4499(%rip), %rdi      # 0xe17f
leaq 0x4ae1(%rip), %rsi      # 0xe7ce
leaq 0x4b39(%rip), %rcx      # 0xe82d
movl $0xb1, %edx
callq 0x20f0
leaq 0x44df(%rip), %rdi      # 0xe1e4
leaq 0x4ac2(%rip), %rsi      # 0xe7ce
leaq 0x4b1a(%rip), %rcx      # 0xe82d
movl $0xb2, %edx
callq 0x20f0
nopl (%rax)