// File: /workspace/github_c_makefile_star2_O3/markpeerdeman[P]mbemm/src/emath/solve.c
// Function: solve_4
void solve_4(matrix Ma, int m, int n,const vector F,intvector nietnul)
{ 
  int    a,i,j,k,l;
  double c,tmp;
  double fabsc, fabsM;
  vector derij; 
  int	 tel1, tel2, tel3, tel4;
  matrix hulpma;
  vector hulpve;
 
  evaluate_level(GREEN,38);
  check_dimensions(Ma, F);
  
  derij=(vector)newvector(m+n);
  
  
  /* First we sweep the upper part of the matrix */
  /* into an upper triangular matrix */

  for (k=0; k<m; k++)
    { 

fprintf(stderr, "processing row: %ld\n", (long)m);

      /* divide row by diagonal coefficient: */
      
      c = M(Ma,k,k); l=k;
      
      /* find largest coefficient for elimination: */
      
      for (i=k; i<m; i++)
	{ 
	  fabsc=fabs(c);
	  fabsM=fabs(M(Ma,i,k));
	  
	  /* if (fabs(M(A,m,k))>=fabs(c)) */
	  
	  if (fabsM>=fabsc)
	    { 
	      c = M(Ma,i,k);
	      l = i;
	    }
	}
      
      /* if (fabs(c)<=TOL) 
	 evaluate_level(RED,33); */
      /* change rows if nescessary: */
      
      if (k!=l)
	{ 
	  tmp = V(F,l);
	  V(F,l) = V(F,k);
	  V(F,k) = tmp;	
	  
	  for (j=k; j<m+n; j++)
	    {
	      if (j<m) 
		{
		  tmp = M(Ma,l,j);
		  M(Ma,l,j)=M(Ma,k,j);
		  M(Ma,k,j)=tmp;	    
		}
	      else 
		if (V(nietnul,j)!=0)
		  {
		    tmp = M(Ma,l,m+V(nietnul,j)-1);
		    M(Ma,l,m+V(nietnul,j)-1)=M(Ma,k,m+V(nietnul,j)-1); 	    
		    M(Ma,k,m+V(nietnul,j)-1) = tmp;
		}
	    }
	}
      
      /* divide current row by diagonal coefficient: */
      
      M(Ma,k,k)=1.0;
      V(F,k) = V(F,k) / c;
      
      for (j=k+1; j<m+n; j++)
	{  
	  if (j<m)
	    M(Ma,k,j) = M(Ma,k,j) / c;
	  else
	    if (V(nietnul,j)!=0) M(Ma,k,m+V(nietnul,i)-1)=M(Ma,k,m+V(nietnul,j)-1)/c;
	}
      
      /* divide remaining rows by diagonal coefficient and eliminate: */
      
      for (i=k+1; i<m; i++)
	{ 
	  c = M(Ma,i,k);
	  M(Ma,i,k)=0.0;
	  V(F,i)=V(F,i)-c*V(F,k);
	  for (j=k+1; j<m+n; j++)
	    { 
	      if (j<m)
		{
		  M(Ma,i,j)=M(Ma,i,j)-c*M(Ma,k,j);
		}
	      else
		if (V(nietnul,j)!=0)
		  {
		    M(Ma,i,m+V(nietnul,j)-1) = M(Ma,i,m+V(nietnul,j)-1) -
		      c*M(Ma,i,m+V(nietnul,j)-1);
		  }
	    }
	}
      /* tijdvalkje(k+2,A); */
    }
  
  
  /* We now sweep the lower part of the matrix */
  /* into a lower triangular matrix */
  
  for (k=m+n-1; k>=m; k--)
    { 
      /* divide row by diagonal coefficient: */
      
      c = M(Ma,k,k); l=k;
      
      /* find largest coefficient for elimination: */
      
      for (i=k; i>=m; i--)
	{ 
	  fabsc=fabs(c);
	  fabsM=fabs(M(Ma,i,k));
	  
	  /* if (fabs(M(A,m,k))>=fabs(c)) */
	  
	  if (fabsM>=fabsc)
	    { 
	      c = M(Ma,i,k);
	      l = m;
	    }
	}
            
      /* if (fabs(c)<=TOL) 
	 evaluate_level(RED,33); */
      /* change rows if nescessary: */
      
      if (k!=l)
	{ 
	  tmp = V(F,l);
	  V(F,l) = V(F,k);
	  V(F,k) = tmp;
	  
	  for (j=k; j>=0; j--)
	    {
	      if (j>=m) 
		{
		  tmp = M(Ma,l,j);
		  M(Ma,l,j)=M(Ma,k,j);  
		  M(Ma,k,j)=tmp;      
		}
	      else 
		if (V(nietnul,j)!=0)
		  {
		    a=m-p+V(nietnul,j)-1;
		    tmp = M(Ma,l,a);
		    M(Ma,l,a)=M(Ma,k,a);	 
		    M(Ma,k,a) = tmp;
		  }
	    }
	}
            
      /* divide current row by diagonal coefficient: */
      
      M(Ma,k,k)=1.0;
      V(F,k) = V(F,k) / c; 
      
      for (j=k-1; j>=0; j--)
	{  
	  if (j>=m)
	    M(Ma,k,j) = M(Ma,k,j) / c;
	  else
	    if (V(nietnul,j)!=0) 
	      M(Ma,k,m-p+V(nietnul,j)-1)=M(Ma,k,m-p+V(nietnul,j)-1)/c;
	}
      
      /* divide remaining rows by diagonal coefficient and eliminate: */
      
      for (i=k-1; i>=m; i--)
	{ 
	  c = M(Ma,i,k);
	  M(Ma,i,k)=0.0;
	  V(F,i)=V(F,i)-c*V(F,k);

	  for (j=k-1; j>=0; j--)
	    { 
	      if (j>=m)
		{
		  M(Ma,i,j)=M(Ma,i,j)-c*M(Ma,k,j);
		}
	      else
		if (V(nietnul,j)!=0)
		  {
		    M(Ma,i,m-p+V(nietnul,j)-1) = M(Ma,i,m-p+V(nietnul,i)-1) -
		      c*M(Ma,i,m-p+V(nietnul,i)-1);
		  }
	    }
	}
      /*	tijdvalkje(k+2,A);  */
    }


  /* We are now left with an upper triangular matrix and a lower triangular matrix */
  /* We now try to find an approximation for the variables corresponding to the */
  /* non-zero colums */

  hulpma=(matrix)newmatrix(p+q,p+q);
  hulpve=(vector)newvector(p+q);

  /* We now fill the additional small matrix with the right elements */  

  tel1=0;
  for (i=0; i<m+n; i++)
    if (V(nietnul,i)!=0)
      {
	V(hulpve,tel1)=V(F,i);
	
	tel2=0;
	for (j=0; j<m+n;j++)
	  if (V(nietnul,j)!=0)
	    {
	      M(hulpma,tel1,tel2)=hetelement(m,i,j,nietnul,Ma);
	      tel2++;
	    }
	tel1++;
      }

  printmatrix(hulpma,"hulpmatrix");
  printvector(hulpve,"hulpvector");
 
  printf("No 4\n");

  /* We now solve the small additonal matrix */
  
  for (k=0; k<p+q; k++)    
    {
      c=M(hulpma,k,k);
      M(hulpma,k,k)=1.0;
      V(hulpve,k) = V(hulpve,k) / c; 
      
      for (j=k-1; j>=0; j--)
	M(hulpma,k,j)=M(Ma,k,j)/c;
	
      
      /* divide remaining rows by diagonal coefficient and eliminate: */
      
      for (i=k+1; i<p+q; i++)
	{ 
	  c = M(hulpma,i,k);
	  M(hulpma,i,k)=0.0;
	  V(hulpve,i)=V(hulpve,i)-c*V(hulpve,k);

	  for (j=k+1; j<p+q; j++)
	    M(hulpma,i,j)=M(hulpma,i,j)-c*M(hulpma,k,j);
	}
    }	

  /* apply backsubstitution to compute the unknowns in hulpve */
  
  for (i=p+q-2;i>=0;i--)
    {
      for (j=p+q-1;j>i;j--)
	{
	  V(hulpve,i) = V(hulpve,i) - M(hulpma,i,j)*V(hulpve,i);
	}
    }

  /* We put the approximation into the vector F */  

  for (i=0; i<m+n; i++)
    {
      tel1=0;
      if (V(nietnul,i)!=0)
	{
	  for (j=0; j<m+n; j++)
	    if (V(nietnul,j)!=0)
	      {
		V(F,i)=V(F,i)-hetelement(m,i,j,nietnul,Ma)*V(hulpve,tel1);
		tel1++;
	      }
	}
    }

  /* Now we can solve the upper triangular matrix */
	  
  for (i=m-2;i>=0;i--)
    {
       for (j=m-1;j>i;j--)
	 {
	   V(F,i) = V(F,i) - M(Ma,i,j)*V(F,i);
	 }
    }

  /* And then we solve the lower triangular matrix */
	  
  for (i=m+1;i<m+n;i++)
    {
      for (j=m;j<i;j++)
	{
	  V(F,i) = V(F,i) - M(Ma,i,j)*V(F,i);
	}
    }
  
  /*  tijdvalkje(-1,A); */
  
  
  evaluate_level(GREEN,34);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x68, %rsp
movq %r8, %r13
movq %rcx, %rbx
movl %edx, %r12d
movl %esi, %r14d
movq %rdi, %r15
pushq $0x26
popq %rsi
xorl %edi, %edi
callq 0xb980
movq %r15, 0x28(%rsp)
movq %r15, %rdi
movq %rbx, 0x38(%rsp)
movq %rbx, %rsi
callq 0x11e0
movq %r12, 0x50(%rsp)
leal (%r12,%r14), %edi
movq %rdi, 0x18(%rsp)
callq 0xbb17
movslq %r14d, %rax
movq %rax, 0x30(%rsp)
testl %r14d, %r14d
jle 0x14bc
movl %r14d, %ebx
movslq 0x18(%rsp), %r12
pushq $0x1
popq %rbp
xorl %esi, %esi
movq %rsi, %r15
movq 0x13d61(%rip), %rax     # 0x14ff8
movq (%rax), %rdi
leaq 0xce05(%rip), %rsi      # 0xe0a6
movq %rbx, %rdx
xorl %eax, %eax
callq 0x1100
movapd 0xcd5d(%rip), %xmm4     # 0xe010
movq 0x28(%rsp), %rax
movq 0x8(%rax), %rax
movq (%rax,%r15,8), %rcx
movsd (%rcx,%r15,8), %xmm0
movq %r15, %rdx
movl %r15d, %esi
movapd %xmm0, %xmm1
andpd %xmm4, %xmm1
movq (%rax,%rdx,8), %rdi
movsd (%rdi,%r15,8), %xmm2
movapd %xmm2, %xmm3
andpd %xmm4, %xmm3
ucomisd %xmm1, %xmm3
cmpnlesd %xmm3, %xmm1
andpd %xmm1, %xmm0
andnpd %xmm2, %xmm1
orpd %xmm1, %xmm0
cmovael %edx, %esi
incq %rdx
cmpq %rdx, %rbx
jne 0x12cc
movl %esi, %edi
movq 0x38(%rsp), %rdx
movq 0x8(%rdx), %rdx
cmpq %rdi, %r15
movapd 0xcd04(%rip), %xmm3     # 0xe020
je 0x139b
movslq %esi, %rsi
movsd (%rdx,%rsi,8), %xmm1
movsd (%rdx,%r15,8), %xmm2
movsd %xmm2, (%rdx,%rsi,8)
movsd %xmm1, (%rdx,%r15,8)
cmpq %r12, %r15
jge 0x139b
movq %r15, %rdi
cmpq %rbx, %rdi
jae 0x1360
movq (%rax,%rsi,8), %r8
movsd (%r8,%rdi,8), %xmm1
movsd (%rcx,%rdi,8), %xmm2
movsd %xmm2, (%r8,%rdi,8)
movsd %xmm1, (%rcx,%rdi,8)
jmp 0x1393
movq 0x8(%r13), %r8
movl (%r8,%rdi,4), %r8d
testl %r8d, %r8d
je 0x1393
movq (%rax,%rsi,8), %r9
addl %r14d, %r8d
movslq %r8d, %r8
movsd -0x8(%r9,%r8,8), %xmm1
movsd -0x8(%rcx,%r8,8), %xmm2
movsd %xmm2, -0x8(%r9,%r8,8)
movsd %xmm1, -0x8(%rcx,%r8,8)
incq %rdi
cmpq %r12, %rdi
jl 0x133f
movabsq $0x3ff0000000000000, %rsi # imm = 0x3FF0000000000000
movq %rsi, (%rcx,%r15,8)
movsd (%rdx,%r15,8), %xmm1
divsd %xmm0, %xmm1
movsd %xmm1, (%rdx,%r15,8)
leaq 0x1(%r15), %rsi
movq %rbp, %rdi
cmpq %r12, %rsi
jge 0x1411
cmpq %rbx, %rdi
jae 0x13da
movsd (%rcx,%rdi,8), %xmm1
divsd %xmm0, %xmm1
movsd %xmm1, (%rcx,%rdi,8)
jmp 0x1408
movq 0x8(%r13), %r8
movl (%r8,%rdi,4), %r9d
testl %r9d, %r9d
je 0x1408
addl %r14d, %r9d
movslq %r9d, %r9
movsd -0x8(%rcx,%r9,8), %xmm1
divsd %xmm0, %xmm1
movslq (%r8,%rbx,4), %r8
addq 0x30(%rsp), %r8
movsd %xmm1, -0x8(%rcx,%r8,8)
incq %rdi
cmpl %edi, 0x18(%rsp)
jg 0x13c5
movq %rbp, %rdi
cmpq %rbx, %rsi
jae 0x14b0
movq (%rax,%rdi,8), %r8
movsd (%r8,%r15,8), %xmm0
andq $0x0, (%r8,%r15,8)
movsd (%rdx,%rdi,8), %xmm1
movsd (%rdx,%r15,8), %xmm2
mulsd %xmm0, %xmm2
subsd %xmm2, %xmm1
movsd %xmm1, (%rdx,%rdi,8)
cmpq %r12, %rsi
jge 0x14a4
xorpd %xmm3, %xmm0
movq %rbp, %r9
cmpq %rbx, %r9
jae 0x146d
movsd (%rcx,%r9,8), %xmm1
mulsd %xmm0, %xmm1
addsd (%r8,%r9,8), %xmm1
movsd %xmm1, (%r8,%r9,8)
jmp 0x149a
movq 0x8(%r13), %r10
movl (%r10,%r9,4), %r10d
testl %r10d, %r10d
je 0x149a
addl %r14d, %r10d
movslq %r10d, %r10
movsd -0x8(%r8,%r10,8), %xmm1
movapd %xmm0, %xmm2
mulsd %xmm1, %xmm2
addsd %xmm1, %xmm2
movsd %xmm2, -0x8(%r8,%r10,8)
incq %r9
cmpl %r9d, 0x18(%rsp)
jg 0x1450
incq %rdi
cmpq %rbx, %rdi
jne 0x141d
incq %rbp
cmpq %rbx, %rsi
jne 0x128d
movslq 0x18(%rsp), %rax
movq %rax, 0x48(%rsp)
cmpl $0x0, 0x50(%rsp)
movq 0x30(%rsp), %r15
jle 0x181b
movq 0x28(%rsp), %rax
movq 0x8(%rax), %r10
leal -0x1(%r14), %eax
movl %eax, 0x44(%rsp)
movq 0x38(%rsp), %rax
movq 0x8(%rax), %rdx
movq 0x18(%rsp), %rcx
leal -0x1(%rcx), %eax
movslq %eax, %rsi
leal -0x2(%rcx), %edi
leaq 0x2(%r15), %rax
movq %rax, 0x58(%rsp)
movq %rsi, 0x10(%rsp)
leaq 0x1(%rsi), %r9
movapd 0xcaf8(%rip), %xmm0     # 0xe010
movapd 0xcb00(%rip), %xmm1     # 0xe020
movq 0x48(%rsp), %r11
movq %r10, 0x8(%rsp)
movq %r11, %rbp
movl %edi, %ebx
leaq (,%rbx,8), %rax
movq %rax, 0x20(%rsp)
decq %r11
movq -0x8(%r10,%rbp,8), %r12
movsd -0x8(%r12,%rbp,8), %xmm2
cmpq %r15, %rbp
movq %r9, 0x60(%rsp)
jle 0x164b
movl 0x10(%rsp), %r8d
leaq (,%r8,8), %rcx
movl %r11d, %esi
movapd %xmm2, %xmm3
andpd %xmm0, %xmm3
movq -0x8(%r10,%r9,8), %rax
movsd (%rax,%r11,8), %xmm4
movapd %xmm4, %xmm5
andpd %xmm0, %xmm5
ucomisd %xmm3, %xmm5
cmpnlesd %xmm5, %xmm3
andpd %xmm3, %xmm2
andnpd %xmm4, %xmm3
orpd %xmm3, %xmm2
cmovael %r14d, %esi
decq %r9
cmpq %r15, %r9
jg 0x1569
cmpl %r11d, %esi
je 0x164b
movslq %esi, %r10
movsd (%rdx,%r10,8), %xmm3
movsd (%rdx,%r11,8), %xmm4
movsd %xmm4, (%rdx,%r10,8)
movsd %xmm3, (%rdx,%r11,8)
testq %rbp, %rbp
jle 0x17f1
leaq 0x14317(%rip), %rax     # 0x158f0
movl (%rax), %eax
movl %eax, (%rsp)
movl %ebp, %r9d
cmpl %r14d, %r8d
jge 0x160f
movq 0x8(%r13), %rax
movl (%rax,%r8,4), %esi
testl %esi, %esi
je 0x1637
addl 0x44(%rsp), %esi
subl (%rsp), %esi
movslq %esi, %r15
leaq (,%r15,8), %rsi
movq 0x8(%rsp), %rax
addq (%rax,%r10,8), %rsi
jmp 0x161e
movq 0x8(%rsp), %rax
movq (%rax,%r10,8), %rsi
addq %rcx, %rsi
movq %r8, %r15
movsd (%rsi), %xmm3
movsd (%r12,%r15,8), %xmm4
movsd %xmm4, (%rsi)
movsd %xmm3, (%r12,%r15,8)
movq 0x30(%rsp), %r15
decq %r8
decl %r9d
addq $-0x8, %rcx
testl %r9d, %r9d
jg 0x15e1
movq 0x8(%rsp), %r10
movabsq $0x3ff0000000000000, %rax # imm = 0x3FF0000000000000
movq %rax, (%r12,%r11,8)
movsd (%rdx,%r11,8), %xmm3
divsd %xmm2, %xmm3
movsd %xmm3, (%rdx,%r11,8)
cmpq $0x2, %rbp
setge %al
movl %eax, (%rsp)
jl 0x16e2
leaq -0x8(%r12), %rcx
leaq 0x1426f(%rip), %rax     # 0x158f0
movl (%rax), %esi
movq 0x20(%rsp), %rax
leaq (%r12,%rax), %r9
movq 0x10(%rsp), %rax
movl %eax, %r10d
movq %rbx, %r15
decl %r10d
movq %r9, %r8
cmpl %r14d, %r10d
jge 0x16bc
movq 0x8(%r13), %rax
movl (%rax,%r15,4), %r8d
testl %r8d, %r8d
je 0x16ca
addl %r14d, %r8d
subl %esi, %r8d
movslq %r8d, %rax
leaq (%rcx,%rax,8), %r8
movsd (%r8), %xmm3
divsd %xmm2, %xmm3
movsd %xmm3, (%r8)
decq %r15
addq $-0x8, %r9
testl %r10d, %r10d
jg 0x1697
movq 0x30(%rsp), %r15
movq 0x8(%rsp), %r10
jmp 0x16e9
movl $0x0, (%rsp)
cmpq 0x58(%rsp), %rbp
jl 0x17d6
shlq $0x2, %rbx
movslq %edi, %rcx
addq 0x20(%rsp), %r12
movq (%r10,%rcx,8), %rbp
movsd (%rbp,%r11,8), %xmm2
andq $0x0, (%rbp,%r11,8)
movsd (%rdx,%rcx,8), %xmm3
movsd (%rdx,%r11,8), %xmm4
mulsd %xmm2, %xmm4
subsd %xmm4, %xmm3
movsd %xmm3, (%rdx,%rcx,8)
cmpb $0x0, (%rsp)
je 0x17bc
xorpd %xmm1, %xmm2
leaq 0x141b2(%rip), %rax     # 0x158f0
movslq (%rax), %rax
movq 0x30(%rsp), %r10
subq %rax, %r10
movq 0x20(%rsp), %rax
leaq (%rax,%rbp), %r9
xorl %r15d, %r15d
movq 0x10(%rsp), %rax
movl %eax, %esi
leal (%rdi,%r15), %eax
cmpl %r14d, %eax
jge 0x179d
movq 0x8(%r13), %rax
leaq (%rax,%rbx), %r8
movl (%r8,%r15,4), %r8d
testl %r8d, %r8d
je 0x17b3
movslq (%rax,%rcx,4), %rax
addq %r10, %rax
movsd -0x8(%rbp,%rax,8), %xmm3
movapd %xmm2, %xmm4
mulsd %xmm3, %xmm4
addsd %xmm3, %xmm4
addl %r10d, %r8d
movslq %r8d, %rax
movsd %xmm4, -0x8(%rbp,%rax,8)
jmp 0x17b3
movsd (%r12,%r15,8), %xmm3
mulsd %xmm2, %xmm3
addsd (%r9,%r15,8), %xmm3
movsd %xmm3, (%r9,%r15,8)
decl %esi
decq %r15
testl %esi, %esi
jg 0x175c
leaq -0x1(%rcx), %rax
movq 0x30(%rsp), %r15
cmpq %r15, %rcx
movq %rax, %rcx
movq 0x8(%rsp), %r10
jg 0x1700
decq 0x10(%rsp)
decq %rdi
movq 0x60(%rsp), %r9
decq %r9
cmpq %r15, %r11
jg 0x152a
jmp 0x181b
movabsq $0x3ff0000000000000, %rax # imm = 0x3FF0000000000000
movq %rax, (%r12,%r11,8)
movsd (%rdx,%r11,8), %xmm3
divsd %xmm2, %xmm3
movsd %xmm3, (%rdx,%r11,8)
movl $0x0, (%rsp)
jmp 0x16db
leaq 0x140ce(%rip), %rbx     # 0x158f0
leaq 0x140cb(%rip), %r15     # 0x158f4
movl (%r15), %edi
addl (%rbx), %edi
movl %edi, %esi
callq 0xb9e0
movq %rax, 0x10(%rsp)
movl (%r15), %edi
addl (%rbx), %edi
callq 0xbb17
movq %rax, (%rsp)
movq 0x18(%rsp), %rax
movl %eax, %ecx
movq %rcx, 0x8(%rsp)
testl %eax, %eax
movq 0x38(%rsp), %r12
jle 0x18f0
xorl %r15d, %r15d
xorl %ecx, %ecx
movq 0x8(%r13), %rax
cmpl $0x0, (%rax,%r15,4)
je 0x18e2
movq 0x8(%r12), %rax
movsd (%rax,%r15,8), %xmm0
movq (%rsp), %rax
movq 0x8(%rax), %rax
movl %ecx, 0x20(%rsp)
movslq %ecx, %rbx
movsd %xmm0, (%rax,%rbx,8)
xorl %ebp, %ebp
xorl %r12d, %r12d
movq 0x8(%r13), %rax
cmpl $0x0, (%rax,%rbp,4)
je 0x18cd
movl %r14d, %edi
movl %r15d, %esi
movl %ebp, %edx
movq %r13, %rcx
movq 0x28(%rsp), %r8
callq 0xcac2
movq 0x10(%rsp), %rax
movq 0x8(%rax), %rax
movq (%rax,%rbx,8), %rax
movslq %r12d, %r12
movsd %xmm0, (%rax,%r12,8)
incl %r12d
incq %rbp
cmpq %rbp, 0x8(%rsp)
jne 0x1895
movl 0x20(%rsp), %ecx
incl %ecx
movq 0x38(%rsp), %r12
incq %r15
cmpq 0x8(%rsp), %r15
jne 0x1866
leaq 0xc7c4(%rip), %rsi      # 0xe0bb
movq 0x10(%rsp), %rdi
callq 0xc528
leaq 0xc7be(%rip), %rsi      # 0xe0c6
movq (%rsp), %rdi
callq 0xc610
leaq 0xc7b9(%rip), %rdi      # 0xe0d1
callq 0x1060
leaq 0x13fd0(%rip), %rax     # 0x158f4
movl (%rax), %eax
leaq 0x13fc3(%rip), %rcx     # 0x158f0
addl (%rcx), %eax
jle 0x1a93
movq 0x10(%rsp), %rcx
movq 0x8(%rcx), %rcx
movq (%rsp), %rdx
movq 0x8(%rdx), %rdx
movl %eax, %ebp
pushq $0x1
popq %rdi
xorl %r8d, %r8d
movapd 0xc6ca(%rip), %xmm0     # 0xe020
movq %rdi, %r9
xorl %r15d, %r15d
movq %rbp, 0x20(%rsp)
movq %r15, %r10
movq (%rcx,%r15,8), %r11
movsd (%r11,%r15,8), %xmm1
movabsq $0x3ff0000000000000, %rsi # imm = 0x3FF0000000000000
movq %rsi, (%r11,%r15,8)
movsd (%rdx,%r15,8), %xmm2
divsd %xmm1, %xmm2
movsd %xmm2, (%rdx,%r15,8)
testq %r15, %r15
je 0x19be
movq 0x28(%rsp), %rsi
movq 0x8(%rsi), %rbx
movq (%rbx,%r10,8), %rbx
movq %r9, %r15
leal -0x2(%r15), %r12d
movsd (%rbx,%r12,8), %xmm2
divsd %xmm1, %xmm2
movsd %xmm2, (%r11,%r12,8)
decq %r15
cmpq $0x1, %r15
jg 0x19a1
leaq 0x1(%r10), %r15
cmpq %rbp, %r15
jae 0x1a26
addq %r8, %r11
movq %r9, %rbp
movq (%rcx,%rbp,8), %rbx
movsd (%rbx,%r10,8), %xmm2
andq $0x0, (%rbx,%r10,8)
movsd (%rdx,%rbp,8), %xmm3
movapd %xmm2, %xmm1
xorpd %xmm0, %xmm1
mulsd (%rdx,%r10,8), %xmm2
subsd %xmm2, %xmm3
movsd %xmm3, (%rdx,%rbp,8)
addq %r8, %rbx
movq %rdi, %r12
movsd (%r11,%r12,8), %xmm2
mulsd %xmm1, %xmm2
addsd (%rbx,%r12,8), %xmm2
movsd %xmm2, (%rbx,%r12,8)
incq %r12
leal (%r10,%r12), %esi
cmpl %esi, %eax
jg 0x19fe
incq %rbp
cmpl %ebp, %eax
jg 0x19cd
incq %r9
addq $0x8, %r8
movq 0x20(%rsp), %rbp
cmpq %rbp, %r15
jne 0x1961
cmpl $0x1, %eax
movq 0x38(%rsp), %r12
je 0x1a93
leal -0x2(%rax), %ecx
decl %eax
cmpl %ecx, %eax
jle 0x1a86
movq (%rsp), %rdx
movq 0x8(%rdx), %rdx
movq 0x10(%rsp), %rsi
movq 0x8(%rsi), %rsi
movq (%rsi,%rcx,8), %rsi
movsd (%rdx,%rcx,8), %xmm0
movl %eax, %edi
movl %edi, %r8d
movsd (%rsi,%r8,8), %xmm1
mulsd %xmm0, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rdx,%rcx,8)
decl %edi
cmpl %ecx, %edi
jg 0x1a6a
movq %rcx, %rdx
decq %rdx
testl %ecx, %ecx
movq %rdx, %rcx
jg 0x1a4a
cmpl $0x0, 0x18(%rsp)
jle 0x1b15
xorl %ebx, %ebx
movq 0x8(%r13), %rax
cmpl $0x0, (%rax,%rbx,4)
je 0x1b0b
xorl %r15d, %r15d
xorl %ebp, %ebp
movq 0x8(%r13), %rax
cmpl $0x0, (%rax,%r15,4)
je 0x1b01
movq 0x8(%r12), %rax
movsd (%rax,%rbx,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movl %r14d, %edi
movl %ebx, %esi
movl %r15d, %edx
movq %r13, %rcx
movq 0x28(%rsp), %r8
callq 0xcac2
movq (%rsp), %rax
movq 0x8(%rax), %rax
movslq %ebp, %rbp
mulsd (%rax,%rbp,8), %xmm0
movsd 0x18(%rsp), %xmm1
subsd %xmm0, %xmm1
movq 0x8(%r12), %rax
movsd %xmm1, (%rax,%rbx,8)
incl %ebp
incq %r15
cmpq %r15, 0x8(%rsp)
jne 0x1aab
incq %rbx
cmpq 0x8(%rsp), %rbx
jne 0x1a9c
cmpl $0x2, %r14d
movq 0x30(%rsp), %r10
movq 0x48(%rsp), %r11
jl 0x1b7a
leal -0x2(%r14), %eax
leal -0x1(%r14), %ecx
leaq -0x1(%r10), %rdx
cmpq %rcx, %rax
jge 0x1b6e
movq 0x8(%r12), %rsi
movq 0x28(%rsp), %rdi
movq 0x8(%rdi), %rdi
movq (%rdi,%rax,8), %rdi
movsd (%rsi,%rax,8), %xmm0
movq %rdx, %r8
movl %r8d, %r9d
movsd (%rdi,%r9,8), %xmm1
mulsd %xmm0, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rsi,%rax,8)
decq %r8
cmpq %r8, %rax
jl 0x1b50
leaq -0x1(%rax), %rsi
testq %rax, %rax
movq %rsi, %rax
jg 0x1b31
cmpl $0x2, 0x50(%rsp)
jl 0x1bd9
pushq $0x1
popq %rax
leal 0x1(%r14), %ecx
movslq %ecx, %rcx
shlq $0x3, %r10
movl %r14d, %edx
cmpl %r14d, %edx
jl 0x1bcc
movq 0x8(%r12), %rdx
movq 0x28(%rsp), %rsi
movq 0x8(%rsi), %rsi
movsd (%rdx,%rcx,8), %xmm0
movq (%rsi,%rcx,8), %rsi
addq %r10, %rsi
xorl %edi, %edi
movsd (%rsi,%rdi,8), %xmm1
mulsd %xmm0, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, (%rdx,%rcx,8)
incq %rdi
cmpl %edi, %eax
jne 0x1bb3
movl %ecx, %edx
incq %rcx
incq %rax
cmpq %r11, %rcx
jl 0x1b92
pushq $0x22
popq %rsi
xorl %edi, %edi
addq $0x68, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xb980