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

```asm
<francis_iter>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  movslq %esi,%rbx
  mov    %ebx,%r12d
  imul   $0x1e,%ebx,%edi
  lea    -0x1(%rbx),%r14d
  imul   %ebx,%r12d
  sub    $0x88,%rsp
  mov    %edi,0x54(%rsp)
  lea    0x0(,%r12,8),%edi
  mov    %ebx,0x20(%rsp)
  mov    %rdx,(%rsp)
  call   <malloc@plt>
  test   %r12d,%r12d
  mov    (%rsp),%r10
  mov    %rax,%r15
  je     L199
  lea    -0x1(%r12),%eax
  cmp    $0x2,%eax
  jbe    L460
  lea    0x8(%rbp),%rax
  cmp    %rax,%r10
  je     L460
  mov    %r12d,%edx
  xor    %eax,%eax
  shr    %edx
  shl    $0x4,%rdx
  nopl   0x0(%rax)
L502:
  movupd 0x0(%rbp,%rax,1),%xmm5
  movups %xmm5,(%r10,%rax,1)
  add    $0x10,%rax
  cmp    %rdx,%rax
  jne    L502
  mov    %r12d,%eax
  and    $0xfffffffe,%eax
  and    $0x1,%r12d
  je     L199
  movsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,(%r10,%rax,8)
L199:
  mov    0x20(%rsp),%ebp
  mov    %r10,%rdi
  mov    %r10,(%rsp)
  mov    %ebp,%esi
  call   <hessenberg@plt>
  lea    0x1(%rbp),%eax
  mov    (%rsp),%rbp
  movl   $0x0,0x50(%rsp)
  mov    %eax,0x70(%rsp)
  lea    0x0(,%rbx,8),%rax
  pxor   %xmm1,%xmm1
  mov    %rax,0x8(%rsp)
  add    $0x8,%rax
  mov    %rax,0x30(%rsp)
L478:
  mov    0x54(%rsp),%edi
  cmp    %edi,0x50(%rsp)
  setl   %dl
  cmpl   $0x1,0x20(%rsp)
  jne    L2
  nopl   0x0(%rax)
L68:
  cmp    $0x1,%r14d
  jle    L407
  test   %dl,%dl
  jne    L243
  jmp    L407
  nopl   0x0(%rax,%rax,1)
L230:
  lea    (%rbx,%r14,1),%eax
  cltq
  ucomisd -0x10(%rbp,%rax,8),%xmm1
  jp     L408
  jne    L408
  sub    $0x2,%r14d
  cmp    $0x1,%r14d
  jle    L74
L243:
  lea    (%r14,%r14,1),%eax
  lea    -0x1(%r14),%ebx
  cltq
  ucomisd -0x8(%rbp,%rax,8),%xmm1
  jp     L230
  jne    L230
  mov    %ebx,%r14d
  cmp    $0x1,%r14d
  jg     L243
L74:
  test   %r14d,%r14d
  je     L68
L121:
  mov    $0x1,%r14d
L244:
  mov    %r14d,%eax
  xor    %ebx,%ebx
  jmp    L261
  nopl   0x0(%rax)
L408:
  mov    $0x2,%eax
L517:
  mov    %eax,%r12d
  imul   %ebx,%eax
  neg    %r12d
  movslq %r12d,%r12
  cltq
  shl    $0x3,%r12
  lea    -0x8(%rbp,%rax,8),%r13
  jmp    L145
  cs nopw 0x0(%rax,%rax,1)
L372:
  add    %r12,%r13
  sub    $0x1,%ebx
  je     L244
L145:
  movsd  0x0(%r13),%xmm0
  call   <fabs@plt>
  test   %eax,%eax
  jne    L372
  mov    %r14d,%eax
  sub    %ebx,%eax
  cmp    %r14d,%ebx
  jg     L102
L261:
  mov    0x20(%rsp),%edx
  add    $0x1,%eax
  movslq %ebx,%rsi
  mov    0x8(%rsp),%rdi
  mov    %eax,(%rsp)
  cltq
  mov    %r15,%r8
  imul   %ebx,%edx
  shl    $0x3,%rax
  mov    %r15,0x68(%rsp)
  sub    %rax,%rdi
  mov    %rax,0x18(%rsp)
  mov    %r14d,%eax
  sub    %ebx,%eax
  mov    %rdi,0x10(%rsp)
  movslq %edx,%r13
  lea    0x1(%rax),%ecx
  mov    %eax,0x2c(%rsp)
  mov    $0x1,%eax
  add    %rsi,%r13
  mov    %edx,0x74(%rsp)
  shl    $0x3,%r13
  mov    %rsi,0x40(%rsp)
  lea    0x0(%rbp,%r13,1),%r11
  sub    %r15,%r13
  add    %rbp,%r13
  cmp    %ebx,%r14d
  mov    %r11,%r9
  mov    %r11,0x58(%rsp)
  cmovge %ecx,%eax
  mov    %r13,0x78(%rsp)
  mov    %ebx,%r11d
  mov    %r9,%r10
  mov    %eax,%edi
  mov    %eax,0x28(%rsp)
  shr    %eax
  sub    $0x1,%eax
  and    $0xfffffffe,%edi
  shl    $0x4,%rax
  mov    %edi,0x24(%rsp)
  lea    (%rdi,%rbx,1),%r12d
  mov    %edx,%edi
  mov    %rax,%rcx
  mov    %rax,0x60(%rsp)
  mov    %rsi,%rax
  xor    %esi,%esi
  add    $0x10,%rcx
  neg    %rax
  mov    %rcx,0x38(%rsp)
  shl    $0x3,%rax
  mov    %rax,0x48(%rsp)
  mov    %r13,%rax
  xchg   %ax,%ax
L117:
  cmpl   $0x2,0x2c(%rsp)
  lea    0x8(%r10),%rcx
  seta   %r9b
  cmp    %ebx,%r14d
  setge  %dl
  and    %edx,%r9d
  mov    %r8,%rdx
  cmp    %rcx,%r8
  je     L434
  test   %r9b,%r9b
  je     L434
  mov    0x38(%rsp),%rcx
  add    %r8,%rcx
  nopl   0x0(%rax,%rax,1)
L366:
  movupd (%rdx,%rax,1),%xmm3
  add    $0x10,%rdx
  movups %xmm3,-0x10(%rdx)
  cmp    %rcx,%rdx
  jne    L366
  mov    0x28(%rsp),%edx
  cmp    %edx,0x24(%rsp)
  je     L233
  lea    (%r12,%rdi,1),%edx
  movslq %edx,%rdx
  movsd  0x0(%rbp,%rdx,8),%xmm0
  lea    (%r12,%rsi,1),%edx
  sub    %ebx,%edx
  movslq %edx,%rdx
  movsd  %xmm0,(%r15,%rdx,8)
L233:
  add    $0x1,%r11d
  add    0x8(%rsp),%r10
  add    0x18(%rsp),%r8
  add    0x20(%rsp),%edi
  add    (%rsp),%esi
  add    0x10(%rsp),%rax
  cmp    %r11d,%r14d
  jge    L117
  mov    0x68(%rsp),%r10
  mov    (%rsp),%esi
  mov    %r15,%rdi
  mov    %r9b,0x2c(%rsp)
  mov    0x74(%rsp),%edx
  mov    0x78(%rsp),%r13
  mov    %r10,0x38(%rsp)
  mov    %edx,0x68(%rsp)
  call   <francisQR@plt>
  mov    0x60(%rsp),%r9
  mov    0x68(%rsp),%edx
  mov    %ebx,%edi
  mov    0x38(%rsp),%r10
  mov    0x58(%rsp),%r8
  xor    %ecx,%ecx
  pxor   %xmm1,%xmm1
  add    $0x10,%r9
  nopl   0x0(%rax)
L510:
  lea    0x8(%r10),%rax
  cmp    %r8,%rax
  je     L164
  cmpb   $0x0,0x2c(%rsp)
  je     L164
  lea    (%r9,%r10,1),%rsi
  mov    %r10,%rax
  nop
L85:
  movupd (%rax),%xmm4
  movups %xmm4,0x0(%r13,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%rsi
  jne    L85
  mov    0x28(%rsp),%esi
  cmp    %esi,0x24(%rsp)
  je     L272
  lea    (%r12,%rcx,1),%eax
  sub    %ebx,%eax
  cltq
  movsd  (%r15,%rax,8),%xmm0
  lea    (%r12,%rdx,1),%eax
  cltq
  movsd  %xmm0,0x0(%rbp,%rax,8)
L272:
  add    $0x1,%edi
  add    0x18(%rsp),%r10
  add    0x8(%rsp),%r8
  add    (%rsp),%ecx
  add    0x20(%rsp),%edx
  add    0x10(%rsp),%r13
  cmp    %edi,%r14d
  jge    L510
L71:
  cmp    %r14d,%ebx
  jge    L47
  mov    0x20(%rsp),%edi
  lea    0x1(%rbx),%edx
  movslq %ebx,%rsi
  mov    %r15,(%rsp)
  mov    %rbp,0x10(%rsp)
  mov    %edx,%r15d
  mov    %edi,%eax
  imul   %edx,%eax
  movslq %eax,%rcx
  sub    %edi,%eax
  add    %rcx,%rsi
  movslq %eax,%r12
  lea    0x0(%rbp,%rsi,8),%r13
  sub    %rcx,%r12
  mov    %r13,%rbp
  jmp    L435
  nopl   0x0(%rax)
L166:
  add    $0x1,%r15d
L435:
  movsd  0x0(%rbp),%xmm0
  call   <fabs@plt>
  movsd  0x0(%rbp,%r12,8),%xmm0
  mov    %eax,%r13d
  call   <fabs@plt>
  movsd  0x8(%rbp),%xmm0
  mov    %eax,%ebx
  call   <fabs@plt>
  pxor   %xmm0,%xmm0
  pxor   %xmm2,%xmm2
  pxor   %xmm1,%xmm1
  add    %eax,%ebx
  cvtsi2sd %r13d,%xmm2
  cvtsi2sd %ebx,%xmm0
  mulsd  D4(%rip),%xmm0
  comisd %xmm2,%xmm0
  jb     L321
  movsd  %xmm1,0x0(%rbp)
L321:
  add    0x30(%rsp),%rbp
  cmp    %r14d,%r15d
  jne    L166
  mov    (%rsp),%r15
  mov    0x10(%rsp),%rbp
L47:
  addl   $0x1,0x50(%rsp)
  jmp    L478
L407:
  mov    %r15,%rdi
  call   <free@plt>
  mov    0x54(%rsp),%edi
  xor    %eax,%eax
  cmp    %edi,0x50(%rsp)
  setne  %al
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L164:
  mov    0x48(%rsp),%rax
  mov    %r10,%rsi
  lea    (%r8,%rax,1),%r11
  mov    0x40(%rsp),%rax
L456:
  movsd  (%rsi),%xmm0
  add    $0x8,%rsi
  movsd  %xmm0,(%r11,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r14d
  jge    L456
  jmp    L272
L434:
  mov    0x48(%rsp),%rcx
  mov    0x40(%rsp),%rdx
  lea    (%r10,%rcx,1),%r13
  mov    %r8,%rcx
  nopl   (%rax)
L520:
  movsd  0x0(%r13,%rdx,8),%xmm0
  add    $0x1,%rdx
  add    $0x8,%rcx
  movsd  %xmm0,-0x8(%rcx)
  cmp    %edx,%r14d
  jge    L520
  jmp    L233
L2:
  mov    0x70(%rsp),%ecx
  mov    0x20(%rsp),%esi
  cs nopw 0x0(%rax,%rax,1)
L497:
  cmp    $0x1,%r14d
  jle    L407
  test   %dl,%dl
  jne    L177
  jmp    L407
  nopl   0x0(%rax,%rax,1)
L10:
  mov    %esi,%eax
  imul   %ebx,%eax
  add    %r14d,%eax
  cltq
  ucomisd -0x10(%rbp,%rax,8),%xmm1
  jp     L337
  jne    L337
  sub    $0x2,%r14d
  cmp    $0x1,%r14d
  jle    L421
L177:
  mov    %r14d,%eax
  lea    -0x1(%r14),%ebx
  imul   %ecx,%eax
  cltq
  ucomisd -0x8(%rbp,%rax,8),%xmm1
  jp     L10
  jne    L10
  mov    %ebx,%r14d
  cmp    $0x1,%r14d
  jg     L177
L421:
  test   %r14d,%r14d
  je     L497
  jmp    L121
  nop
L337:
  mov    0x70(%rsp),%eax
  jmp    L517
L460:
  xor    %eax,%eax
L63:
  movsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,(%r10,%rax,8)
  add    $0x1,%rax
  cmp    %r12,%rax
  jne    L63
  jmp    L199
L102:
  lea    0x1(%rax),%esi
  mov    %r15,%rdi
  call   <francisQR@plt>
  pxor   %xmm1,%xmm1
  jmp    L71
```

