// File: /workspace/github_cpp_makefile_star2_O3/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, %ebx
movq %r8, 0x60(%rsp)
movq %rcx, 0x30(%rsp)
movl %edx, %r15d
movq %rsi, 0x38(%rsp)
movq %rdi, 0x68(%rsp)
cmpb $0x0, 0x128(%rsp)
je 0x1c4c8
movq 0x30(%rsp), %rax
imull %r15d, %eax
movslq %eax, %rdx
shlq $0x3, %rdx
movq 0x120(%rsp), %rdi
xorl %esi, %esi
callq 0x3190
movl %ebx, %eax
andl $0x80000001, %eax       # imm = 0x80000001
cmpl $0x1, %eax
jne 0x1cf8c
movl %ebx, %r14d
shrl %r14d
movq 0x30(%rsp), %rax
movl %eax, %ebp
subl %r14d, %ebp
callq 0x4288
movq %rax, 0x48(%rsp)
movl %ebx, 0x28(%rsp)
movl %ebx, %ebx
movl %r14d, %eax
movq %rax, 0x58(%rsp)
movq %r14, 0x8(%rsp)
movl %ebp, 0x1c(%rsp)
cmpl %ebp, %r14d
jge 0x1c5e4
movl %r15d, %eax
movl 0x28(%rsp), %edi
subl %edi, %eax
incl %eax
movq 0x68(%rsp), %rcx
movq 0x40(%rcx), %rcx
movq %rcx, 0x10(%rsp)
movq 0x8(%rsp), %r8
movl %r8d, %ecx
negl %ecx
movslq %ecx, %rdx
leal 0x1(%r8), %esi
movq 0x30(%rsp), %rcx
leal (%r8,%rcx), %r9d
movl %edi, %ecx
andl $0x7ffffffe, %ecx       # imm = 0x7FFFFFFE
subl %ecx, %r9d
movl %r9d, 0x20(%rsp)
testl %eax, %eax
jle 0x1c5d6
movq 0x10(%rsp), %rcx
leaq (%rcx,%r8,4), %r9
movslq (%r9), %rcx
addq 0x58(%rsp), %rcx
movq 0x120(%rsp), %rdi
leaq (%rdi,%rcx,8), %r14
movq 0x38(%rsp), %r10
xorl %r11d, %r11d
xorpd %xmm0, %xmm0
movq %rdx, %r12
movq 0x60(%rsp), %r13
movslq (%r9,%r12,4), %rcx
leaq (%r10,%rcx,8), %rcx
xorl %ebp, %ebp
xorl %edi, %edi
movsd (%rcx,%rdi,8), %xmm1
mulsd (%r13,%rdi,8), %xmm1
addsd %xmm1, %xmm0
incq %rdi
addq $-0x8, %rbp
cmpl %edi, %ebx
jne 0x1c596
incq %r12
subq %rbp, %r13
cmpl %r12d, %esi
jne 0x1c58a
addsd (%r14), %xmm0
movsd %xmm0, (%r14)
addq $0x8, %r14
incq %r11
addq $0x8, %r10
cmpq %rax, %r11
jne 0x1c57e
incq %r8
cmpl %r8d, 0x20(%rsp)
jne 0x1c555
movl %r15d, %eax
subl 0x8(%rsp), %eax
movq %rax, 0x78(%rsp)
callq 0x4288
subq 0x48(%rsp), %rax
addq %rax, 0x1d59f(%rip)     # 0x39ba0
callq 0x4288
movq %rax, 0xa8(%rsp)
cmpl $0x2, 0x28(%rsp)
jb 0x1c844
movq 0x68(%rsp), %rax
movq 0x40(%rax), %rdx
movq 0x8(%rsp), %rax
movl %eax, %ecx
negl %ecx
movl %ecx, 0x20(%rsp)
movl 0x28(%rsp), %esi
andl $0x7ffffffe, %esi       # imm = 0x7FFFFFFE
movl %r15d, %ecx
movl %esi, 0x40(%rsp)
subl %esi, %ecx
movl %ecx, 0x10(%rsp)
leal 0x1(,%rax,2), %ecx
xorl %r8d, %r8d
movq 0x78(%rsp), %rax
cmpl %eax, 0x8(%rsp)
jge 0x1c713
movslq (%rdx,%r8,4), %rax
movq 0x58(%rsp), %rsi
addq %rsi, %rax
movq 0x120(%rsp), %rsi
leaq (%rsi,%rax,8), %r9
xorl %r10d, %r10d
movslq %r10d, %r11
movq 0x60(%rsp), %r14
movl 0x20(%rsp), %eax
movl %eax, %r12d
movl %r12d, %eax
addl %r8d, %eax
js 0x1c69c
cmpl 0x30(%rsp), %eax
jl 0x1c6a5
jmp 0x1cf4e
addl 0x30(%rsp), %eax
cmpl $-0x1, %eax
je 0x1c6fb
movslq %eax, %r13
movq %r11, %rbp
movl %ecx, %edi
testl %r10d, %r10d
js 0x1ce37
cmpl %r15d, %ebp
jge 0x1ce37
movslq (%rdx,%r13,4), %rax
addq %rbp, %rax
movq 0x38(%rsp), %rsi
movsd (%rsi,%rax,8), %xmm0
mulsd (%r14), %xmm0
addsd (%r9), %xmm0
addq $0x8, %r14
movsd %xmm0, (%r9)
incq %rbp
decl %edi
jne 0x1c6ad
leal 0x1(%r12), %eax
cmpl 0x8(%rsp), %r12d
movl %eax, %r12d
jne 0x1c689
jmp 0x1c701
leaq (%r14,%rbx,8), %r14
jmp 0x1c6ea
addq $0x8, %r9
incl %r10d
cmpl 0x10(%rsp), %r10d
jne 0x1c67a
incq %r8
cmpq 0x58(%rsp), %r8
jne 0x1c650
cmpl $0x2, 0x28(%rsp)
jb 0x1c844
movq 0x68(%rsp), %rax
movq 0x40(%rax), %rdx
movq 0x8(%rsp), %rsi
movl %esi, %eax
negl %eax
movl %eax, 0x48(%rsp)
movq 0x30(%rsp), %rax
movl %eax, %edi
subl 0x40(%rsp), %edi
movslq 0x1c(%rsp), %rcx
cltq
movq %rax, 0x40(%rsp)
leal 0x1(,%rsi,2), %eax
movl %eax, 0x20(%rsp)
movq 0x78(%rsp), %rax
cmpl %eax, 0x8(%rsp)
jge 0x1c834
movslq (%rdx,%rcx,4), %rsi
addq 0x58(%rsp), %rsi
movq 0x120(%rsp), %rax
leaq (%rax,%rsi,8), %r11
xorl %r14d, %r14d
movslq %r14d, %rsi
movq 0x60(%rsp), %r13
movl 0x48(%rsp), %eax
movl %eax, %ebp
testl %edi, %edi
js 0x1cf6d
leal (%rcx,%rbp), %r8d
movq 0x30(%rsp), %r9
cmpl %r9d, %r8d
movl $0x0, %eax
cmovll %eax, %r9d
subl %r9d, %r8d
cmpl $-0x1, %r8d
je 0x1c80c
movslq %r8d, %r8
movq %rsi, %r12
movl 0x20(%rsp), %r9d
testl %r14d, %r14d
js 0x1cef1
cmpl %r15d, %r12d
jge 0x1cef1
movslq (%rdx,%r8,4), %r10
addq %r12, %r10
movq 0x38(%rsp), %rax
movsd (%rax,%r10,8), %xmm0
mulsd (%r13), %xmm0
addsd (%r11), %xmm0
addq $0x8, %r13
movsd %xmm0, (%r11)
incq %r12
decl %r9d
jne 0x1c7ca
jmp 0x1c811
leaq (%r13,%rbx,8), %r13
leal 0x1(%rbp), %r8d
cmpl 0x8(%rsp), %ebp
movl %r8d, %ebp
jne 0x1c799
addq $0x8, %r11
incl %r14d
cmpl 0x10(%rsp), %r14d
jne 0x1c78b
incq %rcx
incl %edi
cmpq 0x40(%rsp), %rcx
jl 0x1c764
movl 0x1c(%rsp), %eax
cmpl %eax, 0x8(%rsp)
jge 0x1caa3
movl %r15d, %eax
movq %rax, 0x20(%rsp)
movq 0x8(%rsp), %r13
movl %r13d, %eax
negl %eax
movslq 0x78(%rsp), %rcx
movq %rcx, 0x88(%rsp)
leal 0x1(%r13), %ecx
movl %ecx, 0x40(%rsp)
movl %eax, 0x48(%rsp)
cltq
movq %rax, 0x70(%rsp)
movl 0x28(%rsp), %ecx
andl $0x7ffffffe, %ecx       # imm = 0x7FFFFFFE
movl %r15d, %eax
subl %ecx, %eax
movl %eax, 0x54(%rsp)
movslq %r15d, %rax
movq %rax, 0x90(%rsp)
movq 0x30(%rsp), %rax
addl %r13d, %eax
subl %ecx, %eax
movl %eax, 0x98(%rsp)
orq $0x1, %rbx
movl %r15d, %eax
subl %r13d, %eax
subl %r13d, %eax
movq %rax, 0x80(%rsp)
leal -0x2(,%r13,2), %eax
incq %rax
movq %rax, 0xa0(%rsp)
cmpl $0x2, 0x28(%rsp)
jb 0x1ca92
movq %r13, %rcx
imulq 0x90(%rsp), %rcx
movq 0x120(%rsp), %rax
leaq (%rax,%rcx,8), %r11
xorl %r8d, %r8d
movq 0x70(%rsp), %rdx
movq %r11, 0x10(%rsp)
movq 0x20(%rsp), %rax
leaq (%rax,%rdx), %r12
movq 0x70(%rsp), %rbp
movq 0x60(%rsp), %r9
leaq (,%rbp,4), %r14
movq 0x68(%rsp), %rax
addq 0x40(%rax), %r14
xorl %r10d, %r10d
xorl %edi, %edi
movl 0x48(%rsp), %eax
movl %eax, %ecx
leaq (%rdx,%rdi), %rsi
testl %esi, %esi
js 0x1c94c
cmpl %r15d, %esi
jl 0x1c95c
jmp 0x1ce56
leal (%r12,%rdi), %esi
cmpl $-0x1, %esi
je 0x1c980
leal (%rcx,%r8), %esi
addl %r15d, %esi
movslq (%r14,%r13,4), %rax
movslq %esi, %rsi
addq %rax, %rsi
movq 0x38(%rsp), %rax
movsd (%rax,%rsi,8), %xmm0
mulsd (%r9,%rdi,8), %xmm0
addsd (%r11), %xmm0
movsd %xmm0, (%r11)
incl %ecx
incq %rdi
addq $-0x8, %r10
cmpl %edi, %ebx
jne 0x1c93a
incq %rbp
subq %r10, %r9
cmpl %ebp, 0x40(%rsp)
jne 0x1c91e
incl %r8d
addq $0x8, %r11
incq %rdx
cmpl 0x8(%rsp), %r8d
jne 0x1c90b
cmpl $0x2, 0x28(%rsp)
jb 0x1ca92
movq 0x88(%rsp), %rax
movq 0x10(%rsp), %rcx
leaq (%rcx,%rax,8), %r14
movq 0xa0(%rsp), %r12
movq 0x80(%rsp), %r11
movl 0x54(%rsp), %r8d
movq 0x78(%rsp), %rax
movq %rax, 0x10(%rsp)
movq 0x70(%rsp), %r9
movq 0x60(%rsp), %rbp
leaq (,%r9,4), %rdi
movq 0x68(%rsp), %rax
addq 0x40(%rax), %rdi
xorl %r10d, %r10d
xorl %ecx, %ecx
movl 0x48(%rsp), %eax
movl %eax, %esi
testl %r8d, %r8d
js 0x1ced2
leaq (%r11,%rcx), %rdx
cmpl 0x20(%rsp), %edx
jl 0x1ca32
cmpl %ecx, %r12d
je 0x1ca56
movq 0x10(%rsp), %rax
leal (%rsi,%rax), %edx
subl %r15d, %edx
movslq (%rdi,%r13,4), %rax
movslq %edx, %rdx
addq %rax, %rdx
movq 0x38(%rsp), %rax
movsd (%rax,%rdx,8), %xmm0
mulsd (%rbp,%rcx,8), %xmm0
addsd (%r14), %xmm0
movsd %xmm0, (%r14)
incl %esi
incq %rcx
addq $-0x8, %r10
cmpl %ecx, %ebx
jne 0x1ca0f
incq %r9
subq %r10, %rbp
cmpl %r9d, 0x40(%rsp)
jne 0x1c9f3
movq 0x10(%rsp), %rcx
incl %ecx
addq $0x8, %r14
incl %r8d
incq %r11
decq %r12
movq %rcx, 0x10(%rsp)
cmpl %r15d, %ecx
jl 0x1c9e9
incq %r13
cmpl %r13d, 0x98(%rsp)
jne 0x1c8db
cmpl $0x2, 0x28(%rsp)
jae 0x1cad0
callq 0x4288
subq 0xa8(%rsp), %rax
addq %rax, 0x1d0ea(%rip)     # 0x39ba8
addq $0xe8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq 0x8(%rsp), %rcx
movl %ecx, %edx
negl %edx
movslq %r15d, %rax
movq %rax, 0xd0(%rsp)
movslq 0x1c(%rsp), %rax
movq %rax, 0xc8(%rsp)
leal (%rcx,%rcx), %eax
leaq 0x1(%rax), %r13
movl %r15d, %esi
subl %eax, %esi
movq %rsi, 0xc0(%rsp)
movl %r15d, %eax
subl %ecx, %eax
movq %rax, 0xb8(%rsp)
movl %edx, 0x1c(%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, 0x90(%rsp)
leaq (%rdx,%rcx,8), %rax
movq %rax, 0x88(%rsp)
movq 0xb0(%rsp), %rax
movq %rax, 0x28(%rsp)
movq 0xb8(%rsp), %rax
movq %rax, 0x80(%rsp)
movq 0xc0(%rsp), %rcx
xorl %eax, %eax
movq %rsi, 0xd8(%rsp)
movq 0x78(%rsp), %rdx
addl %eax, %edx
movslq %edx, %rdx
movq %rdx, 0x98(%rsp)
movq 0x60(%rsp), %r11
movl 0x1c(%rsp), %edx
movl %edx, %r9d
movl %r9d, %edx
addl %esi, %edx
sets %sil
movq 0x30(%rsp), %r8
cmpl %r8d, %edx
setl %r10b
movq 0x68(%rsp), %rdi
movq 0x40(%rdi), %r14
movq 0xe0(%rsp), %rdi
movq %r9, 0xa0(%rsp)
leal (%r9,%rdi), %ebx
cmpl %r8d, %ebx
movl %r8d, %edi
movl $0x0, %r9d
cmovll %r9d, %edi
movl %ebx, 0x40(%rsp)
subl %edi, %ebx
movl %ebx, 0x10(%rsp)
movslq %ebx, %rdi
movq %rdi, 0x70(%rsp)
orb %sil, %r10b
movb %r10b, 0x48(%rsp)
movl %edx, %esi
sarl $0x1f, %esi
andl %r8d, %esi
addl %edx, %esi
movl %esi, 0x54(%rsp)
movslq %esi, %rdx
movq %rdx, 0x20(%rsp)
xorl %edi, %edi
xorl %r8d, %r8d
movl 0x1c(%rsp), %edx
movl %edx, %r12d
leal (%r12,%rax), %r10d
movq 0x28(%rsp), %rdx
addq %r8, %rdx
testl %edx, %edx
js 0x1cc51
movl %r10d, %r9d
cmpl %r15d, %edx
jl 0x1cc5d
jmp 0x1cf2f
movq 0x80(%rsp), %rsi
leal (%rsi,%r8), %r9d
cmpb $0x0, 0x48(%rsp)
je 0x1ceb3
cmpl $-0x1, %r9d
je 0x1cca7
cmpl $-0x1, 0x54(%rsp)
je 0x1cca7
movq 0x20(%rsp), %rsi
movslq (%r14,%rsi,4), %rbp
movslq %r9d, %r9
addq %rbp, %r9
movq 0x38(%rsp), %rbx
movsd (%rbx,%r9,8), %xmm0
mulsd (%r11,%r8,8), %xmm0
movq 0x90(%rsp), %rsi
addsd (%rsi,%rax,8), %xmm0
movsd %xmm0, (%rsi,%rax,8)
leaq (%rcx,%r8), %r9
testl %r9d, %r9d
js 0x1ce75
cmpl %r15d, %r9d
movl %r15d, %ebp
movl $0x0, %esi
cmovll %esi, %ebp
leal (%rcx,%r8), %r9d
subl %ebp, %r9d
cmpl $-0x1, %r9d
je 0x1cd21
cmpl $-0x1, 0x20(%rsp)
je 0x1cd21
movq 0x20(%rsp), %rsi
movl (%r14,%rsi,4), %ebx
subl %ebp, %ebx
movl %r15d, %esi
movq %r13, %r15
addl %ecx, %ebx
addl %r8d, %ebx
movslq %ebx, %rbx
movq 0x38(%rsp), %r13
movsd (%r13,%rbx,8), %xmm0
movq %r15, %r13
movl %esi, %r15d
mulsd (%r11,%r8,8), %xmm0
movq 0x90(%rsp), %rsi
movq 0x98(%rsp), %rbx
addsd (%rsi,%rbx,8), %xmm0
movsd %xmm0, (%rsi,%rbx,8)
testl %edx, %edx
js 0x1cd2f
cmpl %r15d, %edx
jl 0x1cd3b
jmp 0x1cf10
movq 0x80(%rsp), %rdx
leal (%rdx,%r8), %r10d
cmpl $0x0, 0x40(%rsp)
js 0x1ce94
cmpl $-0x1, %r10d
je 0x1cd85
cmpl $-0x1, 0x10(%rsp)
je 0x1cd85
movq 0x70(%rsp), %rdx
movslq (%r14,%rdx,4), %rdx
movslq %r10d, %r10
addq %rdx, %r10
movq 0x38(%rsp), %rdx
movsd (%rdx,%r10,8), %xmm0
mulsd (%r11,%r8,8), %xmm0
movq 0x88(%rsp), %rdx
addsd (%rdx,%rax,8), %xmm0
movsd %xmm0, (%rdx,%rax,8)
cmpl $-0x1, %r9d
je 0x1cdd0
cmpl $-0x1, 0x10(%rsp)
je 0x1cdd0
movq 0x70(%rsp), %rdx
movl (%r14,%rdx,4), %edx
subl %ebp, %edx
addl %ecx, %edx
addl %r8d, %edx
movslq %edx, %rdx
movq 0x38(%rsp), %r9
movsd (%r9,%rdx,8), %xmm0
mulsd (%r11,%r8,8), %xmm0
movq 0x88(%rsp), %rdx
movq 0x98(%rsp), %rsi
addsd (%rdx,%rsi,8), %xmm0
movsd %xmm0, (%rdx,%rsi,8)
incl %r12d
incq %r8
addq $-0x8, %rdi
cmpl %r8d, %r13d
jne 0x1cc34
movq 0xa0(%rsp), %rsi
leal 0x1(%rsi), %edx
subq %rdi, %r11
cmpl 0x8(%rsp), %esi
movl %edx, %r9d
movq 0xd8(%rsp), %rsi
jne 0x1cbb6
incq %rax
incq %rcx
incq 0x80(%rsp)
incq 0x28(%rsp)
cmpq 0x58(%rsp), %rax
jne 0x1cb98
incq %rsi
cmpq 0x58(%rsp), %rsi
jne 0x1cb23
jmp 0x1caaa
leaq 0x10836(%rip), %rdi     # 0x2d674
leaq 0x104dd(%rip), %rsi     # 0x2d322
leaq 0x107c7(%rip), %rcx     # 0x2d613
movl $0x3cb, %edx            # imm = 0x3CB
callq 0x31c0
leaq 0x10741(%rip), %rdi     # 0x2d59e
leaq 0x104be(%rip), %rsi     # 0x2d322
leaq 0x107a8(%rip), %rcx     # 0x2d613
movl $0x40b, %edx            # imm = 0x40B
callq 0x31c0
leaq 0x10722(%rip), %rdi     # 0x2d59e
leaq 0x1049f(%rip), %rsi     # 0x2d322
leaq 0x10789(%rip), %rcx     # 0x2d613
movl $0x455, %edx            # imm = 0x455
callq 0x31c0
leaq 0x10703(%rip), %rdi     # 0x2d59e
leaq 0x10480(%rip), %rsi     # 0x2d322
leaq 0x1076a(%rip), %rcx     # 0x2d613
movl $0x46b, %edx            # imm = 0x46B
callq 0x31c0
leaq 0x106e4(%rip), %rdi     # 0x2d59e
leaq 0x10461(%rip), %rsi     # 0x2d322
leaq 0x1074b(%rip), %rcx     # 0x2d613
movl $0x44b, %edx            # imm = 0x44B
callq 0x31c0
leaq 0x106c5(%rip), %rdi     # 0x2d59e
leaq 0x10442(%rip), %rsi     # 0x2d322
leaq 0x1072c(%rip), %rcx     # 0x2d613
movl $0x426, %edx            # imm = 0x426
callq 0x31c0
leaq 0x1077c(%rip), %rdi     # 0x2d674
leaq 0x10423(%rip), %rsi     # 0x2d322
leaq 0x1070d(%rip), %rcx     # 0x2d613
movl $0x3f1, %edx            # imm = 0x3F1
callq 0x31c0
leaq 0x10687(%rip), %rdi     # 0x2d59e
leaq 0x10404(%rip), %rsi     # 0x2d322
leaq 0x106ee(%rip), %rcx     # 0x2d613
movl $0x468, %edx            # imm = 0x468
callq 0x31c0
leaq 0x10668(%rip), %rdi     # 0x2d59e
leaq 0x103e5(%rip), %rsi     # 0x2d322
leaq 0x106cf(%rip), %rcx     # 0x2d613
movl $0x446, %edx            # imm = 0x446
callq 0x31c0
leaq 0x10649(%rip), %rdi     # 0x2d59e
leaq 0x103c6(%rip), %rsi     # 0x2d322
leaq 0x106b0(%rip), %rcx     # 0x2d613
movl $0x3bf, %edx            # imm = 0x3BF
callq 0x31c0
leaq 0x1062a(%rip), %rdi     # 0x2d59e
leaq 0x103a7(%rip), %rsi     # 0x2d322
leaq 0x10691(%rip), %rcx     # 0x2d613
movl $0x3e3, %edx            # imm = 0x3E3
callq 0x31c0
leaq 0x1066d(%rip), %rdi     # 0x2d600
leaq 0x10388(%rip), %rsi     # 0x2d322
leaq 0x10672(%rip), %rcx     # 0x2d613
movl $0x385, %edx            # imm = 0x385
callq 0x31c0
nop