// File: /workspace/github_c_makefile_star2_O0/PoulinV[P]AxiCLASS/build/../tools/arrays.c
// Function: array_spline_table_one_column
int array_spline_table_one_column(
		       double * x, /* vector of size x_size */
		       int x_size,
		       double * y_array, /* array of size x_size*y_size with elements
					  y_array[index_y*x_size+index_x] */
		       int y_size,
		       int index_y,
		       double * ddy_array, /* array of size x_size*y_size */
		       short spline_mode,
		       ErrorMsg errmsg
		       ) {

  double p;
  double qn;
  double un;
  double * u;
  double sig;
  int index_x;
  double dy_first;
  double dy_last;

  u = (double*)malloc((x_size-1) * sizeof(double));
  if (u == NULL) {
    class_sprintf(errmsg,"%s(L:%d) Cannot allocate u",__func__,__LINE__);
    return _FAILURE_;
  }

  if (x_size==2) spline_mode = _SPLINE_NATURAL_; // in the case of only 2 x-values, only the natural spline method is appropriate, for _SPLINE_EST_DERIV_ at least 3 x-values are needed.

  /************************************************/

  index_x=0;

  if (spline_mode == _SPLINE_NATURAL_) {
    ddy_array[index_y*x_size+index_x] = 0.0;
    u[index_x] = 0.0;
  }
  else {
    if (spline_mode == _SPLINE_EST_DERIV_) {

      dy_first =
	((x[2]-x[0])*(x[2]-x[0])*
	 (y_array[index_y*x_size+1]-y_array[index_y*x_size+0])-
	 (x[1]-x[0])*(x[1]-x[0])*
	 (y_array[index_y*x_size+2]-y_array[index_y*x_size+0]))/
	((x[2]-x[0])*(x[1]-x[0])*(x[2]-x[1]));

      ddy_array[index_y*x_size+index_x] = -0.5;

      u[index_x] =
	(3./(x[1] -  x[0]))*
	((y_array[index_y*x_size+1]-y_array[index_y*x_size+0])/
	 (x[1] - x[0])-dy_first);

    }
    else {
      class_sprintf(errmsg,"%s(L:%d) Spline mode not identified: %d",__func__,__LINE__,spline_mode);
      return _FAILURE_;
    }
  }

  /************************************************/

  for (index_x=1; index_x < x_size-1; index_x++) {

    sig = (x[index_x] - x[index_x-1])/(x[index_x+1] - x[index_x-1]);

    p = sig * ddy_array[index_y*x_size+(index_x-1)] + 2.0;

    ddy_array[index_y*x_size+index_x] = (sig-1.0)/p;

    u[index_x] =
      (y_array[index_y*x_size+(index_x+1)] - y_array[index_y*x_size+index_x])
      / (x[index_x+1] - x[index_x])
      - (y_array[index_y*x_size+index_x] - y_array[index_y*x_size+(index_x-1)])
      / (x[index_x] - x[index_x-1]);

    u[index_x] = (6.0 * u[index_x] /
		  (x[index_x+1] - x[index_x-1])
		  - sig * u[index_x-1]) / p;

  }

  /************************************************/

  if (spline_mode == _SPLINE_NATURAL_) {

      qn=un=0.0;

  }
  else {
    if (spline_mode == _SPLINE_EST_DERIV_) {

      dy_last =
	((x[x_size-3]-x[x_size-1])*(x[x_size-3]-x[x_size-1])*
	 (y_array[index_y*x_size+(x_size-2)]-y_array[index_y*x_size+(x_size-1)])-
	 (x[x_size-2]-x[x_size-1])*(x[x_size-2]-x[x_size-1])*
	 (y_array[index_y*x_size+(x_size-3)]-y_array[index_y*x_size+(x_size-1)]))/
	((x[x_size-3]-x[x_size-1])*(x[x_size-2]-x[x_size-1])*(x[x_size-3]-x[x_size-2]));

      qn=0.5;

      un=
	(3./(x[x_size-1] - x[x_size-2]))*
	(dy_last-(y_array[index_y*x_size+(x_size-1)] - y_array[index_y*x_size+(x_size-2)])/
	 (x[x_size-1] - x[x_size-2]));

    }
    else {
      class_sprintf(errmsg,"%s(L:%d) Spline mode not identified: %d",__func__,__LINE__,spline_mode);
      return _FAILURE_;
    }
  }

  /************************************************/

  index_x=x_size-1;

  ddy_array[index_y*x_size+index_x] =
    (un - qn * u[index_x-1]) /
    (qn * ddy_array[index_y*x_size+(index_x-1)] + 1.0);

  for (index_x=x_size-2; index_x >= 0; index_x--) {

    ddy_array[index_y*x_size+index_x] = ddy_array[index_y*x_size+index_x] *
      ddy_array[index_y*x_size+(index_x+1)] + u[index_x];

  }

  free(u);

  return _SUCCESS_;
}