tool_calls: name: read_data
arguments: data_label: D4
type_infered: byte[8]

tool: {"D4": [0x11,0xea,0x2d,0x81,0x99,0x97,0x71,0x3d]}

assistant: ```c
int francis_iter(double *A, int N, double *H) {
  int success, brkpoint;
  int i, j, it, p, q, t, u;
  double *temp;
  success = 0;
  brkpoint = 30 * N;
  it = 0;
  p = N - 1;
  temp = (double *)malloc(sizeof(double) * N * N);
  for (i = 0; i < N * N; ++i) {
    H[i] = A[i];
  }
  hessenberg(H, N);
  while (p > 1 && it < brkpoint) {
    while (p > 1 && (H[N * p + p - 1] == 0 || H[N * (p - 1) + p - 2] == 0)) {
      if (H[N * p + p - 1] == 0) {
        p--;
      } else if (H[N * (p - 1) + p - 2] == 0) {
        p = p - 2;
      }
    }
    if (p > 0) {
      q = p - 1;
      while (q > 0 && fabs(H[N * q + q - 1]) != 0) {
        q--;
      }
      for (i = q; i <= p; i++) {
        t = i * N;
        u = (i - q) * (p - q + 1);
        for (j = q; j <= p; j++) {
          temp[u + j - q] = H[t + j];
        }
      }
      francisQR(temp, p - q + 1);
      for (i = q; i <= p; i++) {
        t = i * N;
        u = (i - q) * (p - q + 1);
        for (j = q; j <= p; j++) {
          H[t + j] = temp[u + j - q];
        }
      }
      for (i = q; i <= p - 1; ++i) {
        if (fabs(H[(i + 1) * N + i]) <= 1e-12 * (fabs(H[i * N + i]) + fabs(H[(i + 1) * N + i + 1]))) {
          H[(i + 1) * N + i] = 0.;
        }
      }
      it++;
    }
  }
  if (it == brkpoint) {
    success = 0;
  } else {
    success = 1;
  }
  free(temp);
  return success;
}
```