// File: /workspace/github_cpp_makefile_star2_O1/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 $0xe8, %rsp
movl %r9d, 0x28(%rsp)
movq %r8, 0x58(%rsp)
movq %rcx, 0x30(%rsp)
movl %edx, %r15d
movq %rsi, %r12
movq %rdi, 0x60(%rsp)
cmpb $0x0, 0x128(%rsp)
je 0x1c9c0
movq 0x30(%rsp), %rax
imull %r15d, %eax
movslq %eax, %rdx
shlq $0x3, %rdx
movq 0x120(%rsp), %rdi
xorl %esi, %esi
callq 0x3190
movl 0x28(%rsp), %ecx
movl %ecx, %eax
andl $0x80000001, %eax       # imm = 0x80000001
cmpl $0x1, %eax
jne 0x1d492
movl %ecx, %ebx
shrl %ebx
movq 0x30(%rsp), %rax
movl %eax, %ebp
subl %ebx, %ebp
callq 0x4288
movl 0x28(%rsp), %edi
movq %rax, 0x40(%rsp)
movl %edi, %r13d
movl %ebx, %eax
movq %rax, 0x50(%rsp)
movq %rbx, 0x20(%rsp)
movl %ebp, 0x70(%rsp)
cmpl %ebp, %ebx
jge 0x1caff
movl %r15d, %r8d
subl %edi, %r8d
incl %r8d
movq 0x60(%rsp), %rax
movq 0x40(%rax), %rax
movq %rax, 0x38(%rsp)
movq 0x20(%rsp), %rdx
movl %edx, %eax
negl %eax
movslq %eax, %rcx
leal 0x1(%rdx), %esi
movq 0x30(%rsp), %rax
leal (%rdx,%rax), %r9d
movl %edi, %eax
andl $0x7ffffffe, %eax       # imm = 0x7FFFFFFE
subl %eax, %r9d
movl %r9d, 0xc(%rsp)
movl %r8d, 0x18(%rsp)
movl %r8d, %r8d
movq %rdx, 0x10(%rsp)
cmpl $0x0, 0x18(%rsp)
jle 0x1cae8
movq 0x38(%rsp), %rax
movq 0x10(%rsp), %rdx
leaq (%rax,%rdx,4), %r10
movslq (%r10), %rax
addq 0x50(%rsp), %rax
movq 0x120(%rsp), %rdx
leaq (%rdx,%rax,8), %r14
movq %r12, %r11
xorl %ebx, %ebx
xorpd %xmm0, %xmm0
movq %rcx, %rdx
movq 0x58(%rsp), %rax
testl %edi, %edi
jle 0x1cac7
movslq (%r10,%rdx,4), %rdi
leaq (%r11,%rdi,8), %rdi
xorl %ebp, %ebp
xorl %r9d, %r9d
movsd (%rdi,%r9,8), %xmm1
mulsd (%rax,%r9,8), %xmm1
addsd %xmm1, %xmm0
incq %r9
addq $-0x8, %rbp
cmpl %r9d, %r13d
jne 0x1caa4
subq %rbp, %rax
movl 0x28(%rsp), %edi
incq %rdx
cmpl %edx, %esi
jne 0x1ca93
addsd (%r14), %xmm0
movsd %xmm0, (%r14)
addq $0x8, %r14
incq %rbx
addq $0x8, %r11
cmpq %r8, %rbx
jne 0x1ca87
movq 0x10(%rsp), %rdx
incq %rdx
movq %rdx, 0x10(%rsp)
cmpl %edx, 0xc(%rsp)
jne 0x1ca55
movl %r15d, %eax
subl 0x20(%rsp), %eax
movq %rax, 0x68(%rsp)
movl %edi, %ebx
callq 0x4288
subq 0x40(%rsp), %rax
addq %rax, 0x1d082(%rip)     # 0x39ba0
callq 0x4288
movq %rax, 0xa8(%rsp)
cmpl $0x2, %ebx
jb 0x1cc2a
movq 0x60(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x20(%rsp), %rdx
movl %edx, %ecx
negl %ecx
movl %ecx, 0x10(%rsp)
movl 0x28(%rsp), %ecx
andl $0x7ffffffe, %ecx       # imm = 0x7FFFFFFE
movl %r15d, %esi
subl %ecx, %esi
movl %esi, 0x18(%rsp)
leal 0x1(,%rdx,2), %ecx
xorl %edi, %edi
movq 0x68(%rsp), %rdx
cmpl %edx, 0x20(%rsp)
jge 0x1cc1c
movslq (%rax,%rdi,4), %rsi
movq 0x50(%rsp), %rdx
addq %rdx, %rsi
movq 0x120(%rsp), %rdx
leaq (%rdx,%rsi,8), %r8
xorl %r9d, %r9d
movslq %r9d, %r10
movq 0x58(%rsp), %r11
movl 0x10(%rsp), %edx
movl %edx, %ebx
movl %ebx, %esi
addl %edi, %esi
js 0x1cbaf
cmpl 0x30(%rsp), %esi
jl 0x1cbb3
jmp 0x1d473
addl 0x30(%rsp), %esi
cmpl $-0x1, %esi
je 0x1cbfb
movslq %esi, %r14
movq %r10, %rbp
movl %ecx, %esi
testl %r9d, %r9d
js 0x1d31e
cmpl %r15d, %ebp
jge 0x1d31e
movslq (%rax,%r14,4), %rdx
addq %rbp, %rdx
movsd (%r12,%rdx,8), %xmm0
mulsd (%r11), %xmm0
addsd (%r8), %xmm0
addq $0x8, %r11
movsd %xmm0, (%r8)
incq %rbp
decl %esi
jne 0x1cbc0
jmp 0x1cbff
leaq (%r11,%r13,8), %r11
leal 0x1(%rbx), %edx
cmpl 0x20(%rsp), %ebx
movl %edx, %ebx
jne 0x1cb9e
addq $0x8, %r8
incl %r9d
cmpl 0x18(%rsp), %r9d
jne 0x1cb90
incq %rdi
cmpq 0x50(%rsp), %rdi
jne 0x1cb66
cmpl $0x2, 0x28(%rsp)
jb 0x1cd58
movq 0x60(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x20(%rsp), %r8
movl %r8d, %ecx
negl %ecx
movl %ecx, 0x18(%rsp)
movl 0x28(%rsp), %ecx
andl $0x7ffffffe, %ecx       # imm = 0x7FFFFFFE
movq 0x30(%rsp), %rsi
movl %esi, %edx
subl %ecx, %edx
movl %r15d, %edi
subl %ecx, %edi
movl %edi, 0xc(%rsp)
movslq 0x70(%rsp), %rdi
movslq %esi, %rcx
movq %rcx, 0x38(%rsp)
leal 0x1(,%r8,2), %ecx
movl %ecx, 0x10(%rsp)
movq 0x68(%rsp), %rcx
cmpl %ecx, 0x20(%rsp)
jge 0x1cd48
movslq (%rax,%rdi,4), %rcx
addq 0x50(%rsp), %rcx
movq 0x120(%rsp), %rsi
leaq (%rsi,%rcx,8), %r11
xorl %ebx, %ebx
movslq %ebx, %rsi
movq 0x58(%rsp), %r8
movl 0x18(%rsp), %ecx
testl %edx, %edx
js 0x1d454
leal (%rcx,%rdi), %r9d
movq 0x30(%rsp), %r10
cmpl %r10d, %r9d
movl %r10d, %ebp
movl $0x0, %r10d
cmovll %r10d, %ebp
subl %ebp, %r9d
cmpl $-0x1, %r9d
je 0x1cd23
movslq %r9d, %rbp
movq %rsi, %r14
movl 0x10(%rsp), %r9d
testl %ebx, %ebx
js 0x1d35c
cmpl %r15d, %r14d
jge 0x1d35c
movslq (%rax,%rbp,4), %r10
addq %r14, %r10
movsd (%r12,%r10,8), %xmm0
mulsd (%r8), %xmm0
addsd (%r11), %xmm0
addq $0x8, %r8
movsd %xmm0, (%r11)
incq %r14
decl %r9d
jne 0x1cce8
jmp 0x1cd27
leaq (%r8,%r13,8), %r8
leal 0x1(%rcx), %r9d
cmpl 0x20(%rsp), %ecx
movl %r9d, %ecx
jne 0x1ccb3
addq $0x8, %r11
incl %ebx
cmpl 0xc(%rsp), %ebx
jne 0x1cca7
incq %rdi
incl %edx
cmpq 0x38(%rsp), %rdi
jl 0x1cc81
movl 0x70(%rsp), %eax
cmpl %eax, 0x20(%rsp)
movl 0x28(%rsp), %edx
jge 0x1cf82
movq 0x20(%rsp), %r14
movl %r14d, %eax
negl %eax
movslq 0x68(%rsp), %rcx
movq %rcx, 0x98(%rsp)
leal 0x1(%r14), %ecx
movl %ecx, 0x10(%rsp)
movl %eax, 0xc(%rsp)
cltq
movq %rax, 0x40(%rsp)
movl %edx, %eax
andl $0x7ffffffe, %eax       # imm = 0x7FFFFFFE
movl %r15d, %ecx
subl %eax, %ecx
movl %ecx, 0x4c(%rsp)
movslq %r15d, %rcx
movq %rcx, 0x90(%rsp)
movq 0x30(%rsp), %rcx
addl %r14d, %ecx
subl %eax, %ecx
movl %ecx, 0xa0(%rsp)
leal 0x1(%r14,%r14), %ecx
movl %ecx, 0x38(%rsp)
orq $0x1, %r13
movl %r15d, %ecx
subl %r14d, %ecx
subl %r14d, %ecx
movq %rcx, 0x78(%rsp)
xorl %ebx, %ebx
movq %r14, %rax
imulq 0x90(%rsp), %rax
movq 0x120(%rsp), %rcx
leaq (%rcx,%rax,8), %rax
movq %rax, 0x80(%rsp)
cmpl $0x2, %edx
jb 0x1ceab
xorl %r11d, %r11d
movl 0xc(%rsp), %eax
movl %eax, 0x18(%rsp)
movq 0x80(%rsp), %r8
movq 0x40(%rsp), %rsi
movq 0x58(%rsp), %rcx
leaq (,%rsi,4), %r9
movq 0x60(%rsp), %rax
addq 0x40(%rax), %r9
movl 0x18(%rsp), %eax
movl %eax, %edi
movl 0x38(%rsp), %r10d
movl 0xc(%rsp), %eax
movl %eax, %ebp
testl %edi, %edi
js 0x1ce58
cmpl %r15d, %edi
jge 0x1d37b
leal (%r11,%rbp), %edx
jmp 0x1ce5c
leal (%r15,%rdi), %edx
cmpl $-0x1, %edx
je 0x1ce7f
movslq (%r9,%r14,4), %rax
movslq %edx, %rdx
addq %rax, %rdx
movsd (%r12,%rdx,8), %xmm0
mulsd (%rcx), %xmm0
addsd (%r8), %xmm0
movsd %xmm0, (%r8)
addq $0x8, %rcx
incl %ebp
incl %edi
decl %r10d
jne 0x1ce45
incq %rsi
cmpl %esi, 0x10(%rsp)
jne 0x1ce23
incl %r11d
addq $0x8, %r8
incl 0x18(%rsp)
cmpl 0x20(%rsp), %r11d
jne 0x1ce19
cmpl $0x2, 0x28(%rsp)
jb 0x1cf6d
movq 0x98(%rsp), %rax
movq 0x80(%rsp), %rcx
leaq (%rcx,%rax,8), %rcx
movq 0x78(%rsp), %r11
movl 0x4c(%rsp), %eax
movq 0x68(%rsp), %rdx
movl %edx, 0x18(%rsp)
movq 0x40(%rsp), %rsi
movq 0x58(%rsp), %r9
leaq (,%rsi,4), %r10
movq 0x60(%rsp), %rdx
addq 0x40(%rdx), %r10
xorl %edi, %edi
xorl %ebp, %ebp
testl %eax, %eax
js 0x1d33d
leal (%r11,%rbp), %r8d
cmpl %r15d, %r8d
movl %r15d, %edx
cmovll %ebx, %edx
subl %edx, %r8d
cmpl $-0x1, %r8d
je 0x1cf3d
movl (%r10,%r14,4), %r8d
subl %edx, %r8d
addl %r11d, %r8d
addl %ebp, %r8d
movslq %r8d, %rdx
movsd (%r12,%rdx,8), %xmm0
mulsd (%r9,%rbp,8), %xmm0
addsd (%rcx), %xmm0
movsd %xmm0, (%rcx)
incq %rbp
addq $-0x8, %rdi
cmpl %ebp, %r13d
jne 0x1cefb
incq %rsi
subq %rdi, %r9
cmpl %esi, 0x10(%rsp)
jne 0x1cee6
movl 0x18(%rsp), %edx
incl %edx
addq $0x8, %rcx
incl %eax
incq %r11
cmpl %r15d, %edx
jl 0x1ced8
incq %r14
cmpl %r14d, 0xa0(%rsp)
movl 0x28(%rsp), %edx
jne 0x1cddd
cmpl $0x2, %edx
jae 0x1cfad
callq 0x4288
subq 0xa8(%rsp), %rax
addq %rax, 0x1cc0d(%rip)     # 0x39ba8
addq $0xe8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq 0x20(%rsp), %rcx
movl %ecx, %edx
negl %edx
movslq %r15d, %rax
movq %rax, 0xd0(%rsp)
movslq 0x70(%rsp), %rax
movq %rax, 0xc8(%rsp)
leal (%rcx,%rcx), %eax
leaq 0x1(%rax), %rsi
movq %rsi, 0xa0(%rsp)
movl %r15d, %esi
subl %eax, %esi
movq %rsi, 0xc0(%rsp)
movl %r15d, %eax
subl %ecx, %eax
movq %rax, 0xb8(%rsp)
movl %edx, 0x8c(%rsp)
movslq %edx, %rax
movq %rax, 0xb0(%rsp)
xorl %esi, %esi
movq %rsi, %rax
movq 0xd0(%rsp), %rdx
imulq %rdx, %rax
movq 0xc8(%rsp), %rcx
addq %rsi, %rcx
movq %rcx, 0xe0(%rsp)
imulq %rdx, %rcx
movq 0x120(%rsp), %rdx
leaq (%rdx,%rax,8), %rax
movq %rax, 0x80(%rsp)
leaq (%rdx,%rcx,8), %rax
movq %rax, 0x90(%rsp)
movq 0xb0(%rsp), %rax
movq %rax, 0x38(%rsp)
movq 0xb8(%rsp), %rax
movq %rax, 0x78(%rsp)
movq 0xc0(%rsp), %rdi
xorl %ecx, %ecx
movq %rsi, 0xd8(%rsp)
movq 0x68(%rsp), %rax
addl %ecx, %eax
cltq
movq %rax, 0x40(%rsp)
movq 0x58(%rsp), %r11
movl 0x8c(%rsp), %eax
movl %eax, %r9d
movl %r9d, %eax
addl %esi, %eax
sets %dl
movq 0x30(%rsp), %r8
cmpl %r8d, %eax
setl %bl
movq 0x60(%rsp), %rsi
movq 0x40(%rsi), %rsi
movq %rsi, 0x10(%rsp)
movq 0xe0(%rsp), %rsi
movq %r9, 0x70(%rsp)
leal (%r9,%rsi), %r10d
cmpl %r8d, %r10d
movl %r8d, %esi
movl $0x0, %r9d
cmovll %r9d, %esi
movl %r10d, 0xc(%rsp)
subl %esi, %r10d
movl %r10d, 0x4c(%rsp)
movslq %r10d, %rsi
movq %rsi, 0x18(%rsp)
orb %dl, %bl
movl %eax, %edx
sarl $0x1f, %edx
andl %r8d, %edx
addl %eax, %edx
movl %edx, 0x98(%rsp)
movslq %edx, %rax
movq %rax, 0x28(%rsp)
xorl %r9d, %r9d
xorl %r10d, %r10d
movl 0x8c(%rsp), %eax
movl %eax, %r13d
leal (%rcx,%r13), %ebp
movq 0x38(%rsp), %rax
leaq (%rax,%r10), %rdx
testl %edx, %edx
js 0x1d139
movl %ebp, %eax
cmpl %r15d, %edx
jl 0x1d141
jmp 0x1d416
movq 0x78(%rsp), %rax
addl %r10d, %eax
testb %bl, %bl
je 0x1d3f7
cmpl $-0x1, %eax
je 0x1d186
cmpl $-0x1, 0x28(%rsp)
je 0x1d186
movq 0x10(%rsp), %rsi
movq 0x28(%rsp), %r8
movslq (%rsi,%r8,4), %rsi
cltq
addq %rsi, %rax
movsd (%r12,%rax,8), %xmm0
mulsd (%r11,%r10,8), %xmm0
movq 0x80(%rsp), %rax
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
leaq (%rdi,%r10), %r14
testl %r14d, %r14d
js 0x1d3b9
cmpl %r15d, %r14d
movl %r15d, %esi
movl $0x0, %eax
cmovll %eax, %esi
testb %bl, %bl
je 0x1d3d8
movl %esi, %r8d
negl %r8d
movq %r14, %rax
addq %r8, %rax
cmpl $-0x1, %eax
je 0x1d1fe
cmpl $-0x1, 0x98(%rsp)
je 0x1d1fe
movq 0x10(%rsp), %rax
movq 0x28(%rsp), %r8
movl (%rax,%r8,4), %eax
subl %esi, %eax
addl %edi, %eax
addl %r10d, %eax
cltq
movsd (%r12,%rax,8), %xmm0
mulsd (%r11,%r10,8), %xmm0
movq 0x80(%rsp), %rax
movq 0x40(%rsp), %rsi
addsd (%rax,%rsi,8), %xmm0
movsd %xmm0, (%rax,%rsi,8)
testl %edx, %edx
js 0x1d20c
cmpl %r15d, %edx
jl 0x1d215
jmp 0x1d435
movq 0x78(%rsp), %rax
leal (%rax,%r10), %ebp
cmpl $0x0, 0xc(%rsp)
js 0x1d39a
cmpl $-0x1, %ebp
je 0x1d25e
cmpl $-0x1, 0x18(%rsp)
je 0x1d25e
movq 0x10(%rsp), %rax
movq 0x18(%rsp), %rdx
movslq (%rax,%rdx,4), %rax
movslq %ebp, %rdx
addq %rax, %rdx
movsd (%r12,%rdx,8), %xmm0
mulsd (%r11,%r10,8), %xmm0
movq 0x90(%rsp), %rax
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
cmpl %r15d, %r14d
movl %r15d, %eax
movl $0x0, %edx
cmovll %edx, %eax
leal (%rdi,%r10), %edx
subl %eax, %edx
cmpl $-0x1, %edx
je 0x1d2b8
cmpl $-0x1, 0x4c(%rsp)
je 0x1d2b8
movq 0x10(%rsp), %rdx
movq 0x18(%rsp), %rsi
movl (%rdx,%rsi,4), %edx
subl %eax, %edx
addl %edi, %edx
addl %r10d, %edx
movslq %edx, %rax
movsd (%r12,%rax,8), %xmm0
mulsd (%r11,%r10,8), %xmm0
movq 0x90(%rsp), %rax
movq 0x40(%rsp), %rdx
addsd (%rax,%rdx,8), %xmm0
movsd %xmm0, (%rax,%rdx,8)
incl %r13d
incq %r10
addq $-0x8, %r9
cmpl %r10d, 0xa0(%rsp)
jne 0x1d11c
movq 0x70(%rsp), %rdx
leal 0x1(%rdx), %eax
subq %r9, %r11
cmpl 0x20(%rsp), %edx
movl %eax, %r9d
movq 0xd8(%rsp), %rsi
jne 0x1d09a
incq %rcx
incq %rdi
incq 0x78(%rsp)
incq 0x38(%rsp)
cmpq 0x50(%rsp), %rcx
jne 0x1d07d
incq %rsi
cmpq 0x50(%rsp), %rsi
jne 0x1d00b
jmp 0x1cf87
leaq 0x112e7(%rip), %rdi     # 0x2e60c
leaq 0x10f8e(%rip), %rsi     # 0x2e2ba
leaq 0x11278(%rip), %rcx     # 0x2e5ab
movl $0x3cb, %edx            # imm = 0x3CB
callq 0x31c0
leaq 0x111f2(%rip), %rdi     # 0x2e536
leaq 0x10f6f(%rip), %rsi     # 0x2e2ba
leaq 0x11259(%rip), %rcx     # 0x2e5ab
movl $0x426, %edx            # imm = 0x426
callq 0x31c0
leaq 0x112a9(%rip), %rdi     # 0x2e60c
leaq 0x10f50(%rip), %rsi     # 0x2e2ba
leaq 0x1123a(%rip), %rcx     # 0x2e5ab
movl $0x3f1, %edx            # imm = 0x3F1
callq 0x31c0
leaq 0x111b4(%rip), %rdi     # 0x2e536
leaq 0x10f31(%rip), %rsi     # 0x2e2ba
leaq 0x1121b(%rip), %rcx     # 0x2e5ab
movl $0x40b, %edx            # imm = 0x40B
callq 0x31c0
leaq 0x11195(%rip), %rdi     # 0x2e536
leaq 0x10f12(%rip), %rsi     # 0x2e2ba
leaq 0x111fc(%rip), %rcx     # 0x2e5ab
movl $0x46b, %edx            # imm = 0x46B
callq 0x31c0
leaq 0x11176(%rip), %rdi     # 0x2e536
leaq 0x10ef3(%rip), %rsi     # 0x2e2ba
leaq 0x111dd(%rip), %rcx     # 0x2e5ab
movl $0x455, %edx            # imm = 0x455
callq 0x31c0
leaq 0x11157(%rip), %rdi     # 0x2e536
leaq 0x10ed4(%rip), %rsi     # 0x2e2ba
leaq 0x111be(%rip), %rcx     # 0x2e5ab
movl $0x45c, %edx            # imm = 0x45C
callq 0x31c0
leaq 0x11138(%rip), %rdi     # 0x2e536
leaq 0x10eb5(%rip), %rsi     # 0x2e2ba
leaq 0x1119f(%rip), %rcx     # 0x2e5ab
movl $0x44b, %edx            # imm = 0x44B
callq 0x31c0
leaq 0x11119(%rip), %rdi     # 0x2e536
leaq 0x10e96(%rip), %rsi     # 0x2e2ba
leaq 0x11180(%rip), %rcx     # 0x2e5ab
movl $0x446, %edx            # imm = 0x446
callq 0x31c0
leaq 0x110fa(%rip), %rdi     # 0x2e536
leaq 0x10e77(%rip), %rsi     # 0x2e2ba
leaq 0x11161(%rip), %rcx     # 0x2e5ab
movl $0x468, %edx            # imm = 0x468
callq 0x31c0
leaq 0x110db(%rip), %rdi     # 0x2e536
leaq 0x10e58(%rip), %rsi     # 0x2e2ba
leaq 0x11142(%rip), %rcx     # 0x2e5ab
movl $0x3e3, %edx            # imm = 0x3E3
callq 0x31c0
leaq 0x110bc(%rip), %rdi     # 0x2e536
leaq 0x10e39(%rip), %rsi     # 0x2e2ba
leaq 0x11123(%rip), %rcx     # 0x2e5ab
movl $0x3bf, %edx            # imm = 0x3BF
callq 0x31c0
leaq 0x110ff(%rip), %rdi     # 0x2e598
leaq 0x10e1a(%rip), %rsi     # 0x2e2ba
leaq 0x11104(%rip), %rcx     # 0x2e5ab
movl $0x385, %edx            # imm = 0x385
callq 0x31c0
nop