LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/ImageProcs.cpp
// Function: cielab2xyz(double, double, double, double&, double&, double&)
void cielab2xyz(Float_t l, Float_t a, Float_t b,
				Float_t &x, Float_t &y, Float_t &z)
{
	y = (l + 16.0) / 116.0;
	x = a / 500.0 + y;
	z = y - b / 200.0;
	
	if (pow(y, 3) > .008856)
		y = pow(y, 3);
	else y = (y - 16.0 / 116.0) / 7.787;
	if (pow(x, 3) > .008856)
		x = pow(x, 3);
	else x = (x - 16.0 / 116.0) / 7.787;
	if (pow(z, 3) > .008856)
		z = pow(z, 3);
	else z = (z - 16.0 / 116.0) / 7.787;
	
	x *= 95.047;
	y *= 100.0;
	z *= 108.883;
}

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
movq %rdx, %rbx
movq %rsi, %r14
movq %rdi, %r15
addsd 0x1a61b(%rip), %xmm0    # 0x2cbe0
divsd 0x1a5f3(%rip), %xmm0    # 0x2cbc0
movsd %xmm0, (%rsi)
divsd 0x1a5f7(%rip), %xmm1    # 0x2cbd0
addsd %xmm0, %xmm1
movsd %xmm1, (%rdi)
divsd 0x1a5ff(%rip), %xmm2    # 0x2cbe8
addsd (%rsi), %xmm2
movsd %xmm2, (%rdx)
movsd (%rsi), %xmm0
movsd 0x1a5f3(%rip), %xmm1    # 0x2cbf0
callq 0x34f0
ucomisd 0x1a58e(%rip), %xmm0    # 0x2cb98
movsd (%r14), %xmm0
jbe 0x12620
movsd 0x1a5d7(%rip), %xmm1    # 0x2cbf0
callq 0x34f0
jmp 0x12630
addsd 0x1a5d0(%rip), %xmm0    # 0x2cbf8
divsd 0x1a570(%rip), %xmm0    # 0x2cba0
movsd %xmm0, (%r14)
movsd (%r15), %xmm0
movsd 0x1a5ae(%rip), %xmm1    # 0x2cbf0
callq 0x34f0
ucomisd 0x1a549(%rip), %xmm0    # 0x2cb98
movsd (%r15), %xmm0
jbe 0x12665
movsd 0x1a592(%rip), %xmm1    # 0x2cbf0
callq 0x34f0
jmp 0x12675
addsd 0x1a58b(%rip), %xmm0    # 0x2cbf8
divsd 0x1a52b(%rip), %xmm0    # 0x2cba0
movsd %xmm0, (%r15)
movsd (%rbx), %xmm0
movsd 0x1a56a(%rip), %xmm1    # 0x2cbf0
callq 0x34f0
ucomisd 0x1a505(%rip), %xmm0    # 0x2cb98
movsd (%rbx), %xmm0
jbe 0x126a8
movsd 0x1a54f(%rip), %xmm1    # 0x2cbf0
callq 0x34f0
jmp 0x126b8
addsd 0x1a548(%rip), %xmm0    # 0x2cbf8
divsd 0x1a4e8(%rip), %xmm0    # 0x2cba0
movsd %xmm0, (%rbx)
movsd (%r15), %xmm0
mulsd 0x1a4b7(%rip), %xmm0    # 0x2cb80
movsd %xmm0, (%r15)
movsd (%r14), %xmm0
mulsd 0x1a4ad(%rip), %xmm0    # 0x2cb88
movsd %xmm0, (%r14)
movsd (%rbx), %xmm0
mulsd 0x1a4a4(%rip), %xmm0    # 0x2cb90
movsd %xmm0, (%rbx)
popq %rbx
popq %r14
popq %r15
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/ImageProcs.cpp
// Function: calcNonoutlierMeanPixelValue(ImageReal&, double)
Float_t calcNonoutlierMeanPixelValue(ImageReal &grayImage, Float_t outlierSigmas)
{
	int totalPixels = grayImage.totalPixels();
	
	//Find the mean pixel value
	Float_t meanVal = 0;
	for (int i = 0; i < totalPixels; i++)
		meanVal += grayImage[i];
	meanVal /= totalPixels;
	
	//Find the variance of pixel value
	Float_t variance = 0;
	for (int i = 0; i < totalPixels; i++)
		variance += pow(grayImage[i] - meanVal, 2);
	variance /= (totalPixels * 3);
	
	//Find the standard deviation of pixel value
	Float_t stddev = sqrt(variance);
	
	//Define outliers as some number of standard deviations
	Float_t minOutlier = meanVal - stddev * outlierSigmas;
	Float_t maxOutlier = meanVal + stddev * outlierSigmas;
	
	//Find the mean pixel value of nonoutliers
	Float_t nonOutlierMeanVal = 0;
	int numNonOutliers = 0;
		for (int i = 0; i < totalPixels; i++)
			if (grayImage[i] >= minOutlier && grayImage[i] <= maxOutlier)
			{
				nonOutlierMeanVal += grayImage[i];
				numNonOutliers++;
			}
	nonOutlierMeanVal /= numNonOutliers;
	
	return nonOutlierMeanVal;
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x18, %rsp
movapd %xmm0, %xmm2
movq %rdi, %rbx
movslq 0x30(%rdi), %r14
xorpd %xmm0, %xmm0
xorpd %xmm3, %xmm3
testq %r14, %r14
jle 0x139eb
movq 0x38(%rbx), %rax
xorl %ecx, %ecx
addsd (%rax,%rcx,8), %xmm3
incq %rcx
cmpq %rcx, %r14
jne 0x139de
cvtsi2sd %r14d, %xmm1
divsd %xmm1, %xmm3
testl %r14d, %r14d
jle 0x13a1c
movq 0x38(%rbx), %rax
xorpd %xmm0, %xmm0
xorl %ecx, %ecx
movsd (%rax,%rcx,8), %xmm1
subsd %xmm3, %xmm1
mulsd %xmm1, %xmm1
addsd %xmm1, %xmm0
incq %rcx
cmpq %rcx, %r14
jne 0x13a03
leal (%r14,%r14,2), %eax
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
divsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x13a3b
sqrtsd %xmm0, %xmm0
jmp 0x13a58
movsd %xmm2, 0x10(%rsp)
movsd %xmm3, 0x8(%rsp)
callq 0x3240
movsd 0x8(%rsp), %xmm3
movsd 0x10(%rsp), %xmm2
testl %r14d, %r14d
jle 0x13aba
mulsd %xmm2, %xmm0
movapd %xmm3, %xmm2
subsd %xmm0, %xmm2
addsd %xmm3, %xmm0
movq 0x38(%rbx), %rax
xorpd %xmm1, %xmm1
xorl %ecx, %ecx
xorl %edx, %edx
movsd (%rax,%rcx,8), %xmm3
movapd %xmm2, %xmm4
cmplepd %xmm3, %xmm4
movapd %xmm3, %xmm5
cmplepd %xmm0, %xmm5
andpd %xmm4, %xmm5
movq %xmm5, %rsi
testb $0x1, %sil
je 0x13aa3
addsd %xmm3, %xmm1
subl %esi, %edx
incq %rcx
cmpq %rcx, %r14
jne 0x13a79
xorps %xmm0, %xmm0
cvtsi2sd %edx, %xmm0
divsd %xmm0, %xmm1
jmp 0x13ac2
movsd 0x19156(%rip), %xmm1    # 0x2cc18
movapd %xmm1, %xmm0
addq $0x18, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/ImageProcs.cpp
// Function: calcNonoutlierMeanPixelValue(Image3ChnlReal&, double)
Float_t calcNonoutlierMeanPixelValue(Image3ChnlReal &rgbImage, Float_t outlierSigmas)
{
	int totalPixels = rgbImage.totalPixels();
	
	//Find the mean pixel value
	Float_t meanVal = 0;
	for (int c = 0; c < 3; c++)
		for (int i = 0; i < totalPixels; i++)
			meanVal += (*rgbImage[c])[i];
	meanVal /= (totalPixels * 3);
	
	//Find the variance of pixel value
	Float_t variance = 0;
	for (int c = 0; c < 3; c++)
		for (int i = 0; i < totalPixels; i++)
			variance += pow((*rgbImage[c])[i] - meanVal, 2);
	variance /= (totalPixels * 3);
	
	//Find the standard deviation of pixel value
	Float_t stddev = sqrt(variance);
	
	//Define outliers as some number of standard deviations
	Float_t minOutlier = meanVal - stddev * outlierSigmas;
	Float_t maxOutlier = meanVal + stddev * outlierSigmas;
	
	//Find the mean pixel value of nonoutliers
	Float_t nonOutlierMeanVal = 0;
	int numNonOutliers = 0;
	for (int c = 0; c < 3; c++)
		for (int i = 0; i < totalPixels; i++)
			if ((*rgbImage[c])[i] >= minOutlier && (*rgbImage[c])[i] <= maxOutlier)
			{
				nonOutlierMeanVal += (*rgbImage[c])[i];
				numNonOutliers++;
			}
	nonOutlierMeanVal /= numNonOutliers;
	
	return nonOutlierMeanVal;
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x18, %rsp
movapd %xmm0, %xmm3
movq %rdi, %rbx
movl 0x30(%rdi), %r14d
movq 0x38(%rdi), %rax
xorpd %xmm4, %xmm4
xorl %ecx, %ecx
testl %r14d, %r14d
jle 0x13b06
movq (%rax,%rcx,8), %rdx
movq 0x38(%rdx), %rdx
xorl %esi, %esi
addsd (%rdx,%rsi,8), %xmm4
incq %rsi
cmpq %rsi, %r14
jne 0x13af9
incq %rcx
cmpq $0x3, %rcx
jne 0x13aea
leal (%r14,%r14,2), %eax
cvtsi2sd %eax, %xmm1
divsd %xmm1, %xmm4
movq 0x38(%rbx), %rax
xorpd %xmm0, %xmm0
xorl %ecx, %ecx
testl %r14d, %r14d
jle 0x13b4d
movq (%rax,%rcx,8), %rdx
movq 0x38(%rdx), %rdx
xorl %esi, %esi
movsd (%rdx,%rsi,8), %xmm2
subsd %xmm4, %xmm2
mulsd %xmm2, %xmm2
addsd %xmm2, %xmm0
incq %rsi
cmpq %rsi, %r14
jne 0x13b34
incq %rcx
cmpq $0x3, %rcx
jne 0x13b25
divsd %xmm1, %xmm0
xorpd %xmm5, %xmm5
ucomisd %xmm5, %xmm0
jb 0x13b6a
sqrtsd %xmm0, %xmm0
jmp 0x13b8b
movsd %xmm3, 0x10(%rsp)
movsd %xmm4, 0x8(%rsp)
callq 0x3240
xorpd %xmm5, %xmm5
movsd 0x8(%rsp), %xmm4
movsd 0x10(%rsp), %xmm3
mulsd %xmm3, %xmm0
movapd %xmm4, %xmm1
subsd %xmm0, %xmm1
addsd %xmm4, %xmm0
movq 0x38(%rbx), %rax
xorl %ecx, %ecx
xorl %edx, %edx
testl %r14d, %r14d
jle 0x13be7
movq (%rax,%rcx,8), %rsi
movq 0x38(%rsi), %rsi
xorl %edi, %edi
movsd (%rsi,%rdi,8), %xmm2
movapd %xmm1, %xmm3
cmplepd %xmm2, %xmm3
movapd %xmm2, %xmm4
cmplepd %xmm0, %xmm4
andpd %xmm3, %xmm4
movq %xmm4, %r8
testb $0x1, %r8b
je 0x13bdc
addsd %xmm2, %xmm5
subl %r8d, %edx
incq %rdi
cmpq %rdi, %r14
jne 0x13bb2
incq %rcx
cmpq $0x3, %rcx
jne 0x13ba3
xorps %xmm0, %xmm0
cvtsi2sd %edx, %xmm0
divsd %xmm0, %xmm5
movapd %xmm5, %xmm0
addq $0x18, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: assertValidDims(int, int)
bool assertValidDims(int width, int height)
{
	int numOnes = 0;
	for (int bitPos = 0; bitPos < sizeof(int) * 8; bitPos++)
		if (width & (1 << bitPos))
			++numOnes;
	if (width < 16 || numOnes != 1)
		return false;
	
	numOnes = 0;
	for (int bitPos = 0; bitPos < sizeof(int) * 8; bitPos++)
		if (height & (1 << bitPos))
			++numOnes;
	if (height < 16 || numOnes != 1)
		return false;
	
	return true;
}

// Assembly:
xorl %eax, %eax
xorl %ecx, %ecx
btl %eax, %edi
adcl $0x0, %ecx
incl %eax
cmpl $0x20, %eax
jne 0x192a0
xorl %eax, %eax
cmpl $0x10, %edi
jl 0x192d8
cmpl $0x1, %ecx
jne 0x192d8
xorl %ecx, %ecx
xorl %eax, %eax
btl %ecx, %esi
adcl $0x0, %eax
incl %ecx
cmpl $0x20, %ecx
jne 0x192bd
cmpl $0x10, %esi
setge %cl
cmpl $0x1, %eax
sete %al
andb %cl, %al
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: doFourierAnalysis(double const*, double*, int, int)
bool doFourierAnalysis(const Float_t *grayImage, Float_t *fourierTransform, int width, int height)
{
	if (!assertValidDims(width, height))
		return false;
	
	int err = fft2dInit(log2(height), log2(width));
	if (err != 0)
		return false;
	
	if (!fourierTransform)
		fourierTransform = new Float_t(width * height * sizeof(Float_t));
	if (fourierTransform != grayImage)
		memcpy(fourierTransform, grayImage, width * height * sizeof(Float_t));
	rfft2d(fourierTransform, log2(height), log2(width));
	
	return true;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x10, %rsp
movl %ecx, %ebp
movl %edx, %r15d
movq %rsi, %rbx
movq %rdi, %r14
movl %edx, %edi
movl %ecx, %esi
callq 0x1929c
testb %al, %al
je 0x19a4d
cvtsi2sd %ebp, %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0x3510
cvttsd2si %xmm0, %r12
xorps %xmm0, %xmm0
cvtsi2sd %r15d, %xmm0
movsd %xmm0, (%rsp)
callq 0x3510
cvttsd2si %xmm0, %rsi
movq %r12, %rdi
callq 0x23658
testl %eax, %eax
sete %al
jne 0x19a4f
testq %rbx, %rbx
jne 0x19a04
movl $0x8, %edi
callq 0x3300
movq %rax, %rbx
movl %ebp, %eax
imull %r15d, %eax
cltq
shlq $0x3, %rax
movq %rax, %xmm0
punpckldq 0x11634(%rip), %xmm0 # xmm0 = xmm0[0],mem[0],xmm0[1],mem[1]
subpd 0x1163c(%rip), %xmm0    # 0x2b030
movapd %xmm0, %xmm1
unpckhpd %xmm0, %xmm1            # xmm1 = xmm1[1],xmm0[1]
addsd %xmm0, %xmm1
movsd %xmm1, (%rbx)
cmpq %r14, %rbx
je 0x19a1f
imull %r15d, %ebp
movslq %ebp, %rdx
shlq $0x3, %rdx
movq %rbx, %rdi
movq %r14, %rsi
callq 0x3260
movsd 0x8(%rsp), %xmm0
callq 0x3510
cvttsd2si %xmm0, %r14
movsd (%rsp), %xmm0
callq 0x3510
cvttsd2si %xmm0, %rdx
movq %rbx, %rdi
movq %r14, %rsi
callq 0x2407f
movb $0x1, %al
jmp 0x19a4f
xorl %eax, %eax
addq $0x10, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: doFourierAnalysis(ImageReal const&, ImageReal&)
bool doFourierAnalysis(const ImageReal &grayImage, ImageReal &fourierTransform)
{
	int width = grayImage.width(), height = grayImage.height();
	
	if (!assertValidDims(width, height))
		return false;
	
	int err = fft2dInit(log2(height), log2(width));
	if (err != 0)
		return false;
	
	fourierTransform = grayImage;
	rfft2d(fourierTransform.img(), log2(height), log2(width));
	
	return true;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x18, %rsp
movq %rsi, %rbx
movq %rdi, %r14
movl 0x28(%rdi), %ebp
movl 0x2c(%rdi), %r15d
movl %ebp, %edi
movl %r15d, %esi
callq 0x1929c
testb %al, %al
je 0x19afa
cvtsi2sd %r15d, %xmm0
movsd %xmm0, 0x10(%rsp)
callq 0x3510
cvttsd2si %xmm0, %r15
xorps %xmm0, %xmm0
cvtsi2sd %ebp, %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0x3510
cvttsd2si %xmm0, %rsi
movq %r15, %rdi
callq 0x23658
testl %eax, %eax
sete %al
jne 0x19afc
movq %rbx, %rdi
movq %r14, %rsi
callq 0xd036
movq 0x38(%rbx), %rbx
movsd 0x10(%rsp), %xmm0
callq 0x3510
cvttsd2si %xmm0, %r14
movsd 0x8(%rsp), %xmm0
callq 0x3510
cvttsd2si %xmm0, %rdx
movq %rbx, %rdi
movq %r14, %rsi
callq 0x2407f
movb $0x1, %al
jmp 0x19afc
xorl %eax, %eax
addq $0x18, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: multiplyFourierTransforms(double*, double*, int, int, double*)
bool multiplyFourierTransforms(Float_t *ft1, Float_t *ft2,
	int width, int height, Float_t *prod)
{
	if (!assertValidDims(width, height))
		return false;
	
	if (!prod)
		prod = new Float_t[width * height];
	rspect2dprod(ft1, ft2, prod, height, width);
	
	return true;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %r8, %r13
movl %ecx, %r12d
movl %edx, %ebp
movq %rsi, %r14
movq %rdi, %r15
movl %edx, %edi
movl %ecx, %esi
callq 0x1929c
movl %eax, %ebx
testb %al, %al
je 0x19b6d
testq %r13, %r13
jne 0x19b59
movl %r12d, %eax
imull %ebp, %eax
cltq
leaq (,%rax,8), %rcx
testl %eax, %eax
movq $-0x1, %rdi
cmovnsq %rcx, %rdi
callq 0x3030
movq %rax, %r13
movslq %r12d, %rcx
movslq %ebp, %r8
movq %r15, %rdi
movq %r14, %rsi
movq %r13, %rdx
callq 0x246f3
movl %ebx, %eax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: multiplyFourierTransformsAndInvert(double*, double*, int, int, double*)
bool multiplyFourierTransformsAndInvert(Float_t *ft1, Float_t *ft2,
	int width, int height, Float_t *convolvedImg)
{
	if (!multiplyFourierTransforms(ft1, ft2, width, height, convolvedImg))
		return false;
	
	int err = fft2dInit(log2(height), log2(width));
	if (err != 0)
		return false;
	
	//Invert the Fourier Transform
	rifft2d(convolvedImg, log2(height), log2(width));
	
	return true;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x10, %rsp
movq %r8, %rbx
movl %ecx, %r14d
movl %edx, %ebp
callq 0x19b07
testb %al, %al
je 0x19bff
cvtsi2sd %r14d, %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0x3510
cvttsd2si %xmm0, %r14
xorps %xmm0, %xmm0
cvtsi2sd %ebp, %xmm0
movsd %xmm0, (%rsp)
callq 0x3510
cvttsd2si %xmm0, %rsi
movq %r14, %rdi
callq 0x23658
testl %eax, %eax
sete %al
jne 0x19c01
movsd 0x8(%rsp), %xmm0
callq 0x3510
cvttsd2si %xmm0, %r14
movsd (%rsp), %xmm0
callq 0x3510
cvttsd2si %xmm0, %rdx
movq %rbx, %rdi
movq %r14, %rsi
callq 0x243ac
movb $0x1, %al
jmp 0x19c01
xorl %eax, %eax
addq $0x10, %rsp
popq %rbx
popq %r14
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: convolveImageWithFilter(double*, double*, int, int, double*)
bool convolveImageWithFilter(Float_t *img, Float_t *filterFT,
	int width, int height, Float_t *imgConv)
{
	if (!assertValidDims(width, height))
		return false;
	
	//Get the Fourier Transform of the image
	Float_t *imgFT = new Float_t[width * height];
	assert(doFourierAnalysis(img, imgFT, width, height));
	
	//Multiplty the Fourier Transforms and invert the FT
	return multiplyFourierTransformsAndInvert(imgFT, filterFT, width, height, imgConv);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %r8, %rbx
movl %ecx, %ebp
movl %edx, %r14d
movq %rsi, %r15
movq %rdi, %r12
movl %edx, %edi
movl %ecx, %esi
callq 0x1929c
testb %al, %al
je 0x19c8a
movl %ebp, %eax
imull %r14d, %eax
cltq
leaq (,%rax,8), %rcx
testl %eax, %eax
movq $-0x1, %rdi
cmovnsq %rcx, %rdi
callq 0x3030
movq %rax, %r13
movq %r12, %rdi
movq %rax, %rsi
movl %r14d, %edx
movl %ebp, %ecx
callq 0x1995b
testb %al, %al
je 0x19c9b
movq %r13, %rdi
movq %r15, %rsi
movl %r14d, %edx
movl %ebp, %ecx
movq %rbx, %r8
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x19b7e
xorl %eax, %eax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x13fa9(%rip), %rdi     # 0x2dc4b
leaq 0x13f0a(%rip), %rsi     # 0x2dbb3
leaq 0x13fc8(%rip), %rcx     # 0x2dc78
movl $0xfa, %edx
callq 0x31c0
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: embedImageInPow2Bounds(ImageReal const&, int, int, ImageReal&, int&, int&)
bool embedImageInPow2Bounds(const ImageReal &img, int newWidth, int newHeight, ImageReal &imgOut, int &offsetX, int &offsetY)
{
	if (!assertValidDims(newWidth, newHeight))
		return false;
	
	imgOut.resize(newWidth, newHeight);
	int width = img.width(), height = img.height();
	offsetX = (newWidth - width) / 2;
	offsetY = (newHeight - height) / 2;
	
	//By default, an embedded odd-length signal will be positioned such that the center pixel
	//resides at lower right, but I prefer that it be at upper left, so shift it by one.
	if (width % 2 == 1)
		offsetX++;
	if (height % 2 == 1)
		offsetY++;
	
	for (int y = 0; y < height; y++)
		for (int x = 0; x < width; x++)
			imgOut(x + offsetX, y + offsetY) = img(x, y);
	
	return true;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %r9, %r15
movq %r8, %rbx
movq %rcx, %r12
movl %edx, %ebp
movl %esi, %r14d
movq %rdi, %r13
movl %esi, %edi
movl %edx, %esi
callq 0x1929c
movb %al, 0x7(%rsp)
testb %al, %al
je 0x19d9c
movq (%r12), %rax
movq %r12, %rdi
movl %r14d, %esi
movl %ebp, %edx
callq *0x28(%rax)
movl 0x28(%r13), %eax
movl 0x2c(%r13), %ecx
subl %eax, %r14d
movl %r14d, %edx
shrl $0x1f, %edx
addl %r14d, %edx
sarl %edx
movl %edx, (%rbx)
subl %ecx, %ebp
movl %ebp, %edx
shrl $0x1f, %edx
addl %ebp, %edx
sarl %edx
movl %edx, (%r15)
movl %eax, %edx
andl $0x80000001, %edx       # imm = 0x80000001
cmpl $0x1, %edx
jne 0x19d2f
incl (%rbx)
movl %ecx, %edx
andl $0x80000001, %edx       # imm = 0x80000001
cmpl $0x1, %edx
jne 0x19d3f
incl (%r15)
testl %ecx, %ecx
jle 0x19d9c
movq 0x38(%r13), %rdx
movl 0x28(%r13), %esi
movq 0x38(%r12), %rdi
movl 0x28(%r12), %r8d
movl (%r15), %r9d
imull %r8d, %r9d
addl (%rbx), %r9d
xorl %r10d, %r10d
xorl %r11d, %r11d
movq %rax, %r14
movl %r10d, %ebp
movl %r9d, %r15d
testl %eax, %eax
jle 0x19d8e
movl %ebp, %r12d
movsd (%rdx,%r12,8), %xmm0
movl %r15d, %r12d
movsd %xmm0, (%rdi,%r12,8)
incl %r15d
incl %ebp
decq %r14
jne 0x19d72
incl %r11d
addl %r8d, %r9d
addl %esi, %r10d
cmpl %ecx, %r11d
jne 0x19d65
movb 0x7(%rsp), %al
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: embedImageInPow2Bounds(Image3ChnlReal const&, int, int, Image3ChnlReal&, int&, int&)
bool embedImageInPow2Bounds(const Image3ChnlReal &img, int newWidth, int newHeight, Image3ChnlReal &imgOut, int &offsetX, int &offsetY)
{
	if (!assertValidDims(newWidth, newHeight))
		return false;
	
	imgOut.resize(newWidth, newHeight);
	int width = img.width(), height = img.height();
	offsetX = (newWidth - width) / 2;
	offsetY = (newHeight - height) / 2;
	
	//By default, an embedded odd-length signal will be positioned such that the center pixel
	//resides at lower right, but I prefer that it be at upper left, so shift it by one.
	if (width % 2 == 1)
		offsetX++;
	if (height % 2 == 1)
		offsetY++;
	
	for (int c = 0; c < 3; c++)
		for (int y = 0; y < height; y++)
			for (int x = 0; x < width; x++)
				(*imgOut[c])(x + offsetX, y + offsetY) = (*img[c])(x, y);
	
	return true;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq %r9, 0x18(%rsp)
movq %r8, 0x10(%rsp)
movq %rcx, 0x20(%rsp)
movl %edx, %r13d
movl %esi, %ebp
movq %rdi, %r12
movl %esi, %edi
movl %edx, %esi
callq 0x1929c
movb %al, 0xf(%rsp)
testb %al, %al
je 0x19ed2
movq 0x20(%rsp), %rdi
movq (%rdi), %rax
movl %ebp, %esi
movl %r13d, %edx
callq *0x28(%rax)
movl 0x28(%r12), %r14d
movl 0x2c(%r12), %ecx
subl %r14d, %ebp
movl %ebp, %eax
shrl $0x1f, %eax
addl %ebp, %eax
sarl %eax
movq 0x10(%rsp), %rdx
movl %eax, (%rdx)
subl %ecx, %r13d
movl %r13d, %eax
shrl $0x1f, %eax
addl %r13d, %eax
sarl %eax
movq 0x18(%rsp), %rdx
movl %eax, (%rdx)
movl %r14d, %eax
andl $0x80000001, %eax       # imm = 0x80000001
cmpl $0x1, %eax
jne 0x19e3f
movq 0x10(%rsp), %rax
incl (%rax)
movl %ecx, %eax
andl $0x80000001, %eax       # imm = 0x80000001
cmpl $0x1, %eax
jne 0x19e52
movq 0x18(%rsp), %rax
incl (%rax)
xorl %edx, %edx
testl %ecx, %ecx
jle 0x19ec9
movq 0x38(%r12), %rsi
movq 0x20(%rsp), %rax
movq 0x38(%rax), %rdi
movq 0x10(%rsp), %rax
movl (%rax), %r8d
movq 0x18(%rsp), %rax
movl (%rax), %r9d
xorl %r10d, %r10d
testl %r14d, %r14d
jle 0x19ebe
movq (%rsi,%rdx,8), %rax
movq 0x38(%rax), %r11
movq (%rdi,%rdx,8), %rbx
movq 0x38(%rbx), %r13
movl 0x28(%rbx), %ebp
imull %r9d, %ebp
addl %r8d, %ebp
movl 0x28(%rax), %ebx
imull %r10d, %ebx
movq %r14, %rax
movl %ebx, %r15d
movsd (%r11,%r15,8), %xmm0
movl %ebp, %r15d
movsd %xmm0, (%r13,%r15,8)
incl %ebp
incl %ebx
decq %rax
jne 0x19ea2
incl %r10d
incl %r9d
cmpl %ecx, %r10d
jne 0x19e79
incq %rdx
cmpq $0x3, %rdx
jne 0x19e54
movb 0xf(%rsp), %al
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: cropImage(ImageReal const&, int, int, int, int, ImageReal&)
void cropImage(const ImageReal &img, int newWidth, int newHeight, int offsetX, int offsetY, ImageReal &imgOut)
{
	int width = img.width(), height = img.height();
	assert(offsetX + newWidth <= width && offsetY + newHeight <= height);
	
	imgOut.resize(newWidth, newHeight);
	
	for (int y = 0; y < newHeight; y++)
		for (int x = 0; x < newWidth; x++)
			imgOut(x, y) = img(x + offsetX, y + offsetY);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movl %ecx, %ebx
movl %esi, %r15d
leal (%rbx,%r15), %eax
cmpl 0x28(%rdi), %eax
jg 0x19f8f
movl %r8d, %r12d
movl %edx, %r14d
movq %rdi, %rbp
leal (%r12,%r14), %eax
cmpl 0x2c(%rdi), %eax
jg 0x19f8f
movq %r9, %r13
movq (%r9), %rax
movq %r9, %rdi
movl %r15d, %esi
movl %r14d, %edx
callq *0x28(%rax)
testl %r14d, %r14d
jle 0x19f80
movq 0x38(%rbp), %rax
movl 0x28(%rbp), %ecx
movq 0x38(%r13), %rdx
movl 0x28(%r13), %esi
movl %r15d, %edi
imull %ecx, %r12d
addl %r12d, %ebx
xorl %r8d, %r8d
xorl %r9d, %r9d
movq %rdi, %r10
movl %r8d, %r11d
movl %ebx, %ebp
testl %r15d, %r15d
jle 0x19f73
movl %ebp, %r12d
movsd (%rax,%r12,8), %xmm0
movl %r11d, %r12d
movsd %xmm0, (%rdx,%r12,8)
incl %ebp
incl %r11d
decq %r10
jne 0x19f57
incl %r9d
addl %ecx, %ebx
addl %esi, %r8d
cmpl %r14d, %r9d
jne 0x19f4a
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x13d2a(%rip), %rdi     # 0x2dcc0
leaq 0x13c16(%rip), %rsi     # 0x2dbb3
leaq 0x13d59(%rip), %rcx     # 0x2dcfd
movl $0x134, %edx            # imm = 0x134
callq 0x31c0
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: cropImage(Image3ChnlReal const&, int, int, int, int, Image3ChnlReal&)
void cropImage(const Image3ChnlReal &img, int newWidth, int newHeight, int offsetX, int offsetY, Image3ChnlReal &imgOut)
{
	int width = img.width(), height = img.height();
	assert(offsetX + newWidth <= width && offsetY + newHeight <= height);
	
	imgOut.resize(newWidth, newHeight);
	
	for (int c = 0; c < 3; c++)
		for (int y = 0; y < newHeight; y++)
			for (int x = 0; x < newWidth; x++)
				(*imgOut[c])(x, y) = (*img[c])(x + offsetX, y + offsetY);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %r9, (%rsp)
movl %ecx, %r15d
movl %esi, %r13d
leal (%r15,%r13), %eax
cmpl 0x28(%rdi), %eax
jg 0x1a08c
movl %edx, %r12d
leal (%r8,%r12), %eax
cmpl 0x2c(%rdi), %eax
jg 0x1a08c
movq %rdi, 0x8(%rsp)
movq %r8, 0x10(%rsp)
movq (%rsp), %rdi
movq (%rdi), %rax
movl %r13d, %esi
movl %r12d, %edx
callq *0x28(%rax)
movl %r13d, %ebx
xorl %ecx, %ecx
testl %r12d, %r12d
jle 0x1a074
movq 0x8(%rsp), %rax
movq 0x38(%rax), %rdx
movq (%rsp), %rax
movq 0x38(%rax), %rsi
xorl %edi, %edi
movq 0x10(%rsp), %rax
movl %eax, %r8d
testl %r13d, %r13d
jle 0x1a06a
movq (%rdx,%rcx,8), %rax
movq 0x38(%rax), %r9
movq (%rsi,%rcx,8), %r11
movq 0x38(%r11), %r10
movl 0x28(%r11), %r11d
imull %edi, %r11d
movl 0x28(%rax), %r14d
imull %r8d, %r14d
addl %r15d, %r14d
movq %rbx, %rax
movl %r14d, %ebp
movsd (%r9,%rbp,8), %xmm0
movl %r11d, %ebp
movsd %xmm0, (%r10,%rbp,8)
incl %r14d
incl %r11d
decq %rax
jne 0x1a04d
incl %edi
incl %r8d
cmpl %r12d, %edi
jne 0x1a022
incq %rcx
cmpq $0x3, %rcx
jne 0x1a002
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x13c2d(%rip), %rdi     # 0x2dcc0
leaq 0x13b19(%rip), %rsi     # 0x2dbb3
leaq 0x13c9f(%rip), %rcx     # 0x2dd40
movl $0x140, %edx            # imm = 0x140
callq 0x31c0
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: shiftImage(ImageReal const&, int, int, ImageReal&)
void shiftImage(const ImageReal &img, int shiftX, int shiftY, ImageReal &imgOut)
{
	assert(&img != &imgOut);
	
	if (shiftX < 0)
		shiftX += img.width();
	shiftX %= img.width();
	if (shiftY < 0)
		shiftY += img.height();
	shiftY %= img.height();
	
	imgOut.resize(img.width(), img.height());
	int width = img.width(), height = img.height();
	
	int xx, yy;
	for (int y = 0; y < height; y++)
	{
		yy = (y + shiftY) % height;
		for (int x = 0; x < width; x++)
		{
			xx = (x + shiftX) % width;
			imgOut(x, y) = img(xx, yy);
		}
	}
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
cmpq %rcx, %rdi
je 0x1a160
movq %rcx, %r14
movl %esi, %ebp
movq %rdi, %r15
movl 0x28(%rdi), %r12d
movl 0x2c(%rdi), %ecx
movl %edx, %eax
sarl $0x1f, %eax
andl %ecx, %eax
addl %edx, %eax
cltd
idivl %ecx
movl %edx, %ebx
movq (%r14), %rax
movq %r14, %rdi
movl %r12d, %esi
movl %ecx, %edx
callq *0x28(%rax)
movl 0x2c(%r15), %edi
testl %edi, %edi
jle 0x1a157
movl 0x28(%r15), %r8d
movl %ebp, %eax
sarl $0x1f, %eax
andl %r12d, %eax
addl %ebp, %eax
cltd
idivl %r12d
movl %edx, %ecx
movq 0x38(%r15), %r9
movq 0x38(%r14), %r10
movl 0x28(%r14), %r11d
xorl %r14d, %r14d
xorl %r15d, %r15d
testl %r8d, %r8d
jle 0x1a14c
leal (%r15,%rbx), %eax
cltd
idivl %edi
movl %edx, %esi
imull %r8d, %esi
xorl %r12d, %r12d
leal (%rcx,%r12), %eax
cltd
idivl %r8d
addl %esi, %edx
movsd (%r9,%rdx,8), %xmm0
leal (%r14,%r12), %eax
movsd %xmm0, (%r10,%rax,8)
incq %r12
cmpq %r12, %r8
jne 0x1a12a
incl %r15d
addq %r11, %r14
cmpl %edi, %r15d
jne 0x1a115
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
leaq 0x13c26(%rip), %rdi     # 0x2dd8d
leaq 0x13a45(%rip), %rsi     # 0x2dbb3
leaq 0x13c28(%rip), %rcx     # 0x2dd9d
movl $0x14c, %edx            # imm = 0x14C
callq 0x31c0
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FourierTransform.cpp
// Function: shiftImage(Image3ChnlReal const&, int, int, Image3ChnlReal&)
void shiftImage(const Image3ChnlReal &img, int shiftX, int shiftY, Image3ChnlReal &imgOut)
{
	assert(&img != &imgOut);
	
	if (shiftX < 0)
		shiftX += img.width();
	shiftX %= img.width();
	if (shiftY < 0)
		shiftY += img.height();
	shiftY %= img.height();
	
	imgOut.resize(img.width(), img.height());
	int width = img.width(), height = img.height();
	
	int xx, yy;
	for (int c = 0; c < 3; c++)
		for (int y = 0; y < height; y++)
		{
			yy = (y + shiftY) % height;
			for (int x = 0; x < width; x++)
			{
				xx = (x + shiftX) % width;
				(*imgOut[c])(x, y) = (*img[c])(xx, yy);
			}
		}
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
cmpq %rcx, %rdi
je 0x1a270
movl %edx, %r10d
movq %rdi, %rbx
movl %esi, %eax
sarl $0x1f, %eax
movl 0x28(%rdi), %r8d
movl 0x2c(%rdi), %r9d
andl %r8d, %eax
addl %esi, %eax
cltd
idivl %r8d
movl %edx, %r15d
movl %r10d, %eax
sarl $0x1f, %eax
andl %r9d, %eax
addl %r10d, %eax
cltd
idivl %r9d
movl %edx, %r12d
movq (%rcx), %rax
movq %rcx, 0x10(%rsp)
movq %rcx, %rdi
movl %r8d, %esi
movl %r9d, %edx
callq *0x28(%rax)
movl 0x28(%rbx), %esi
movq %rbx, 0x8(%rsp)
movl 0x2c(%rbx), %edi
xorl %r8d, %r8d
testl %edi, %edi
jle 0x1a258
movq 0x8(%rsp), %rax
movq 0x38(%rax), %r9
movq 0x10(%rsp), %rax
movq 0x38(%rax), %r10
xorl %r11d, %r11d
testl %esi, %esi
jle 0x1a250
leal (%r11,%r12), %eax
cltd
idivl %edi
movl %edx, %ecx
movq (%r9,%r8,8), %rax
movq 0x38(%rax), %r13
imull 0x28(%rax), %ecx
movq (%r10,%r8,8), %rax
movq 0x38(%rax), %rbp
movl 0x28(%rax), %ebx
imull %r11d, %ebx
xorl %r14d, %r14d
leal (%r15,%r14), %eax
cltd
idivl %esi
addl %ecx, %edx
movsd (%r13,%rdx,8), %xmm0
leal (%rbx,%r14), %eax
movsd %xmm0, (%rbp,%rax,8)
incq %r14
cmpq %r14, %rsi
jne 0x1a22e
incl %r11d
cmpl %edi, %r11d
jne 0x1a203
incq %r8
cmpq $0x3, %r8
jne 0x1a1ea
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x13b16(%rip), %rdi     # 0x2dd8d
leaq 0x13935(%rip), %rsi     # 0x2dbb3
leaq 0x13b52(%rip), %rcx     # 0x2ddd7
movl $0x166, %edx            # imm = 0x166
callq 0x31c0
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::CreateScratch()
void SteerablePyramid::CreateScratch()
{
	//Note that scratch_2D_A_ will be memory leaked when the last instance of SteerablePyramid is deleted!!!
	//Reasonably simple to fix, but this is just a dinky little experimental program.  Tsk tsk.
	if (scratchSize_ < width_ * height_)
	{
		if (scratch_2D_A_)
			delete [] scratch_2D_A_;
		if (scratch_2D_B_)
			delete [] scratch_2D_B_;
		
		scratchSize_ = width_ * height_;
		
		scratch_2D_A_ = new Float_t[scratchSize_];
		assert(scratch_2D_A_);
		scratch_2D_B_ = new Float_t[scratchSize_];
		assert(scratch_2D_B_);
	}
}

// Assembly:
movl 0x34(%rdi), %eax
imull 0x30(%rdi), %eax
cmpl %eax, 0x1e9a1(%rip)     # 0x39af0
jge 0x1b1bc
pushq %rbx
movq %rdi, %rbx
movq 0x1e984(%rip), %rdi     # 0x39ae0
testq %rdi, %rdi
je 0x1b166
callq 0x33c0
movq 0x1e97b(%rip), %rdi     # 0x39ae8
testq %rdi, %rdi
je 0x1b177
callq 0x33c0
movslq 0x30(%rbx), %rax
movslq 0x34(%rbx), %rcx
imulq %rax, %rcx
movl %ecx, 0x1e967(%rip)     # 0x39af0
movq %rcx, %rax
shlq $0x3, %rax
testl %ecx, %ecx
movq $-0x1, %rbx
cmovnsq %rax, %rbx
movq %rbx, %rdi
callq 0x3030
movq %rax, 0x1e934(%rip)     # 0x39ae0
movq %rbx, %rdi
callq 0x3030
movq %rax, 0x1e92d(%rip)     # 0x39ae8
popq %rbx
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::operator=(SteerablePyramid const&)
const SteerablePyramid& SteerablePyramid::operator=(const SteerablePyramid& sp)
{
	if (&sp == this)
		return *this;
	
	Init(sp.width_, sp.height_);
	
	useFTconvolution_ = sp.useFTconvolution_;
	
	int totalPixels = width_ * height_;
	memcpy(h0Coeffs_, sp.h0Coeffs_, totalPixels * sizeof(Float_t));
	
	for (int i = 0; i < numLevels_; i++)
	{
		int totalPixels = (width_ * height_) / (int)(pow(4, i));
		memcpy(b1Coeffs_[i], sp.b1Coeffs_[i], totalPixels * sizeof(Float_t));
		memcpy(b2Coeffs_[i], sp.b2Coeffs_[i], totalPixels * sizeof(Float_t));
		memcpy(b3Coeffs_[i], sp.b3Coeffs_[i], totalPixels * sizeof(Float_t));
		memcpy(b4Coeffs_[i], sp.b4Coeffs_[i], totalPixels * sizeof(Float_t));
		memcpy(l1Coeffs_[i], sp.l1Coeffs_[i], (totalPixels / 4) * sizeof(Float_t));
	}
	
	return *this;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
cmpq %rdi, %rsi
je 0x1b327
movq %rsi, %r14
movl 0x30(%rsi), %esi
movl 0x34(%r14), %edx
movq %rbx, %rdi
callq 0x1b4ce
movb 0x48(%r14), %al
movb %al, 0x48(%rbx)
movslq 0x30(%rbx), %rax
movslq 0x34(%rbx), %rdx
imulq %rax, %rdx
movq (%rbx), %rdi
movq (%r14), %rsi
shlq $0x3, %rdx
callq 0x3260
cmpl $0x0, 0x38(%rbx)
jle 0x1b327
xorl %r12d, %r12d
xorl %r13d, %r13d
movl 0x34(%rbx), %ebp
imull 0x30(%rbx), %ebp
movl %r12d, %eax
xorps %xmm0, %xmm0
cvtsi2sd %rax, %xmm0
callq 0x32a0
cvttsd2si %xmm0, %ecx
movl %ebp, %eax
cltd
idivl %ecx
movq 0x8(%rbx), %rcx
movq (%rcx,%r13,8), %rdi
movq 0x8(%r14), %rcx
movq (%rcx,%r13,8), %rsi
movslq %eax, %rbp
leaq (,%rbp,8), %r15
movq %r15, %rdx
callq 0x3260
movq 0x10(%rbx), %rax
movq (%rax,%r13,8), %rdi
movq 0x10(%r14), %rax
movq (%rax,%r13,8), %rsi
movq %r15, %rdx
callq 0x3260
movq 0x18(%rbx), %rax
movq (%rax,%r13,8), %rdi
movq 0x18(%r14), %rax
movq (%rax,%r13,8), %rsi
movq %r15, %rdx
callq 0x3260
movq 0x20(%rbx), %rax
movq (%rax,%r13,8), %rdi
movq 0x20(%r14), %rax
movq (%rax,%r13,8), %rsi
movq %r15, %rdx
callq 0x3260
movq 0x28(%rbx), %rax
movq (%rax,%r13,8), %rdi
movq 0x28(%r14), %rax
leal 0x3(%rbp), %ecx
testl %ebp, %ebp
cmovnsl %ebp, %ecx
movq (%rax,%r13,8), %rsi
sarl $0x2, %ecx
movslq %ecx, %rdx
shlq $0x3, %rdx
callq 0x3260
incq %r13
movslq 0x38(%rbx), %rax
addl $0x2, %r12d
cmpq %rax, %r13
jl 0x1b261
movq %rbx, %rax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::Destroy()
void SteerablePyramid::Destroy()
{
	if (yLookup_)
		delete [] yLookup_;
	yLookup_ = NULL;
	
	if (h0Coeffs_)
		delete [] h0Coeffs_;
	h0Coeffs_ = NULL;
	
	if (b1Coeffs_)
	{
		for (int i = 0; i < numLevels_; i++)
			delete [] b1Coeffs_[i];
		delete [] b1Coeffs_;
	}
	b1Coeffs_ = NULL;
	
	if (b2Coeffs_)
	{
		for (int i = 0; i < numLevels_; i++)
			delete [] b2Coeffs_[i];
		delete [] b2Coeffs_;
	}
	b2Coeffs_ = NULL;
	
	if (b3Coeffs_)
	{
		for (int i = 0; i < numLevels_; i++)
			delete [] b3Coeffs_[i];
		delete [] b3Coeffs_;
	}
	b3Coeffs_ = NULL;
	
	if (b4Coeffs_)
	{
		for (int i = 0; i < numLevels_; i++)
			delete [] b4Coeffs_[i];
		delete [] b4Coeffs_;
	}
	b4Coeffs_ = NULL;
	
	if (l1Coeffs_)
	{
		for (int i = 0; i < numLevels_; i++)
			delete [] l1Coeffs_[i];
		delete [] l1Coeffs_;
	}
	l1Coeffs_ = NULL;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
movq 0x40(%rdi), %rdi
testq %rdi, %rdi
je 0x1b355
callq 0x33c0
movq $0x0, 0x40(%rbx)
movq (%rbx), %rdi
testq %rdi, %rdi
je 0x1b36a
callq 0x33c0
movq $0x0, (%rbx)
cmpq $0x0, 0x8(%rbx)
je 0x1b3ad
cmpl $0x0, 0x38(%rbx)
jle 0x1b39f
xorl %r14d, %r14d
movq 0x8(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x1b393
callq 0x33c0
incq %r14
movslq 0x38(%rbx), %rax
cmpq %rax, %r14
jl 0x1b381
movq 0x8(%rbx), %rdi
testq %rdi, %rdi
je 0x1b3ad
callq 0x33c0
movq $0x0, 0x8(%rbx)
cmpq $0x0, 0x10(%rbx)
je 0x1b3f1
cmpl $0x0, 0x38(%rbx)
jle 0x1b3e3
xorl %r14d, %r14d
movq 0x10(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x1b3d7
callq 0x33c0
incq %r14
movslq 0x38(%rbx), %rax
cmpq %rax, %r14
jl 0x1b3c5
movq 0x10(%rbx), %rdi
testq %rdi, %rdi
je 0x1b3f1
callq 0x33c0
movq $0x0, 0x10(%rbx)
cmpq $0x0, 0x18(%rbx)
je 0x1b435
cmpl $0x0, 0x38(%rbx)
jle 0x1b427
xorl %r14d, %r14d
movq 0x18(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x1b41b
callq 0x33c0
incq %r14
movslq 0x38(%rbx), %rax
cmpq %rax, %r14
jl 0x1b409
movq 0x18(%rbx), %rdi
testq %rdi, %rdi
je 0x1b435
callq 0x33c0
movq $0x0, 0x18(%rbx)
cmpq $0x0, 0x20(%rbx)
je 0x1b479
cmpl $0x0, 0x38(%rbx)
jle 0x1b46b
xorl %r14d, %r14d
movq 0x20(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x1b45f
callq 0x33c0
incq %r14
movslq 0x38(%rbx), %rax
cmpq %rax, %r14
jl 0x1b44d
movq 0x20(%rbx), %rdi
testq %rdi, %rdi
je 0x1b479
callq 0x33c0
movq $0x0, 0x20(%rbx)
cmpq $0x0, 0x28(%rbx)
je 0x1b4bd
cmpl $0x0, 0x38(%rbx)
jle 0x1b4af
xorl %r14d, %r14d
movq 0x28(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x1b4a3
callq 0x33c0
incq %r14
movslq 0x38(%rbx), %rax
cmpq %rax, %r14
jl 0x1b491
movq 0x28(%rbx), %rdi
testq %rdi, %rdi
je 0x1b4bd
callq 0x33c0
movq $0x0, 0x28(%rbx)
addq $0x8, %rsp
popq %rbx
popq %r14
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::GetFilters(std::vector<std::vector<double, std::allocator<double>>, std::allocator<std::vector<double, std::allocator<double>>>>&)
void SteerablePyramid::GetFilters(vector<vector<Float_t> > &filters)
{
	InitFilters();
	
	filters.clear();
	filters.resize(7);
	
	for (int i = 0; i < 7; i++)
	{
		filters[i].resize(skFilterDims[i] * skFilterDims[i]);
		for (int j = 0; j < skFilterDims[i] * skFilterDims[i]; j++)
			filters[i][j] = skFilters[i][j];
	}
}

// Assembly:
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
movq %rdi, %rbx
callq 0x1aff0
movq (%rbx), %rsi
movq %rbx, %rdi
callq 0x22564
movl $0x7, %esi
movq %rbx, %rdi
callq 0x21c70
xorl %r14d, %r14d
leaq 0x1d355(%rip), %r15     # 0x39b00
leaq 0x1d36e(%rip), %r12     # 0x39b20
leaq (%r14,%r14,2), %r13
leaq (,%r13,8), %rdi
addq (%rbx), %rdi
movl (%r15,%r14,4), %esi
imull %esi, %esi
callq 0xa084
movl (%r15,%r14,4), %eax
testl %eax, %eax
je 0x1c7fd
imull %eax, %eax
movq (%r12,%r14,8), %rcx
movq (%rbx), %rdx
movq (%rdx,%r13,8), %rdx
cmpl $0x1, %eax
adcl $0x0, %eax
xorl %esi, %esi
movsd (%rcx,%rsi,8), %xmm0
movsd %xmm0, (%rdx,%rsi,8)
incq %rsi
cmpq %rsi, %rax
jne 0x1c7eb
incq %r14
cmpq $0x7, %r14
jne 0x1c7b2
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::DetermineNumLevels()
void SteerablePyramid::DetermineNumLevels()
{
	numLevels_ = 0;
	int shorterDim = min(width_, height_);
	while (shorterDim > skFilterDims[2])	//Limit is determined by the dimensions of the L1 filter (probably 17x17)
	{
		numLevels_++;
		shorterDim /= 2;
	}
}

// Assembly:
movl $0x0, 0x38(%rdi)
movl 0x30(%rdi), %eax
movl 0x34(%rdi), %ecx
cmpl %eax, %ecx
cmovll %ecx, %eax
movl 0x1d2e0(%rip), %ecx     # 0x39b08
cmpl %ecx, %eax
jle 0x1c843
movl 0x38(%rdi), %edx
movl %eax, %esi
incl %edx
shrl $0x1f, %esi
addl %eax, %esi
sarl %esi
movl %esi, %eax
cmpl %ecx, %esi
jg 0x1c831
movl %edx, 0x38(%rdi)
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::SetAllZero()
void SteerablePyramid::SetAllZero()
{
	for (int i = 0; i < numLevels_; i++)
	{
		int horDim = width_ / pow(2, i);
		int verDim = height_ / pow(2, i);
		int totalPixels = horDim * verDim;
		int totalBytes = totalPixels * sizeof(Float_t);
		
		if (i == 0)
			memset(h0Coeffs_, 0, totalBytes);
		memset(b1Coeffs_[i], 0, totalBytes);
		memset(b2Coeffs_[i], 0, totalBytes);
		memset(b3Coeffs_[i], 0, totalBytes);
		memset(b4Coeffs_[i], 0, totalBytes);
		memset(l1Coeffs_[i], 0, totalBytes / 4);
	}
}

// Assembly:
cmpl $0x0, 0x38(%rdi)
jle 0x1c94b
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x18, %rsp
movq %rdi, %rbx
xorl %r15d, %r15d
xorps %xmm0, %xmm0
cvtsi2sdl 0x30(%rbx), %xmm0
movsd %xmm0, 0x8(%rsp)
xorps %xmm0, %xmm0
cvtsi2sd %r15d, %xmm0
movsd %xmm0, 0x10(%rsp)
callq 0x32a0
movsd 0x8(%rsp), %xmm1
divsd %xmm0, %xmm1
cvttsd2si %xmm1, %ebp
xorps %xmm0, %xmm0
cvtsi2sdl 0x34(%rbx), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x10(%rsp), %xmm0
callq 0x32a0
movsd 0x8(%rsp), %xmm1
divsd %xmm0, %xmm1
cvttsd2si %xmm1, %eax
imull %ebp, %eax
shll $0x3, %eax
movslq %eax, %r14
testq %r15, %r15
jne 0x1c8d3
movq (%rbx), %rdi
xorl %esi, %esi
movq %r14, %rdx
callq 0x3190
movq 0x8(%rbx), %rax
movq (%rax,%r15,8), %rdi
xorl %esi, %esi
movq %r14, %rdx
callq 0x3190
movq 0x10(%rbx), %rax
movq (%rax,%r15,8), %rdi
xorl %esi, %esi
movq %r14, %rdx
callq 0x3190
movq 0x18(%rbx), %rax
movq (%rax,%r15,8), %rdi
xorl %esi, %esi
movq %r14, %rdx
callq 0x3190
movq 0x20(%rbx), %rax
movq (%rax,%r15,8), %rdi
xorl %esi, %esi
movq %r14, %rdx
callq 0x3190
movq 0x28(%rbx), %rax
movq (%rax,%r15,8), %rdi
sarl $0x2, %r14d
movslq %r14d, %rdx
xorl %esi, %esi
callq 0x3190
incq %r15
movslq 0x38(%rbx), %rax
cmpq %rax, %r15
jl 0x1c862
addq $0x18, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::AccumulateHistograms(double*, int, int, std::vector<int, std::allocator<int>>&, double)
void SteerablePyramid::AccumulateHistograms(Float_t *coeffs, int horDim, int verDim, vector<int> &histogram, Float_t binSpan)
{
	int numBins = histogram.size();
	
	int totalPixels = horDim * verDim;
	int bin;
	for (int i = 0; i < totalPixels; i++)
	{
		bin = fabs(coeffs[i]) / binSpan;
		if (bin >= numBins)
		{
			cout << "Out of bounds histogram: "
					<< fabs(coeffs[i]) << "   "
					<< binSpan << "   "
					<< fabs(coeffs[i]) / binSpan
					<< "   "
					<< bin << "   "
					<< numBins << endl;
			bin = numBins - 1;
		}
		histogram[bin]++;
	}
}

// Assembly:
imull %ecx, %edx
testl %edx, %edx
jle 0x21aff
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq %rsi, %r14
movq 0x8(%r8), %r15
subq (%r8), %r15
shrq $0x2, %r15
leal -0x1(%r15), %eax
movl %eax, 0x14(%rsp)
movl %edx, %edx
xorl %ebx, %ebx
movapd 0x966f(%rip), %xmm2     # 0x2b010
movsd %xmm0, 0x8(%rsp)
movq %rdx, 0x20(%rsp)
movq %r8, 0x18(%rsp)
movsd (%r14,%rbx,8), %xmm1
andpd %xmm2, %xmm1
divsd %xmm0, %xmm1
cvttsd2si %xmm1, %r13d
cmpl %r15d, %r13d
jl 0x21adc
movl $0x19, %edx
movq 0x155f7(%rip), %r12     # 0x36fd0
movq %r12, %rdi
leaq 0xda06(%rip), %rsi      # 0x2f3e9
callq 0x3360
movsd (%r14,%rbx,8), %xmm0
movapd 0x961a(%rip), %xmm1     # 0x2b010
andpd %xmm1, %xmm0
movq %r12, %rdi
callq 0x3480
movq %rax, %r12
movl $0x3, %edx
movq %rax, %rdi
leaq 0x9efb(%rip), %rbp      # 0x2b90f
movq %rbp, %rsi
callq 0x3360
movq %r12, %rdi
movsd 0x8(%rsp), %xmm0
callq 0x3480
movq %rax, %r12
movl $0x3, %edx
movq %rax, %rdi
movq %rbp, %rsi
callq 0x3360
movsd (%r14,%rbx,8), %xmm0
andpd 0x95c5(%rip), %xmm0     # 0x2b010
divsd 0x8(%rsp), %xmm0
movq %r12, %rdi
callq 0x3480
movq %rax, %r12
movl $0x3, %edx
movq %rax, %rdi
movq %rbp, %rsi
callq 0x3360
movq %r12, %rdi
movl %r13d, %esi
callq 0x34e0
movq %rax, %r12
movl $0x3, %edx
movq %rax, %rdi
movq %rbp, %rsi
callq 0x3360
movq %r12, %rdi
movl %r15d, %esi
callq 0x34e0
movq %rax, %r12
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %r12, %rdi
movl $0xa, %esi
callq 0x32e0
movsbl %al, %esi
movq %r12, %rdi
callq 0x3040
movq %rax, %rdi
callq 0x31f0
movapd 0x9549(%rip), %xmm2     # 0x2b010
movq 0x20(%rsp), %rdx
movq 0x18(%rsp), %r8
movsd 0x8(%rsp), %xmm0
movl 0x14(%rsp), %r13d
movslq %r13d, %rax
movq (%r8), %rcx
incl (%rcx,%rax,4)
incq %rbx
cmpq %rbx, %rdx
jne 0x219b1
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: FilterFT::~FilterFT()
~FilterFT()
	
	{
		for (int i = 0; i < 7; i++)
			for (int j = 0; j < filtersFTs[i].size(); j++)
			{
				assert(filtersFTs[i][j]);
				delete [] filtersFTs[i][j];
				delete [] filtersInvFTs[i][j];
			}
	}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
leaq 0x8(%rdi), %r14
leaq 0xb0(%rdi), %rax
movq %rax, (%rsp)
xorl %r12d, %r12d
leaq (%r12,%r12,2), %rcx
movq (%r14,%rcx,8), %rax
cmpq %rax, 0x8(%r14,%rcx,8)
je 0x21bed
leaq (%r14,%rcx,8), %r13
movq (%rsp), %rdx
leaq (%rdx,%rcx,8), %rbp
xorl %r15d, %r15d
movq (%rax,%r15,8), %rdi
testq %rdi, %rdi
je 0x21c50
callq 0x33c0
movq (%rbp), %rax
movq (%rax,%r15,8), %rdi
testq %rdi, %rdi
je 0x21bd6
callq 0x33c0
incq %r15
movq (%r13), %rax
movq 0x8(%r13), %rcx
subq %rax, %rcx
sarq $0x3, %rcx
cmpq %r15, %rcx
ja 0x21bb2
incq %r12
cmpq $0x7, %r12
jne 0x21b94
movl $0x150, %r14d           # imm = 0x150
movq -0x10(%rbx,%r14), %rdi
testq %rdi, %rdi
je 0x21c12
movq (%rbx,%r14), %rsi
subq %rdi, %rsi
callq 0x3310
addq $-0x18, %r14
cmpq $0xa8, %r14
jne 0x21bfc
movl $0xa8, %r14d
movq -0x10(%rbx,%r14), %rdi
testq %rdi, %rdi
je 0x21c3b
movq (%rbx,%r14), %rsi
subq %rdi, %rsi
callq 0x3310
addq $-0x18, %r14
jne 0x21c25
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0xd7dc(%rip), %rdi      # 0x2f433
leaq 0xc65c(%rip), %rsi      # 0x2e2ba
leaq 0xd7df(%rip), %rcx      # 0x2f444
movl $0x37, %edx
callq 0x31c0
nop
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: FilterFT::operator=(FilterFT const&)
const FilterFT& operator=(const FilterFT &fft)
	{
		if (&fft == this)
			return *this;
		
		width_ = fft.width_;
		height_ = fft.height_;
		
		for (int i = 0; i < 7; i++)
		{
			filtersFTs[i].resize(fft.filtersFTs[i].size());
			filtersInvFTs[i].resize(fft.filtersInvFTs[i].size());
			for (int j = 0; j < filtersFTs[i].size(); j++)
			{
				int totalPixels = (width_ * height_) / pow(4, j);
				filtersFTs[i][j] = new Float_t[totalPixels];
				filtersInvFTs[i][j] = new Float_t[totalPixels];
				memcpy(filtersFTs[i][j], fft.filtersFTs[i][j], totalPixels * sizeof(Float_t));
				memcpy(filtersInvFTs[i][j], fft.filtersInvFTs[i][j], totalPixels * sizeof(Float_t));
			}
		}
		
		return *this;
	}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x48, %rsp
movq %rdi, (%rsp)
cmpq %rdi, %rsi
je 0x2244a
movl (%rsi), %eax
movq (%rsp), %rcx
movl %eax, (%rcx)
movl 0x4(%rsi), %eax
movl %eax, 0x4(%rcx)
leaq 0x8(%rcx), %rax
movq %rax, 0x20(%rsp)
leaq 0x8(%rsi), %rax
movq %rax, 0x10(%rsp)
movl $0xb0, %eax
addq %rax, %rcx
movq %rcx, 0x18(%rsp)
addq %rsi, %rax
movq %rax, 0x8(%rsp)
xorl %eax, %eax
movq %rax, 0x28(%rsp)
leaq (,%rax,8), %rax
leaq (%rax,%rax,2), %r12
movq 0x20(%rsp), %rbx
leaq (%rbx,%r12), %r14
movq 0x10(%rsp), %rax
movq 0x8(%rax,%r12), %rsi
subq (%rax,%r12), %rsi
sarq $0x3, %rsi
movq %r14, %rdi
callq 0x21f42
movq 0x18(%rsp), %rax
leaq (%rax,%r12), %r15
movq 0x8(%rsp), %rax
movq 0x8(%rax,%r12), %rsi
subq (%rax,%r12), %rsi
sarq $0x3, %rsi
movq %r15, %rdi
callq 0x21f42
movq 0x8(%rbx,%r12), %rax
cmpq (%rbx,%r12), %rax
je 0x22438
movq %r12, %rcx
movq 0x10(%rsp), %rax
addq %r12, %rax
movq %rax, 0x30(%rsp)
addq 0x8(%rsp), %rcx
movq %rcx, 0x38(%rsp)
xorl %ebp, %ebp
xorl %ebx, %ebx
movq (%rsp), %rcx
movl 0x4(%rcx), %eax
imull (%rcx), %eax
xorps %xmm0, %xmm0
cvtsi2sd %eax, %xmm0
movsd %xmm0, 0x40(%rsp)
movl %ebp, %eax
xorps %xmm0, %xmm0
cvtsi2sd %rax, %xmm0
callq 0x32a0
movsd 0x40(%rsp), %xmm1
divsd %xmm0, %xmm1
cvttsd2si %xmm1, %eax
cltq
leaq (,%rax,8), %r12
movq $-0x1, %r13
testl %eax, %eax
js 0x223ca
movq %r12, %r13
movq %r13, %rdi
callq 0x3030
movq (%r14), %rcx
movq %rax, (%rcx,%rbx,8)
movq %r13, %rdi
callq 0x3030
movq (%r15), %rcx
movq %rax, (%rcx,%rbx,8)
movq (%r14), %rax
movq (%rax,%rbx,8), %rdi
movq 0x30(%rsp), %rax
movq (%rax), %rax
movq (%rax,%rbx,8), %rsi
movq %r12, %rdx
callq 0x3260
movq (%r15), %rax
movq (%rax,%rbx,8), %rdi
movq 0x38(%rsp), %rax
movq (%rax), %rax
movq (%rax,%rbx,8), %rsi
movq %r12, %rdx
callq 0x3260
movq 0x8(%r14), %rax
subq (%r14), %rax
incq %rbx
sarq $0x3, %rax
addl $0x2, %ebp
cmpq %rbx, %rax
ja 0x2237e
movq 0x28(%rsp), %rax
incq %rax
cmpq $0x7, %rax
jne 0x222fa
movq (%rsp), %rax
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/ImagePairStats.cpp
// Function: calcImgLabDeltaE(Image3ChnlReal const&, Image3ChnlReal const&)
Float_t calcImgLabDeltaE(const Image3ChnlReal &img1, const Image3ChnlReal &img2)
{
	assert(img1.width() > 0 && img1.height() > 0);
	assert(img1.width() == img2.width() && img1.height() == img2.height());
	
	Float_t deltaEmean = 0;
	int width = img1.width(), height = img1.height();
	int numNonborderPixels = 0;
	for (int y = skEdgeSkip; y < height - skEdgeSkip; y++)
		for (int x = skEdgeSkip; x < width - skEdgeSkip; x++)
		{
			deltaEmean += calcLabDeltaE((*img1[0])(x, y), (*img1[1])(x, y), (*img1[2])(x, y),
										(*img2[0])(x, y), (*img2[1])(x, y), (*img2[2])(x, y));
			numNonborderPixels++;
		}
	
	deltaEmean /= (Float_t)numNonborderPixels;
		
	return deltaEmean;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movl 0x28(%rdi), %ebp
testl %ebp, %ebp
jle 0x22b37
movq %rdi, %r14
movl 0x2c(%rdi), %r15d
testl %r15d, %r15d
jle 0x22b37
movq %rsi, %rbx
cmpl 0x28(%rsi), %ebp
jne 0x22b56
cmpl 0x2c(%rbx), %r15d
jne 0x22b56
xorpd %xmm4, %xmm4
xorl %r12d, %r12d
xorpd %xmm3, %xmm3
xorl %r13d, %r13d
movq 0x38(%r14), %rax
movq (%rax), %rcx
movq 0x8(%rax), %rdx
movq 0x38(%rcx), %rsi
movl 0x28(%rcx), %ecx
imull %r12d, %ecx
addl %r13d, %ecx
movsd (%rsi,%rcx,8), %xmm1
movq 0x38(%rdx), %rcx
movl 0x28(%rdx), %edx
imull %r12d, %edx
addl %r13d, %edx
movsd (%rcx,%rdx,8), %xmm2
movq 0x10(%rax), %rax
movq 0x38(%rax), %rcx
movl 0x28(%rax), %eax
imull %r12d, %eax
addl %r13d, %eax
movsd (%rcx,%rax,8), %xmm0
movq 0x38(%rbx), %rax
movq (%rax), %rcx
movq 0x8(%rax), %rdx
movq 0x38(%rcx), %rsi
movl 0x28(%rcx), %ecx
imull %r12d, %ecx
addl %r13d, %ecx
movq 0x38(%rdx), %rdi
movl 0x28(%rdx), %edx
imull %r12d, %edx
addl %r13d, %edx
movq 0x10(%rax), %rax
movq 0x38(%rax), %r8
movl 0x28(%rax), %eax
imull %r12d, %eax
subsd (%rsi,%rcx,8), %xmm1
addl %r13d, %eax
subsd (%rdi,%rdx,8), %xmm2
mulsd %xmm2, %xmm2
mulsd %xmm1, %xmm1
addsd %xmm2, %xmm1
subsd (%r8,%rax,8), %xmm0
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
ucomisd %xmm4, %xmm0
jb 0x22ae6
sqrtsd %xmm0, %xmm0
jmp 0x22af9
movsd %xmm3, (%rsp)
callq 0x3240
xorpd %xmm4, %xmm4
movsd (%rsp), %xmm3
addsd %xmm0, %xmm3
incl %r13d
cmpl %r13d, %ebp
jne 0x22a35
incl %r12d
cmpl %r15d, %r12d
jne 0x22a32
imull %r15d, %ebp
xorps %xmm0, %xmm0
cvtsi2sd %ebp, %xmm0
divsd %xmm0, %xmm3
movapd %xmm3, %xmm0
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0xc91c(%rip), %rdi      # 0x2f45a
leaq 0xc93b(%rip), %rsi      # 0x2f480
leaq 0xc947(%rip), %rcx      # 0x2f493
movl $0x22, %edx
callq 0x31c0
leaq 0xaf42(%rip), %rdi      # 0x2da9f
leaq 0xc91c(%rip), %rsi      # 0x2f480
leaq 0xc928(%rip), %rcx      # 0x2f493
movl $0x23, %edx
callq 0x31c0