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

```asm
<corth>:
  endbr64
  push   %r15
  mov    %edx,%eax
  mov    %ecx,%r15d
  add    $0x1,%edx
  push   %r14
  sub    $0x1,%ecx
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %edi,0x14(%rsp)
  mov    0xc0(%rsp),%rbx
  mov    %esi,0x78(%rsp)
  mov    0xc8(%rsp),%rbp
  mov    %ecx,0x7c(%rsp)
  cmp    %edx,%ecx
  jl     L69
  movslq %edi,%rsi
  mov    %r9,%r13
  mov    %edi,%r9d
  mov    %r8,%r12
  mov    %rsi,%rdi
  movslq %edx,%r8
  pxor   %xmm4,%xmm4
  neg    %rsi
  imul   %edi,%ecx
  movq   D22(%rip),%xmm5
  shl    $0x3,%rsi
  movslq %ecx,%rdx
  mov    %r9d,%ecx
  imul   %eax,%ecx
  movslq %ecx,%rdi
  mov    %ecx,0x10(%rsp)
  sub    %rdi,%rdx
  lea    0x0(,%rdx,8),%rdi
  mov    %r15d,%edx
  sub    %eax,%edx
  mov    %rdi,0x8(%rsp)
  lea    -0x1(%rdx),%edi
  lea    -0x1(%rax),%edx
  imul   %r9d,%edx
  mov    %edi,0x18(%rsp)
  lea    0x1(%r9),%edi
  neg    %r9d
  mov    %edi,0x24(%rsp)
  add    %edx,%eax
  mov    %eax,0x1c(%rsp)
  movslq %r9d,%rax
  lea    0x0(,%rax,8),%r11
  nopl   0x0(%rax,%rax,1)
L53:
  movq   $0x0,-0x8(%rbx,%r8,8)
  mov    %r8d,%r9d
  mov    %r8d,%edi
  movq   $0x0,-0x8(%rbp,%r8,8)
  mov    %r8d,0x20(%rsp)
  cmp    %r8d,%r15d
  jl     L114
  mov    0x18(%rsp),%r10d
  mov    0x10(%rsp),%ecx
  movapd %xmm4,%xmm3
  sub    0x14(%rsp),%ecx
  movslq %ecx,%rcx
  add    %r8,%r10
  lea    -0x1(%r8,%rcx,1),%rax
  lea    (%rcx,%r10,1),%rdx
  shl    $0x3,%rax
  shl    $0x3,%rdx
  nopw   0x0(%rax,%rax,1)
L95:
  movsd  (%r12,%rax,1),%xmm0
  movsd  0x0(%r13,%rax,1),%xmm1
  add    $0x8,%rax
  andpd  %xmm5,%xmm0
  andpd  %xmm5,%xmm1
  addsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm3
  cmp    %rax,%rdx
  jne    L95
  ucomisd %xmm4,%xmm3
  jp     L16
  jne    L16
L114:
  mov    0x14(%rsp),%edi
  add    %rsi,0x8(%rsp)
  add    $0x1,%r8
  add    %edi,0x10(%rsp)
  mov    0x24(%rsp),%edi
  add    %edi,0x1c(%rsp)
  subl   $0x1,0x18(%rsp)
  jne    L53
L69:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L16:
  lea    (%r15,%r9,1),%eax
  movapd %xmm4,%xmm2
  cltq
  mov    %rax,%rdx
  sub    %r10,%rax
  sub    %r8,%rdx
  lea    -0x10(,%rax,8),%r14
  lea    -0x8(,%rdx,8),%r9
  lea    0x0(,%rcx,8),%rdx
  lea    (%r12,%rdx,1),%r10
  mov    %r9,%rax
  add    %r13,%rdx
  nopl   (%rax)
L48:
  movsd  (%r10,%rax,1),%xmm0
  divsd  %xmm3,%xmm0
  movsd  %xmm0,(%rbx,%rax,1)
  movsd  (%rdx,%rax,1),%xmm0
  divsd  %xmm3,%xmm0
  movsd  %xmm0,0x0(%rbp,%rax,1)
  movsd  (%rbx,%rax,1),%xmm1
  mulsd  %xmm0,%xmm0
  sub    $0x8,%rax
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm2
  addsd  %xmm0,%xmm2
  cmp    %rax,%r14
  jne    L48
  ucomisd %xmm2,%xmm4
  ja     L169
  movapd %xmm2,%xmm6
  sqrtsd %xmm6,%xmm6
L21:
  movsd  -0x8(%rbx,%r8,8),%xmm0
  movsd  -0x8(%rbp,%r8,8),%xmm1
  mov    $0x2,%eax
  mov    %rsi,0x70(%rsp)
  mov    %r11,0x68(%rsp)
  mov    %r9,0x60(%rsp)
  mov    %edi,0x50(%rsp)
  mov    %r8,0x38(%rsp)
  movsd  %xmm3,0x58(%rsp)
  movsd  %xmm2,0x48(%rsp)
  movsd  %xmm6,0x40(%rsp)
  call   <pythag@plt>
  movslq 0x1c(%rsp),%rdx
  mov    0x38(%rsp),%r8
  pxor   %xmm4,%xmm4
  movsd  0x40(%rsp),%xmm6
  movsd  0x48(%rsp),%xmm2
  lea    0x0(,%rdx,8),%rdi
  movsd  0x58(%rsp),%xmm3
  mov    0x60(%rsp),%r9
  mov    0x68(%rsp),%r11
  mov    0x70(%rsp),%rsi
  mov    %rdi,0x30(%rsp)
  add    %r12,%rdi
  test   %eax,%eax
  mov    %rdi,0x28(%rsp)
  movq   D22(%rip),%xmm5
  mov    0x50(%rsp),%edi
  je     L177
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm6,%xmm1
  divsd  %xmm0,%xmm6
  movsd  D17(%rip),%xmm0
  addsd  %xmm1,%xmm2
  movsd  -0x8(%rbx,%r8,8),%xmm1
  addsd  %xmm6,%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,-0x8(%rbx,%r8,8)
  mulsd  -0x8(%rbp,%r8,8),%xmm0
  movsd  %xmm0,-0x8(%rbp,%r8,8)
L0:
  mov    0x78(%rsp),%eax
  mov    0x20(%rsp),%ecx
  cmp    %ecx,%eax
  jl     L202
  movslq 0x10(%rsp),%rdx
  lea    0x1(%rax),%r10d
  shl    $0x3,%rdx
  lea    (%r12,%rdx,1),%rcx
  add    %r13,%rdx
L193:
  mov    %r9,%rax
  movapd %xmm4,%xmm7
  movapd %xmm4,%xmm1
  nopl   0x0(%rax,%rax,1)
L54:
  movsd  (%rbx,%rax,1),%xmm0
  movsd  (%rcx,%rax,1),%xmm9
  movsd  (%rdx,%rax,1),%xmm10
  movapd %xmm0,%xmm8
  mulsd  %xmm9,%xmm8
  mulsd  %xmm10,%xmm0
  addsd  %xmm8,%xmm7
  movsd  0x0(%rbp,%rax,1),%xmm8
  sub    $0x8,%rax
  addsd  %xmm0,%xmm1
  mulsd  %xmm8,%xmm9
  movapd %xmm8,%xmm11
  mulsd  %xmm10,%xmm11
  subsd  %xmm9,%xmm1
  addsd  %xmm11,%xmm7
  cmp    %rax,%r14
  jne    L54
  divsd  %xmm2,%xmm7
  mov    %r8,%rax
  divsd  %xmm2,%xmm1
  nopl   (%rax)
L46:
  movsd  -0x8(%rbx,%rax,8),%xmm8
  movsd  -0x8(%rcx,%rax,8),%xmm0
  mulsd  %xmm7,%xmm8
  subsd  %xmm8,%xmm0
  movsd  -0x8(%rbp,%rax,8),%xmm8
  mulsd  %xmm1,%xmm8
  addsd  %xmm8,%xmm0
  movsd  %xmm0,-0x8(%rcx,%rax,8)
  movsd  -0x8(%rbp,%rax,8),%xmm8
  movsd  -0x8(%rdx,%rax,8),%xmm0
  mulsd  %xmm7,%xmm8
  subsd  %xmm8,%xmm0
  movsd  -0x8(%rbx,%rax,8),%xmm8
  mulsd  %xmm1,%xmm8
  subsd  %xmm8,%xmm0
  movsd  %xmm0,-0x8(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r15d
  jge    L46
  add    $0x1,%edi
  sub    %rsi,%rcx
  sub    %rsi,%rdx
  cmp    %r10d,%edi
  jne    L193
L202:
  test   %r15d,%r15d
  jle    L14
  movslq 0x10(%rsp),%rcx
  mov    0x7c(%rsp),%r10d
  xor    %edi,%edi
  shl    $0x3,%rcx
L38:
  mov    0x8(%rsp),%rax
  movapd %xmm4,%xmm7
  movapd %xmm4,%xmm1
  lea    (%rax,%rcx,1),%rdx
  mov    %r9,%rax
  nopl   0x0(%rax)
L175:
  movsd  (%rbx,%rax,1),%xmm0
  movsd  (%r12,%rdx,1),%xmm9
  movsd  0x0(%r13,%rdx,1),%xmm10
  add    %r11,%rdx
  movapd %xmm0,%xmm8
  mulsd  %xmm9,%xmm8
  mulsd  %xmm10,%xmm0
  addsd  %xmm8,%xmm7
  movsd  0x0(%rbp,%rax,1),%xmm8
  sub    $0x8,%rax
  addsd  %xmm0,%xmm1
  mulsd  %xmm8,%xmm9
  movapd %xmm8,%xmm11
  mulsd  %xmm10,%xmm11
  addsd  %xmm9,%xmm1
  subsd  %xmm11,%xmm7
  cmp    %rax,%r14
  jne    L175
  divsd  %xmm2,%xmm7
  mov    %rcx,%rdx
  mov    %r8,%rax
  divsd  %xmm2,%xmm1
  nopl   0x0(%rax)
L36:
  movsd  -0x8(%rbx,%rax,8),%xmm8
  movsd  (%r12,%rdx,1),%xmm0
  mulsd  %xmm7,%xmm8
  subsd  %xmm8,%xmm0
  movsd  -0x8(%rbp,%rax,8),%xmm8
  mulsd  %xmm1,%xmm8
  subsd  %xmm8,%xmm0
  movsd  %xmm0,(%r12,%rdx,1)
  movsd  -0x8(%rbx,%rax,8),%xmm8
  movsd  -0x8(%rbp,%rax,8),%xmm0
  add    $0x1,%rax
  mulsd  %xmm1,%xmm8
  mulsd  %xmm7,%xmm0
  addsd  0x0(%r13,%rdx,1),%xmm0
  subsd  %xmm8,%xmm0
  movsd  %xmm0,0x0(%r13,%rdx,1)
  sub    %rsi,%rdx
  cmp    %eax,%r15d
  jge    L36
  lea    0x1(%rdi),%rax
  add    $0x8,%rcx
  cmp    %rdi,%r10
  je     L14
  mov    %rax,%rdi
  jmp    L38
  nopl   0x0(%rax,%rax,1)
L14:
  movsd  -0x8(%rbx,%r8,8),%xmm0
  mov    0x28(%rsp),%rax
  xorpd  D4(%rip),%xmm6
  mulsd  %xmm3,%xmm0
  movsd  %xmm0,-0x8(%rbx,%r8,8)
  mulsd  -0x8(%rbp,%r8,8),%xmm3
  movsd  %xmm3,-0x8(%rbp,%r8,8)
  movsd  (%rax),%xmm0
  mulsd  %xmm6,%xmm0
  movsd  %xmm0,(%rax)
  mov    0x30(%rsp),%rax
  add    %r13,%rax
  mulsd  (%rax),%xmm6
  movsd  %xmm6,(%rax)
  jmp    L114
L177:
  mov    0x28(%rsp),%rax
  movsd  %xmm6,-0x8(%rbx,%r8,8)
  movsd  %xmm3,(%rax)
  jmp    L0
L169:
  movapd %xmm2,%xmm0
  mov    %rsi,0x58(%rsp)
  mov    %r11,0x50(%rsp)
  mov    %r8,0x48(%rsp)
  mov    %r9,0x40(%rsp)
  mov    %edi,0x30(%rsp)
  movsd  %xmm3,0x38(%rsp)
  movsd  %xmm2,0x28(%rsp)
  call   <sqrt@plt>
  mov    0x58(%rsp),%rsi
  mov    0x50(%rsp),%r11
  mov    0x48(%rsp),%r8
  mov    0x40(%rsp),%r9
  movapd %xmm0,%xmm6
  movsd  0x38(%rsp),%xmm3
  mov    0x30(%rsp),%edi
  movsd  0x28(%rsp),%xmm2
  jmp    L21
```

