user: Translate the below assembly code into C programming language:

```asm
<BackSubstitute_Complex_Vector>:
  endbr64
  push   %r15
  lea    -0x1(%rcx),%eax
  mov    %rdi,%r15
  mov    %rdx,%r9
  push   %r14
  mov    %rsi,%r11
  mov    %r9,%r10
  push   %r13
  push   %r12
  movslq %eax,%r12
  push   %rbp
  mov    %r12d,%r14d
  push   %rbx
  mov    %ecx,%ebx
  imul   %r8d,%ebx
  sub    $0xd8,%rsp
  movq   D12(%rip),%xmm8
  mov    %ecx,0x40(%rsp)
  movslq %ebx,%rax
  sub    %r8d,%ecx
  lea    (%r15,%rax,8),%rax
  mov    %rdi,0x18(%rsp)
  mov    %r12,%rdi
  mov    %rdx,0x20(%rsp)
  movslq %ecx,%rdx
  sub    %r8d,%edi
  mov    %r12d,0x44(%rsp)
  shl    $0x3,%r12
  mov    %rsi,0xc0(%rsp)
  lea    (%rax,%r12,1),%rbp
  lea    (%rax,%rdx,8),%rsi
  movslq %edi,%rdx
  movsd  0x0(%rbp),%xmm1
  lea    0x8(%r12),%rcx
  lea    (%rax,%rdx,8),%rdi
  movsd  %xmm0,0xc8(%rsp)
  movsd  (%rsi),%xmm0
  add    %rcx,%r11
  add    %rcx,%r10
  mov    %rcx,0x38(%rsp)
  movapd %xmm1,%xmm3
  mov    %r11,0x10(%rsp)
  movsd  (%r11),%xmm4
  lea    (%rax,%rcx,1),%r13
  movapd %xmm0,%xmm2
  andpd  %xmm8,%xmm3
  andpd  %xmm8,%xmm2
  comisd %xmm2,%xmm3
  jbe    L130
  movsd  0x0(%r13),%xmm0
  movq   D41(%rip),%xmm7
  subsd  %xmm4,%xmm0
  movaps %xmm7,0x90(%rsp)
  xorpd  %xmm7,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdi)
  movsd  (%r10),%xmm0
  xorpd  %xmm7,%xmm0
  divsd  0x0(%rbp),%xmm0
  movsd  %xmm0,(%rsi)
L125:
  mov    D6(%rip),%rax
  mov    0x18(%rsp),%rcx
  mov    %rax,0x0(%rbp)
  lea    (%r8,%r8,1),%eax
  sub    %eax,%ebx
  movq   $0x0,0x0(%r13)
  movslq %ebx,%rax
  lea    (%rcx,%rax,8),%r13
  mov    0x40(%rsp),%ecx
  mov    %ecx,%ebp
  sub    $0x2,%ebp
  js     L70
  movslq %r8d,%rbx
  lea    -0x8(%r13,%r12,1),%r15
  mov    %r8d,0x8(%rsp)
  movslq %ebp,%rbp
  shl    $0x3,%rbx
  pxor   %xmm7,%xmm7
  mov    %rbx,%rax
  neg    %rax
  mov    %rax,0x28(%rsp)
  mov    $0xfffffffffffffff8,%rax
  sub    %rbx,%rax
  mov    %rax,0x30(%rsp)
  mov    0x38(%rsp),%rax
  lea    0x0(%r13,%rax,1),%r12
  lea    (%r8,%rcx,1),%eax
  mov    %r15,%r8
  mov    %r13,%r15
  cltq
  lea    -0x8(%r13,%rax,8),%rcx
  mov    %rcx,%r13
  nopw   0x0(%rax,%rax,1)
L105:
  mov    0x10(%rsp),%rax
  movsd  (%r8),%xmm6
  mov    %ebp,%edi
  mov    %ebp,%r9d
  movsd  (%r12),%xmm4
  movsd  (%rax),%xmm3
  mov    0x8(%rsp),%eax
  imul   %r14d,%eax
  subsd  %xmm3,%xmm6
  cltq
  shl    $0x3,%rax
  cmp    %r14d,0x40(%rsp)
  jle    L13
  mov    0x44(%rsp),%esi
  movslq %r14d,%r11
  add    0x38(%rsp),%rax
  movapd %xmm7,%xmm5
  lea    (%r15,%r11,8),%rdx
  add    0x18(%rsp),%rax
  sub    %r14d,%esi
  add    %r11,%rsi
  lea    0x8(%r15,%rsi,8),%rsi
  nopw   0x0(%rax,%rax,1)
L36:
  movsd  (%rdx),%xmm0
  movsd  -0x8(%rax),%xmm1
  add    $0x8,%rdx
  mulsd  %xmm0,%xmm1
  mulsd  (%rax),%xmm0
  add    %rbx,%rax
  addsd  %xmm1,%xmm4
  addsd  %xmm0,%xmm5
  cmp    %rdx,%rsi
  jne    L36
L134:
  mov    0x20(%rsp),%rax
  movsd  (%rax,%rbp,8),%xmm1
  comisd %xmm1,%xmm7
  ja     L76
  ucomisd %xmm7,%xmm1
  movsd  (%r10),%xmm0
  lea    -0x8(%r12),%r11
  jp     L57
  je     L2
L57:
  mov    0xc0(%rsp),%rax
  mulsd  %xmm1,%xmm1
  add    0x8(%rsp),%edi
  mov    $0x0,%ecx
  movsd  0x8(%r8),%xmm9
  movapd %xmm10,%xmm13
  movslq %edi,%rdi
  movsd  (%rax,%rbp,8),%xmm2
  andpd  %xmm8,%xmm13
  movsd  (%r15,%rdi,8),%xmm11
  movapd %xmm9,%xmm12
  movsd  %xmm13,0x58(%rsp)
  subsd  %xmm3,%xmm2
  andpd  %xmm8,%xmm12
  movapd %xmm2,%xmm3
  addsd  %xmm2,%xmm3
  mulsd  %xmm2,%xmm2
  mulsd  %xmm0,%xmm3
  addsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  subsd  %xmm1,%xmm2
  ucomisd %xmm7,%xmm2
  setnp  %al
  cmovne %ecx,%eax
  test   %al,%al
  je     L7
  ucomisd %xmm7,%xmm3
  setnp  %al
  cmovne %ecx,%eax
  test   %al,%al
  je     L7
  movapd %xmm6,%xmm2
  movapd %xmm11,%xmm1
  andpd  %xmm8,%xmm2
  andpd  %xmm8,%xmm1
  addsd  %xmm12,%xmm2
  addsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  andpd  %xmm8,%xmm1
  addsd  %xmm13,%xmm2
  addsd  %xmm1,%xmm2
  mulsd  0xc8(%rsp),%xmm2
L7:
  movapd %xmm10,%xmm14
  movapd %xmm10,%xmm13
  mov    %r12,%rsi
  mov    %r11,%rdi
  movsd  0x50(%rsp),%xmm1
  mulsd  %xmm5,%xmm14
  mov    $0x4,%eax
  movsd  0x48(%rsp),%xmm15
  mulsd  %xmm4,%xmm13
  mov    %r10,0x88(%rsp)
  mulsd  %xmm9,%xmm1
  mov    %r8,0x80(%rsp)
  mulsd  %xmm9,%xmm15
  mov    %r9d,0x70(%rsp)
  movsd  %xmm12,0x78(%rsp)
  movsd  %xmm11,0x68(%rsp)
  movsd  %xmm10,0x60(%rsp)
  movsd  %xmm6,0xb8(%rsp)
  subsd  %xmm14,%xmm1
  movapd %xmm4,%xmm14
  mulsd  %xmm0,%xmm14
  subsd  %xmm13,%xmm15
  movsd  %xmm4,0xb0(%rsp)
  movsd  %xmm9,0xa8(%rsp)
  mulsd  %xmm5,%xmm0
  movsd  %xmm5,0xa0(%rsp)
  subsd  %xmm14,%xmm1
  addsd  %xmm15,%xmm0
  call   <Complex_Division@plt>
  movsd  0x78(%rsp),%xmm12
  lea    0x8(%r13),%rsi
  mov    0x88(%rsp),%r10
  movsd  0x60(%rsp),%xmm10
  movsd  0x68(%rsp),%xmm11
  pxor   %xmm7,%xmm7
  movsd  (%r10),%xmm3
  mov    0x70(%rsp),%r9d
  movq   D12(%rip),%xmm8
  mov    0x80(%rsp),%r8
  movapd %xmm3,%xmm0
  andpd  %xmm8,%xmm0
  addsd  0x58(%rsp),%xmm0
  comisd %xmm0,%xmm12
  jbe    L84
  mulsd  (%r12),%xmm3
  movsd  -0x8(%r12),%xmm0
  mov    %r9d,%r14d
  movsd  0xa8(%rsp),%xmm9
  movsd  0xb8(%rsp),%xmm6
  movapd 0x90(%rsp),%xmm2
  movsd  0xb0(%rsp),%xmm4
  movsd  0xa0(%rsp),%xmm5
  mulsd  %xmm6,%xmm0
  xorpd  %xmm2,%xmm4
  divsd  %xmm9,%xmm3
  xorpd  %xmm2,%xmm5
  subsd  %xmm0,%xmm4
  addsd  %xmm3,%xmm4
  movsd  %xmm4,0x0(%r13)
  movsd  (%r10),%xmm0
  mulsd  -0x8(%r12),%xmm0
  mulsd  (%r12),%xmm6
  divsd  %xmm9,%xmm0
  subsd  %xmm6,%xmm5
  subsd  %xmm0,%xmm5
  movsd  %xmm5,0x8(%r13)
L26:
  sub    $0x1,%rbp
  add    0x28(%rsp),%r15
  sub    %rbx,%r12
  sub    %rbx,%r13
  add    0x30(%rsp),%r8
  test   %ebp,%ebp
  jns    L105
L70:
  add    $0xd8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L76:
  movsd  %xmm5,0x50(%rsp)
  movapd %xmm6,%xmm10
  movsd  %xmm4,0x48(%rsp)
  jmp    L26
  nopl   0x0(%rax,%rax,1)
L2:
  movapd %xmm0,%xmm3
  mov    %r12,%rsi
  mov    %r11,%rdi
  mov    $0x4,%eax
  movapd 0x90(%rsp),%xmm7
  movapd %xmm6,%xmm2
  mov    %r10,0x70(%rsp)
  mov    %r8,0x68(%rsp)
  xorpd  %xmm7,%xmm5
  xorpd  %xmm7,%xmm4
  mov    %r9d,0x58(%rsp)
  movapd %xmm5,%xmm1
  movapd %xmm4,%xmm0
  movsd  %xmm10,0x60(%rsp)
  call   <Complex_Division@plt>
  mov    0x58(%rsp),%r14d
  mov    0x68(%rsp),%r8
  pxor   %xmm7,%xmm7
  movq   D12(%rip),%xmm8
  movsd  0x60(%rsp),%xmm10
  mov    0x70(%rsp),%r10
  jmp    L26
  nopw   0x0(%rax,%rax,1)
L84:
  movsd  (%r12),%xmm4
  movsd  0x50(%rsp),%xmm1
  movapd %xmm10,%xmm2
  mov    %r13,%rdi
  movsd  0x48(%rsp),%xmm0
  mov    $0x4,%eax
  movapd 0x90(%rsp),%xmm7
  mov    %r10,0x70(%rsp)
  mulsd  %xmm11,%xmm4
  mov    %r8,0x68(%rsp)
  mulsd  -0x8(%r12),%xmm11
  xorpd  %xmm7,%xmm1
  xorpd  %xmm7,%xmm0
  movsd  %xmm10,0x58(%rsp)
  mov    %r9d,0x60(%rsp)
  subsd  %xmm4,%xmm1
  subsd  %xmm11,%xmm0
  call   <Complex_Division@plt>
  mov    0x60(%rsp),%r14d
  mov    0x68(%rsp),%r8
  pxor   %xmm7,%xmm7
  movsd  0x58(%rsp),%xmm10
  mov    0x70(%rsp),%r10
  movq   D12(%rip),%xmm8
  jmp    L26
  nopl   0x0(%rax,%rax,1)
L13:
  movapd %xmm7,%xmm5
  jmp    L134
L130:
  movsd  (%rdi),%xmm2
  movq   D41(%rip),%xmm7
  pxor   %xmm1,%xmm1
  mov    $0x4,%eax
  movsd  (%r10),%xmm3
  mov    %r8d,0x30(%rsp)
  subsd  %xmm4,%xmm2
  xorpd  %xmm7,%xmm0
  mov    %r10,0x8(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movaps %xmm7,0x90(%rsp)
  call   <Complex_Division@plt>
  mov    0x30(%rsp),%r8d
  mov    0x8(%rsp),%r10
  movq   D12(%rip),%xmm8
  movsd  0x28(%rsp),%xmm10
  jmp    L125
```

