// File: /workspace/github_cpp_makefile_star2_O2/kebwi[P]Trifid/src/SteerablePyramid.cpp
// Function: SteerablePyramid::ConvolveImage(double const*, int, int, double*, int, double*, bool)
void SteerablePyramid::ConvolveImage(const Float_t* image, int width, int height, Float_t kernel[], int kernelDim, Float_t* result, bool initValues)
{
	Float_t* resultPtr = result;
	if (initValues)
		memset(result, 0, width * height * sizeof(Float_t));
	
	assert(kernelDim % 2 == 1);	//Kernel dim should be odd
	int kernelHalfDim = kernelDim / 2;
	int widthMinusKernelHalfDim = width - kernelHalfDim, widthMinusKernelDimMinusOne = width - (kernelDim - 1);
	int heightMinusKernelHalfDim = height - kernelHalfDim;
	int y, x;
	int xx, yy;
	int ky, kx;
	
	Float_t* kernelPtr;
	const Float_t* imagePtr;
	Float_t sum;
	
	//Convolve the interior
	unsigned long t = getCurrentTime_ms();
	for (y = kernelHalfDim; y < heightMinusKernelHalfDim; y++)
	{
		resultPtr = &result[yLookup_[y] + kernelHalfDim];
		for (x = 0; x < widthMinusKernelDimMinusOne; x++)
		{
			//For each component in the kernel
			kernelPtr = kernel;
			sum = 0;
			for (ky = -kernelHalfDim; ky <= kernelHalfDim; ky++)
			{
				imagePtr = image + yLookup_[y + ky] + x;
				for (kx = 0; kx < kernelDim; kx++)
				{
					sum += *imagePtr++ * *kernelPtr++;
					
					//if (kernelDim == 17 && width == 32 && y == kernelHalfDim && x == 0)
					//	cout << setw(6) << x + kx << setw(4) << y + ky;
				}
				//if (kernelDim == 17 && width == 32 && y == kernelHalfDim && x == 0)
				//	cout << endl;
			}
			*resultPtr++ += sum;
		}
	}
	gT2 += getCurrentTime_ms() - t;
	
	//Convolve the border
	t = getCurrentTime_ms();
	
	//Horizontal stripes top and bottom
	for (y = 0; y < kernelHalfDim; y++)
	{
		resultPtr = &result[yLookup_[y] + kernelHalfDim];
		for (x = kernelHalfDim; x < widthMinusKernelHalfDim; x++, resultPtr++)
		{
			//For each component in the kernel
			kernelPtr = kernel;
			for (ky = -kernelHalfDim; ky <= kernelHalfDim; ky++)
			{
				yy = y + ky;
				
				if (yy < 0)
					yy = handleNegEdge(yy, height);
				else if (yy >= height)
					assert(false);
				
				if (yy == -1)
				{
					kernelPtr += kernelDim;
					continue;
				}
				
				for (kx = -kernelHalfDim; kx <= kernelHalfDim; kx++)
				{
					//Find the offset pixel location
					xx = x + kx;
					assert(xx >= 0 && xx < width);
					*resultPtr += image[yLookup_[yy] + xx] * *kernelPtr++;
					
					//if (kernelDim == 17 && width == 32 && y == 0 && x == kernelHalfDim)
					//	cout << setw(6) << xx << setw(4) << yy;
				}
				//if (kernelDim == 17 && width == 32 && y == 0 && x == kernelHalfDim)
				//	cout << endl;
			}
		}
	}
				
	for (y = heightMinusKernelHalfDim; y < height; y++)
	{
		resultPtr = &result[yLookup_[y] + kernelHalfDim];
		for (x = kernelHalfDim; x < widthMinusKernelHalfDim; x++, resultPtr++)
		{
			//For each component in the kernel
			kernelPtr = kernel;
			for (ky = -kernelHalfDim; ky <= kernelHalfDim; ky++)
			{
				yy = y + ky;
				
				if (yy < 0)
					assert(false);
				else if (yy >= height)
					yy = handlePosEdge(yy, height);
				
				if (yy == -1)
				{
					kernelPtr += kernelDim;
					continue;
				}
				
				for (kx = -kernelHalfDim; kx <= kernelHalfDim; kx++)
				{
					//Find the offset pixel location
					xx = x + kx;
					assert(xx >= 0 && xx < width);
					*resultPtr += image[yLookup_[yy] + xx] * *kernelPtr++;
				}
			}
		}
	}
	
	//Vertical stripes left and right
	for (y = kernelHalfDim; y < heightMinusKernelHalfDim; y++)
	{
		resultPtr = &result[y * width];
		for (x = 0; x < kernelHalfDim; x++, resultPtr++)
		{
			//For each component in the kernel
			kernelPtr = kernel;
			for (ky = -kernelHalfDim; ky <= kernelHalfDim; ky++)
			{
				yy = y + ky;
				for (kx = -kernelHalfDim; kx <= kernelHalfDim; kx++)
				{
					//Find the offset pixel location
					xx = x + kx;
					
					if (xx < 0)
						xx = handleNegEdge(xx, width);
					else if (xx >= width)
						assert(false);
					
					if (xx == -1)
					{
						kernelPtr++;
						continue;
					}
					
					*resultPtr += image[yLookup_[yy] + xx] * *kernelPtr++;
				}
			}
		}
				
		resultPtr = &result[y * width] + widthMinusKernelHalfDim;
		for (x = widthMinusKernelHalfDim; x < width; x++, resultPtr++)
		{
			//For each component in the kernel
			kernelPtr = kernel;
			for (ky = -kernelHalfDim; ky <= kernelHalfDim; ky++)
			{
				yy = y + ky;
				for (kx = -kernelHalfDim; kx <= kernelHalfDim; kx++)
				{
					//Find the offset pixel location
					xx = x + kx;
					
					if (xx < 0)
						assert(false);
					else if (xx >= width)
						xx = handlePosEdge(xx, width);
					
					if (xx == -1)
					{
						kernelPtr++;
						continue;
					}
					
					*resultPtr += image[yLookup_[yy] + xx] * *kernelPtr++;
				}
			}
		}
	}
	
	//Corners
	for (y = 0; y < kernelHalfDim; y++)
		for (x = 0; x < kernelHalfDim; x++)
		{
			//For each component in the kernel
			kernelPtr = kernel;
			for (ky = -kernelHalfDim; ky <= kernelHalfDim; ky++)
				for (kx = -kernelHalfDim; kx <= kernelHalfDim; kx++, kernelPtr++)
				{
					//Find the offset pixel location
					xx = x + kx;
					yy = y + ky;
					
					if (xx < 0)
						xx = handleNegEdge(xx, width);
					else if (xx >= width)
						assert(false);
					
					if (yy < 0)
						yy = handleNegEdge(yy, height);
					else if (yy >= height)
						assert(false);
					
					if (xx != -1 && yy != -1)
						result[y * width + x] += image[yLookup_[yy] + xx] * *kernelPtr;
					
					//Find the offset pixel location
					xx = widthMinusKernelHalfDim + x + kx;
					yy = y + ky;
					
					if (xx < 0)
						assert(false);
					else if (xx >= width)
						xx = handlePosEdge(xx, width);
					
					if (yy < 0)
						yy = handleNegEdge(yy, height);
					else if (yy >= height)
						assert(false);
					
					if (xx != -1 && yy != -1)
						result[y * width + (widthMinusKernelHalfDim + x)] += image[yLookup_[yy] + xx] * *kernelPtr;
					
					//Find the offset pixel location
					xx = x + kx;
					yy = heightMinusKernelHalfDim + y + ky;
					
					if (xx < 0)
						xx = handleNegEdge(xx, width);
					else if (xx >= width)
						assert(false);
					
					if (yy < 0)
						assert(false);
					else if (yy >= height)
						yy = handlePosEdge(yy, height);
					
					if (xx != -1 && yy != -1)
						result[(heightMinusKernelHalfDim + y) * width + x] += image[yLookup_[yy] + xx] * *kernelPtr;
					
					//Find the offset pixel location
					xx = widthMinusKernelHalfDim + x + kx;
					yy = heightMinusKernelHalfDim + y + ky;
					
					if (xx < 0)
						assert(false);
					else if (xx >= width)
						xx = handlePosEdge(xx, width);
					
					if (yy < 0)
						assert(false);
					else if (yy >= height)
						yy = handlePosEdge(yy, height);
					
					if (xx != -1 && yy != -1)
						result[(heightMinusKernelHalfDim + y) * width + (widthMinusKernelHalfDim + x)] += image[yLookup_[yy] + xx] * *kernelPtr;
				}
		}
	gT3 += getCurrentTime_ms() - t;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0xd8, %rsp
movl %r9d, %r12d
movq %r8, 0x28(%rsp)
movl %ecx, 0x8(%rsp)
movl %edx, %r15d
movq %rsi, 0x20(%rsp)
movq %rdi, 0x30(%rsp)
cmpb $0x0, 0x118(%rsp)
je 0x168d4
movl 0x8(%rsp), %eax
imull %r15d, %eax
movslq %eax, %rdx
shlq $0x3, %rdx
movq 0x110(%rsp), %rdi
xorl %esi, %esi
callq 0x3170
movl %r12d, %eax
andl $0x80000001, %eax       # imm = 0x80000001
cmpl $0x1, %eax
jne 0x17274
movl %r12d, %r14d
shrl %r14d
movl %r15d, %ebx
subl %r12d, %ebx
incl %ebx
movl 0x8(%rsp), %ebp
subl %r14d, %ebp
callq 0x3c71
movq %rax, %r13
movq 0x30(%rsp), %rax
movq 0x40(%rax), %rax
movl %r14d, %ecx
negl %ecx
movl %ecx, 0x4(%rsp)
movslq %ecx, %rcx
movq %rcx, 0x48(%rsp)
xorl %ecx, %ecx
testl %ebx, %ebx
cmovlel %ecx, %ebx
leal 0x1(%r14), %ecx
movq %rcx, 0x50(%rsp)
movl %r12d, %r12d
cmpl %ebp, %r14d
movl %ebp, 0xc(%rsp)
movl %ebp, %ecx
cmovgl %r14d, %ecx
movq %rcx, 0x60(%rsp)
movl %r14d, %edi
movq %r14, %rcx
movq %r14, 0x10(%rsp)
movq %rdi, 0x58(%rsp)
cmpq 0x60(%rsp), %rcx
je 0x169e4
leaq (%rax,%rcx,4), %rdx
movslq (%rdx), %rsi
addq %rdi, %rsi
movq 0x110(%rsp), %rdi
leaq (%rdi,%rsi,8), %r8
movq 0x20(%rsp), %rsi
xorl %edi, %edi
cmpq %rbx, %rdi
je 0x169d7
xorpd %xmm0, %xmm0
movq 0x48(%rsp), %r9
movq 0x28(%rsp), %r10
cmpq 0x50(%rsp), %r9
je 0x169c0
movslq (%rdx,%r9,4), %r11
leaq (%rsi,%r11,8), %r11
xorl %r14d, %r14d
cmpl %r14d, %r12d
je 0x169bb
movsd (%r11,%r14,8), %xmm1
mulsd (%r10), %xmm1
addq $0x8, %r10
addsd %xmm1, %xmm0
incq %r14
jmp 0x1699e
incq %r9
jmp 0x1698c
addsd (%r8), %xmm0
movsd %xmm0, (%r8)
addq $0x8, %r8
incq %rdi
addq $0x8, %rsi
jmp 0x16979
incq %rcx
movq 0x58(%rsp), %rdi
jmp 0x16951
movl %r15d, %ebp
movq 0x10(%rsp), %rbx
subl %ebx, %ebp
callq 0x3c71
subq %r13, %rax
addq %rax, 0x1c193(%rip)     # 0x32b90
callq 0x3c71
movq %rax, 0x98(%rsp)
movq 0x30(%rsp), %rax
movq 0x40(%rax), %rcx
cmpl %ebp, %ebx
movl %ebp, 0x38(%rsp)
cmovgl %ebx, %ebp
movl %ebp, 0x40(%rsp)
leal 0x1(,%rbx,2), %edx
xorl %esi, %esi
movq 0x58(%rsp), %rdi
cmpq %rdi, %rsi
je 0x16af6
movslq (%rcx,%rsi,4), %rax
addq %rdi, %rax
movq 0x110(%rsp), %rdi
leaq (%rdi,%rax,8), %rdi
xorl %r8d, %r8d
movq 0x10(%rsp), %rax
movl %eax, %r9d
cmpl 0x40(%rsp), %r9d
je 0x16aee
movl %r8d, 0x18(%rsp)
movslq %r8d, %r10
movl 0x4(%rsp), %r11d
movq 0x28(%rsp), %rbx
cmpl 0x50(%rsp), %r11d
je 0x16ada
movl %r11d, %eax
addl %esi, %eax
js 0x16a8b
cmpl 0x8(%rsp), %eax
jl 0x16a94
jmp 0x17255
addl 0x8(%rsp), %eax
cmpl $-0x1, %eax
je 0x16ad1
movslq %eax, %r14
movq %r10, %r13
movl %edx, %eax
subl $0x1, %eax
jb 0x16ad5
cmpl %r15d, %r13d
jge 0x16bf7
movslq (%rcx,%r14,4), %rbp
addq %r13, %rbp
movq 0x20(%rsp), %r8
movsd (%r8,%rbp,8), %xmm0
mulsd (%rbx), %xmm0
addq $0x8, %rbx
addsd (%rdi), %xmm0
movsd %xmm0, (%rdi)
incq %r13
jmp 0x16a9c
leaq (%rbx,%r12,8), %rbx
incl %r11d
jmp 0x16a72
incl %r9d
addq $0x8, %rdi
movl 0x18(%rsp), %r8d
incl %r8d
jmp 0x16a55
incq %rsi
jmp 0x16a29
movq 0x30(%rsp), %rax
movq 0x40(%rax), %rcx
movslq 0xc(%rsp), %r8
movslq 0x8(%rsp), %rax
movq %rax, 0x68(%rsp)
movq 0x10(%rsp), %rax
leal 0x1(,%rax,2), %r10d
xorl %edi, %edi
movq %r8, 0xb8(%rsp)
movq 0x20(%rsp), %r11
cmpq 0x68(%rsp), %r8
jge 0x16c16
movslq (%rcx,%r8,4), %rax
addq 0x58(%rsp), %rax
movq 0x110(%rsp), %rdx
leaq (%rdx,%rax,8), %r9
xorl %edx, %edx
movq 0x10(%rsp), %rax
cmpl 0x40(%rsp), %eax
je 0x16bef
movl %eax, 0xc(%rsp)
movl %edx, 0x18(%rsp)
movslq %edx, %rax
movl 0x4(%rsp), %r14d
movq 0x28(%rsp), %r13
cmpl 0x50(%rsp), %r14d
je 0x16bda
movl %r14d, %edx
addl %r8d, %edx
js 0x17236
movl 0x8(%rsp), %esi
cmpl %esi, %edx
cmovll %edi, %esi
subl %esi, %edx
cmpl $-0x1, %edx
je 0x16bd0
movslq %edx, %rdx
movq %rax, %rbx
movl %r10d, %esi
subl $0x1, %esi
jb 0x16bd5
cmpl %r15d, %ebx
jge 0x16ddb
movslq (%rcx,%rdx,4), %rbp
addq %rbx, %rbp
movsd (%r11,%rbp,8), %xmm0
mulsd (%r13), %xmm0
addq $0x8, %r13
addsd (%r9), %xmm0
movsd %xmm0, (%r9)
incq %rbx
jmp 0x16b9c
leaq (%r13,%r12,8), %r13
incl %r14d
jmp 0x16b70
movl 0xc(%rsp), %eax
incl %eax
addq $0x8, %r9
movl 0x18(%rsp), %edx
incl %edx
jmp 0x16b51
incq %r8
jmp 0x16b2a
leaq 0xea16(%rip), %rdi      # 0x25614
leaq 0xe6bd(%rip), %rsi      # 0x252c2
leaq 0xe9a7(%rip), %rcx      # 0x255b3
movl $0x3cb, %edx            # imm = 0x3CB
callq 0x31a0
movslq 0x38(%rsp), %rax
movq %rax, 0x88(%rsp)
movslq %r15d, %rax
movq %rax, 0x80(%rsp)
movq 0x10(%rsp), %r8
leal (%r8,%r8), %eax
leal 0x1(%r8,%r8), %edi
movl %r15d, %edx
notl %edx
movl %eax, 0xc(%rsp)
negl %eax
movl %eax, 0x18(%rsp)
cmpq 0x60(%rsp), %r8
je 0x16dfa
movq %r8, %rax
imulq 0x80(%rsp), %rax
movq 0x110(%rsp), %rcx
leaq (%rcx,%rax,8), %rbx
xorl %r10d, %r10d
movl 0x4(%rsp), %eax
movq %rbx, 0x40(%rsp)
cmpl 0x10(%rsp), %r10d
je 0x16d18
movq 0x48(%rsp), %r14
movq 0x28(%rsp), %r12
cmpq 0x50(%rsp), %r14
je 0x16d0a
leaq (,%r14,4), %r13
movq 0x30(%rsp), %rcx
addq 0x40(%rcx), %r13
movl %eax, %r11d
movl %edi, %ecx
movl 0x4(%rsp), %esi
movl %esi, %r9d
subl $0x1, %ecx
jb 0x16d05
testl %r11d, %r11d
js 0x16cca
movl %r11d, %ebp
cmpl %r15d, %r11d
jl 0x16cd6
jmp 0x17156
cmpl %r11d, %edx
je 0x16cf9
leal (%r9,%r10), %ebp
addl %r15d, %ebp
movslq (%r13,%r8,4), %rsi
movslq %ebp, %rbp
addq %rsi, %rbp
movq 0x20(%rsp), %rsi
movsd (%rsi,%rbp,8), %xmm0
mulsd (%r12), %xmm0
addsd (%rbx), %xmm0
movsd %xmm0, (%rbx)
addq $0x8, %r12
incl %r9d
incl %r11d
jmp 0x16cb3
incq %r14
jmp 0x16c8f
incl %r10d
addq $0x8, %rbx
incl %eax
jmp 0x16c7a
movq 0x88(%rsp), %rax
movq 0x40(%rsp), %rcx
leaq (%rcx,%rax,8), %r9
pushq $-0x1
popq %r10
movl %r15d, %ebp
movl 0x38(%rsp), %eax
movl %eax, %ebx
cmpl %r15d, %ebx
jge 0x16dd3
movq 0x48(%rsp), %r14
movq 0x28(%rsp), %r12
cmpq 0x50(%rsp), %r14
je 0x16dc3
leaq (,%r14,4), %rax
movq 0x30(%rsp), %rcx
addq 0x40(%rcx), %rax
movl 0x18(%rsp), %ecx
movl 0x4(%rsp), %esi
movl %esi, %r13d
cmpl $0x1, %ecx
je 0x16dbe
movl %ebp, %r11d
addl %ecx, %r11d
js 0x17137
cmpl %r15d, %r11d
jl 0x16d8e
cmpl %ecx, %r10d
je 0x16db3
leal (%rbx,%r13), %r11d
subl %r15d, %r11d
movslq (%rax,%r8,4), %rsi
movslq %r11d, %r11
addq %rsi, %r11
movq 0x20(%rsp), %rsi
movsd (%rsi,%r11,8), %xmm0
mulsd (%r12), %xmm0
addsd (%r9), %xmm0
movsd %xmm0, (%r9)
addq $0x8, %r12
incl %r13d
incl %ecx
jmp 0x16d6c
incq %r14
jmp 0x16d49
incl %ebx
addq $0x8, %r9
incl %ebp
decl %r10d
jmp 0x16d36
incq %r8
jmp 0x16c4b
leaq 0xe832(%rip), %rdi      # 0x25614
leaq 0xe4d9(%rip), %rsi      # 0x252c2
leaq 0xe7c3(%rip), %rcx      # 0x255b3
movl $0x3f1, %edx            # imm = 0x3F1
callq 0x31a0
movq 0x10(%rsp), %rax
leaq 0x1(,%rax,2), %rbp
movl 0x38(%rsp), %eax
movq %rax, 0xb0(%rsp)
movslq 0x4(%rsp), %rax
movq %rax, 0xa8(%rsp)
movl %r15d, %eax
subl 0xc(%rsp), %eax
movq %rax, 0xa0(%rsp)
xorl %r8d, %r8d
movq 0x58(%rsp), %rdi
cmpq %rdi, %r8
je 0x17175
movq %r8, %rax
movq 0x80(%rsp), %rdx
imulq %rdx, %rax
movq 0xb8(%rsp), %rcx
addq %r8, %rcx
movq %rcx, 0xc8(%rsp)
imulq %rdx, %rcx
movq 0x110(%rsp), %rdx
leaq (%rdx,%rax,8), %rax
movq %rax, 0x38(%rsp)
leaq (%rdx,%rcx,8), %rax
movq %rax, 0x90(%rsp)
movq 0xa0(%rsp), %r9
movq 0xa8(%rsp), %rax
movq %rax, 0x60(%rsp)
movq 0xb0(%rsp), %rax
movq %rax, 0x70(%rsp)
xorl %ecx, %ecx
movq %r8, 0xc0(%rsp)
cmpq %rdi, %rcx
je 0x1712f
movq 0x88(%rsp), %rax
addq %rcx, %rax
movq %rax, 0x68(%rsp)
movl 0x4(%rsp), %eax
movl %eax, %r11d
movq 0x28(%rsp), %rsi
cmpl 0x50(%rsp), %r11d
je 0x1711a
movl %r11d, %eax
addl %r8d, %eax
sets %dl
movl 0x8(%rsp), %r10d
cmpl %r10d, %eax
setl %bl
movq 0x30(%rsp), %rdi
movq 0x40(%rdi), %rdi
movq 0xc8(%rsp), %r8
movq %r11, 0xd0(%rsp)
leal (%r11,%r8), %r14d
cmpl %r10d, %r14d
movl %r10d, %r8d
movl $0x0, %r11d
cmovll %r11d, %r8d
movl %r14d, 0xc(%rsp)
subl %r8d, %r14d
movl %r14d, 0x48(%rsp)
movslq %r14d, %r8
movq %r8, 0x10(%rsp)
orb %dl, %bl
movb %bl, 0x40(%rsp)
movl %eax, %edx
sarl $0x1f, %edx
andl %r10d, %edx
addl %eax, %edx
movl %edx, 0x7c(%rsp)
movslq %edx, %rax
movq %rax, 0x18(%rsp)
xorl %edx, %edx
movl 0x4(%rsp), %eax
movl %eax, %r12d
movq %rsi, %r13
cmpl %edx, %ebp
je 0x170fa
leal (%r12,%rcx), %ebx
movq 0x60(%rsp), %rax
addq %rdx, %rax
testl %eax, %eax
js 0x16f82
movl %ebx, %r10d
cmpl %r15d, %eax
jl 0x16f8b
jmp 0x17217
movq 0x70(%rsp), %r8
leal (%r8,%rdx), %r10d
cmpb $0x0, 0x40(%rsp)
je 0x171ba
cmpl $-0x1, %r10d
je 0x16fd3
cmpl $-0x1, 0x7c(%rsp)
je 0x16fd3
movq 0x18(%rsp), %r8
movslq (%rdi,%r8,4), %r14
movslq %r10d, %r10
addq %r14, %r10
movq 0x20(%rsp), %r11
movsd (%r11,%r10,8), %xmm0
mulsd (%rsi,%rdx,8), %xmm0
movq 0x38(%rsp), %r8
addsd (%r8,%rcx,8), %xmm0
movsd %xmm0, (%r8,%rcx,8)
leaq (%r9,%rdx), %r10
testl %r10d, %r10d
js 0x1719b
cmpl %r15d, %r10d
movl %r15d, %r14d
movl $0x0, %r8d
cmovll %r8d, %r14d
leal (%r9,%rdx), %r10d
subl %r14d, %r10d
cmpl $-0x1, %r10d
je 0x17045
cmpl $-0x1, 0x18(%rsp)
je 0x17045
movq %rbp, %r8
movq 0x18(%rsp), %r11
movl (%rdi,%r11,4), %ebp
subl %r14d, %ebp
addl %r9d, %ebp
addl %edx, %ebp
movslq %ebp, %r11
movq 0x20(%rsp), %rbp
movsd (%rbp,%r11,8), %xmm0
movq %r8, %rbp
mulsd (%rsi,%rdx,8), %xmm0
movq 0x38(%rsp), %r8
movq 0x68(%rsp), %r11
addsd (%r8,%r11,8), %xmm0
movsd %xmm0, (%r8,%r11,8)
testl %eax, %eax
js 0x17053
cmpl %r15d, %eax
jl 0x1705b
jmp 0x171f8
movq 0x70(%rsp), %rax
leal (%rax,%rdx), %ebx
cmpl $0x0, 0xc(%rsp)
js 0x171d9
cmpl $-0x1, %ebx
je 0x170a3
cmpl $-0x1, 0x48(%rsp)
je 0x170a3
movq 0x10(%rsp), %rax
movslq (%rdi,%rax,4), %rax
movslq %ebx, %r11
addq %rax, %r11
movq 0x20(%rsp), %rax
movsd (%rax,%r11,8), %xmm0
mulsd (%rsi,%rdx,8), %xmm0
movq 0x90(%rsp), %rax
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
cmpl $-0x1, %r10d
je 0x170eb
cmpl $-0x1, 0x48(%rsp)
je 0x170eb
movq 0x10(%rsp), %rax
movl (%rdi,%rax,4), %eax
subl %r14d, %eax
addl %r9d, %eax
addl %edx, %eax
cltq
movq 0x20(%rsp), %r10
movsd (%r10,%rax,8), %xmm0
mulsd (%rsi,%rdx,8), %xmm0
movq 0x90(%rsp), %rax
movq 0x68(%rsp), %r8
addsd (%rax,%r8,8), %xmm0
movsd %xmm0, (%rax,%r8,8)
incl %r12d
addq $0x8, %r13
incq %rdx
jmp 0x16f5d
movq 0xd0(%rsp), %r11
incl %r11d
movq %r13, %rsi
movq 0x58(%rsp), %rdi
movq 0xc0(%rsp), %r8
jmp 0x16ed4
incq %rcx
incq 0x70(%rsp)
incq 0x60(%rsp)
incq %r9
jmp 0x16eaf
incq %r8
jmp 0x16e37
leaq 0xe400(%rip), %rdi      # 0x2553e
leaq 0xe17d(%rip), %rsi      # 0x252c2
leaq 0xe467(%rip), %rcx      # 0x255b3
movl $0x426, %edx            # imm = 0x426
callq 0x31a0
leaq 0xe3e1(%rip), %rdi      # 0x2553e
leaq 0xe15e(%rip), %rsi      # 0x252c2
leaq 0xe448(%rip), %rcx      # 0x255b3
movl $0x40b, %edx            # imm = 0x40B
callq 0x31a0
callq 0x3c71
subq 0x98(%rsp), %rax
addq %rax, 0x1ba0f(%rip)     # 0x32b98
addq $0xd8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0xe39c(%rip), %rdi      # 0x2553e
leaq 0xe119(%rip), %rsi      # 0x252c2
leaq 0xe403(%rip), %rcx      # 0x255b3
movl $0x455, %edx            # imm = 0x455
callq 0x31a0
leaq 0xe37d(%rip), %rdi      # 0x2553e
leaq 0xe0fa(%rip), %rsi      # 0x252c2
leaq 0xe3e4(%rip), %rcx      # 0x255b3
movl $0x44b, %edx            # imm = 0x44B
callq 0x31a0
leaq 0xe35e(%rip), %rdi      # 0x2553e
leaq 0xe0db(%rip), %rsi      # 0x252c2
leaq 0xe3c5(%rip), %rcx      # 0x255b3
movl $0x46b, %edx            # imm = 0x46B
callq 0x31a0
leaq 0xe33f(%rip), %rdi      # 0x2553e
leaq 0xe0bc(%rip), %rsi      # 0x252c2
leaq 0xe3a6(%rip), %rcx      # 0x255b3
movl $0x468, %edx            # imm = 0x468
callq 0x31a0
leaq 0xe320(%rip), %rdi      # 0x2553e
leaq 0xe09d(%rip), %rsi      # 0x252c2
leaq 0xe387(%rip), %rcx      # 0x255b3
movl $0x446, %edx            # imm = 0x446
callq 0x31a0
leaq 0xe301(%rip), %rdi      # 0x2553e
leaq 0xe07e(%rip), %rsi      # 0x252c2
leaq 0xe368(%rip), %rcx      # 0x255b3
movl $0x3e3, %edx            # imm = 0x3E3
callq 0x31a0
leaq 0xe2e2(%rip), %rdi      # 0x2553e
leaq 0xe05f(%rip), %rsi      # 0x252c2
leaq 0xe349(%rip), %rcx      # 0x255b3
movl $0x3bf, %edx            # imm = 0x3BF
callq 0x31a0
leaq 0xe325(%rip), %rdi      # 0x255a0
leaq 0xe040(%rip), %rsi      # 0x252c2
leaq 0xe32a(%rip), %rcx      # 0x255b3
movl $0x385, %edx            # imm = 0x385
callq 0x31a0
nop