user: Converte this assembly snippet into its equivalent C code:

```asm
<Hessenberg_Form_Orthogonal>:
  endbr64
  push   %r15
  xor    %eax,%eax
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xd8,%rsp
  mov    %rdi,0x98(%rsp)
  mov    %rsi,0xb8(%rsp)
  mov    %edx,0x48(%rsp)
  cmp    $0x2,%edx
  jg     L200
L139:
  add    $0xd8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L200:
  mov    %rdi,%r15
  lea    0x0(,%rdx,8),%edi
  mov    %edx,%ebx
  movslq %edx,%r14
  call   <malloc@plt>
  mov    %rax,0x88(%rsp)
  test   %rax,%rax
  je     L405
  shl    $0x3,%r14
  lea    -0x1(%rbx),%esi
  mov    %rax,0x60(%rsp)
  lea    -0x2(%rbx),%r12d
  lea    (%r15,%r14,1),%rdi
  lea    0x8(%rax),%r15
  add    $0x18,%rax
  mov    %esi,0x24(%rsp)
  mov    %ebx,0x74(%rsp)
  movsd  D26(%rip),%xmm4
  pxor   %xmm3,%xmm3
  mov    %rdi,0x30(%rsp)
  lea    0x8(%r14),%rdi
  movq   D38(%rip),%xmm5
  mov    %rdi,0x78(%rsp)
  mov    %rsi,0x80(%rsp)
  mov    %ebx,0x28(%rsp)
  movq   $0x8,0x40(%rsp)
  movq   $0x1,0x18(%rsp)
  movq   $0x0,0x58(%rsp)
  mov    %rax,0xc8(%rsp)
  mov    %r15,0x8(%rsp)
  nopl   0x0(%rax)
L324:
  mov    0x18(%rsp),%rax
  mov    %eax,0x70(%rsp)
  mov    %eax,0x68(%rsp)
  mov    %eax,0xc0(%rsp)
  cmp    %eax,0x48(%rsp)
  jle    L327
  mov    0x30(%rsp),%rax
  mov    0x18(%rsp),%ebp
  movapd %xmm3,%xmm1
  mov    0x48(%rsp),%edx
  nopw   0x0(%rax,%rax,1)
L321:
  movsd  (%rax),%xmm0
  mov    %ebp,%ebx
  add    $0x1,%ebp
  add    %r14,%rax
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  cmp    %ebp,%edx
  jne    L321
  ucomisd %xmm3,%xmm1
  jp     L52
  mov    0x8(%rsp),%rax
  lea    0x8(%rax),%rax
  mov    %rax,0x10(%rsp)
  lea    -0x1(%r12),%eax
  mov    %eax,0x4c(%rsp)
  jne    L52
L360:
  mov    0x10(%rsp),%rdi
  addq   $0x1,0x18(%rsp)
  mov    0x78(%rsp),%rbx
  mov    0x48(%rsp),%esi
  add    %rbx,0x30(%rsp)
  mov    0x18(%rsp),%rax
  addq   $0x1,0x58(%rsp)
  mov    0x4c(%rsp),%r12d
  subl   $0x1,0x28(%rsp)
  addq   $0x8,0x40(%rsp)
  addq   $0x8,0x60(%rsp)
  subl   $0x1,0x24(%rsp)
  add    %esi,0x74(%rsp)
  mov    %rdi,0x8(%rsp)
  cmp    0x80(%rsp),%rax
  jne    L324
  mov    0x88(%rsp),%rdi
  call   <free@plt>
  xor    %eax,%eax
  jmp    L139
  nopl   (%rax)
L52:
  ucomisd %xmm1,%xmm3
  ja     L306
  sqrtsd %xmm1,%xmm1
L60:
  mov    0x30(%rsp),%rax
  movsd  (%rax),%xmm0
  comisd %xmm3,%xmm0
  jb     L131
  xorpd  %xmm5,%xmm1
L131:
  subsd  %xmm1,%xmm0
  mov    0x8(%rsp),%rax
  mov    0x70(%rsp),%esi
  lea    0x1(%rsi),%edx
  movsd  %xmm0,(%rax)
  mov    0x30(%rsp),%rax
  movsd  %xmm1,(%rax)
  add    %r14,%rax
  cmp    %edx,%ebp
  jle    L202
  mov    0x8(%rsp),%rsi
  lea    -0x1(%r12),%ecx
  mov    %ecx,0x4c(%rsp)
  add    0x58(%rsp),%rcx
  add    $0x8,%rsi
  mov    %rsi,0x10(%rsp)
  mov    %rsi,%rdx
  mov    0xc8(%rsp),%rsi
  lea    (%rsi,%rcx,8),%rcx
  nopl   0x0(%rax)
L314:
  movsd  (%rax),%xmm0
  add    $0x8,%rdx
  movsd  %xmm0,-0x8(%rdx)
  movq   $0x0,(%rax)
  add    %r14,%rax
  cmp    %rcx,%rdx
  jne    L314
L13:
  mov    0x8(%rsp),%r9
  movapd %xmm4,%xmm6
  mov    0x30(%rsp),%rcx
  mov    0x98(%rsp),%rdi
  movslq 0x74(%rsp),%rax
  mulsd  (%r9),%xmm1
  add    0x18(%rsp),%rax
  lea    0x8(%rcx),%rsi
  lea    (%rdi,%rax,8),%rdi
  lea    0x0(,%r12,8),%rax
  add    0x58(%rsp),%r12
  lea    0x10(%rcx,%rax,1),%r8
  mov    %rax,0x90(%rsp)
  mov    0x88(%rsp),%rax
  lea    0x10(%rax,%r12,8),%rcx
  divsd  %xmm1,%xmm6
  nopl   0x0(%rax)
L18:
  mov    %rdi,%rdx
  mov    %r9,%rax
  movapd %xmm3,%xmm1
  nopw   0x0(%rax,%rax,1)
L390:
  movsd  (%rax),%xmm0
  mulsd  (%rdx),%xmm0
  add    $0x8,%rax
  add    %r14,%rdx
  addsd  %xmm0,%xmm1
  cmp    %rcx,%rax
  jne    L390
  mulsd  %xmm6,%xmm1
  mov    %rsi,%rdx
  mov    %r9,%rax
  nopw   0x0(%rax,%rax,1)
L419:
  movsd  (%rax),%xmm2
  movsd  (%rdx),%xmm0
  add    $0x8,%rax
  mulsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rdx)
  add    %r14,%rdx
  cmp    %rcx,%rax
  jne    L419
  add    $0x8,%rsi
  add    $0x8,%rdi
  cmp    %r8,%rsi
  jne    L18
  mov    0x98(%rsp),%r15
  mov    0x60(%rsp),%rdx
  mov    %ebx,0x38(%rsp)
  mov    0x40(%rsp),%rax
  mov    0x24(%rsp),%edi
  mov    %ebp,0xc4(%rsp)
  mov    %r15,%r11
  mov    %r15,%rsi
  mov    %rdx,%rcx
  mov    0x28(%rsp),%ebx
  add    %rax,%r11
  sub    %rdx,%rax
  add    $0x18,%rdx
  mov    0x88(%rsp),%r13
  add    %rax,%rsi
  mov    %rax,0xb0(%rsp)
  mov    %edi,%eax
  and    $0xfffffffe,%edi
  shr    %eax
  mov    %edi,0x2c(%rsp)
  add    0x70(%rsp),%edi
  sub    %r11,%rcx
  sub    $0x1,%eax
  mov    %rdx,0xa8(%rsp)
  shl    $0x4,%rax
  cmpl   $0x2,0x28(%rsp)
  lea    (%rdx,%rax,1),%r9
  movslq %edi,%rdx
  movslq 0x68(%rsp),%rdi
  lea    0x0(,%rdx,8),%r10
  cmovne %rdx,%rdi
  add    $0x10,%rax
  lea    0x0(%r13,%r10,1),%r12
  mov    %r15,%rdx
  mov    %rax,0x50(%rsp)
  mov    0x90(%rsp),%rax
  shl    $0x3,%rdi
  add    $0x8,%rax
  lea    0x0(%r13,%rdi,1),%r8
  xor    %r13d,%r13d
  mov    %rax,0xa0(%rsp)
  cmp    $0x2,%ebx
  je     L225
L174:
  mov    0x8(%rsp),%rax
  movapd %xmm3,%xmm0
  nopl   0x0(%rax)
L361:
  movupd (%rax),%xmm2
  movupd -0x8(%rsi,%rax,1),%xmm7
  add    $0x10,%rax
  mulpd  %xmm7,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %r9,%rax
  jne    L361
  mov    0x24(%rsp),%r15d
  cmp    %r15d,0x2c(%rsp)
  je     L367
L322:
  movsd  (%rdx,%rdi,1),%xmm1
  mulsd  (%r8),%xmm1
  addsd  %xmm1,%xmm0
L367:
  mulsd  %xmm6,%xmm0
  cmp    $0x2,%ebx
  setne  %bpl
  cmp    %r11,0x10(%rsp)
  je     L101
  mov    0x50(%rsp),%rax
  movapd %xmm0,%xmm7
  unpcklpd %xmm7,%xmm7
  lea    (%rax,%r11,1),%r15
  mov    %r11,%rax
  test   %bpl,%bpl
  je     L101
  nopl   0x0(%rax)
L137:
  movupd 0x8(%rax,%rcx,1),%xmm2
  movupd (%rax),%xmm1
  add    $0x10,%rax
  mulpd  %xmm7,%xmm2
  subpd  %xmm2,%xmm1
  movups %xmm1,-0x10(%rax)
  cmp    %r15,%rax
  jne    L137
  mov    0x24(%rsp),%r15d
  cmp    %r15d,0x2c(%rsp)
  je     L392
  mulsd  (%r12),%xmm0
  lea    (%rdx,%r10,1),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax)
L392:
  add    %r14,%rdx
  lea    0x1(%r13),%eax
  add    %r14,%r11
  add    %r14,%rsi
  sub    %r14,%rcx
  cmp    %r13d,0x38(%rsp)
  je     L27
  mov    %eax,%r13d
  cmp    $0x2,%ebx
  jne    L174
L225:
  movapd %xmm3,%xmm0
  jmp    L322
  nopl   (%rax)
L101:
  mov    0x40(%rsp),%rax
  mov    0xa0(%rsp),%r15
  add    %rdx,%rax
  add    %r11,%r15
  nopl   0x0(%rax,%rax,1)
L302:
  movsd  0x8(%rax,%rcx,1),%xmm2
  movsd  (%rax),%xmm1
  add    $0x8,%rax
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  %xmm1,-0x8(%rax)
  cmp    %r15,%rax
  jne    L302
  jmp    L392
  nopw   0x0(%rax,%rax,1)
L27:
  mov    0x40(%rsp),%r11
  mov    0x60(%rsp),%rcx
  mov    %bpl,0x50(%rsp)
  mov    $0x1,%r13d
  mov    0xb8(%rsp),%r15
  mov    0xb0(%rsp),%rsi
  mov    0xc4(%rsp),%ebp
  mov    0x70(%rsp),%edx
  add    %r15,%r11
  add    %r15,%rsi
  mov    0xa8(%rsp),%r9
  mov    0x38(%rsp),%ebx
  sub    %r11,%rcx
  cmp    0xc0(%rsp),%ebp
  cmovg  0x24(%rsp),%r13d
  mov    0x88(%rsp),%rbp
  mov    %ebx,0x2c(%rsp)
  mov    %r13d,%eax
  mov    %r13d,%edi
  shr    %eax
  and    $0xfffffffe,%edi
  sub    $0x1,%eax
  add    %edi,%edx
  mov    %edi,0x70(%rsp)
  mov    0x68(%rsp),%edi
  shl    $0x4,%rax
  mov    0x70(%rsp),%ebx
  add    %rax,%r9
  cmpl   $0x2,0x28(%rsp)
  cmovne %edx,%edi
  add    $0x10,%rax
  movslq %edx,%rdx
  xor    %r12d,%r12d
  mov    %rax,0x38(%rsp)
  mov    0x90(%rsp),%rax
  lea    0x0(,%rdx,8),%r10
  mov    %r15,%rdx
  movslq %edi,%rdi
  shl    $0x3,%rdi
  add    $0x8,%rax
  lea    0x0(%rbp,%rdi,1),%r8
  add    %r10,%rbp
  cmpl   $0x2,0x28(%rsp)
  mov    %rax,0x68(%rsp)
  je     L371
  xchg   %ax,%ax
L253:
  mov    0x8(%rsp),%rax
  movapd %xmm3,%xmm0
  nopl   0x0(%rax)
L94:
  movupd -0x8(%rsi,%rax,1),%xmm2
  movupd (%rax),%xmm7
  add    $0x10,%rax
  mulpd  %xmm7,%xmm2
  addsd  %xmm2,%xmm0
  unpckhpd %xmm2,%xmm2
  addsd  %xmm2,%xmm0
  cmp    %r9,%rax
  jne    L94
  cmp    %r13d,%ebx
  je     L212
L70:
  movsd  (%rdx,%rdi,1),%xmm1
  mulsd  (%r8),%xmm1
  addsd  %xmm1,%xmm0
L212:
  mulsd  %xmm6,%xmm0
  cmp    %r11,0x10(%rsp)
  je     L142
  mov    0x38(%rsp),%rax
  movapd %xmm0,%xmm7
  cmpb   $0x0,0x50(%rsp)
  unpcklpd %xmm7,%xmm7
  lea    (%rax,%r11,1),%r15
  mov    %r11,%rax
  je     L142
L224:
  movupd 0x8(%rcx,%rax,1),%xmm2
  movupd (%rax),%xmm1
  add    $0x10,%rax
  mulpd  %xmm7,%xmm2
  subpd  %xmm2,%xmm1
  movups %xmm1,-0x10(%rax)
  cmp    %rax,%r15
  jne    L224
  cmp    %r13d,%ebx
  je     L252
  mulsd  0x0(%rbp),%xmm0
  lea    (%rdx,%r10,1),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rax)
L252:
  add    %r14,%rdx
  lea    0x1(%r12),%eax
  add    %r14,%r11
  add    %r14,%rsi
  sub    %r14,%rcx
  cmp    %r12d,0x2c(%rsp)
  je     L360
  cmpl   $0x2,0x28(%rsp)
  mov    %eax,%r12d
  jne    L253
L371:
  movapd %xmm3,%xmm0
  jmp    L70
  nopl   0x0(%rax)
L142:
  mov    0x40(%rsp),%rax
  mov    0x68(%rsp),%r15
  add    %rdx,%rax
  add    %r11,%r15
L113:
  movsd  0x8(%rcx,%rax,1),%xmm2
  movsd  (%rax),%xmm1
  add    $0x8,%rax
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movsd  %xmm1,-0x8(%rax)
  cmp    %rax,%r15
  jne    L113
  jmp    L252
L327:
  mov    0x8(%rsp),%rax
  add    $0x8,%rax
  mov    %rax,0x10(%rsp)
  lea    -0x1(%r12),%eax
  mov    %eax,0x4c(%rsp)
  jmp    L360
L202:
  mov    0x8(%rsp),%rax
  add    $0x8,%rax
  mov    %rax,0x10(%rsp)
  lea    -0x1(%r12),%eax
  mov    %eax,0x4c(%rsp)
  jmp    L13
L405:
  mov    $0xffffffff,%eax
  jmp    L139
L306:
  movapd %xmm1,%xmm0
  call   <sqrt@plt>
  mov    D26(%rip),%rax
  pxor   %xmm3,%xmm3
  movq   D38(%rip),%xmm5
  movapd %xmm0,%xmm1
  movq   %rax,%xmm4
  jmp    L60
```