tool_calls: name: read_data
arguments: label: D12
type: qword
name: read_data
arguments: label: D41
type: f32[2]
name: read_data
arguments: label: D6
type: f64

tool: D12: 9223372036854775807

tool: {"D41": [0.0,-0.0]}

tool: D6: 1

assistant: ```c
void BackSubstitute_Complex_Vector(double *H, double eigen_real[], double eigen_imag[], int row, double zero_tolerance, int n) {
  double *pH;
  double *pV;
  double x, y;
  double u[4];
  double v[2];
  double w[2];
  int i, j, k;
  k = row - 1;
  pH = H + n * row;
  if (fabs(pH[k]) > fabs(pH[row - n])) {
    pH[k - n] = -(pH[row] - eigen_real[row]) / pH[k];
    pH[row - n] = -eigen_imag[row] / pH[k];
  } else
    Complex_Division(-pH[row - n], 0.0, pH[k - n] - eigen_real[row], eigen_imag[row], &pH[k - n], &pH[row - n]);
  pH[k] = 1.0;
  pH[row] = 0.0;
  for (i = row - 2, pH = H + n * i; i >= 0; pH -= n, i--) {
    u[0] = pH[i] - eigen_real[row];
    w[0] = pH[row];
    w[1] = 0.0;
    pV = H + k * n;
    for (j = k; j < row; j++, pV += n) {
      w[0] += pH[j] * pV[row - 1];
      w[1] += pH[j] * pV[row];
    }
    if (eigen_imag[i] < 0.0) {
      u[3] = u[0];
      v[0] = w[0];
      v[1] = w[1];
    } else {
      k = i;
      if (eigen_imag[i] == 0.0) {
        Complex_Division(-w[0], -w[1], u[0], eigen_imag[row], &pH[row - 1], &pH[row]);
      } else {
        u[1] = pH[i + 1];
        u[2] = pH[n + i];
        x = eigen_real[i] - eigen_real[row];
        y = 2.0 * x * eigen_imag[row];
        x = x * x + eigen_imag[i] * eigen_imag[i] - eigen_imag[row] * eigen_imag[row];
        if (x == 0.0 && y == 0.0)
          x = zero_tolerance * (fabs(u[0]) + fabs(u[1]) + fabs(u[2]) + fabs(u[3]) + fabs(eigen_imag[row]));
        Complex_Division(u[1] * v[0] - u[3] * w[0] + w[1] * eigen_imag[row], u[1] * v[1] - u[3] * w[1] - w[0] * eigen_imag[row], x, y, &pH[row - 1], &pH[row]);
        if (fabs(u[1]) > (fabs(u[3]) + fabs(eigen_imag[row]))) {
          pH[n + row - 1] = -w[0] - u[0] * pH[row - 1] + eigen_imag[row] * pH[row] / u[1];
          pH[n + row] = -w[1] - u[0] * pH[row] - eigen_imag[row] * pH[row - 1] / u[1];
        } else {
          Complex_Division(-v[0] - u[2] * pH[row - 1], -v[1] - u[2] * pH[row], u[3], eigen_imag[row], &pH[n + row - 1], &pH[n + row]);
        }
      }
    }
  }
}
```