// Assembly:
subq $0x88, %rsp
movq 0x98(%rsp), %rax
movw 0x90(%rsp), %ax
movq %rdi, 0x78(%rsp)
movl %esi, 0x74(%rsp)
movq %rdx, 0x68(%rsp)
movl %ecx, 0x64(%rsp)
movl %r8d, 0x60(%rsp)
movq %r9, 0x58(%rsp)
movl 0x74(%rsp), %eax
subl $0x1, %eax
movslq %eax, %rdi
shlq $0x3, %rdi
callq 0x8460
movq %rax, 0x38(%rsp)
cmpq $0x0, 0x38(%rsp)
jne 0x20fc8
jmp 0x20f65
xorl %eax, %eax
movl %eax, %esi
leaq 0x1b2e94(%rip), %rdx    # 0x1d3e04
leaq 0x1b308f(%rip), %rcx    # 0x1d4006
movl $0x4be, %r8d            # imm = 0x4BE
movq %rsi, %rdi
movb $0x0, %al
callq 0x8540
movl %eax, 0x14(%rsp)
movq 0x98(%rsp), %rdi
movl 0x14(%rsp), %eax
addl $0x1, %eax
movslq %eax, %rsi
leaq 0x1b2e60(%rip), %rdx    # 0x1d3e04
leaq 0x1b305b(%rip), %rcx    # 0x1d4006
movl $0x4be, %r8d            # imm = 0x4BE
movb $0x0, %al
callq 0x8540
movl $0x1, 0x84(%rsp)
jmp 0x21945
cmpl $0x2, 0x74(%rsp)
jne 0x20fd9
movw $0x0, 0x90(%rsp)
movl $0x0, 0x2c(%rsp)
movswl 0x90(%rsp), %eax
cmpl $0x0, %eax
jne 0x21022
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq 0x38(%rsp), %rax
movslq 0x2c(%rsp), %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
jmp 0x21278
movswl 0x90(%rsp), %eax
cmpl $0x1, %eax
jne 0x211ff
movq 0x78(%rsp), %rax
movsd 0x10(%rax), %xmm0
movq 0x78(%rsp), %rax
subsd (%rax), %xmm0
movq 0x78(%rsp), %rax
movsd 0x10(%rax), %xmm1
movq 0x78(%rsp), %rax
subsd (%rax), %xmm1
mulsd %xmm1, %xmm0
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl $0x0, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm2
movq 0x78(%rsp), %rax
movsd 0x8(%rax), %xmm1
movq 0x78(%rsp), %rax
subsd (%rax), %xmm1
movq 0x78(%rsp), %rax
movsd 0x8(%rax), %xmm3
movq 0x78(%rsp), %rax
subsd (%rax), %xmm3
mulsd %xmm3, %xmm1
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl $0x2, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm3
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl $0x0, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm3
mulsd %xmm3, %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x78(%rsp), %rax
movsd 0x10(%rax), %xmm1
movq 0x78(%rsp), %rax
subsd (%rax), %xmm1
movq 0x78(%rsp), %rax
movsd 0x8(%rax), %xmm2
movq 0x78(%rsp), %rax
subsd (%rax), %xmm2
mulsd %xmm2, %xmm1
movq 0x78(%rsp), %rax
movsd 0x10(%rax), %xmm2
movq 0x78(%rsp), %rax
subsd 0x8(%rax), %xmm2
mulsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, 0x20(%rsp)
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
movsd 0x1b2529(%rip), %xmm0   # 0x1d36a0
movsd %xmm0, (%rax,%rcx,8)
movq 0x78(%rsp), %rax
movsd 0x8(%rax), %xmm1
movq 0x78(%rsp), %rax
subsd (%rax), %xmm1
movsd 0x1b2501(%rip), %xmm0   # 0x1d3698
divsd %xmm1, %xmm0
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl $0x0, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movsd 0x8(%rax), %xmm2
movq 0x78(%rsp), %rax
subsd (%rax), %xmm2
divsd %xmm2, %xmm1
subsd 0x20(%rsp), %xmm1
mulsd %xmm1, %xmm0
movq 0x38(%rsp), %rax
movslq 0x2c(%rsp), %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0x21276
jmp 0x21201
movswl 0x90(%rsp), %r9d
xorl %eax, %eax
movl %eax, %esi
leaq 0x1b2c0a(%rip), %rdx    # 0x1d3e1f
leaq 0x1b2dea(%rip), %rcx    # 0x1d4006
movl $0x4df, %r8d            # imm = 0x4DF
movq %rsi, %rdi
movb $0x0, %al
callq 0x8540
movl %eax, 0x10(%rsp)
movq 0x98(%rsp), %rdi
movl 0x10(%rsp), %eax
addl $0x1, %eax
movslq %eax, %rsi
movswl 0x90(%rsp), %r9d
leaq 0x1b2bcd(%rip), %rdx    # 0x1d3e1f
leaq 0x1b2dad(%rip), %rcx    # 0x1d4006
movl $0x4df, %r8d            # imm = 0x4DF
movb $0x0, %al
callq 0x8540
movl $0x1, 0x84(%rsp)
jmp 0x21945
jmp 0x21278
movl $0x1, 0x2c(%rsp)
movl 0x2c(%rsp), %eax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jge 0x214c9
movq 0x78(%rsp), %rax
movslq 0x2c(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x78(%rsp), %rax
movl 0x2c(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm0
movq 0x78(%rsp), %rax
movl 0x2c(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x2c(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, 0x30(%rsp)
movsd 0x30(%rsp), %xmm0
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x2c(%rsp), %edx
subl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
mulsd %xmm1, %xmm0
movsd 0x1b1967(%rip), %xmm1   # 0x1d2c80
addsd %xmm1, %xmm0
movsd %xmm0, 0x50(%rsp)
movsd 0x30(%rsp), %xmm0
movsd 0x1b0f17(%rip), %xmm1   # 0x1d2248
subsd %xmm1, %xmm0
divsd 0x50(%rsp), %xmm0
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x2c(%rsp), %edx
addl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm0
movq 0x78(%rsp), %rax
movl 0x2c(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movslq 0x2c(%rsp), %rcx
subsd (%rax,%rcx,8), %xmm1
divsd %xmm1, %xmm0
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x2c(%rsp), %edx
subl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movslq 0x2c(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x78(%rsp), %rax
movl 0x2c(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm2
divsd %xmm2, %xmm1
subsd %xmm1, %xmm0
movq 0x38(%rsp), %rax
movslq 0x2c(%rsp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x38(%rsp), %rax
movslq 0x2c(%rsp), %rcx
movsd 0x1b224e(%rip), %xmm1   # 0x1d3688
mulsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x2c(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x78(%rsp), %rax
movl 0x2c(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm0
divsd %xmm0, %xmm1
movsd 0x30(%rsp), %xmm0
movq 0x38(%rsp), %rax
movl 0x2c(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
divsd 0x50(%rsp), %xmm0
movq 0x38(%rsp), %rax
movslq 0x2c(%rsp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl 0x2c(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x2c(%rsp)
jmp 0x21280
movswl 0x90(%rsp), %eax
cmpl $0x0, %eax
jne 0x214ed
xorps %xmm0, %xmm0
movsd %xmm0, 0x40(%rsp)
xorps %xmm0, %xmm0
movsd %xmm0, 0x48(%rsp)
jmp 0x2180a
movswl 0x90(%rsp), %eax
cmpl $0x1, %eax
jne 0x21791
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x3, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm0
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x3, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm1
mulsd %xmm1, %xmm0
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x74(%rsp), %edx
subl $0x2, %edx
addl %edx, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x74(%rsp), %edx
subl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm2
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm3
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm3
mulsd %xmm3, %xmm1
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x74(%rsp), %edx
subl $0x3, %edx
addl %edx, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm3
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x74(%rsp), %edx
subl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm3
mulsd %xmm3, %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x3, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm2
mulsd %xmm2, %xmm1
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x3, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm2
mulsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, 0x18(%rsp)
movsd 0x1b15b1(%rip), %xmm0   # 0x1d2c88
movsd %xmm0, 0x48(%rsp)
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm1
movsd 0x1b1f8b(%rip), %xmm0   # 0x1d3698
divsd %xmm1, %xmm0
movsd 0x18(%rsp), %xmm1
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x74(%rsp), %edx
subl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x68(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x74(%rsp), %edx
subl $0x2, %edx
addl %edx, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm2
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm3
movq 0x78(%rsp), %rax
movl 0x74(%rsp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
subsd (%rax,%rcx,8), %xmm3
divsd %xmm3, %xmm2
subsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x40(%rsp)
jmp 0x21808
jmp 0x21793
movswl 0x90(%rsp), %r9d
xorl %eax, %eax
movl %eax, %esi
leaq 0x1b2678(%rip), %rdx    # 0x1d3e1f
leaq 0x1b2858(%rip), %rcx    # 0x1d4006
movl $0x514, %r8d            # imm = 0x514
movq %rsi, %rdi
movb $0x0, %al
callq 0x8540
movl %eax, 0xc(%rsp)
movq 0x98(%rsp), %rdi
movl 0xc(%rsp), %eax
addl $0x1, %eax
movslq %eax, %rsi
movswl 0x90(%rsp), %r9d
leaq 0x1b263b(%rip), %rdx    # 0x1d3e1f
leaq 0x1b281b(%rip), %rcx    # 0x1d4006
movl $0x514, %r8d            # imm = 0x514
movb $0x0, %al
callq 0x8540
movl $0x1, 0x84(%rsp)
jmp 0x21945
jmp 0x2180a
movl 0x74(%rsp), %eax
subl $0x1, %eax
movl %eax, 0x2c(%rsp)
movsd 0x40(%rsp), %xmm1
movsd 0x48(%rsp), %xmm0
movq 0x38(%rsp), %rax
movl 0x2c(%rsp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x48(%rsp), %xmm1
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x2c(%rsp), %edx
subl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
mulsd %xmm2, %xmm1
movsd 0x1b09c3(%rip), %xmm2   # 0x1d2248
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movl 0x74(%rsp), %eax
subl $0x2, %eax
movl %eax, 0x2c(%rsp)
cmpl $0x0, 0x2c(%rsp)
jl 0x21930
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
movl 0x2c(%rsp), %edx
addl $0x1, %edx
addl %edx, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x38(%rsp), %rax
movslq 0x2c(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x58(%rsp), %rax
movl 0x60(%rsp), %ecx
imull 0x74(%rsp), %ecx
addl 0x2c(%rsp), %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movl 0x2c(%rsp), %eax
addl $-0x1, %eax
movl %eax, 0x2c(%rsp)
jmp 0x218b2
movq 0x38(%rsp), %rdi
callq 0x8420
movl $0x0, 0x84(%rsp)
movl 0x84(%rsp), %eax
addq $0x88, %rsp
retq
nopw %cs:(%rax,%rax)