tool_calls: name: disassemble_data
arguments: label: D22
type_guessed: qword
name: disassemble_data
arguments: label: D17
type_guessed: f64
name: disassemble_data
arguments: label: D4
type_guessed: f32[4]

tool: D22: 9223372036854775807

tool: D17: 1

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

assistant: ```c
void corth(int nm, int n, int low, int igh, double *ar, double *ai, double *ortr, double *orti) {
  int i, j, m, ii, jj, la, mp, kp1;
  double f, g, h, fi, fr, scale;
  la = igh - 1;
  kp1 = low + 1;
  if (la < kp1)
    return;
  for (m = kp1; m <= la; m++) {
    h = 0.0;
    ortr[(m)-1] = 0.0;
    orti[(m)-1] = 0.0;
    scale = 0.0;
    for (i = m; i <= igh; i++)
      scale += fabs(ar[(i)-1 + ((m - 1) - 1) * (nm)]) + fabs(ai[(i)-1 + ((m - 1) - 1) * (nm)]);
    if (scale == 0.0)
      continue;
    mp = m + igh;
    for (ii = m; ii <= igh; ii++) {
      i = mp - ii;
      ortr[(i)-1] = ar[(i)-1 + ((m - 1) - 1) * (nm)] / scale;
      orti[(i)-1] = ai[(i)-1 + ((m - 1) - 1) * (nm)] / scale;
      h = h + ortr[(i)-1] * ortr[(i)-1] + orti[(i)-1] * orti[(i)-1];
    }
    g = sqrt(h);
    f = pythag(ortr[(m)-1], orti[(m)-1]);
    if (f != 0.0) {
      h += f * g;
      g /= f;
      ortr[(m)-1] = (1.0 + g) * ortr[(m)-1];
      orti[(m)-1] = (1.0 + g) * orti[(m)-1];
    } else {
      ortr[(m)-1] = g;
      ar[(m)-1 + ((m - 1) - 1) * (nm)] = scale;
    }
    for (j = m; j <= n; j++) {
      fr = 0.0;
      fi = 0.0;
      for (ii = m; ii <= igh; ii++) {
        i = mp - ii;
        fr = fr + ortr[(i)-1] * ar[(i)-1 + ((j)-1) * (nm)] + orti[(i)-1] * ai[(i)-1 + ((j)-1) * (nm)];
        fi = fi + ortr[(i)-1] * ai[(i)-1 + ((j)-1) * (nm)] - orti[(i)-1] * ar[(i)-1 + ((j)-1) * (nm)];
      }
      fr /= h;
      fi /= h;
      for (i = m; i <= igh; i++) {
        ar[(i)-1 + ((j)-1) * (nm)] = ar[(i)-1 + ((j)-1) * (nm)] - fr * ortr[(i)-1] + fi * orti[(i)-1];
        ai[(i)-1 + ((j)-1) * (nm)] = ai[(i)-1 + ((j)-1) * (nm)] - fr * orti[(i)-1] - fi * ortr[(i)-1];
      }
    }
    for (i = 1; i <= igh; i++) {
      fr = 0.0;
      fi = 0.0;
      for (jj = m; jj <= igh; jj++) {
        j = mp - jj;
        fr = fr + ortr[(j)-1] * ar[(i)-1 + ((j)-1) * (nm)] - orti[(j)-1] * ai[(i)-1 + ((j)-1) * (nm)];
        fi = fi + ortr[(j)-1] * ai[(i)-1 + ((j)-1) * (nm)] + orti[(j)-1] * ar[(i)-1 + ((j)-1) * (nm)];
      }
      fr /= h;
      fi /= h;
      for (j = m; j <= igh; j++) {
        ar[(i)-1 + ((j)-1) * (nm)] = ar[(i)-1 + ((j)-1) * (nm)] - fr * ortr[(j)-1] - fi * orti[(j)-1];
        ai[(i)-1 + ((j)-1) * (nm)] = ai[(i)-1 + ((j)-1) * (nm)] + fr * orti[(j)-1] - fi * ortr[(j)-1];
      }
    }
    ortr[(m)-1] = scale * ortr[(m)-1];
    orti[(m)-1] = scale * orti[(m)-1];
    ar[(m)-1 + ((m - 1) - 1) * (nm)] = -g * ar[(m)-1 + ((m - 1) - 1) * (nm)];
    ai[(m)-1 + ((m - 1) - 1) * (nm)] = -g * ai[(m)-1 + ((m - 1) - 1) * (nm)];
  }
}
```