// File: /workspace/github_cpp_makefile_star2_O1/ashlall[P]StreamStats/ChiSquareExperiments/VarySizeTwoSampleBoth.cpp
// Function: main
int main(int argc, char* argv[])
{
  if (argc < 12)
    {
      cout << "usage: VarySizeTwoSample lower-size upper-size num-stream(>0) [N|U|P|E] location1 scale1 location2 scale2 [0|1] memory-percent num-buckets" << endl;
      throw ParameterError();
    }

  double lower = atof (argv[1]);
  double upper = atof (argv[2]);
  int data_repeats = atoi (argv[3]);
  char distribution_type = argv[4][0];
  double location1 = atof (argv[5]);
  double scale1 = atof (argv[6]);
  double location2 = atof (argv[7]);
  double scale2 = atof (argv[8]);
  int all_quantiles = atoi (argv[9]);
  double memory_percent = atof (argv[10]);
  int num_buckets = atoi (argv[11]);
  int seed1 = 1;
  int seed2 = 2;

  // ensures that the parameters will not create error
  if (data_repeats <= 0)
    {
      cout << "The number of streams must be greater than zero.\n";
      throw ParameterError();
    }
  if (distribution_type != 'N' && distribution_type != 'U' && distribution_type != 'P' && distribution_type != 'E')
    {
      cout << "The distribution is either N, U, P, or E.\n";
      throw ParameterError();
    }
  if (lower <= 0)
    {
      cout << "The lower stream size must be greater than 0.\n";
      throw ParameterError();
    }
  if (memory_percent <= 0)
    {
      cout << "The memory percent must be greater than 0.\n";
      throw ParameterError();
    }
  if (num_buckets <= 0)
    {
      cout << "The number of bins must be greater than 0.\n";
      throw ParameterError();
    }

  // finds the number of different sizes the experiment will run on
  int stream_size;
  int num_sizes = 0, size = lower;
  while (size <= upper)
    {
      num_sizes++;
      size *= 10;
    }

  // creates the arrays that will hold the calculated statistics 
  double actual_values[data_repeats][num_sizes];
  double GK_values[data_repeats][num_sizes];
  double QD_values[data_repeats][num_sizes];
  double RS_values[data_repeats][num_sizes];
  int sizes[num_sizes];

  // creates and intializes the log file
  ofstream data_file;
  char str[150];
  name_file(str, argv, 0);
  data_file.open(str);
  data_file << "Distribution: " << distribution_type << endl;
  data_file << "location_1 = " << location1 << " scale_1 = " << scale1 << endl;
  data_file << "location_2 = " << location2 << " scale_2 = " << scale2 << endl;
  data_file << "number of bins = " << num_buckets << endl;
  data_file << "memory percent = " << memory_percent << endl;

  for (int i = 0; i < data_repeats; i++) // runs all test data_repeats times
    {
      int j = 0;
      stream_size = lower;
      data_file << "Pass " << i+1 << ":" << endl;
      while (stream_size <= upper)  // runs tests for every stream size
	{ 
	  int sample_size = memory_percent * stream_size;
	  if (!i)
	    sizes[j] = stream_size;
	  data_file << "stream_size = " << stream_size << endl;

	  // generates the data based on parameters
	  DataGenerator data1(distribution_type, stream_size, seed1, location1, scale1);
	  double *stream1 = data1.get_stream();
	  DataGenerator data2(distribution_type, stream_size, seed2, location2, scale2);
	  double *stream2 = data2.get_stream();

	  // computes GK estimate
	  ChiSquareContinuous GK_sketch1(sample_size, 1);
	  for (int k = 0; k < stream_size; k++)
	    GK_sketch1.insert(stream1[k]);
	  ChiSquareContinuous GK_sketch2(sample_size, 1);
	  for (int k = 0; k < stream_size; k++)
	    GK_sketch2.insert(stream2[k]);
	  double GK_stat = GK_sketch1.two_sample_statistic(GK_sketch2, num_buckets);
	  //double GK_stat = GK_sketch2.two_sample_statistic(GK_sketch1, num_buckets);
	  GK_values[i][j] = GK_stat;
	  data_file << "GK = " << GK_stat << endl;

	  if (all_quantiles)
	    {
	      // computes QDigest estimate
	      ChiSquareContinuous QD_sketch1(sample_size, 2);
	      for (int k = 0; k < stream_size; k++)
		QD_sketch1.insert(stream1[k]);
	      ChiSquareContinuous QD_sketch2(sample_size, 2);
	      for (int k = 0; k < stream_size; k++)
		QD_sketch2.insert(stream2[k]);
	      double QD_stat = QD_sketch1.two_sample_statistic(QD_sketch2, num_buckets);
	      QD_values[i][j] = QD_stat;
	      data_file << "QD = " << QD_stat << endl;

	      // computes Reservoir Sampling estimate
	      ChiSquareContinuous RS_sketch1(sample_size, 3);
	      for (int k = 0; k < stream_size; k++)
		RS_sketch1.insert(stream1[k]);
	      ChiSquareContinuous RS_sketch2(sample_size, 3);
	      for (int k = 0; k < stream_size; k++)
		RS_sketch2.insert(stream2[k]);
	      double RS_stat = RS_sketch1.two_sample_statistic(RS_sketch2, num_buckets);
	      RS_values[i][j]= RS_stat;
	      data_file << "RS = " << RS_stat<< endl;
	    }
       
	  // computes actual statistic
	  double *upper_interval = GK_sketch1.get_upper();
	  double *lower_interval = GK_sketch1.get_lower();
	  double actual = data2.get_stat_two_sample(data1, num_buckets, upper_interval, lower_interval);
	  actual_values[i][j] = actual;
	  data_file << "Real = " << actual << endl;

	  stream_size *= 10;
	  j++;
	}     
      seed1++;
      seed2++;
    }
  data_file.close();

  // creates pvalue file
  ofstream data2_file;
  char st[150];
  name_file(st, argv, 3);
  data2_file.open(st);

  // creates table file
  name_file(str, argv, 1);
  data_file.open(str);

  int deg_freedom = num_buckets;
  for (int i = 0; i < num_sizes; i++)
  {
    data_file << sizes[i] << "\t";
    double error = 0;

    // adds pvalue error from the GK sketch to the table file
    for (int j = 0; j < data_repeats; j++)
      error += abs(pochisq(GK_values[j][i], deg_freedom) - pochisq(actual_values[j][i], deg_freedom));
    data_file << error / data_repeats;

    if (all_quantiles)
    {
      // adds pvalue error from the QDigest sketch to the table file
      error = 0;
      for (int j = 0; j < data_repeats; j++)
        error += abs(pochisq(QD_values[j][i], deg_freedom) - pochisq(actual_values[j][i], deg_freedom));
      data_file << "\t" << error / data_repeats << "\t";

      // adds pvalue error from the Reservoir Samplig sketch to the table file
      error = 0;
      for (int j = 0; j < data_repeats; j++)
        error += abs(pochisq(RS_values[j][i], deg_freedom) - pochisq(actual_values[j][i], deg_freedom));
      data_file << error / data_repeats;
    }
    data_file << endl;

    // adds all pvalues to the pvalue file
    data2_file << "stream_size = " << sizes[i] << endl;
    for (int j = 0; j < data_repeats; j++)
      {
        data2_file << pochisq(actual_values[j][i], deg_freedom) << " actual" << endl;
        data2_file << pochisq(GK_values[j][i], deg_freedom) << " GK" << endl;
        if (all_quantiles)
        {
          data2_file << pochisq(QD_values[j][i], deg_freedom) << " QD" << endl;
          data2_file << pochisq(RS_values[j][i], deg_freedom) << " RS" << endl;
        }
      }
  }
  data2_file.close();
  data_file.close();

  // creates extra file
  name_file(str, argv, 2);
  data_file.open(str);

  // adds actual statistics to the extra file
  for (int i = 0; i < num_sizes; i++)
    {
      for (int j = 0; j < data_repeats; j++)
	data_file << actual_values[j][i] << "\t";
    }
  data_file << endl;

  // adds GK statistics to the extra file
  for (int i = 0; i < num_sizes; i++)
    {
      for(int j = 0; j <data_repeats; j++)
	data_file << GK_values[j][i] << "\t";
    }
  data_file << endl;
  if (all_quantiles)
    {
      // adds QDigest statistics to the extra file
      for (int i = 0; i < num_sizes; i++)
	{
	  for(int j = 0; j <data_repeats; j++)
	    data_file << QD_values[j][i] << "\t";
	}
      data_file << endl;
      
      // adds Reservoir Sampling statistics to the extra file
      for (int i = 0; i < num_sizes; i++)
	{
	  for(int j = 0; j <data_repeats; j++)
	    data_file << RS_values[j][i] << "\t";
	}
      data_file << endl;
    }

  data_file << endl;
  data_file.close();
  return 0;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x718, %rsp            # imm = 0x718
cmpl $0xb, %edi
jle 0x65f2
movq %rsi, %r14
movq 0x8(%rsi), %rdi
callq 0x2240
movapd %xmm0, -0x60(%rbp)
movq 0x10(%r14), %rdi
callq 0x2240
movsd %xmm0, -0x98(%rbp)
movq 0x18(%r14), %rdi
callq 0x22b0
movl %eax, %ebx
movq 0x20(%r14), %rax
movq 0x28(%r14), %rdi
movb (%rax), %r15b
callq 0x2240
movsd %xmm0, -0x128(%rbp)
movq 0x30(%r14), %rdi
callq 0x2240
movsd %xmm0, -0x120(%rbp)
movq 0x38(%r14), %rdi
callq 0x2240
movsd %xmm0, -0x118(%rbp)
movq 0x40(%r14), %rdi
callq 0x2240
movsd %xmm0, -0x110(%rbp)
movq 0x48(%r14), %rdi
callq 0x22b0
movl %eax, -0x64(%rbp)
movq 0x50(%r14), %rdi
callq 0x2240
movsd %xmm0, -0xe8(%rbp)
movq 0x58(%r14), %rdi
callq 0x22b0
movl %ebx, -0x44(%rbp)
testl %ebx, %ebx
jle 0x660f
movzbl %r15b, %edx
addl $-0x45, %edx
cmpl $0x10, %edx
ja 0x664f
movl $0x10a01, %ecx          # imm = 0x10A01
btl %edx, %ecx
jae 0x664f
xorpd %xmm0, %xmm0
movapd -0x60(%rbp), %xmm1
ucomisd %xmm1, %xmm0
jae 0x661f
ucomisd -0xe8(%rbp), %xmm0
jae 0x662f
testl %eax, %eax
jle 0x663f
movb %r15b, -0x29(%rbp)
movl %eax, -0x38(%rbp)
movq %r14, -0xe0(%rbp)
cvttsd2si %xmm1, %eax
movl %eax, -0xa4(%rbp)
cvttpd2dq %xmm1, %xmm0
cvtdq2pd %xmm0, %xmm0
xorl %eax, %eax
movsd -0x98(%rbp), %xmm1
ucomisd %xmm0, %xmm1
jb 0x53af
xorl %eax, %eax
movl -0xa4(%rbp), %ecx
incl %eax
addl %ecx, %ecx
leal (%rcx,%rcx,4), %ecx
xorps %xmm0, %xmm0
cvtsi2sd %ecx, %xmm0
ucomisd %xmm0, %xmm1
jae 0x539b
movl -0x44(%rbp), %ecx
movl %eax, %edx
imulq %rdx, %rcx
leaq 0xf(,%rcx,8), %rax
andq $-0x10, %rax
movq %rsp, %rcx
subq %rax, %rcx
movq %rcx, -0xb8(%rbp)
movq %rcx, %rsp
movq %rsp, %r13
subq %rax, %r13
movq %r13, %rsp
movq %rsp, %r12
subq %rax, %r12
movq %r12, %rsp
movq %rsp, %r15
subq %rax, %r15
movq %r15, %rsp
movq %rsp, %rcx
movq %rdx, -0x40(%rbp)
leaq 0xf(,%rdx,4), %rax
andq $-0x10, %rax
subq %rax, %rcx
movq %rcx, -0x108(%rbp)
movq %rcx, %rsp
leaq -0x400(%rbp), %rbx
movq %rbx, %rdi
callq 0x2310
leaq -0x740(%rbp), %r14
movq %r14, %rdi
movq -0xe0(%rbp), %rsi
xorl %edx, %edx
callq 0x67b2
movq %rbx, %rdi
movq %r14, %rsi
movl $0x10, %edx
callq 0x2060
leaq 0x4013(%rip), %rsi      # 0x9460
leaq -0x400(%rbp), %rdi
movl $0xe, %edx
callq 0x2220
leaq -0x600(%rbp), %rsi
movb -0x29(%rbp), %al
movb %al, (%rsi)
leaq -0x400(%rbp), %rdi
movl $0x1, %edx
callq 0x2220
movq %rax, %rbx
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq 0x3fc3(%rip), %rsi      # 0x946f
leaq -0x400(%rbp), %rdi
movl $0xd, %edx
callq 0x2220
leaq -0x400(%rbp), %rdi
movsd -0x128(%rbp), %xmm0
callq 0x22c0
movq %rax, %rbx
leaq 0x3fa2(%rip), %rsi      # 0x947d
movl $0xb, %edx
movq %rax, %rdi
callq 0x2220
movq %rbx, %rdi
movsd -0x120(%rbp), %xmm0
callq 0x22c0
movq %rax, %rbx
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq 0x3f60(%rip), %rsi      # 0x9489
leaq -0x400(%rbp), %rdi
movl $0xd, %edx
callq 0x2220
leaq -0x400(%rbp), %rdi
movsd -0x118(%rbp), %xmm0
callq 0x22c0
movq %rax, %rbx
leaq 0x3f3f(%rip), %rsi      # 0x9497
movl $0xb, %edx
movq %rax, %rdi
callq 0x2220
movq %rbx, %rdi
movsd -0x110(%rbp), %xmm0
callq 0x22c0
movq %rax, %rbx
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq 0x3efd(%rip), %rsi      # 0x94a3
leaq -0x400(%rbp), %rdi
movl $0x11, %edx
callq 0x2220
leaq -0x400(%rbp), %rdi
movl -0x38(%rbp), %esi
callq 0x2330
movq %rax, %rbx
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq 0x3ebe(%rip), %rsi      # 0x94b5
leaq -0x400(%rbp), %rdi
movl $0x11, %edx
callq 0x2220
leaq -0x400(%rbp), %rdi
movsd -0xe8(%rbp), %xmm0
callq 0x22c0
movq %rax, %rbx
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
xorps %xmm0, %xmm0
cvtsi2sdl -0xa4(%rbp), %xmm0
movsd %xmm0, -0x100(%rbp)
movl -0x44(%rbp), %eax
cmpl $0x2, %eax
movl $0x1, -0x90(%rbp)
movl $0x1, %ecx
cmovgel %eax, %ecx
movq %rcx, -0xc0(%rbp)
movl $0x2, -0xb0(%rbp)
movq $0x0, -0x88(%rbp)
leaq -0x400(%rbp), %rbx
movq %r15, -0xc8(%rbp)
movq %r12, -0xd0(%rbp)
movq %r13, -0xd8(%rbp)
movl $0x5, %edx
movq %rbx, %rdi
leaq 0x3e0f(%rip), %rsi      # 0x94c7
callq 0x2220
movq -0x88(%rbp), %rax
leaq 0x1(%rax), %rsi
leaq -0x400(%rbp), %rdi
movq %rsi, -0x130(%rbp)
callq 0x2330
movq %rax, %rbx
movl $0x1, %edx
movq %rax, %rdi
leaq 0x3de0(%rip), %rsi      # 0x94cd
callq 0x2220
movq (%rbx), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
movsd -0x98(%rbp), %xmm0
ucomisd -0x100(%rbp), %xmm0
leaq -0x400(%rbp), %rbx
jb 0x5ca0
movq -0x88(%rbp), %rax
imulq -0x40(%rbp), %rax
leaq (,%rax,8), %rcx
addq %r13, %rcx
movq %rcx, -0xf8(%rbp)
leaq (%r12,%rax,8), %rcx
movq %rcx, -0x140(%rbp)
leaq (%r15,%rax,8), %rcx
movq %rcx, -0x138(%rbp)
movq -0xb8(%rbp), %rcx
leaq (%rcx,%rax,8), %rax
movq %rax, -0xf0(%rbp)
movq $0x0, -0x50(%rbp)
movsd -0x100(%rbp), %xmm0
movsd %xmm0, -0x80(%rbp)
movl -0xa4(%rbp), %eax
movq %rax, -0x60(%rbp)
cmpq $0x0, -0x88(%rbp)
movq -0x60(%rbp), %r12
jne 0x57b8
movq -0x108(%rbp), %rax
movq -0x50(%rbp), %rcx
movl %r12d, (%rax,%rcx,4)
movl $0xe, %edx
movq %rbx, %rdi
leaq 0x3d08(%rip), %rsi      # 0x94cf
callq 0x2220
movb -0x29(%rbp), %bl
leaq -0x400(%rbp), %rdi
movl %r12d, %esi
callq 0x2330
movq %rax, %r14
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %r14, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %r14, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
movsbl %bl, %r15d
leaq -0x160(%rbp), %rdi
movl %r15d, %esi
movl %r12d, %edx
movl -0x90(%rbp), %ecx
movsd -0x128(%rbp), %xmm0
movsd -0x120(%rbp), %xmm1
callq 0x250c
movq -0x160(%rbp), %r14
leaq -0x150(%rbp), %rdi
movl %r15d, %esi
movl %r12d, %edx
movl -0xb0(%rbp), %ecx
movsd -0x118(%rbp), %xmm0
movsd -0x110(%rbp), %xmm1
callq 0x250c
movsd -0x80(%rbp), %xmm0
mulsd -0xe8(%rbp), %xmm0
cvttsd2si %xmm0, %eax
movq -0x150(%rbp), %r13
xorps %xmm0, %xmm0
cvtsi2sd %eax, %xmm0
leaq -0x600(%rbp), %rdi
movsd %xmm0, -0x80(%rbp)
movl $0x1, %esi
callq 0x46a4
movq -0x60(%rbp), %rax
movl %eax, %r12d
testl %eax, %eax
jle 0x58c2
movq -0x5e0(%rbp), %r15
xorl %ebx, %ebx
movsd (%r14,%rbx,8), %xmm0
movq (%r15), %rax
movq %r15, %rdi
callq *(%rax)
incq %rbx
cmpq %rbx, %r12
jne 0x58ac
leaq -0x6a0(%rbp), %rdi
movsd -0x80(%rbp), %xmm0
movl $0x1, %esi
callq 0x46a4
cmpl $0x0, -0x60(%rbp)
jle 0x58fe
movq -0x680(%rbp), %r15
xorl %ebx, %ebx
movsd (%r13,%rbx,8), %xmm0
movq (%r15), %rax
movq %r15, %rdi
callq *(%rax)
incq %rbx
cmpq %rbx, %r12
jne 0x58e7
leaq -0x600(%rbp), %rdi
leaq -0x6a0(%rbp), %rsi
movl -0x38(%rbp), %edx
callq 0x4f80
leaq -0x400(%rbp), %rbx
movq -0xf8(%rbp), %rax
movq -0x50(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl $0x5, %edx
movq %rbx, %rdi
leaq 0x3ba4(%rip), %rsi      # 0x94de
movsd %xmm0, -0xa0(%rbp)
callq 0x2220
movsd -0xa0(%rbp), %xmm0
movq %rbx, %rdi
callq 0x22c0
movq %rax, %r15
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %r15, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %r15, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
cmpl $0x0, -0x64(%rbp)
je 0x5ba9
leaq -0x200(%rbp), %rdi
movsd -0x80(%rbp), %xmm0
movl $0x2, %esi
callq 0x46a4
cmpl $0x0, -0x60(%rbp)
jle 0x59c6
movq -0x1e0(%rbp), %r15
xorl %ebx, %ebx
movsd (%r14,%rbx,8), %xmm0
movq (%r15), %rax
movq %r15, %rdi
callq *(%rax)
incq %rbx
cmpq %rbx, %r12
jne 0x59b0
leaq -0x1d8(%rbp), %rdi
movsd -0x80(%rbp), %xmm0
movl $0x2, %esi
callq 0x46a4
cmpl $0x0, -0x60(%rbp)
jle 0x5a02
movq -0x1b8(%rbp), %r15
xorl %ebx, %ebx
movsd (%r13,%rbx,8), %xmm0
movq (%r15), %rax
movq %r15, %rdi
callq *(%rax)
incq %rbx
cmpq %rbx, %r12
jne 0x59eb
leaq -0x200(%rbp), %rdi
leaq -0x1d8(%rbp), %rsi
movl -0x38(%rbp), %edx
callq 0x4f80
leaq -0x400(%rbp), %rbx
movq -0x140(%rbp), %rax
movq -0x50(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl $0x5, %edx
movq %rbx, %rdi
leaq 0x3aa6(%rip), %rsi      # 0x94e4
movsd %xmm0, -0xa0(%rbp)
callq 0x2220
movsd -0xa0(%rbp), %xmm0
movq %rbx, %rdi
callq 0x22c0
movq %rax, %r15
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %r15, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %r15, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq -0x1b0(%rbp), %rdi
movsd -0x80(%rbp), %xmm0
movl $0x3, %esi
callq 0x46a4
cmpl $0x0, -0x60(%rbp)
jle 0x5ac0
movq -0x190(%rbp), %r15
xorl %ebx, %ebx
movsd (%r14,%rbx,8), %xmm0
movq (%r15), %rax
movq %r15, %rdi
callq *(%rax)
incq %rbx
cmpq %rbx, %r12
jne 0x5aaa
leaq -0x188(%rbp), %rdi
movsd -0x80(%rbp), %xmm0
movl $0x3, %esi
callq 0x46a4
cmpl $0x0, -0x60(%rbp)
jle 0x5afc
movq -0x168(%rbp), %r15
xorl %ebx, %ebx
movsd (%r13,%rbx,8), %xmm0
movq (%r15), %rax
movq %r15, %rdi
callq *(%rax)
incq %rbx
cmpq %rbx, %r12
jne 0x5ae5
leaq -0x1b0(%rbp), %rdi
leaq -0x188(%rbp), %rsi
movl -0x38(%rbp), %edx
callq 0x4f80
leaq -0x400(%rbp), %rbx
movq -0x138(%rbp), %rax
movq -0x50(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl $0x5, %edx
movq %rbx, %rdi
leaq 0x39b2(%rip), %rsi      # 0x94ea
movsd %xmm0, -0x80(%rbp)
callq 0x2220
movsd -0x80(%rbp), %xmm0
movq %rbx, %rdi
callq 0x22c0
movq %rax, %r15
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %r15, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %r15, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq -0x188(%rbp), %rdi
callq 0x4844
leaq -0x1b0(%rbp), %rdi
callq 0x4844
leaq -0x1d8(%rbp), %rdi
callq 0x4844
leaq -0x200(%rbp), %rdi
callq 0x4844
movq -0x5f0(%rbp), %rcx
movq -0x5e8(%rbp), %r8
leaq -0x150(%rbp), %rdi
leaq -0x160(%rbp), %rsi
movl -0x38(%rbp), %edx
callq 0x27e6
movq -0xf0(%rbp), %rax
movq -0x50(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl $0x7, %edx
movq %rbx, %rdi
leaq 0x3904(%rip), %rsi      # 0x94f0
movsd %xmm0, -0x80(%rbp)
callq 0x2220
movsd -0x80(%rbp), %xmm0
movq %rbx, %rdi
callq 0x22c0
movq %rax, %r15
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %r15, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %r15, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq -0x6a0(%rbp), %rdi
callq 0x4844
leaq -0x600(%rbp), %rdi
callq 0x4844
testq %r13, %r13
je 0x5c52
movq %r13, %rdi
callq 0x2250
testq %r14, %r14
je 0x5c5f
movq %r14, %rdi
callq 0x2250
movq -0x60(%rbp), %rax
addl %eax, %eax
leal (%rax,%rax,4), %eax
incq -0x50(%rbp)
movq %rax, -0x60(%rbp)
cvtsi2sd %eax, %xmm1
movsd -0x98(%rbp), %xmm0
movsd %xmm1, -0x80(%rbp)
ucomisd %xmm1, %xmm0
movq -0xc8(%rbp), %r15
movq -0xd0(%rbp), %r12
movq -0xd8(%rbp), %r13
jae 0x579b
incl -0x90(%rbp)
incl -0xb0(%rbp)
movq -0x130(%rbp), %rcx
movq %rcx, -0x88(%rbp)
cmpq -0xc0(%rbp), %rcx
jne 0x56a9
leaq -0x400(%rbp), %rdi
callq 0x2270
leaq -0x600(%rbp), %rdi
callq 0x2310
leaq -0x6a0(%rbp), %rbx
movq %rbx, %rdi
movq -0xe0(%rbp), %r14
movq %r14, %rsi
movl $0x3, %edx
callq 0x67b2
leaq -0x600(%rbp), %rdi
movq %rbx, %rsi
movl $0x10, %edx
callq 0x2060
leaq -0x740(%rbp), %rbx
movq %rbx, %rdi
movq %r14, %rsi
movl $0x1, %edx
callq 0x67b2
leaq -0x400(%rbp), %rdi
movq %rbx, %rsi
movl $0x10, %edx
callq 0x2060
cmpq $0x0, -0x40(%rbp)
je 0x61f2
xorps %xmm0, %xmm0
cvtsi2sdl -0x44(%rbp), %xmm0
movsd %xmm0, -0xb0(%rbp)
movq -0x40(%rbp), %rax
leaq (,%rax,8), %r12
xorl %ecx, %ecx
movq -0xc8(%rbp), %rax
movq -0xd0(%rbp), %rsi
movq -0xd8(%rbp), %rdx
movq -0xb8(%rbp), %r14
movq %r12, -0x60(%rbp)
movq %rsi, -0x88(%rbp)
movq %r14, -0x90(%rbp)
movq %rdx, -0xa0(%rbp)
movq %rax, -0x98(%rbp)
movq -0x108(%rbp), %rax
movq %rcx, -0xf0(%rbp)
movl (%rax,%rcx,4), %esi
leaq -0x400(%rbp), %rdi
movl %esi, -0xf8(%rbp)
callq 0x2330
movl $0x1, %edx
movq %rax, %rdi
leaq 0x3725(%rip), %rsi      # 0x94f8
callq 0x2220
xorpd %xmm1, %xmm1
movq -0xc0(%rbp), %r15
xorl %r13d, %r13d
movl -0x38(%rbp), %ebx
movq -0xa0(%rbp), %r14
movq -0x90(%rbp), %r12
movsd %xmm1, -0x50(%rbp)
movsd (%r14,%r13), %xmm0
movl %ebx, %edi
callq 0x421e
movaps %xmm0, -0x80(%rbp)
movsd (%r12,%r13), %xmm0
movl %ebx, %edi
callq 0x421e
movsd -0x50(%rbp), %xmm1
movapd -0x80(%rbp), %xmm2
subsd %xmm0, %xmm2
andpd 0x3370(%rip), %xmm2     # 0x91a0
addsd %xmm2, %xmm1
addq -0x60(%rbp), %r13
decq %r15
jne 0x5df7
divsd -0xb0(%rbp), %xmm1
leaq -0x400(%rbp), %rdi
movapd %xmm1, %xmm0
callq 0x22c0
movq -0x60(%rbp), %r12
cmpl $0x0, -0x64(%rbp)
je 0x5f93
xorpd %xmm1, %xmm1
movq -0xc0(%rbp), %r15
xorl %r13d, %r13d
movl -0x38(%rbp), %ebx
movq -0x90(%rbp), %r14
movq -0x88(%rbp), %r12
movsd %xmm1, -0x50(%rbp)
movsd (%r12,%r13), %xmm0
movl %ebx, %edi
callq 0x421e
movaps %xmm0, -0x80(%rbp)
movsd (%r14,%r13), %xmm0
movl %ebx, %edi
callq 0x421e
movsd -0x50(%rbp), %xmm1
movapd -0x80(%rbp), %xmm2
subsd %xmm0, %xmm2
andpd 0x32e5(%rip), %xmm2     # 0x91a0
addsd %xmm2, %xmm1
addq -0x60(%rbp), %r13
decq %r15
jne 0x5e82
movsd %xmm1, -0x50(%rbp)
movl $0x1, %edx
leaq -0x400(%rbp), %rdi
leaq 0x3618(%rip), %rsi      # 0x94f8
callq 0x2220
movsd -0x50(%rbp), %xmm0
divsd -0xb0(%rbp), %xmm0
leaq -0x400(%rbp), %rdi
callq 0x22c0
movl $0x1, %edx
movq %rax, %rdi
leaq 0x35eb(%rip), %rsi      # 0x94f8
callq 0x2220
xorpd %xmm1, %xmm1
movq -0xc0(%rbp), %r15
xorl %r13d, %r13d
movl -0x38(%rbp), %ebx
movq -0x98(%rbp), %r12
movq -0x90(%rbp), %r14
movsd %xmm1, -0x50(%rbp)
movsd (%r12,%r13), %xmm0
movl %ebx, %edi
callq 0x421e
movaps %xmm0, -0x80(%rbp)
movsd (%r14,%r13), %xmm0
movl %ebx, %edi
callq 0x421e
movsd -0x50(%rbp), %xmm1
movapd -0x80(%rbp), %xmm2
subsd %xmm0, %xmm2
andpd 0x3236(%rip), %xmm2     # 0x91a0
addsd %xmm2, %xmm1
addq -0x60(%rbp), %r13
decq %r15
jne 0x5f31
divsd -0xb0(%rbp), %xmm1
leaq -0x400(%rbp), %rdi
movapd %xmm1, %xmm0
callq 0x22c0
movq -0x60(%rbp), %r12
movq -0x400(%rbp), %rax
movq -0x18(%rax), %rdi
leaq -0x400(%rbp), %rax
addq %rax, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
leaq -0x400(%rbp), %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
movl $0xe, %edx
leaq -0x600(%rbp), %rdi
leaq 0x34f3(%rip), %rsi      # 0x94cf
callq 0x2220
leaq -0x600(%rbp), %rdi
movl -0xf8(%rbp), %esi
callq 0x2330
movq %rax, %rbx
movq (%rax), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
movq -0xc0(%rbp), %r13
xorl %r15d, %r15d
movq -0x90(%rbp), %r14
movsd (%r14,%r15), %xmm0
movl -0x38(%rbp), %edi
callq 0x421e
leaq -0x600(%rbp), %rdi
callq 0x22c0
movq %rax, %rbx
movl $0x7, %edx
movq %rax, %rdi
leaq 0x34a0(%rip), %rsi      # 0x94fa
callq 0x2220
movq (%rbx), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
movq -0xa0(%rbp), %rax
movsd (%rax,%r15), %xmm0
movl -0x38(%rbp), %edi
callq 0x421e
leaq -0x600(%rbp), %rdi
callq 0x22c0
movq %rax, %rbx
movl $0x3, %edx
movq %rax, %rdi
leaq 0x3449(%rip), %rsi      # 0x9502
callq 0x2220
movq (%rbx), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
cmpl $0x0, -0x64(%rbp)
je 0x61ad
movq -0x88(%rbp), %rax
movsd (%rax,%r15), %xmm0
movl -0x38(%rbp), %edi
callq 0x421e
leaq -0x600(%rbp), %rdi
callq 0x22c0
movq %rax, %rbx
movl $0x3, %edx
movq %rax, %rdi
leaq 0x33e4(%rip), %rsi      # 0x9506
callq 0x2220
movq (%rbx), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
movq -0x98(%rbp), %rax
movsd (%rax,%r15), %xmm0
movl -0x38(%rbp), %edi
callq 0x421e
leaq -0x600(%rbp), %rdi
callq 0x22c0
movq %rax, %rbx
movl $0x3, %edx
movq %rax, %rdi
leaq 0x3389(%rip), %rsi      # 0x950a
callq 0x2220
movq (%rbx), %rax
movq -0x18(%rax), %rdi
addq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
addq %r12, %r15
decq %r13
jne 0x602e
movq -0xf0(%rbp), %rcx
incq %rcx
addq $0x8, %r14
movq -0xa0(%rbp), %rdx
addq $0x8, %rdx
movq -0x88(%rbp), %rsi
addq $0x8, %rsi
movq -0x98(%rbp), %rax
addq $0x8, %rax
cmpq -0x40(%rbp), %rcx
jne 0x5d85
leaq -0x600(%rbp), %rdi
callq 0x2270
leaq -0x400(%rbp), %rdi
callq 0x2270
leaq -0x740(%rbp), %r15
movq %r15, %rdi
movq -0xe0(%rbp), %rsi
movl $0x2, %edx
callq 0x67b2
leaq -0x400(%rbp), %rdi
movq %r15, %rsi
movl $0x10, %edx
callq 0x2060
cmpq $0x0, -0x40(%rbp)
movq -0xb8(%rbp), %r14
je 0x62c1
movl -0x44(%rbp), %eax
cmpl $0x2, %eax
movl $0x1, %ecx
cmovgel %eax, %ecx
movq %rcx, -0x38(%rbp)
movq -0x40(%rbp), %rax
leaq (,%rax,8), %r13
xorl %eax, %eax
leaq -0x400(%rbp), %r15
leaq 0x3283(%rip), %r12      # 0x94f8
movq %rax, -0x60(%rbp)
movq -0x38(%rbp), %rbx
movq %r14, -0xb8(%rbp)
movsd (%r14), %xmm0
movq %r15, %rdi
callq 0x22c0
movl $0x1, %edx
movq %rax, %rdi
movq %r12, %rsi
callq 0x2220
addq %r13, %r14
decq %rbx
jne 0x6284
movq -0x60(%rbp), %rax
incq %rax
movq -0xb8(%rbp), %r14
addq $0x8, %r14
cmpq -0x40(%rbp), %rax
jne 0x6275
movq -0x400(%rbp), %rax
leaq -0x400(%rbp), %rdi
addq -0x18(%rax), %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
leaq -0x400(%rbp), %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
cmpq $0x0, -0x40(%rbp)
movq -0xd8(%rbp), %r14
je 0x637c
movl -0x44(%rbp), %eax
cmpl $0x2, %eax
movl $0x1, %ecx
cmovgel %eax, %ecx
movq %rcx, -0x38(%rbp)
movq -0x40(%rbp), %rax
leaq (,%rax,8), %rbx
xorl %eax, %eax
leaq -0x400(%rbp), %r15
leaq 0x31c8(%rip), %r12      # 0x94f8
movq %rax, -0x60(%rbp)
movq -0x38(%rbp), %r13
movq %r14, -0xd8(%rbp)
movsd (%r14), %xmm0
movq %r15, %rdi
callq 0x22c0
movl $0x1, %edx
movq %rax, %rdi
movq %r12, %rsi
callq 0x2220
addq %rbx, %r14
decq %r13
jne 0x633f
movq -0x60(%rbp), %rax
incq %rax
movq -0xd8(%rbp), %r14
addq $0x8, %r14
cmpq -0x40(%rbp), %rax
jne 0x6330
movq -0x400(%rbp), %rax
leaq -0x400(%rbp), %rdi
addq -0x18(%rax), %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
leaq -0x400(%rbp), %rbx
movq %rbx, %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
cmpl $0x0, -0x64(%rbp)
je 0x653d
cmpq $0x0, -0x40(%rbp)
movq -0xd0(%rbp), %r13
je 0x6445
movl -0x44(%rbp), %eax
cmpl $0x2, %eax
movl $0x1, %ecx
cmovgel %eax, %ecx
movq %rcx, -0x38(%rbp)
movq -0x40(%rbp), %rax
leaq (,%rax,8), %r14
xorl %eax, %eax
leaq -0x400(%rbp), %r15
leaq 0x3100(%rip), %r12      # 0x94f8
movq %rax, -0x60(%rbp)
movq -0x38(%rbp), %rbx
movq %r13, -0xd0(%rbp)
movsd (%r13), %xmm0
movq %r15, %rdi
callq 0x22c0
movl $0x1, %edx
movq %rax, %rdi
movq %r12, %rsi
callq 0x2220
addq %r14, %r13
decq %rbx
jne 0x6407
movq -0x60(%rbp), %rax
incq %rax
movq -0xd0(%rbp), %r13
addq $0x8, %r13
cmpq -0x40(%rbp), %rax
jne 0x63f8
movq -0x400(%rbp), %rax
movq -0x18(%rax), %rdi
leaq -0x400(%rbp), %rax
addq %rax, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
leaq -0x400(%rbp), %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
cmpq $0x0, -0x40(%rbp)
movq -0xc8(%rbp), %r15
je 0x6503
movl -0x44(%rbp), %eax
cmpl $0x2, %eax
movl $0x1, %ecx
cmovgel %eax, %ecx
movq %rcx, -0x38(%rbp)
movq -0x40(%rbp), %rax
leaq (,%rax,8), %r14
xorl %eax, %eax
leaq -0x400(%rbp), %r12
leaq 0x3041(%rip), %r13      # 0x94f8
movq %rax, -0x60(%rbp)
movq -0x38(%rbp), %rbx
movq %r15, -0xc8(%rbp)
movsd (%r15), %xmm0
movq %r12, %rdi
callq 0x22c0
movl $0x1, %edx
movq %rax, %rdi
movq %r13, %rsi
callq 0x2220
addq %r14, %r15
decq %rbx
jne 0x64c6
movq -0x60(%rbp), %rax
incq %rax
movq -0xc8(%rbp), %r15
addq $0x8, %r15
cmpq -0x40(%rbp), %rax
jne 0x64b7
movq -0x400(%rbp), %rax
leaq -0x400(%rbp), %rdi
addq -0x18(%rax), %rdi
movl $0xa, %esi
callq 0x21f0
leaq -0x400(%rbp), %rbx
movsbl %al, %esi
leaq -0x400(%rbp), %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
movq -0x400(%rbp), %rax
addq -0x18(%rax), %rbx
movq %rbx, %rdi
movl $0xa, %esi
callq 0x21f0
movsbl %al, %esi
leaq -0x400(%rbp), %rdi
callq 0x2040
movq %rax, %rdi
callq 0x2140
leaq -0x400(%rbp), %r14
movq %r14, %rdi
callq 0x2270
movq 0x6a46(%rip), %rax      # 0xcfc8
movq (%rax), %rbx
movq 0x18(%rax), %r15
leaq -0x5f8(%rbp), %rdi
movq %rbx, -0x8(%rdi)
movq -0x18(%rbx), %rax
movq %r15, -0x600(%rbp,%rax)
callq 0x21b0
movl $0xf8, %r12d
leaq -0x600(%rbp), %rdi
addq %r12, %rdi
callq 0x20b0
leaq -0x3f8(%rbp), %rdi
movq %rbx, -0x8(%rdi)
movq -0x18(%rbx), %rax
movq %r15, -0x400(%rbp,%rax)
callq 0x21b0
addq %r12, %r14
movq %r14, %rdi
callq 0x20b0
xorl %eax, %eax
leaq -0x28(%rbp), %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq 0x69d7(%rip), %rdi      # 0xcfd0
leaq 0x2cf1(%rip), %rsi      # 0x92f1
callq 0x2200
movq %rax, %rdi
callq 0x20f0
jmp 0x6662
movq 0x69ba(%rip), %rdi      # 0xcfd0
leaq 0x2d5f(%rip), %rsi      # 0x937c
jmp 0x665d
movq 0x69aa(%rip), %rdi      # 0xcfd0
leaq 0x2dac(%rip), %rsi      # 0x93d9
jmp 0x665d
movq 0x699a(%rip), %rdi      # 0xcfd0
leaq 0x2dcb(%rip), %rsi      # 0x9408
jmp 0x665d
movq 0x698a(%rip), %rdi      # 0xcfd0
leaq 0x2de7(%rip), %rsi      # 0x9434
jmp 0x665d
movq 0x697a(%rip), %rdi      # 0xcfd0
leaq 0x2d51(%rip), %rsi      # 0x93ae
callq 0x2200
movl $0x1, %edi
callq 0x2100
leaq 0x6605(%rip), %rsi      # 0xcc78
movq %rax, %rdi
xorl %edx, %edx
callq 0x2320
jmp 0x66d1
jmp 0x6765
jmp 0x6765
jmp 0x66d1
jmp 0x6765
jmp 0x66d1
jmp 0x6765
jmp 0x66ee
jmp 0x66ff
jmp 0x6721
jmp 0x6710
jmp 0x6760
movq %rax, %rbx
jmp 0x674d
movq %rax, %rbx
jmp 0x673c
jmp 0x66d1
jmp 0x6765
jmp 0x6765
jmp 0x6765
jmp 0x6765
jmp 0x66dd
jmp 0x66ff
movq %rax, %rbx
jmp 0x6774
jmp 0x6721
jmp 0x6721
movq %rax, %rbx
leaq -0x188(%rbp), %rdi
callq 0x4844
jmp 0x66f1
movq %rax, %rbx
leaq -0x1b0(%rbp), %rdi
callq 0x4844
jmp 0x6702
movq %rax, %rbx
leaq -0x1d8(%rbp), %rdi
callq 0x4844
jmp 0x6713
movq %rax, %rbx
leaq -0x200(%rbp), %rdi
callq 0x4844
jmp 0x6724
movq %rax, %rbx
leaq -0x6a0(%rbp), %rdi
callq 0x4844
leaq -0x600(%rbp), %rdi
callq 0x4844
movq -0x150(%rbp), %rdi
testq %rdi, %rdi
je 0x674d
callq 0x2250
movq -0x160(%rbp), %rdi
testq %rdi, %rdi
je 0x6774
callq 0x2250
jmp 0x6774
movq %rax, %rbx
jmp 0x6730
movq %rax, %rbx
leaq -0x600(%rbp), %rdi
callq 0x2290
movq 0x684d(%rip), %rax      # 0xcfc8
movq (%rax), %rcx
movq 0x18(%rax), %rax
leaq -0x3f8(%rbp), %rdi
movq %rcx, -0x8(%rdi)
movq -0x18(%rcx), %rcx
movq %rax, -0x400(%rbp,%rcx)
callq 0x21b0
leaq -0x308(%rbp), %rdi
callq 0x20b0
movq %rbx, %rdi
callq 0x2350