tool_calls: name: extract_symbolic_data
arguments: tag: D26
type: f32[2]
name: extract_symbolic_data
arguments: tag: D38
type: f32[2]

tool: {"D26": [0.0,-1.875]}

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

assistant: ```c
int Hessenberg_Form_Orthogonal(double *A, double *U, int n) {
  int i, k, col;
  double *u;
  double *p_row, *psubdiag;
  double *pA, *pU;
  double sss;
  double scale;
  double innerproduct;
  if (n <= 2)
    return 0;
  u = (double *)malloc(n * sizeof(double));
  if (u == NULL)
    return -1;
  for (psubdiag = A + n, col = 0; col < (n - 2); psubdiag += (n + 1), col++) {
    for (pA = psubdiag, sss = 0.0, i = col + 1; i < n; pA += n, i++)
      sss += *pA * *pA;
    if (sss == 0.0)
      continue;
    sss = sqrt(sss);
    if (*psubdiag >= 0.0)
      sss = -sss;
    u[col + 1] = *psubdiag - sss;
    *psubdiag = sss;
    for (pA = psubdiag + n, i = col + 2; i < n; pA += n, i++) {
      u[i] = *pA;
      *pA = 0.0;
    }
    scale = -1.0 / (sss * u[col + 1]);
    for (p_row = psubdiag - col, i = col + 1; i < n; i++) {
      pA = A + n * (col + 1) + i;
      for (innerproduct = 0.0, k = col + 1; k < n; pA += n, k++)
        innerproduct += u[k] * *pA;
      innerproduct *= scale;
      for (pA = p_row + i, k = col + 1; k < n; pA += n, k++)
        *pA -= u[k] * innerproduct;
    }
    for (p_row = A, i = 0; i < n; p_row += n, i++) {
      for (innerproduct = 0.0, k = col + 1; k < n; k++)
        innerproduct += u[k] * *(p_row + k);
      innerproduct *= scale;
      for (k = col + 1; k < n; k++)
        *(p_row + k) -= u[k] * innerproduct;
    }
    for (i = 0, pU = U; i < n; pU += n, i++) {
      for (innerproduct = 0.0, k = col + 1; k < n; k++)
        innerproduct += u[k] * *(pU + k);
      innerproduct *= scale;
      for (k = col + 1; k < n; k++)
        *(pU + k) -= u[k] * innerproduct;
    }
  }
  free(u);
  return 0;
}
```