// File: /workspace/github_cpp_makefile_star2_O3/PeterDenton[P]Nu-Pert/src/GF.cpp
// Function: GF::D(int, int, double, double, int)
double D(int alpha, int Beta, double a, double delta, int order)
{
	assert (order >= 0 && order <= 1); // order valid for 0, 1
	if (alpha == Beta)
		return 0;
	double sdelta = sin(delta);
	double phi = Hat::phi(a);
	double cphi = cos(phi);
	double sphi = sin(phi);
	double cphisq = pow(cphi, 2);
	double sphisq = pow(sphi, 2);
	double psi = Check::psi(a);
	double cpsi = cos(psi);
	double spsi = sin(psi);
	double cpsisq = pow(cpsi, 2);
	double spsisq = pow(spsi, 2);
	double D0 = sdelta * Check::Jrm(a);
	double D1 = sdelta * c23 * s23 * cphi * (epsp(a) * Dmsqee * pow(spsi, 2) * (sphisq - cphisq * cpsisq) / Check::Dl31(a)
				- epsp(a) * Dmsqee * pow(cpsi, 2) * (sphisq - cphisq * spsisq) / Check::Dl32(a));
	double D = D0 + order * D1;
	if ((alpha + Beta) % 2 == 0)
		D *= -1;
	if (alpha > Beta)
		D *= -1;
	return D;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0xa8, %rsp
cmpl $0x2, %edx
jae 0x6645
movapd %xmm0, %xmm2
movl %esi, %ebx
movl %edi, %r14d
xorpd %xmm0, %xmm0
cmpl %esi, %edi
je 0x6637
movl %edx, %ebp
movaps %xmm1, %xmm0
movsd %xmm2, 0x8(%rsp)
callq 0x2070
movsd %xmm0, 0x18(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x6e4a
movsd %xmm0, 0x40(%rsp)
callq 0x22c0
movaps %xmm0, 0x20(%rsp)
movsd 0x40(%rsp), %xmm0
callq 0x2070
mulsd %xmm0, %xmm0
movapd %xmm0, 0x40(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x29c5
movsd %xmm0, 0x10(%rsp)
callq 0x22c0
movaps %xmm0, 0x90(%rsp)
movsd 0x10(%rsp), %xmm0
callq 0x2070
movaps %xmm0, 0x80(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x2cc9
movsd 0x18(%rsp), %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x10(%rsp)
leaq 0x7dff(%rip), %rax      # 0xe2f8
mulsd (%rax), %xmm1
leaq 0x7ddc(%rip), %rax      # 0xe2e0
mulsd (%rax), %xmm1
mulsd 0x20(%rsp), %xmm1
movsd %xmm1, 0x18(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x39d8
movaps %xmm0, 0x60(%rsp)
leaq 0xcc5d(%rip), %r15      # 0x13188
movapd 0x20(%rsp), %xmm0
movapd %xmm0, %xmm1
xorpd 0x7ae3(%rip), %xmm1     # 0xe020
mulsd %xmm0, %xmm1
movapd %xmm1, 0x70(%rsp)
movsd (%r15), %xmm0
movaps %xmm0, 0x30(%rsp)
movsd 0x8(%rsp), %xmm0
xorl %edi, %edi
callq 0x2a9f
movaps %xmm0, 0x20(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x39d8
movaps %xmm0, 0x50(%rsp)
movsd (%r15), %xmm0
unpcklpd 0x30(%rsp), %xmm0       # xmm0 = xmm0[0],mem[0]
movaps %xmm0, 0x30(%rsp)
movsd 0x8(%rsp), %xmm0
xorl %edi, %edi
callq 0x2b3e
movapd 0x80(%rsp), %xmm2
unpcklpd 0x90(%rsp), %xmm2       # xmm2 = xmm2[0],mem[0]
mulpd %xmm2, %xmm2
movapd 0x50(%rsp), %xmm1
unpcklpd 0x60(%rsp), %xmm1       # xmm1 = xmm1[0],mem[0]
mulpd 0x30(%rsp), %xmm1
movapd %xmm1, %xmm4
movapd %xmm2, %xmm1
movapd %xmm2, %xmm3
shufpd $0x1, %xmm2, %xmm1      # xmm1 = xmm1[1],xmm2[0]
mulpd %xmm4, %xmm1
movapd 0x70(%rsp), %xmm4
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
movapd 0x40(%rsp), %xmm2
unpcklpd %xmm2, %xmm2            # xmm2 = xmm2[0,0]
mulpd %xmm3, %xmm4
addpd %xmm2, %xmm4
mulpd %xmm1, %xmm4
unpcklpd 0x20(%rsp), %xmm0       # xmm0 = xmm0[0],mem[0]
divpd %xmm0, %xmm4
movapd %xmm4, %xmm1
unpckhpd %xmm4, %xmm1            # xmm1 = xmm1[1],xmm4[1]
subsd %xmm4, %xmm1
mulsd 0x18(%rsp), %xmm1
xorps %xmm0, %xmm0
cvtsi2sd %ebp, %xmm0
mulsd %xmm1, %xmm0
addsd 0x10(%rsp), %xmm0
leal (%rbx,%r14), %eax
movapd 0x79fa(%rip), %xmm1     # 0xe020
testb $0x1, %al
jne 0x662e
xorpd %xmm1, %xmm0
cmpl %ebx, %r14d
jle 0x6637
xorpd %xmm1, %xmm0
addq $0xa8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
leaq 0x7f62(%rip), %rdi      # 0xe5ae
leaq 0x7ecd(%rip), %rsi      # 0xe520
leaq 0x7ff7(%rip), %rcx      # 0xe651
movl $0x65, %edx
callq 0x20f0