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

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

	double la = Tilde::la(a);
	double lb = Tilde::lb(a);
	double lc = Tilde::lc(a);

	// 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, order) * std::complex<double>(cos(la * L2E), -sin(la * L2E));
	A += Talphabetai(alpha, beta, 1, a, delta, order) * std::complex<double>(cos(lb * L2E), -sin(lb * L2E));
	A += Talphabetai(alpha, beta, 2, a, delta, order) * std::complex<double>(cos(lc * L2E), -sin(lc * L2E));

	double P = std::norm(A);
	ProbabilityWarning W(&P);
	return P;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0xa0, %rsp
cmpl $0x3, %edi
jae 0xa9cb
movl %esi, %ebp
cmpl $0x3, %esi
jae 0xa9e4
movl %edx, %ebx
cmpl $0x2, %edx
jae 0xa9fd
movl %edi, %r14d
leaq 0x1bdc(%rip), %rax      # 0xc328
mulsd (%rax), %xmm1
mulsd 0x18b0(%rip), %xmm1     # 0xc008
movsd %xmm1, 0x20(%rsp)
leaq 0x1b43(%rip), %rax      # 0xc2a8
leaq 0x6a14(%rip), %rcx      # 0x11180
movsd (%rcx), %xmm1
leaq 0x1b21(%rip), %rcx      # 0xc298
movsd (%rcx), %xmm4
mulsd %xmm1, %xmm4
movsd (%rax), %xmm5
addsd %xmm4, %xmm5
leaq 0x69ea(%rip), %rax      # 0x11178
movsd (%rax), %xmm3
mulsd %xmm3, %xmm5
addsd %xmm0, %xmm5
movsd %xmm5, 0x10(%rsp)
leaq 0x1b19(%rip), %rax      # 0xc2c0
mulsd (%rax), %xmm1
mulsd %xmm3, %xmm1
movsd %xmm1, 0x18(%rsp)
leaq 0x1b14(%rip), %rax      # 0xc2d0
addsd (%rax), %xmm4
mulsd %xmm3, %xmm4
movsd %xmm4, 0x28(%rsp)
movl %ebp, %esi
xorl %edx, %edx
movapd %xmm2, %xmm1
movl %ebx, %ecx
movsd %xmm2, 0x38(%rsp)
movsd %xmm0, 0x30(%rsp)
callq 0xa5e0
leaq 0x58(%rsp), %r15
movsd %xmm0, (%r15)
movsd %xmm1, 0x8(%r15)
movsd 0x10(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
callq 0x2290
movsd %xmm0, 0x8(%rsp)
movsd 0x10(%rsp), %xmm0
callq 0x2060
xorpd 0x17fb(%rip), %xmm0     # 0xc020
leaq 0x48(%rsp), %r12
movsd 0x8(%rsp), %xmm1
movsd %xmm1, (%r12)
movlpd %xmm0, 0x8(%r12)
movq %r15, %rdi
movq %r12, %rsi
callq 0x3719
xorpd %xmm2, %xmm2
addsd %xmm2, %xmm0
movsd %xmm0, 0x8(%rsp)
addsd %xmm2, %xmm1
movsd %xmm1, 0x10(%rsp)
pushq $0x1
popq %rdx
movl %r14d, %edi
movl %ebp, %esi
movsd 0x30(%rsp), %xmm0
movsd 0x38(%rsp), %xmm1
movl %ebx, %ecx
callq 0xa5e0
movsd %xmm0, (%r15)
movsd %xmm1, 0x8(%r15)
movsd 0x18(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
callq 0x2290
movsd %xmm0, 0x40(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x2060
xorpd 0x176a(%rip), %xmm0     # 0xc020
movsd 0x40(%rsp), %xmm1
movsd %xmm1, (%r12)
movlpd %xmm0, 0x8(%r12)
leaq 0x58(%rsp), %r12
leaq 0x48(%rsp), %r15
movq %r12, %rdi
movq %r15, %rsi
callq 0x3719
addsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
addsd 0x10(%rsp), %xmm1
movsd %xmm1, 0x18(%rsp)
pushq $0x2
popq %rdx
movl %r14d, %edi
movl %ebp, %esi
movsd 0x30(%rsp), %xmm0
movsd 0x38(%rsp), %xmm1
movl %ebx, %ecx
callq 0xa5e0
movsd %xmm0, (%r12)
movsd %xmm1, 0x8(%r12)
movsd 0x28(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x28(%rsp)
callq 0x2290
movsd %xmm0, 0x20(%rsp)
movsd 0x28(%rsp), %xmm0
callq 0x2060
xorpd 0x16d2(%rip), %xmm0     # 0xc020
movsd 0x20(%rsp), %xmm1
movsd %xmm1, (%r15)
movlpd %xmm0, 0x8(%r15)
leaq 0x58(%rsp), %rdi
leaq 0x48(%rsp), %rbx
movq %rbx, %rsi
callq 0x3719
addsd 0x8(%rsp), %xmm0
addsd 0x18(%rsp), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, (%rbx)
leaq 0x58(%rsp), %rbx
leaq 0x48(%rsp), %r14
movq %rbx, %rdi
movq %r14, %rsi
callq 0xadac
movsd (%r14), %xmm0
movsd %xmm0, 0x28(%rsp)
movq %rbx, %rdi
callq 0xaa7a
movsd 0x28(%rsp), %xmm0
addq $0xa0, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
leaq 0x17ad(%rip), %rdi      # 0xc17f
leaq 0x2268(%rip), %rsi      # 0xcc41
leaq 0x2329(%rip), %rcx      # 0xcd09
pushq $0x55
jmp 0xaa14
leaq 0x17f9(%rip), %rdi      # 0xc1e4
leaq 0x224f(%rip), %rsi      # 0xcc41
leaq 0x2310(%rip), %rcx      # 0xcd09
pushq $0x56
jmp 0xaa14
leaq 0x1baa(%rip), %rdi      # 0xc5ae
leaq 0x2236(%rip), %rsi      # 0xcc41
leaq 0x22f7(%rip), %rcx      # 0xcd09
pushq $0x57
popq %rdx
callq 0x2100
nop