user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<Hessenberg_Form_Orthogonal>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %rsi,-0x60(%rbp)
  mov    %edx,-0x64(%rbp)
  cmpl   $0x2,-0x64(%rbp)
  jg     L191
  mov    $0x0,%eax
  jmp    L45
L191:
  mov    -0x64(%rbp),%eax
  cltq
  shl    $0x3,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0x10(%rbp)
  cmpq   $0x0,-0x10(%rbp)
  jne    L62
  mov    $0xffffffff,%eax
  jmp    L45
L62:
  mov    -0x64(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x38(%rbp)
  movl   $0x0,-0x44(%rbp)
  jmp    L199
L243:
  mov    -0x38(%rbp),%rax
  mov    %rax,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L245
L143:
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x20(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x64(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x30(%rbp)
  addl   $0x1,-0x4c(%rbp)
L245:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L143
  pxor   %xmm0,%xmm0
  ucomisd -0x20(%rbp),%xmm0
  jp     L160
  pxor   %xmm0,%xmm0
  ucomisd -0x20(%rbp),%xmm0
  je     L306
L160:
  mov    -0x20(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x20(%rbp)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jb     L121
  movsd  -0x20(%rbp),%xmm0
  movq   D1(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
L121:
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x44(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  subsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x64(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x2,%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L22
L126:
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rax,%rdx
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,(%rdx)
  mov    -0x30(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x64(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x30(%rbp)
  addl   $0x1,-0x4c(%rbp)
L22:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L126
  mov    -0x44(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  movsd  D17(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x44(%rbp),%eax
  cltq
  shl    $0x3,%rax
  neg    %rax
  mov    %rax,%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L54
L270:
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  imul   -0x64(%rbp),%eax
  movslq %eax,%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  add    %rdx,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x30(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x48(%rbp)
  jmp    L7
L319:
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x64(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x30(%rbp)
  addl   $0x1,-0x48(%rbp)
L7:
  mov    -0x48(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L319
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x48(%rbp)
  jmp    L17
L39:
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x18(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x30(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x64(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x30(%rbp)
  addl   $0x1,-0x48(%rbp)
L17:
  mov    -0x48(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L39
  addl   $0x1,-0x4c(%rbp)
L54:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L270
  mov    -0x58(%rbp),%rax
  mov    %rax,-0x40(%rbp)
  movl   $0x0,-0x4c(%rbp)
  jmp    L57
L249:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x48(%rbp)
  jmp    L30
L72:
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  addl   $0x1,-0x48(%rbp)
L30:
  mov    -0x48(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L72
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x48(%rbp)
  jmp    L314
L48:
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x18(%rbp),%xmm1
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x48(%rbp)
L314:
  mov    -0x48(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L48
  mov    -0x64(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x40(%rbp)
  addl   $0x1,-0x4c(%rbp)
L57:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L249
  movl   $0x0,-0x4c(%rbp)
  mov    -0x60(%rbp),%rax
  mov    %rax,-0x28(%rbp)
  jmp    L26
L51:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x48(%rbp)
  jmp    L232
L294:
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  addl   $0x1,-0x48(%rbp)
L232:
  mov    -0x48(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L294
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x48(%rbp)
  jmp    L252
L150:
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x18(%rbp),%xmm1
  mov    -0x48(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x48(%rbp)
L252:
  mov    -0x48(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L150
  mov    -0x64(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x28(%rbp)
  addl   $0x1,-0x4c(%rbp)
L26:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L51
  jmp    L299
L306:
  nop
L299:
  mov    -0x64(%rbp),%eax
  cltq
  add    $0x1,%rax
  shl    $0x3,%rax
  add    %rax,-0x38(%rbp)
  addl   $0x1,-0x44(%rbp)
L199:
  mov    -0x64(%rbp),%eax
  sub    $0x2,%eax
  cmp    %eax,-0x44(%rbp)
  jl     L243
  mov    -0x10(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    $0x0,%eax
L45:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D1
type_infered: f32[2]
name: retrieve_data
arguments: data_label: D17
type_infered: f32[2]

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

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

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;
}
```