LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Exact.cpp
// Function: Exact::Talphabetai(int, int, int, double, double)
std::complex<double> Talphabetai(int alpha, int beta, int i, double a, double delta)
{
	assert (alpha >= 0 && alpha <= 2); // alpha valid for 0, 1, 2
	assert (beta >= 0 && beta <= 2); // beta valid for 0, 1, 2
	assert (i >= 0 && i <= 2); // i valid for 0, 1, 2

	Matrix<std::complex<double> > U = UMNSm(a, delta);
	return std::conj(U(alpha, i)) * U(beta, i);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x48, %rsp
cmpl $0x3, %edi
jae 0x53f5
movl %esi, %ebx
cmpl $0x3, %esi
jae 0x5414
movl %edx, %ebp
cmpl $0x3, %edx
jae 0x5433
movl %edi, %r14d
leaq 0x30(%rsp), %r15
movq %r15, %rdi
callq 0x5055
movl %r14d, %esi
movl %ebp, %r14d
movq %r15, %rdi
movq %r14, %rdx
callq 0xa590
movsd (%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x8(%rax), %xmm0
movaps %xmm0, 0x10(%rsp)
movl %ebx, %esi
leaq 0x30(%rsp), %rdi
movq %r14, %rdx
callq 0xa590
movaps 0x10(%rsp), %xmm1
xorps 0x8c7f(%rip), %xmm1     # 0xe020
leaq 0x20(%rsp), %rdi
movsd 0x8(%rsp), %xmm0
movsd %xmm0, (%rdi)
movlps %xmm1, 0x8(%rdi)
movq %rax, %rsi
callq 0x3994
movsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
leaq 0x30(%rsp), %rdi
callq 0xa582
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
addq $0x48, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
leaq 0x8d83(%rip), %rdi      # 0xe17f
leaq 0x8f8d(%rip), %rsi      # 0xe390
leaq 0x8f94(%rip), %rcx      # 0xe39e
movl $0xc9, %edx
callq 0x20f0
leaq 0x8dc9(%rip), %rdi      # 0xe1e4
leaq 0x8f6e(%rip), %rsi      # 0xe390
leaq 0x8f75(%rip), %rcx      # 0xe39e
movl $0xca, %edx
callq 0x20f0
leaq 0x8dc1(%rip), %rdi      # 0xe1fb
leaq 0x8f4f(%rip), %rsi      # 0xe390
leaq 0x8f56(%rip), %rcx      # 0xe39e
movl $0xcb, %edx
callq 0x20f0
movq %rax, %rbx
leaq 0x30(%rsp), %rdi
callq 0xa582
movq %rbx, %rdi
callq 0x22a0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/First_Peak.cpp
// Function: PL(double, double const*)
double PL(double L, const double * p)
{
	int alpha = p[0];
	int beta = p[1];
	double E = p[2];
	double Yrho = p[3];
	double delta = p[4];
	int sign = p[5];
	return sign * Exact::Palphabeta(alpha, beta, Yrho * E * Y_to_a, L / E, delta);
}

// Assembly:
subq $0x18, %rsp
movapd %xmm0, %xmm1
cvttsd2si (%rdi), %eax
cvttsd2si 0x8(%rdi), %esi
movsd 0x28(%rdi), %xmm0
cvttpd2dq %xmm0, %xmm0
cvtdq2pd %xmm0, %xmm0
movaps %xmm0, (%rsp)
movsd 0x10(%rdi), %xmm0
divsd %xmm0, %xmm1
mulsd 0x18(%rdi), %xmm0
movsd 0x20(%rdi), %xmm2
leaq 0x8b0c(%rip), %rcx      # 0xe330
mulsd (%rcx), %xmm0
movl %eax, %edi
callq 0x5467
mulsd (%rsp), %xmm0
addq $0x18, %rsp
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/First_Peak.cpp
// Function: First_Extremum_Exact(int, int, double, double, double, bool)
double First_Extremum_Exact(int alpha, int beta, double E, double Yrho, double delta, bool max)
{
	double L_range_min, L_range_max;
	double p[6];
	p[0] = alpha;
	p[1] = beta;
	p[2] = E;
	p[3] = Yrho;
	p[4] = delta;
	p[5] = (alpha == beta) ? 1 : -1;
	p[5] *= (max) ? 1 : -1; // switch sign for max/min

	// numbers here are selected so that the first max/min are contained in this region, but not the second
	if (max)
	{
		L_range_min = E * 300;
		L_range_max = E * 1000;
	}
	else
	{
		L_range_min = E * 500;
		L_range_max = E * 1800;
	}

	return std::abs(gss_min(PL, L_range_min, L_range_max, tol, p));
}

// Assembly:
subq $0x38, %rsp
cvtsi2sd %edi, %xmm3
cvtsi2sd %esi, %xmm4
movsd %xmm3, (%rsp)
movsd %xmm4, 0x8(%rsp)
movsd %xmm0, 0x10(%rsp)
movsd %xmm1, 0x18(%rsp)
movsd %xmm2, 0x20(%rsp)
xorl %eax, %eax
cmpl %esi, %edi
sete %al
leal -0x1(,%rax,2), %eax
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
testl %edx, %edx
jne 0x5883
xorpd 0x879d(%rip), %xmm1     # 0xe020
movsd %xmm1, 0x28(%rsp)
jne 0x5895
movsd 0x8c45(%rip), %xmm3     # 0xe4d8
jmp 0x589d
movsd 0x8c33(%rip), %xmm3     # 0xe4d0
mulsd %xmm0, %xmm3
jne 0x58ad
movsd 0x8c3d(%rip), %xmm1     # 0xe4e8
jmp 0x58b5
movsd 0x8c2b(%rip), %xmm1     # 0xe4e0
mulsd %xmm0, %xmm1
movsd 0xd8df(%rip), %xmm2     # 0x131a0
leaq -0xe0(%rip), %rdi       # 0x57e8
movq %rsp, %rsi
movapd %xmm3, %xmm0
callq 0xab8c
andpd 0x8be4(%rip), %xmm0     # 0xe4c0
addq $0x38, %rsp
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/First_Peak.cpp
// Function: FEE_x(double, double const*)
double FEE_x(double E, const double *p)
{
	int alpha = p[0];
	int beta = p[1];
	double Yrho = p[2];
	double delta = p[3];
	bool max = p[4] > 0;

	return First_Extremum_Exact(alpha, beta, E, Yrho, delta, max) / E;
}

// Assembly:
pushq %rax
movsd %xmm0, (%rsp)
cvttsd2si (%rdi), %eax
cvttsd2si 0x8(%rdi), %esi
movsd 0x10(%rdi), %xmm1
movsd 0x18(%rdi), %xmm2
movsd 0x20(%rdi), %xmm4
xorpd %xmm3, %xmm3
xorl %edx, %edx
ucomisd %xmm3, %xmm4
seta %dl
movl %eax, %edi
callq 0x5839
divsd (%rsp), %xmm0
popq %rax
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/First_Peak.cpp
// Function: FP_line(int, int, double, double, bool, double*, double*)
void FP_line(int alpha, int beta, double Yrho, double delta, bool max, double *m, double *b)
{
	double p[5];
	p[0] = alpha;
	p[1] = beta;
	p[2] = Yrho;
	p[3] = delta;
	p[4] = 1; // look for a max

	double rsq;
	Linear_Fit(FEE_x, FEE_y, 0.1, 5., 1e2, p, m, b, &rsq);
}

// Assembly:
subq $0x38, %rsp
movq %r8, %r9
movq %rcx, %r8
cvtsi2sd %edi, %xmm2
leaq 0x10(%rsp), %rcx
movsd %xmm2, (%rcx)
xorps %xmm2, %xmm2
cvtsi2sd %esi, %xmm2
movsd %xmm2, 0x8(%rcx)
movsd %xmm0, 0x10(%rcx)
movsd %xmm1, 0x18(%rcx)
movabsq $0x3ff0000000000000, %rax # imm = 0x3FF0000000000000
movq %rax, 0x20(%rcx)
leaq 0x8(%rsp), %rax
movq %rax, (%rsp)
leaq -0xb0(%rip), %rdi       # 0x58e1
leaq -0x7e(%rip), %rsi       # 0x591a
movsd 0x8b50(%rip), %xmm0     # 0xe4f0
movsd 0x8b50(%rip), %xmm1     # 0xe4f8
movl $0x64, %edx
callq 0xa27c
addq $0x38, %rsp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/GF.cpp
// Function: GF::ReUm4(double, double, int, int, int, int, int, int)
double ReUm4(double a, double delta, int alpha, int beta, int w, int x, int y, int z)
{
	assert (alpha >= 0 && alpha <= 2); // alpha valid for 0, 1, 2
	assert (beta >= 0 && beta <= 2); // beta valid for 0, 1, 2
	assert (w >= 0 && w <= 2); // w valid for 0, 1, 2
	assert (x >= 0 && x <= 2); // x valid for 0, 1, 2
	assert (y >= 0 && y <= 2); // y valid for 0, 1, 2
	assert (z >= 0 && z <= 2); // z valid for 0, 1, 2

	Matrix<std::complex<double> > U = Check::UMNSm(a, delta);

	return std::real(U(alpha, w) * std::conj(U(beta, x) * U(alpha, y)) * U(beta, z));
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x68, %rsp
cmpl $0x3, %edi
jae 0x5af6
movl %esi, %r14d
cmpl $0x3, %esi
jae 0x5b15
movl %edx, %r15d
cmpl $0x3, %edx
jae 0x5b34
movl %ecx, %r12d
cmpl $0x3, %ecx
jae 0x5b53
movl %r8d, %ebp
cmpl $0x3, %r8d
jae 0x5b72
movl %r9d, 0x8(%rsp)
cmpl $0x3, %r9d
jae 0x5b91
movl %edi, %r13d
leaq 0x28(%rsp), %rbx
movq %rbx, %rdi
callq 0x32a0
movl %r13d, %r13d
movl %r15d, %edx
movq %rbx, %rdi
movq %r13, %rsi
callq 0xa590
movq %rax, %r15
movl %r14d, %r14d
movl %r12d, %edx
leaq 0x28(%rsp), %rdi
movq %r14, %rsi
callq 0xa590
movq %rax, %r12
movl %ebp, %edx
leaq 0x28(%rsp), %rdi
movq %r13, %rsi
callq 0xa590
movups (%r12), %xmm0
leaq 0x10(%rsp), %rdi
movaps %xmm0, (%rdi)
movq %rax, %rsi
callq 0x3994
leaq 0x10(%rsp), %rdi
movsd (%rdi), %xmm0
movsd 0x8(%rdi), %xmm1
xorps 0x859e(%rip), %xmm1     # 0xe020
leaq 0x58(%rsp), %rsi
movsd %xmm0, (%rsi)
movlps %xmm1, 0x8(%rsi)
movups (%r15), %xmm0
movaps %xmm0, (%rdi)
callq 0x3994
movaps 0x10(%rsp), %xmm0
movaps %xmm0, 0x40(%rsp)
movl 0x8(%rsp), %edx
leaq 0x28(%rsp), %rdi
movq %r14, %rsi
callq 0xa590
leaq 0x10(%rsp), %rdi
movaps 0x40(%rsp), %xmm0
movaps %xmm0, (%rdi)
movq %rax, %rsi
callq 0x3994
movsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
leaq 0x28(%rsp), %rdi
callq 0xa582
movsd 0x8(%rsp), %xmm0
addq $0x68, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x8682(%rip), %rdi      # 0xe17f
leaq 0x8a1c(%rip), %rsi      # 0xe520
leaq 0x8a20(%rip), %rcx      # 0xe52b
movl $0x13, %edx
callq 0x20f0
leaq 0x86c8(%rip), %rdi      # 0xe1e4
leaq 0x89fd(%rip), %rsi      # 0xe520
leaq 0x8a01(%rip), %rcx      # 0xe52b
movl $0x14, %edx
callq 0x20f0
leaq 0x8a2f(%rip), %rdi      # 0xe56a
leaq 0x89de(%rip), %rsi      # 0xe520
leaq 0x89e2(%rip), %rcx      # 0xe52b
movl $0x15, %edx
callq 0x20f0
leaq 0x8a21(%rip), %rdi      # 0xe57b
leaq 0x89bf(%rip), %rsi      # 0xe520
leaq 0x89c3(%rip), %rcx      # 0xe52b
movl $0x16, %edx
callq 0x20f0
leaq 0x8a13(%rip), %rdi      # 0xe58c
leaq 0x89a0(%rip), %rsi      # 0xe520
leaq 0x89a4(%rip), %rcx      # 0xe52b
movl $0x17, %edx
callq 0x20f0
leaq 0x8a05(%rip), %rdi      # 0xe59d
leaq 0x8981(%rip), %rsi      # 0xe520
leaq 0x8985(%rip), %rcx      # 0xe52b
movl $0x18, %edx
callq 0x20f0
jmp 0x5bb2
movq %rax, %rbx
leaq 0x28(%rsp), %rdi
callq 0xa582
movq %rbx, %rdi
callq 0x22a0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/GF.cpp
// Function: GF::G1(int, int, double, double)
double G1(int alpha, int Beta, double a, double delta)
{
	double spsi = sin(Check::psi(a));
	Matrix<std::complex<double> > U = Check::UMNSm(a, delta);
	double deltaalphabeta = alpha == Beta ? 1 : 0;
	return -spsi * std::real((U(alpha, 0) * std::conj(U(Beta, 2)) + U(alpha, 2) * std::conj(U(Beta, 0)))
				* (2. * std::conj(U(alpha, 2)) * U(Beta, 2) - deltaalphabeta));
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0xa8, %rsp
movsd %xmm1, 0x20(%rsp)
movsd %xmm0, 0x8(%rsp)
movl %esi, %ebx
movl %edi, %ebp
callq 0x29f0
callq 0x2070
movaps %xmm0, 0x50(%rsp)
leaq 0x38(%rsp), %r15
movq %r15, %rdi
movsd 0x8(%rsp), %xmm0
movsd 0x20(%rsp), %xmm1
callq 0x32a0
cmpl %ebx, %ebp
je 0x5d34
xorpd %xmm0, %xmm0
jmp 0x5d3c
movsd 0x82dc(%rip), %xmm0     # 0xe018
movsd %xmm0, 0x8(%rsp)
movslq %ebp, %r14
movq %r15, %rdi
movq %r14, %rsi
xorl %edx, %edx
callq 0xa590
movq %rax, %r15
movslq %ebx, %rbx
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %rbx, %rsi
callq 0xa590
movsd (%rax), %xmm0
movsd 0x8(%rax), %xmm1
xorps 0x82a6(%rip), %xmm1     # 0xe020
leaq 0x98(%rsp), %rsi
movsd %xmm0, (%rsi)
movlps %xmm1, 0x8(%rsi)
movupd (%r15), %xmm0
leaq 0x10(%rsp), %rdi
movapd %xmm0, (%rdi)
callq 0x3994
movapd 0x10(%rsp), %xmm0
movapd %xmm0, 0x20(%rsp)
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %r14, %rsi
callq 0xa590
movq %rax, %r15
leaq 0x38(%rsp), %rdi
movq %rbx, %rsi
xorl %edx, %edx
callq 0xa590
movsd (%rax), %xmm0
movsd 0x8(%rax), %xmm1
xorps 0x8243(%rip), %xmm1     # 0xe020
leaq 0x88(%rsp), %rsi
movsd %xmm0, (%rsi)
movlps %xmm1, 0x8(%rsi)
movupd (%r15), %xmm0
leaq 0x10(%rsp), %rdi
movapd %xmm0, (%rdi)
callq 0x3994
movapd 0x20(%rsp), %xmm0
addpd 0x10(%rsp), %xmm0
movapd %xmm0, 0x20(%rsp)
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %r14, %rsi
callq 0xa590
movupd (%rax), %xmm0
movapd %xmm0, 0x60(%rsp)
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %rbx, %rsi
callq 0xa590
movapd 0x60(%rsp), %xmm0
mulpd 0x86b2(%rip), %xmm0     # 0xe500
leaq 0x10(%rsp), %rdi
movapd %xmm0, (%rdi)
movq %rax, %rsi
callq 0x3994
leaq 0x10(%rsp), %rdi
movsd (%rdi), %xmm0
subsd 0x8(%rsp), %xmm0
movsd 0x8(%rdi), %xmm1
leaq 0x78(%rsp), %rsi
movsd %xmm0, (%rsi)
movsd %xmm1, 0x8(%rsi)
movapd 0x20(%rsp), %xmm0
movapd %xmm0, (%rdi)
callq 0x3994
movapd 0x50(%rsp), %xmm0
xorpd 0x8182(%rip), %xmm0     # 0xe020
mulsd 0x10(%rsp), %xmm0
movapd %xmm0, 0x50(%rsp)
leaq 0x38(%rsp), %rdi
callq 0xa582
movaps 0x50(%rsp), %xmm0
addq $0xa8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
jmp 0x5ecf
jmp 0x5ecf
jmp 0x5ecf
jmp 0x5ecf
movq %rax, %rbx
leaq 0x38(%rsp), %rdi
callq 0xa582
movq %rbx, %rdi
callq 0x22a0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/GF.cpp
// Function: GF::G2(int, int, double, double)
double G2(int alpha, int Beta, double a, double delta)
{
	double cpsi = cos(Check::psi(a));
	Matrix<std::complex<double> > U = Check::UMNSm(a, delta);
	double deltaalphabeta = alpha == Beta ? 1 : 0;
	return cpsi * std::real((U(alpha, 1) * std::conj(U(Beta, 2)) + U(alpha, 2) * std::conj(U(Beta, 1)))
				* (2. * std::conj(U(alpha, 2)) * U(Beta, 2) - deltaalphabeta));
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x98, %rsp
movsd %xmm1, 0x20(%rsp)
movsd %xmm0, (%rsp)
movl %esi, %ebx
movl %edi, %ebp
callq 0x29f0
callq 0x22c0
movsd %xmm0, 0x8(%rsp)
leaq 0x38(%rsp), %r15
movq %r15, %rdi
movsd (%rsp), %xmm0
movsd 0x20(%rsp), %xmm1
callq 0x32a0
cmpl %ebx, %ebp
je 0x5f32
xorpd %xmm0, %xmm0
jmp 0x5f3a
movsd 0x80de(%rip), %xmm0     # 0xe018
movsd %xmm0, (%rsp)
movslq %ebp, %r14
movl $0x1, %edx
movq %r15, %rdi
movq %r14, %rsi
callq 0xa590
movq %rax, %r15
movslq %ebx, %rbx
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %rbx, %rsi
callq 0xa590
movsd (%rax), %xmm0
movsd 0x8(%rax), %xmm1
xorps 0x80a6(%rip), %xmm1     # 0xe020
leaq 0x88(%rsp), %rsi
movsd %xmm0, (%rsi)
movlps %xmm1, 0x8(%rsi)
movupd (%r15), %xmm0
leaq 0x10(%rsp), %rdi
movapd %xmm0, (%rdi)
callq 0x3994
movapd 0x10(%rsp), %xmm0
movapd %xmm0, 0x20(%rsp)
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %r14, %rsi
callq 0xa590
movq %rax, %r15
leaq 0x38(%rsp), %rdi
movl $0x1, %edx
movq %rbx, %rsi
callq 0xa590
movsd (%rax), %xmm0
movsd 0x8(%rax), %xmm1
xorps 0x8040(%rip), %xmm1     # 0xe020
leaq 0x78(%rsp), %rsi
movsd %xmm0, (%rsi)
movlps %xmm1, 0x8(%rsi)
movupd (%r15), %xmm0
leaq 0x10(%rsp), %rdi
movapd %xmm0, (%rdi)
callq 0x3994
movapd 0x20(%rsp), %xmm0
addpd 0x10(%rsp), %xmm0
movapd %xmm0, 0x20(%rsp)
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %r14, %rsi
callq 0xa590
movupd (%rax), %xmm0
movapd %xmm0, 0x50(%rsp)
leaq 0x38(%rsp), %rdi
movl $0x2, %edx
movq %rbx, %rsi
callq 0xa590
movapd 0x50(%rsp), %xmm0
mulpd 0x84b2(%rip), %xmm0     # 0xe500
leaq 0x10(%rsp), %rdi
movapd %xmm0, (%rdi)
movq %rax, %rsi
callq 0x3994
leaq 0x10(%rsp), %rdi
movsd (%rdi), %xmm0
subsd (%rsp), %xmm0
movsd 0x8(%rdi), %xmm1
leaq 0x68(%rsp), %rsi
movsd %xmm0, (%rsi)
movsd %xmm1, 0x8(%rsi)
movapd 0x20(%rsp), %xmm0
movapd %xmm0, (%rdi)
callq 0x3994
movsd 0x8(%rsp), %xmm0
mulsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
leaq 0x38(%rsp), %rdi
callq 0xa582
movsd 0x8(%rsp), %xmm0
addq $0x98, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
jmp 0x60c7
jmp 0x60c7
jmp 0x60c7
jmp 0x60c7
movq %rax, %rbx
leaq 0x38(%rsp), %rdi
callq 0xa582
movq %rbx, %rdi
callq 0x22a0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/GF.cpp
// Function: GF::C21(int, int, double, double, int)
double C21(int alpha, int Beta, double a, double delta, int order)
{
	assert (order >= 0 && order <= 1); // order valid for 0, 1
	double C210 = -ReUm4(a, delta, alpha, Beta, 1, 1, 0, 0);
	double C1 = epsp(a) * Dmsqee;
	double F1 = GF::F1(alpha, Beta, a, delta);
	double F2 = GF::F2(alpha, Beta, a, delta);
	double Dl31 = Check::Dl31(a, order);
	double Dl32 = Check::Dl32(a, order);
	double C211 = C1 * (F1 / Dl31 + F2 / Dl32);
	return C210 + order * C211;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x50, %rsp
cmpl $0x2, %edx
jae 0x61d2
movl %edx, %ebx
movl %esi, %ebp
movl %edi, %r14d
movsd %xmm0, 0x8(%rsp)
movsd %xmm1, 0x10(%rsp)
movl $0x1, %edx
movl $0x1, %ecx
xorl %r8d, %r8d
xorl %r9d, %r9d
callq 0x59b8
movsd %xmm0, 0x28(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x3a04
leaq 0xd05b(%rip), %rax      # 0x13188
mulsd (%rax), %xmm0
movsd %xmm0, 0x20(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5bc7
movaps %xmm0, 0x40(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5c5c
movaps %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
movl %ebx, %edi
callq 0x2aca
movaps %xmm0, 0x30(%rsp)
movsd 0x8(%rsp), %xmm0
movl %ebx, %edi
callq 0x2b69
movapd 0x40(%rsp), %xmm2
unpcklpd 0x10(%rsp), %xmm2       # xmm2 = xmm2[0],mem[0]
movapd 0x30(%rsp), %xmm1
unpcklpd %xmm0, %xmm1            # xmm1 = xmm1[0],xmm0[0]
divpd %xmm1, %xmm2
movapd %xmm2, %xmm1
unpckhpd %xmm2, %xmm1            # xmm1 = xmm1[1],xmm2[1]
addsd %xmm2, %xmm1
xorps %xmm0, %xmm0
cvtsi2sd %ebx, %xmm0
mulsd 0x20(%rsp), %xmm1
mulsd %xmm1, %xmm0
subsd 0x28(%rsp), %xmm0
addq $0x50, %rsp
popq %rbx
popq %r14
popq %rbp
retq
leaq 0x83d5(%rip), %rdi      # 0xe5ae
leaq 0x8340(%rip), %rsi      # 0xe520
leaq 0x83e0(%rip), %rcx      # 0xe5c7
movl $0x3d, %edx
callq 0x20f0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/GF.cpp
// Function: GF::C31(int, int, double, double, int)
double C31(int alpha, int Beta, double a, double delta, int order)
{
	assert (order >= 0 && order <= 1); // order valid for 0, 1
	double C310 = -ReUm4(a, delta, alpha, Beta, 2, 2, 0, 0);
	double C1 = epsp(a) * Dmsqee;
	double F1 = GF::F1(alpha, Beta, a, delta);
	double F2 = GF::F2(alpha, Beta, a, delta);
	double G1 = GF::G1(alpha, Beta, a, delta);
	double Dl31 = Check::Dl31(a, order);
	double Dl32 = Check::Dl32(a, order);
	double C311 = C1 * ((F1 + G1) / Dl31 - F2 / Dl32);
	return C310 + order * C311;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x60, %rsp
cmpl $0x2, %edx
jae 0x6309
movl %edx, %ebx
movl %esi, %ebp
movl %edi, %r14d
movsd %xmm0, 0x8(%rsp)
movsd %xmm1, 0x10(%rsp)
movl $0x2, %edx
movl $0x2, %ecx
xorl %r8d, %r8d
xorl %r9d, %r9d
callq 0x59b8
movsd %xmm0, 0x38(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x3a04
leaq 0xcf46(%rip), %rax      # 0x13188
mulsd (%rax), %xmm0
movsd %xmm0, 0x30(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5bc7
movsd %xmm0, 0x28(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5c5c
movaps %xmm0, 0x50(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5ce5
movaps %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
movl %ebx, %edi
callq 0x2aca
movaps %xmm0, 0x40(%rsp)
movsd 0x8(%rsp), %xmm0
movl %ebx, %edi
callq 0x2b69
movapd 0x10(%rsp), %xmm1
addsd 0x28(%rsp), %xmm1
unpcklpd 0x50(%rsp), %xmm1       # xmm1 = xmm1[0],mem[0]
movapd 0x40(%rsp), %xmm2
unpcklpd %xmm0, %xmm2            # xmm2 = xmm2[0],xmm0[0]
divpd %xmm2, %xmm1
movapd %xmm1, %xmm0
unpckhpd %xmm1, %xmm0            # xmm0 = xmm0[1],xmm1[1]
subsd %xmm0, %xmm1
xorps %xmm0, %xmm0
cvtsi2sd %ebx, %xmm0
mulsd 0x30(%rsp), %xmm1
mulsd %xmm1, %xmm0
subsd 0x38(%rsp), %xmm0
addq $0x60, %rsp
popq %rbx
popq %r14
popq %rbp
retq
leaq 0x829e(%rip), %rdi      # 0xe5ae
leaq 0x8209(%rip), %rsi      # 0xe520
leaq 0x82d7(%rip), %rcx      # 0xe5f5
movl $0x4a, %edx
callq 0x20f0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/GF.cpp
// Function: GF::C32(int, int, double, double, int)
double C32(int alpha, int Beta, double a, double delta, int order)
{
	assert (order >= 0 && order <= 1); // order valid for 0, 1
	double C320 = -ReUm4(a, delta, alpha, Beta, 2, 2, 1, 1);
	double C1 = epsp(a) * Dmsqee;
	double F1 = GF::F1(alpha, Beta, a, delta);
	double F2 = GF::F2(alpha, Beta, a, delta);
	double G2 = GF::G2(alpha, Beta, a, delta);
	double Dl31 = Check::Dl31(a, order);
	double Dl32 = Check::Dl32(a, order);
	double C321 = C1 * (-F1 / Dl31 + (F2 + G2) / Dl32);
	return C320 + order * C321;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x60, %rsp
cmpl $0x2, %edx
jae 0x6442
movl %edx, %ebx
movl %esi, %ebp
movl %edi, %r14d
movsd %xmm0, 0x8(%rsp)
movsd %xmm1, 0x10(%rsp)
movl $0x2, %edx
movl $0x2, %ecx
movl $0x1, %r8d
movl $0x1, %r9d
callq 0x59b8
movsd %xmm0, 0x38(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x3a04
leaq 0xce09(%rip), %rax      # 0x13188
mulsd (%rax), %xmm0
movsd %xmm0, 0x30(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5bc7
movaps %xmm0, 0x50(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5c5c
movsd %xmm0, 0x28(%rsp)
movl %r14d, %edi
movl %ebp, %esi
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
callq 0x5ee4
movaps %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
movl %ebx, %edi
callq 0x2aca
movaps %xmm0, 0x40(%rsp)
movsd 0x8(%rsp), %xmm0
movl %ebx, %edi
callq 0x2b69
movapd 0x10(%rsp), %xmm1
addsd 0x28(%rsp), %xmm1
unpcklpd 0x50(%rsp), %xmm1       # xmm1 = xmm1[0],mem[0]
unpcklpd 0x40(%rsp), %xmm0       # xmm0 = xmm0[0],mem[0]
divpd %xmm0, %xmm1
movapd %xmm1, %xmm0
unpckhpd %xmm1, %xmm0            # xmm0 = xmm0[1],xmm1[1]
subsd %xmm0, %xmm1
xorps %xmm0, %xmm0
cvtsi2sd %ebx, %xmm0
mulsd 0x30(%rsp), %xmm1
mulsd %xmm1, %xmm0
subsd 0x38(%rsp), %xmm0
addq $0x60, %rsp
popq %rbx
popq %r14
popq %rbp
retq
leaq 0x8165(%rip), %rdi      # 0xe5ae
leaq 0x80d0(%rip), %rsi      # 0xe520
leaq 0x81cc(%rip), %rcx      # 0xe623
movl $0x58, %edx
callq 0x20f0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Hat.cpp
// Function: Hat::lp(double)
double lp(double a)
{
	double x = Dmsqee + a;
	double y = normal_sign * sqrt(pow(Dmsqee - a, 2) + 4 * s13sq * a * Dmsqee);
	double z = eps * Dmsqee * s12sq;
	return (x + y) / 2. + z;
}

// Assembly:
pushq %rbx
subq $0x10, %rsp
movapd %xmm0, %xmm1
leaq 0xc638(%rip), %rbx      # 0x13188
movsd (%rbx), %xmm2
movapd %xmm2, %xmm3
leaq 0x7749(%rip), %rax      # 0xe2a8
movsd (%rax), %xmm0
mulsd 0x79a5(%rip), %xmm0     # 0xe510
mulsd %xmm1, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm2
leaq 0xc61e(%rip), %rax      # 0x1319c
cvtsi2sdl (%rax), %xmm4
subsd %xmm1, %xmm3
mulsd %xmm3, %xmm3
addsd %xmm3, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x6b9e
sqrtsd %xmm0, %xmm0
jmp 0x6bb9
movsd %xmm2, 0x8(%rsp)
movsd %xmm4, (%rsp)
callq 0x2120
movsd (%rsp), %xmm4
movsd 0x8(%rsp), %xmm2
mulsd %xmm4, %xmm0
leaq 0xc5cc(%rip), %rax      # 0x13190
movsd (%rax), %xmm1
mulsd (%rbx), %xmm1
addsd %xmm0, %xmm2
leaq 0x76c1(%rip), %rax      # 0xe298
mulsd (%rax), %xmm1
mulsd 0x7425(%rip), %xmm2     # 0xe008
addsd %xmm1, %xmm2
movapd %xmm2, %xmm0
addq $0x10, %rsp
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Hat.cpp
// Function: Hat::Pee(double, double)
double Pee(double a, double LE)
{
	double L4E = km_per_GeV_to_per_eV2 * LE / 4; // in eV^-2 now
	double Deltapm = (lp(a) - lm(a)) * L4E;
	double P = 1 - pow(s2phi(a) * sin(Deltapm), 2);
	ProbabilityWarning W(&P);
	return P;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0xc8, %rsp
movapd %xmm0, %xmm3
leaq 0x710c(%rip), %rax      # 0xe328
mulsd (%rax), %xmm1
leaq 0xbf61(%rip), %rbx      # 0x13188
movsd (%rbx), %xmm4
leaq 0xbf6a(%rip), %r15      # 0x1319c
movl (%r15), %ebp
movapd %xmm4, %xmm2
subsd %xmm0, %xmm2
leaq 0x7064(%rip), %rax      # 0xe2a8
movsd (%rax), %xmm8
mulsd 0x72be(%rip), %xmm8     # 0xe510
mulsd %xmm2, %xmm2
mulsd %xmm0, %xmm8
movapd %xmm8, %xmm0
movapd %xmm4, 0x60(%rsp)
mulsd %xmm4, %xmm0
addsd %xmm2, %xmm0
xorpd %xmm2, %xmm2
ucomisd %xmm2, %xmm0
movapd %xmm3, 0x10(%rsp)
jb 0x7284
sqrtsd %xmm0, %xmm5
jmp 0x72b1
movsd %xmm1, 0x8(%rsp)
movsd %xmm8, 0x20(%rsp)
callq 0x2120
xorpd %xmm2, %xmm2
movsd 0x20(%rsp), %xmm8
movapd 0x10(%rsp), %xmm3
movsd 0x8(%rsp), %xmm1
movapd %xmm0, %xmm5
mulsd 0x707f(%rip), %xmm1     # 0xe338
movsd %xmm1, 0x8(%rsp)
leaq 0xbeca(%rip), %r14      # 0x13190
movsd (%r14), %xmm7
movsd (%rbx), %xmm6
leaq 0x6fc2(%rip), %rax      # 0xe298
movsd (%rax), %xmm4
movl (%r15), %r15d
movapd %xmm6, %xmm0
subsd %xmm3, %xmm0
mulsd %xmm0, %xmm0
mulsd %xmm6, %xmm8
addsd %xmm0, %xmm8
ucomisd %xmm2, %xmm8
jb 0x7304
xorps %xmm0, %xmm0
sqrtsd %xmm8, %xmm0
jmp 0x7344
movapd %xmm8, %xmm0
movapd %xmm4, 0x20(%rsp)
movapd %xmm5, 0x50(%rsp)
movapd %xmm6, 0x40(%rsp)
movapd %xmm7, 0x30(%rsp)
callq 0x2120
movapd 0x30(%rsp), %xmm7
movapd 0x40(%rsp), %xmm6
movapd 0x50(%rsp), %xmm5
movapd 0x20(%rsp), %xmm4
movapd 0x10(%rsp), %xmm3
unpcklpd %xmm0, %xmm5            # xmm5 = xmm5[0],xmm0[0]
movapd 0x60(%rsp), %xmm2
unpcklpd %xmm6, %xmm2            # xmm2 = xmm2[0],xmm6[0]
movd %r15d, %xmm0
movd %ebp, %xmm1
punpckldq %xmm0, %xmm1    # xmm1 = xmm1[0],xmm0[0],xmm1[1],xmm0[1]
movapd %xmm3, %xmm0
unpcklpd %xmm3, %xmm0            # xmm0 = xmm0[0],xmm3[0]
cvtdq2pd %xmm1, %xmm1
addpd %xmm2, %xmm0
mulpd %xmm5, %xmm1
movhpd (%r14), %xmm7           # xmm7 = xmm7[0],mem[0]
movhpd (%rbx), %xmm6           # xmm6 = xmm6[0],mem[0]
mulpd %xmm7, %xmm6
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
mulpd %xmm6, %xmm4
movapd %xmm0, %xmm2
addpd %xmm1, %xmm2
subpd %xmm1, %xmm0
movsd %xmm2, %xmm0            # xmm0 = xmm2[0],xmm0[1]
mulpd 0x7320(%rip), %xmm0     # 0xe6c0
addpd %xmm4, %xmm0
movapd %xmm0, %xmm1
unpckhpd %xmm0, %xmm1            # xmm1 = xmm1[1],xmm0[1]
subsd %xmm1, %xmm0
movsd 0x8(%rsp), %xmm1
mulsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
movapd %xmm3, %xmm0
callq 0x6ca2
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x2070
mulsd 0x10(%rsp), %xmm0
mulsd %xmm0, %xmm0
movsd 0x6c2c(%rip), %xmm1     # 0xe018
subsd %xmm0, %xmm1
leaq 0x78(%rsp), %rbx
movsd %xmm1, (%rbx)
leaq 0x80(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rsi
callq 0xd2da
movsd (%rbx), %xmm0
movsd %xmm0, 0x8(%rsp)
movq %r14, %rdi
callq 0xce78
movsd 0x8(%rsp), %xmm0
addq $0xc8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Hat.cpp
// Function: Hat::W(double, int)
Matrix<std::complex<double> > W(double a, int order)
{
	assert (order >= 0 && order <= 1); // order valid for 0, 1

	Matrix<std::complex<double> > W(3, 3);
	double c = cos(phi(a) - t13);
	double s = sin(phi(a) - t13);
	double coef = eps * c12 * s12 * Dmsqee;

	switch (order)
	{
		case 1:
			W(0, 1) += coef * c / Dl0m(a);
			W(1, 0) -= coef * c / Dl0m(a);
			W(1, 2) += coef * s / Dlp0(a);
			W(2, 1) -= coef * s / Dlp0(a);

		case 0:
			// zeroth order
			for (int i = 0; i < 3; i++)
				W(i, i) += 1;
	}
	return W;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x30, %rsp
movsd %xmm0, 0x8(%rsp)
cmpl $0x2, %esi
jae 0x7b28
movl %esi, %ebp
movq %rdi, %rbx
movl $0x3, %esi
movl $0x3, %edx
callq 0xa436
movsd 0x8(%rsp), %xmm0
callq 0x6e76
leaq 0x698d(%rip), %rax      # 0xe318
movsd (%rax), %xmm1
movsd %xmm1, 0x10(%rsp)
subsd %xmm1, %xmm0
callq 0x22c0
movsd %xmm0, 0x18(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x6e76
subsd 0x10(%rsp), %xmm0
callq 0x2070
testb $0x1, %bpl
je 0x7af2
movsd %xmm0, 0x10(%rsp)
leaq 0xb7bf(%rip), %rax      # 0x13190
movsd (%rax), %xmm0
leaq 0x6914(%rip), %rax      # 0xe2f0
mulsd (%rax), %xmm0
leaq 0x68f1(%rip), %rax      # 0xe2d8
mulsd (%rax), %xmm0
leaq 0xb796(%rip), %rax      # 0x13188
mulsd (%rax), %xmm0
movsd %xmm0, 0x28(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x6ed3
movsd %xmm0, 0x20(%rsp)
movl $0x1, %edx
movq %rbx, %rdi
xorl %esi, %esi
callq 0xa590
movsd 0x18(%rsp), %xmm0
mulsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
divsd 0x20(%rsp), %xmm0
addsd (%rax), %xmm0
movsd %xmm0, (%rax)
movsd 0x8(%rsp), %xmm0
callq 0x6ed3
movsd %xmm0, 0x20(%rsp)
movl $0x1, %esi
movq %rbx, %rdi
xorl %edx, %edx
callq 0xa590
movsd 0x18(%rsp), %xmm1
divsd 0x20(%rsp), %xmm1
movsd (%rax), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movsd 0x8(%rsp), %xmm0
callq 0x7138
movsd %xmm0, 0x18(%rsp)
movl $0x1, %esi
movl $0x2, %edx
movq %rbx, %rdi
callq 0xa590
movsd 0x10(%rsp), %xmm0
mulsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
divsd 0x18(%rsp), %xmm0
addsd (%rax), %xmm0
movsd %xmm0, (%rax)
movsd 0x8(%rsp), %xmm0
callq 0x7138
movsd %xmm0, 0x8(%rsp)
movl $0x2, %esi
movl $0x1, %edx
movq %rbx, %rdi
callq 0xa590
movsd 0x10(%rsp), %xmm1
divsd 0x8(%rsp), %xmm1
movsd (%rax), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax)
xorl %r14d, %r14d
movq %rbx, %rdi
movq %r14, %rsi
movq %r14, %rdx
callq 0xa590
movsd (%rax), %xmm0
addsd 0x6509(%rip), %xmm0     # 0xe018
movsd %xmm0, (%rax)
incq %r14
cmpq $0x3, %r14
jne 0x7af5
movq %rbx, %rax
addq $0x30, %rsp
popq %rbx
popq %r14
popq %rbp
retq
leaq 0x6a7f(%rip), %rdi      # 0xe5ae
leaq 0x6b9a(%rip), %rsi      # 0xe6d0
leaq 0x6b9f(%rip), %rcx      # 0xe6dc
movl $0x5e, %edx
callq 0x20f0
jmp 0x7b49
movq %rax, %r14
movq %rbx, %rdi
callq 0xa582
movq %r14, %rdi
callq 0x22a0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Hat.cpp
// Function: Hat::Talphabetai(int, int, int, double, double, int)
std::complex<double> Talphabetai(int alpha, int beta, int i, double a, 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 (i >= 0 && i <= 2); // i valid for 0, 1, 2
	assert (order >= 0 && order <= 1); // order valid for 0, 1

	Matrix<std::complex<double> > Vtmp = V(a, delta, order);
	return std::conj(Vtmp(alpha, i)) * Vtmp(beta, i);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x48, %rsp
cmpl $0x3, %edi
jae 0x7d98
movl %esi, %ebx
cmpl $0x3, %esi
jae 0x7db7
movl %edx, %ebp
cmpl $0x3, %edx
jae 0x7dd6
cmpl $0x2, %ecx
jae 0x7df5
movl %edi, %r14d
leaq 0x30(%rsp), %r15
movq %r15, %rdi
movl %ecx, %esi
callq 0x7bed
movl %r14d, %esi
movl %ebp, %r14d
movq %r15, %rdi
movq %r14, %rdx
callq 0xa590
movsd (%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x8(%rax), %xmm0
movaps %xmm0, 0x10(%rsp)
movl %ebx, %esi
leaq 0x30(%rsp), %rdi
movq %r14, %rdx
callq 0xa590
movaps 0x10(%rsp), %xmm1
xorps 0x62dc(%rip), %xmm1     # 0xe020
leaq 0x20(%rsp), %rdi
movsd 0x8(%rsp), %xmm0
movsd %xmm0, (%rdi)
movlps %xmm1, 0x8(%rdi)
movq %rax, %rsi
callq 0x3994
movsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
leaq 0x30(%rsp), %rdi
callq 0xa582
movsd 0x8(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
addq $0x48, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
leaq 0x63e0(%rip), %rdi      # 0xe17f
leaq 0x692a(%rip), %rsi      # 0xe6d0
leaq 0x6999(%rip), %rcx      # 0xe746
movl $0x81, %edx
callq 0x20f0
leaq 0x6426(%rip), %rdi      # 0xe1e4
leaq 0x690b(%rip), %rsi      # 0xe6d0
leaq 0x697a(%rip), %rcx      # 0xe746
movl $0x82, %edx
callq 0x20f0
leaq 0x641e(%rip), %rdi      # 0xe1fb
leaq 0x68ec(%rip), %rsi      # 0xe6d0
leaq 0x695b(%rip), %rcx      # 0xe746
movl $0x83, %edx
callq 0x20f0
leaq 0x67b2(%rip), %rdi      # 0xe5ae
leaq 0x68cd(%rip), %rsi      # 0xe6d0
leaq 0x693c(%rip), %rcx      # 0xe746
movl $0x84, %edx
callq 0x20f0
movq %rax, %rbx
leaq 0x30(%rsp), %rdi
callq 0xa582
movq %rbx, %rdi
callq 0x22a0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Hyperbolas.cpp
// Function: Hyperbolas::l2(double)
double l2(double a)
{
	double ax = cos(2 * t12) * eps * Dmsqee;
	double ax12 = lb - la(0);
	double ax13 = lc - la(0);
	double atr = (ax13 - ax12) / 2.; // ansatz that simple average is right, perhaps weighted average by the X's is better

	if (a < atr)
		return 0.5 * ((a - ax) + sqrt(pow(a - ax, 2) + pow(X12(a), 2))) + lb;
	else
		return 0.5 * ((la(a) + lc) - sqrt(pow(la(a) - lc, 2) + pow(X23, 2)));

}

// Assembly:
pushq %r14
pushq %rbx
subq $0x18, %rsp
movsd %xmm0, (%rsp)
leaq 0x5f32(%rip), %rax      # 0xe308
movsd (%rax), %xmm0
addsd %xmm0, %xmm0
movsd %xmm0, 0x10(%rsp)
callq 0x22c0
leaq 0xada0(%rip), %rbx      # 0x13190
movsd (%rbx), %xmm2
leaq 0xad8d(%rip), %r14      # 0x13188
movsd (%r14), %xmm3
movsd 0xada8(%rip), %xmm5     # 0x131b0
leaq 0x5e99(%rip), %rax      # 0xe2a8
leaq 0x5e82(%rip), %rcx      # 0xe298
movsd (%rcx), %xmm1
mulsd %xmm2, %xmm1
addsd (%rax), %xmm1
mulsd %xmm3, %xmm1
xorpd %xmm8, %xmm8
movapd %xmm1, %xmm6
addsd %xmm8, %xmm6
subsd %xmm6, %xmm5
movsd 0xad78(%rip), %xmm4     # 0x131b8
movapd %xmm4, %xmm7
subsd %xmm6, %xmm7
subsd %xmm5, %xmm7
movsd (%rsp), %xmm5
mulsd 0x5baf(%rip), %xmm7     # 0xe008
ucomisd %xmm5, %xmm7
jbe 0x84a8
mulsd %xmm2, %xmm0
mulsd %xmm3, %xmm0
subsd %xmm0, %xmm5
movsd %xmm5, (%rsp)
mulsd %xmm5, %xmm5
movsd %xmm5, 0x8(%rsp)
movsd 0x10(%rsp), %xmm0
callq 0x2070
mulsd (%rbx), %xmm0
mulsd (%r14), %xmm0
mulsd %xmm0, %xmm0
addsd 0x8(%rsp), %xmm0
ucomisd 0x5eb8(%rip), %xmm0     # 0xe358
jb 0x84d9
sqrtsd %xmm0, %xmm0
jmp 0x84de
addsd %xmm5, %xmm1
movapd %xmm4, %xmm2
addsd %xmm1, %xmm2
subsd %xmm4, %xmm1
mulsd %xmm1, %xmm1
movsd 0xacfc(%rip), %xmm0     # 0x131c0
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
ucomisd %xmm8, %xmm0
jb 0x84fd
sqrtsd %xmm0, %xmm0
jmp 0x850c
callq 0x2120
movsd (%rsp), %xmm1
addsd %xmm0, %xmm1
mulsd 0x5b19(%rip), %xmm1     # 0xe008
addsd 0xacb9(%rip), %xmm1     # 0x131b0
movapd %xmm1, %xmm0
jmp 0x851c
movsd %xmm2, (%rsp)
callq 0x2120
movsd (%rsp), %xmm2
subsd %xmm0, %xmm2
mulsd 0x5af0(%rip), %xmm2     # 0xe008
movapd %xmm2, %xmm0
addq $0x18, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Matrix.cpp
// Function: Matrix<std::complex<double>>::operator=(Matrix<std::complex<double>>)
Matrix<T> Matrix<T>::operator = (const Matrix other)
{
	rows = other.rows;
	cols = other.cols;
	buffer = new T[rows * cols];
	std::copy(other.buffer, other.buffer + rows * cols, buffer);
	return *this;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rdx, %r15
movq %rsi, %r14
movq %rdi, %rbx
movq (%rdx), %rax
movq %rax, (%rsi)
movq 0x8(%rdx), %rbp
movq %rbp, 0x8(%rsi)
imulq %rax, %rbp
movq %rbp, %rax
shrq $0x3c, %rax
movq %rbp, %r12
shlq $0x4, %r12
xorl %edi, %edi
negq %rax
sbbq %rdi, %rdi
orq %r12, %rdi
callq 0x2030
movq %rax, %r13
testq %rbp, %rbp
je 0xa6b9
movq %r13, %rdi
xorl %esi, %esi
movq %r12, %rdx
callq 0x20d0
movq %r13, 0x10(%r14)
testq %rbp, %rbp
je 0xa6d1
movq 0x10(%r15), %rsi
movq %r13, %rdi
movq %r12, %rdx
callq 0x2270
movq %rbx, %rdi
movq %r14, %rsi
callq 0xa502
movq %rbx, %rax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Matrix.cpp
// Function: Matrix<std::complex<double>>::operator+=(Matrix<std::complex<double>>)
Matrix<T> Matrix<T>::operator += (const Matrix other)
{
	assert(rows == other.rows);
	assert(cols == other.cols);

	for (std::size_t i = 0; i < rows * cols; i++)
		buffer[i] += other.buffer[i];
	return *this;
}

// Assembly:
pushq %rbx
movq (%rsi), %rax
cmpq (%rdx), %rax
jne 0xa753
movq 0x8(%rsi), %rax
cmpq 0x8(%rdx), %rax
jne 0xa772
movq %rdi, %rbx
movq 0x8(%rsi), %rax
imulq (%rsi), %rax
testq %rax, %rax
je 0xa746
xorl %eax, %eax
xorl %ecx, %ecx
movq 0x10(%rdx), %rdi
movq 0x10(%rsi), %r8
movupd (%rdi,%rax), %xmm0
movupd (%r8,%rax), %xmm1
addpd %xmm0, %xmm1
movupd %xmm1, (%r8,%rax)
incq %rcx
movq 0x8(%rsi), %rdi
imulq (%rsi), %rdi
addq $0x10, %rax
cmpq %rdi, %rcx
jb 0xa715
movq %rbx, %rdi
callq 0xa502
movq %rbx, %rax
popq %rbx
retq
leaq 0x422c(%rip), %rdi      # 0xe986
leaq 0x4137(%rip), %rsi      # 0xe898
leaq 0x4231(%rip), %rcx      # 0xe999
movl $0x40, %edx
callq 0x20f0
leaq 0x427f(%rip), %rdi      # 0xe9f8
leaq 0x4118(%rip), %rsi      # 0xe898
leaq 0x4212(%rip), %rcx      # 0xe999
movl $0x41, %edx
callq 0x20f0
nop
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Matrix.cpp
// Function: Matrix<std::complex<double>>::operator+(Matrix<std::complex<double>>) const
Matrix<T> Matrix<T>::operator + (const Matrix other) const
{
	assert(rows == other.rows);
	assert(cols == other.cols);

	Matrix ret(cols, rows);
	for (std::size_t i = 0; i < rows * cols; i++)
		ret.buffer[i] = buffer[i] + other.buffer[i];
	return ret;
}

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
movq %rdx, %rbx
movq (%rsi), %rdx
cmpq (%rbx), %rdx
jne 0xa804
movq %rsi, %r15
movq 0x8(%rsi), %rsi
cmpq 0x8(%rbx), %rsi
jne 0xa823
movq %rdi, %r14
callq 0xa436
movq 0x8(%r15), %rax
imulq (%r15), %rax
testq %rax, %rax
je 0xa7fb
xorl %eax, %eax
xorl %ecx, %ecx
movq 0x10(%r15), %rdx
movq 0x10(%rbx), %rsi
movq 0x10(%r14), %rdi
movupd (%rdx,%rax), %xmm0
movupd (%rsi,%rax), %xmm1
addpd %xmm0, %xmm1
movupd %xmm1, (%rdi,%rax)
incq %rcx
movq 0x8(%r15), %rdx
imulq (%r15), %rdx
addq $0x10, %rax
cmpq %rdx, %rcx
jb 0xa7c8
movq %r14, %rax
popq %rbx
popq %r14
popq %r15
retq
leaq 0x417b(%rip), %rdi      # 0xe986
leaq 0x4086(%rip), %rsi      # 0xe898
leaq 0x41f2(%rip), %rcx      # 0xea0b
movl $0x4b, %edx
callq 0x20f0
leaq 0x41ce(%rip), %rdi      # 0xe9f8
leaq 0x4067(%rip), %rsi      # 0xe898
leaq 0x41d3(%rip), %rcx      # 0xea0b
movl $0x4c, %edx
callq 0x20f0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Matrix.cpp
// Function: Matrix<std::complex<double>>::operator*(Matrix<std::complex<double>>) const
Matrix<T> Matrix<T>::operator * (const Matrix rhs) const
{
	assert(cols == rhs.rows);

	Matrix ret(rows, rhs.cols);
	for (std::size_t i = 0; i < rows; i++)
	{
		for (std::size_t j = 0; j < rhs.cols; j++)
		{
			for (std::size_t k = 0; k < rows; k++)
			{
				ret(i, j) += this->Get(i, k) * rhs(k, j);
			}
		}
	}

	return ret;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x38, %rsp
movq 0x8(%rsi), %rax
cmpq (%rdx), %rax
jne 0xa94d
movq %rdx, %rbx
movq %rsi, %r15
movq %rdi, %r14
movq (%rsi), %rsi
movq 0x8(%rdx), %rdx
callq 0xa436
cmpq $0x0, (%r15)
je 0xa93b
xorl %r12d, %r12d
cmpq $0x0, 0x8(%rbx)
je 0xa92f
xorl %r13d, %r13d
cmpq $0x0, (%r15)
je 0xa922
xorl %ebp, %ebp
movq %r15, %rdi
movq %r12, %rsi
movq %rbp, %rdx
callq 0xa96c
movsd %xmm0, (%rsp)
movsd %xmm1, 0x20(%rsp)
movq %rbx, %rdi
movq %rbp, %rsi
movq %r13, %rdx
callq 0xa5f6
movsd %xmm0, 0x28(%rsp)
movsd %xmm1, 0x30(%rsp)
movsd (%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
leaq 0x10(%rsp), %rdi
leaq 0x28(%rsp), %rsi
callq 0x3994
movaps 0x10(%rsp), %xmm0
movaps %xmm0, (%rsp)
movq %r14, %rdi
movq %r12, %rsi
movq %r13, %rdx
callq 0xa590
movupd (%rax), %xmm0
addpd (%rsp), %xmm0
movupd %xmm0, (%rax)
incq %rbp
cmpq (%r15), %rbp
jb 0xa899
incq %r13
cmpq 0x8(%rbx), %r13
jb 0xa88d
incq %r12
cmpq (%r15), %r12
jb 0xa87f
movq %r14, %rax
addq $0x38, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x411b(%rip), %rdi      # 0xea6f
leaq 0x3f3d(%rip), %rsi      # 0xe898
leaq 0x411e(%rip), %rcx      # 0xea80
movl $0x57, %edx
callq 0x20f0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Matrix.cpp
// Function: Matrix<std::complex<double>>::Transpose()
void Matrix<T>::Transpose()
{
	T * tmp = new T[rows * cols];
	std::copy(buffer, buffer + rows * cols, tmp);

	for (std::size_t i = 0; i < rows; i++)
	{
		for (std::size_t j = 0; j < cols; j++)
			buffer[Get_index(i, j)] = tmp[Get_index(j, i)];
	}
	delete [] tmp;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
movq (%rdi), %r12
imulq 0x8(%rdi), %r12
movq %r12, %rax
shrq $0x3c, %rax
movq %r12, %r15
shlq $0x4, %r15
xorl %edi, %edi
negq %rax
sbbq %rdi, %rdi
orq %r15, %rdi
callq 0x2030
movq %rax, %r14
testq %r12, %r12
je 0xaa1c
movq %r14, %rdi
xorl %esi, %esi
movq %r15, %rdx
callq 0x20d0
testq %r12, %r12
je 0xaa30
movq 0x10(%rbx), %rsi
movq %r14, %rdi
movq %r15, %rdx
callq 0x2270
cmpq $0x0, (%rbx)
je 0xaa8a
xorl %eax, %eax
xorl %ecx, %ecx
movq 0x8(%rbx), %rdx
testq %rdx, %rdx
je 0xaa7e
movq %rcx, %rsi
shlq $0x4, %rsi
addq %r14, %rsi
xorl %edi, %edi
xorl %r8d, %r8d
movq %rax, %r9
imulq %rdx, %r9
imulq %r8, %rdx
shlq $0x4, %rdx
addq 0x10(%rbx), %r9
movups (%rsi,%rdx), %xmm0
movups %xmm0, (%rdi,%r9)
incq %r8
movq 0x8(%rbx), %rdx
addq $0x10, %rdi
cmpq %rdx, %r8
jb 0xaa52
incq %rcx
addq $0x10, %rax
cmpq (%rbx), %rcx
jb 0xaa3a
movq %r14, %rdi
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
jmp 0x21e0
nop
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Matrix.cpp
// Function: Matrix<std::complex<double>>::Hermitian()
void Matrix<T>::Hermitian()
{
	T * tmp = new T[rows * cols];
	std::copy(buffer, buffer + rows * cols, tmp);

	for (std::size_t i = 0; i < rows; i++)
	{
		for (std::size_t j = 0; j < cols; j++)
			buffer[Get_index(i, j)] = std::conj(tmp[Get_index(j, i)]);
	}
	delete [] tmp;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
movq (%rdi), %r12
imulq 0x8(%rdi), %r12
movq %r12, %rax
shrq $0x3c, %rax
movq %r12, %r15
shlq $0x4, %r15
xorl %edi, %edi
negq %rax
sbbq %rdi, %rdi
orq %r15, %rdi
callq 0x2030
movq %rax, %r14
testq %r12, %r12
je 0xaae4
movq %r14, %rdi
xorl %esi, %esi
movq %r15, %rdx
callq 0x20d0
testq %r12, %r12
je 0xaaf8
movq 0x10(%rbx), %rsi
movq %r14, %rdi
movq %r15, %rdx
callq 0x2270
cmpq $0x0, (%rbx)
je 0xab6d
xorl %eax, %eax
movaps 0x3519(%rip), %xmm0     # 0xe020
xorl %ecx, %ecx
movq 0x8(%rbx), %rdi
testq %rdi, %rdi
je 0xab61
movq %rcx, %rdx
shlq $0x4, %rdx
addq %r14, %rdx
movl $0x8, %esi
xorl %r8d, %r8d
movq %rax, %r9
imulq %rdi, %r9
imulq %r8, %rdi
shlq $0x4, %rdi
movsd (%rdx,%rdi), %xmm1
movsd 0x8(%rdx,%rdi), %xmm2
xorps %xmm0, %xmm2
addq 0x10(%rbx), %r9
movsd %xmm1, -0x8(%rsi,%r9)
movlps %xmm2, (%rsi,%r9)
incq %r8
movq 0x8(%rbx), %rdi
addq $0x10, %rsi
cmpq %rdi, %r8
jb 0xab24
incq %rcx
addq $0x10, %rax
cmpq (%rbx), %rcx
jb 0xab09
movq %r14, %rdi
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
jmp 0x21e0
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Minimization.cpp
// Function: double gss_min<double>(double (*)(double, double const*), double, double, double, double*)
T gss_min(T (*f)(T, const T*), T a_initial, T b_initial, T tol, T *p)
{
	T gr = (sqrt(5.) - 1) / 2;
	T a = a_initial;
	T b = b_initial;
	T c = b - gr * (b - a);
	T d = a + gr * (b - a);
	while (std::abs(c - d) > tol)
	{
		if (f(c, p) < f(d, p))
		{
			b = d;
			d = c;
			c = b - gr * (b - a);
		}
		else
		{
			a = c;
			c = d;
			d = a + gr * (b - a);
		}
	}
	T ret = (b + a) / 2.;

	// Warns user in case final minimum is close to an edge
	if (std::abs(ret - a_initial) < 2 * tol or std::abs(b_initial - ret) < 2 * tol)
		EdgeCaseWarning W(a_initial, b_initial, ret);
	return ret;
}

// Assembly:
pushq %r14
pushq %rbx
subq $0xd8, %rsp
movapd %xmm1, %xmm4
subsd %xmm0, %xmm4
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
movapd %xmm1, %xmm3
mulpd 0x3fa2(%rip), %xmm4     # 0xeb50
unpcklpd %xmm0, %xmm3            # xmm3 = xmm3[0],xmm0[0]
addpd %xmm3, %xmm4
movapd %xmm4, %xmm5
unpckhpd %xmm4, %xmm5            # xmm5 = xmm5[1],xmm4[1]
movapd %xmm4, %xmm3
subsd %xmm5, %xmm3
andpd 0x38f2(%rip), %xmm3     # 0xe4c0
ucomisd %xmm2, %xmm3
movapd %xmm1, 0x20(%rsp)
movapd %xmm1, %xmm3
movapd %xmm0, 0x50(%rsp)
movsd %xmm0, 0x8(%rsp)
jbe 0xacbc
movq %rsi, %rbx
movq %rdi, %r14
movapd 0x50(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movapd 0x20(%rsp), %xmm3
movsd %xmm2, 0x18(%rsp)
movapd %xmm5, 0x30(%rsp)
movapd %xmm4, 0x60(%rsp)
movsd %xmm3, 0x10(%rsp)
movapd %xmm4, %xmm0
movq %rbx, %rdi
callq *%r14
movsd %xmm0, 0x48(%rsp)
movapd 0x30(%rsp), %xmm0
movq %rbx, %rdi
callq *%r14
ucomisd 0x48(%rsp), %xmm0
jbe 0xac6c
movapd 0x30(%rsp), %xmm3
movapd %xmm3, %xmm4
subsd 0x8(%rsp), %xmm4
mulsd 0x3f0e(%rip), %xmm4     # 0xeb68
addsd %xmm3, %xmm4
movapd 0x60(%rsp), %xmm0
movsd 0x18(%rsp), %xmm2
jmp 0xac9e
movsd 0x10(%rsp), %xmm3
movapd %xmm3, %xmm0
movapd 0x60(%rsp), %xmm1
subsd %xmm1, %xmm0
mulsd 0x3ed8(%rip), %xmm0     # 0xeb60
addsd %xmm1, %xmm0
movsd %xmm1, 0x8(%rsp)
movsd 0x18(%rsp), %xmm2
movapd 0x30(%rsp), %xmm4
movapd %xmm4, %xmm1
subsd %xmm0, %xmm1
andpd 0x3812(%rip), %xmm1     # 0xe4c0
ucomisd %xmm2, %xmm1
movapd %xmm0, %xmm5
ja 0xac0c
addsd 0x8(%rsp), %xmm3
mulsd 0x333e(%rip), %xmm3     # 0xe008
movapd %xmm3, %xmm1
movapd 0x50(%rsp), %xmm0
subsd %xmm0, %xmm1
andpd 0x37e0(%rip), %xmm1     # 0xe4c0
addsd %xmm2, %xmm2
ucomisd %xmm1, %xmm2
ja 0xad02
movapd 0x20(%rsp), %xmm1
subsd %xmm3, %xmm1
andpd 0x37c4(%rip), %xmm1     # 0xe4c0
ucomisd %xmm1, %xmm2
jbe 0xad2d
leaq 0x78(%rsp), %rbx
movq %rbx, %rdi
movapd 0x20(%rsp), %xmm1
movapd %xmm3, %xmm2
movsd %xmm3, 0x10(%rsp)
callq 0xd150
movq %rbx, %rdi
callq 0xce78
movsd 0x10(%rsp), %xmm3
movapd %xmm3, %xmm0
addq $0xd8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Mixed.cpp
// Function: Mixed::s12msq(double, int)
double s12msq(double D, int order)
{
	double x = (pow(Check::l2(D, order), 2) - alpha() * Check::l2(D, order) + beta()) * Check::Dl31(D, order);
	double s12msq_n = -x;
	double s12msq_d = Check::Dl32(D, order) * (pow(Check::l1(D, order), 2) - alpha() * Check::l1(D, order) + beta()) - x;
	return s12msq_n / s12msq_d;
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x48, %rsp
movl %edi, %ebx
movsd %xmm0, (%rsp)
callq 0x26f6
mulsd %xmm0, %xmm0
movsd %xmm0, 0x20(%rsp)
leaq 0x82e2(%rip), %r14      # 0x13180
leaq 0x342b(%rip), %rax      # 0xe2d0
movsd (%rax), %xmm2
movsd %xmm2, 0x18(%rsp)
leaq 0x33fa(%rip), %rax      # 0xe2b0
movsd (%rax), %xmm0
movsd %xmm0, 0x30(%rsp)
leaq 0x33f9(%rip), %rax      # 0xe2c0
movsd (%rax), %xmm1
movsd %xmm1, 0x38(%rsp)
leaq 0x33d0(%rip), %rax      # 0xe2a8
movapd %xmm2, %xmm3
mulsd %xmm1, %xmm3
addsd (%rax), %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x10(%rsp)
movsd (%r14), %xmm0
mulsd %xmm2, %xmm0
addsd %xmm3, %xmm0
movsd %xmm0, 0x8(%rsp)
movsd (%rsp), %xmm0
movl %ebx, %edi
callq 0x26f6
mulsd 0x8(%rsp), %xmm0
movsd 0x20(%rsp), %xmm1
subsd %xmm0, %xmm1
movsd (%r14), %xmm0
mulsd 0x18(%rsp), %xmm0
mulsd 0x30(%rsp), %xmm0
mulsd 0x38(%rsp), %xmm0
addsd %xmm1, %xmm0
movapd %xmm0, 0x20(%rsp)
movsd (%rsp), %xmm0
movl %ebx, %edi
callq 0x2aca
movapd 0x20(%rsp), %xmm1
xorpd 0x30c8(%rip), %xmm1     # 0xe020
mulsd %xmm1, %xmm0
movsd %xmm0, 0x20(%rsp)
movsd (%rsp), %xmm0
movl %ebx, %edi
callq 0x2b69
movsd %xmm0, 0x8(%rsp)
movsd (%rsp), %xmm0
movl %ebx, %edi
callq 0x252c
mulsd %xmm0, %xmm0
movsd %xmm0, 0x40(%rsp)
movsd (%r14), %xmm0
mulsd 0x18(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd (%rsp), %xmm0
movl %ebx, %edi
callq 0x252c
mulsd 0x10(%rsp), %xmm0
movsd 0x18(%rsp), %xmm1
mulsd (%r14), %xmm1
movsd 0x40(%rsp), %xmm2
subsd %xmm0, %xmm2
mulsd 0x30(%rsp), %xmm1
mulsd 0x38(%rsp), %xmm1
addsd %xmm2, %xmm1
mulsd 0x8(%rsp), %xmm1
movsd 0x20(%rsp), %xmm0
addsd %xmm0, %xmm1
divsd %xmm1, %xmm0
addq $0x48, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O1/PeterDenton[P]Nu-Pert/src/Mixed.cpp
// Function: Mixed::sdeltam(double, double, int)
double sdeltam(double D, double delta, int order)
{
	// tmp
	double Et = E(D, order);
	double Ft = F(D, order);
	// denominator
	double densq1 = pow(Et, 2) * s23sq + pow(Ft, 2) * c23sq + 2 * Et * Ft * c23 * s23 * cos(delta);
	double densq2 = pow(Et, 2) * c23sq + pow(Ft, 2) * s23sq - 2 * Et * Ft * c23 * s23 * cos(delta);
	// numerator c23sq -> c23 typo corrected
	return ((pow(Et, 2) + pow(Ft, 2)) * sin(delta) * s23 * c23) / sqrt(densq1 * densq2);
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x68, %rsp
movl %edi, %ebx
movsd %xmm1, 0x8(%rsp)
movsd %xmm0, 0x20(%rsp)
callq 0xada9
movaps %xmm0, 0x10(%rsp)
movsd 0x20(%rsp), %xmm0
movl %ebx, %edi
callq 0x28c0
leaq 0x7c66(%rip), %rax      # 0x13180
subsd (%rax), %xmm0
leaq 0x2d8b(%rip), %rax      # 0xe2b0
mulsd (%rax), %xmm0
leaq 0x2dc0(%rip), %rax      # 0xe2f0
mulsd (%rax), %xmm0
leaq 0x2d9d(%rip), %rax      # 0xe2d8
mulsd (%rax), %xmm0
leaq 0x2dba(%rip), %rax      # 0xe300
mulsd (%rax), %xmm0
leaq 0x2d4f(%rip), %rbx      # 0xe2a0
movapd 0x10(%rsp), %xmm1
movapd %xmm1, %xmm2
unpcklpd %xmm0, %xmm2            # xmm2 = xmm2[0],xmm0[0]
mulpd %xmm2, %xmm2
movapd %xmm2, 0x50(%rsp)
leaq 0x2d58(%rip), %r14      # 0xe2c8
addsd %xmm1, %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, 0x10(%rsp)
leaq 0x2d73(%rip), %rax      # 0xe2f8
movsd (%rax), %xmm0
movaps %xmm0, 0x20(%rsp)
leaq 0x2d4b(%rip), %rax      # 0xe2e0
movsd (%rax), %xmm0
movaps %xmm0, 0x40(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x22c0
movaps %xmm0, 0x30(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x22c0
movapd 0x2a5f(%rip), %xmm6     # 0xe020
movapd 0x10(%rsp), %xmm7
xorpd %xmm7, %xmm6
movsd (%r14), %xmm1
movsd (%rbx), %xmm2
movapd %xmm2, %xmm3
unpcklpd %xmm1, %xmm3            # xmm3 = xmm3[0],xmm1[0]
movapd 0x50(%rsp), %xmm5
movapd %xmm5, %xmm4
unpckhpd %xmm5, %xmm4            # xmm4 = xmm4[1],xmm5[1]
mulpd %xmm4, %xmm3
movapd %xmm4, %xmm8
movapd %xmm5, %xmm4
unpcklpd %xmm5, %xmm4            # xmm4 = xmm4[0],xmm5[0]
unpcklpd %xmm2, %xmm1            # xmm1 = xmm1[0],xmm2[0]
mulpd %xmm4, %xmm1
addpd %xmm3, %xmm1
unpcklpd %xmm7, %xmm6            # xmm6 = xmm6[0],xmm7[0]
movapd 0x20(%rsp), %xmm2
unpcklpd %xmm2, %xmm2            # xmm2 = xmm2[0,0]
mulpd %xmm6, %xmm2
movapd 0x40(%rsp), %xmm3
unpcklpd %xmm3, %xmm3            # xmm3 = xmm3[0,0]
mulpd %xmm2, %xmm3
unpcklpd 0x30(%rsp), %xmm0       # xmm0 = xmm0[0],mem[0]
mulpd %xmm3, %xmm0
addpd %xmm1, %xmm0
movapd %xmm0, 0x10(%rsp)
addsd %xmm5, %xmm8
movapd %xmm8, 0x30(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x2070
movapd %xmm0, %xmm2
mulsd 0x30(%rsp), %xmm2
mulsd 0x40(%rsp), %xmm2
mulsd 0x20(%rsp), %xmm2
movapd 0x10(%rsp), %xmm1
movapd %xmm1, %xmm0
unpckhpd %xmm1, %xmm0            # xmm0 = xmm0[1],xmm1[1]
mulsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0xb68a
sqrtsd %xmm0, %xmm0
jmp 0xb69b
movsd %xmm2, 0x8(%rsp)
callq 0x2120
movsd 0x8(%rsp), %xmm2
divsd %xmm0, %xmm2
movapd %xmm2, %xmm0
addq $0x68, %rsp
popq %rbx
popq %r14
retq