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

```asm
<lob_bnd_ext>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %rcx,%r13
  push   %r12
  mov    %rdx,%r12
  push   %rbp
  push   %rbx
  mov    0x48(%rsp),%ebx
  mov    0x38(%rsp),%r10d
  mov    0x40(%rsp),%rcx
  mov    %rsi,-0x18(%rsp)
  mov    0x50(%rsp),%rsi
  mov    %rdi,-0x10(%rsp)
  test   %ebx,%ebx
  je     L205
  xor    %eax,%eax
  xor    %edx,%edx
  nopl   0x0(%rax)
L145:
  mov    %eax,%edi
  add    $0x1,%edx
  movq   $0x0,(%rsi,%rdi,8)
  lea    0x1(%rax),%edi
  add    $0x2,%eax
  movq   $0x0,(%rsi,%rdi,8)
  cmp    %edx,%ebx
  jne    L145
L205:
  test   %r9d,%r9d
  je     L227
  lea    -0x1(%rbx),%r15d
  mov    %rcx,%rax
  xor    %ebp,%ebp
  xor    %r14d,%r14d
  mov    %r15d,-0x4(%rsp)
  movsd  D7(%rip),%xmm4
  add    $0x1,%r15
  shl    $0x4,%r15
  nopl   0x0(%rax)
L167:
  mov    %ebp,%edx
  movsd  (%r12,%rdx,8),%xmm2
  lea    0x1(%rbp),%edx
  movsd  (%r12,%rdx,8),%xmm3
  test   %ebx,%ebx
  je     L85
  lea    (%rax,%r15,1),%r11
  xor    %edx,%edx
  nopl   0x0(%rax,%rax,1)
L172:
  movsd  (%rax),%xmm0
  addsd  0x8(%rax),%xmm0
  mov    %edx,%edi
  movapd %xmm2,%xmm1
  lea    (%rsi,%rdi,8),%rdi
  add    $0x10,%rax
  mulsd  %xmm4,%xmm0
  mulsd  %xmm0,%xmm1
  addsd  (%rdi),%xmm1
  mulsd  %xmm3,%xmm0
  movsd  %xmm1,(%rdi)
  lea    0x1(%rdx),%edi
  add    $0x2,%edx
  lea    (%rsi,%rdi,8),%rdi
  addsd  (%rdi),%xmm0
  movsd  %xmm0,(%rdi)
  cmp    %r11,%rax
  jne    L172
L85:
  add    $0x1,%r14d
  add    $0x2,%ebp
  cmp    %r14d,%r9d
  jne    L167
  lea    -0x1(%r10),%edx
  test   %ebx,%ebx
  je     L27
L216:
  lea    -0x1(%r10),%r14d
  mov    -0x10(%rsp),%rax
  xor    %ebp,%ebp
  xor    %r12d,%r12d
  mov    %r14,%rdx
  add    $0x1,%r14
  shl    $0x4,%r14
  nop
L261:
  mov    %ebp,%edi
  movsd  (%rsi,%rdi,8),%xmm1
  lea    0x1(%rbp),%edi
  movsd  (%rsi,%rdi,8),%xmm2
  test   %r10d,%r10d
  je     L229
  mov    %r13,%rdi
  lea    (%rax,%r14,1),%r11
  nopl   0x0(%rax,%rax,1)
L74:
  movsd  (%rdi),%xmm0
  add    $0x10,%rax
  add    $0x8,%rdi
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rax)
  movsd  %xmm0,-0x8(%rax)
  cmp    %r11,%rax
  jne    L74
L229:
  add    $0x1,%r12d
  add    $0x2,%ebp
  cmp    %r12d,%ebx
  jne    L261
  test   %r9d,%r9d
  je     L176
  lea    -0x1(%rbx),%eax
  mov    %eax,-0x4(%rsp)
L27:
  mov    -0x18(%rsp),%rbp
  lea    -0x1(%r9),%eax
  mov    %edx,%r11d
  pxor   %xmm4,%xmm4
  mov    -0x4(%rsp),%r12d
  lea    0x0(,%r10,4),%r13d
  add    $0x1,%r11
  lea    0x8(%rbp,%rax,8),%rax
  shl    $0x3,%r13
  shl    $0x4,%r11
  mov    %rax,-0x18(%rsp)
  add    $0x1,%r12
  shl    $0x4,%r12
  nopw   0x0(%rax,%rax,1)
L41:
  movsd  0x0(%rbp),%xmm5
  test   %ebx,%ebx
  je     L44
  mov    -0x10(%rsp),%rax
  lea    (%rcx,%r12,1),%r9
  xor    %edi,%edi
  nopl   0x0(%rax,%rax,1)
L139:
  lea    0x1(%rdi),%edx
  movsd  (%rcx),%xmm1
  add    $0x10,%rcx
  movsd  -0x8(%rcx),%xmm2
  movsd  (%rsi,%rdx,8),%xmm0
  mov    %edi,%edx
  mulsd  %xmm5,%xmm0
  addsd  (%rsi,%rdx,8),%xmm0
  subsd  %xmm0,%xmm1
  subsd  %xmm0,%xmm2
  comisd %xmm4,%xmm1
  jae    L274
  comisd %xmm2,%xmm4
  jae    L154
  lea    (%rax,%r11,1),%r14
  xor    %edx,%edx
  test   %r10d,%r10d
  je     L83
  nopw   0x0(%rax,%rax,1)
L25:
  lea    0x3(%rdx),%r15d
  add    $0x10,%rax
  movsd  (%r8,%r15,8),%xmm0
  mov    %edx,%r15d
  add    $0x4,%edx
  movsd  (%r8,%r15,8),%xmm6
  movapd %xmm0,%xmm3
  mulsd  %xmm1,%xmm3
  movapd %xmm6,%xmm7
  mulsd  %xmm2,%xmm7
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm6
  addsd  %xmm7,%xmm3
  addsd  -0x10(%rax),%xmm3
  addsd  %xmm6,%xmm0
  movsd  %xmm3,-0x10(%rax)
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %r14,%rax
  jne    L25
L83:
  add    $0x2,%edi
  cmp    %r9,%rcx
  jne    L139
L44:
  add    %r13,%r8
  add    $0x8,%rbp
  cmp    %rbp,-0x18(%rsp)
  jne    L41
L176:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L274:
  test   %r10d,%r10d
  je     L83
  lea    (%rax,%r11,1),%r14
  xor    %edx,%edx
  nopl   0x0(%rax,%rax,1)
L57:
  lea    0x1(%rdx),%r15d
  add    $0x10,%rax
  movsd  (%r8,%r15,8),%xmm0
  mov    %edx,%r15d
  movsd  (%r8,%r15,8),%xmm3
  lea    0x3(%rdx),%r15d
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm0
  addsd  -0x10(%rax),%xmm0
  movsd  %xmm0,-0x10(%rax)
  movsd  (%r8,%r15,8),%xmm0
  lea    0x2(%rdx),%r15d
  add    $0x4,%edx
  movsd  (%r8,%r15,8),%xmm3
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm3
  addsd  %xmm3,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%r14
  jne    L57
  add    $0x2,%edi
  cmp    %r9,%rcx
  jne    L139
  jmp    L44
L154:
  test   %r10d,%r10d
  je     L83
  lea    (%rax,%r11,1),%r15
  xor    %edx,%edx
  nop
L185:
  lea    0x3(%rdx),%r14d
  add    $0x10,%rax
  movsd  (%r8,%r14,8),%xmm0
  lea    0x2(%rdx),%r14d
  movsd  (%r8,%r14,8),%xmm3
  lea    0x1(%rdx),%r14d
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm0
  addsd  -0x10(%rax),%xmm0
  movsd  %xmm0,-0x10(%rax)
  movsd  (%r8,%r14,8),%xmm0
  mov    %edx,%r14d
  add    $0x4,%edx
  movsd  (%r8,%r14,8),%xmm3
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm3
  addsd  %xmm3,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %r15,%rax
  jne    L185
  add    $0x2,%edi
  cmp    %r9,%rcx
  jne    L139
  jmp    L44
L227:
  test   %ebx,%ebx
  jne    L216
  jmp    L176
```

tool_calls: name: parse_data
arguments: tag: D7
type_infered: byte[8]

tool: {"D7": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

assistant: ```c
void lob_bnd_ext(double *restrict b_, const double *restrict z, const double *restrict Q, const double *restrict h, const double *restrict lbnp, unsigned n, unsigned m, const double *restrict br_, unsigned mr, double *restrict a) {
  const double *restrict br = br_;
  double *restrict b = b_;
  unsigned i, j, k;
  for (i = 0; i < mr; ++i)
    a[2 * i + 1] = a[2 * i + 0] = 0;
  for (j = 0; j < n; ++j) {
    double t, q0 = Q[2 * j], q1 = Q[2 * j + 1];
    for (i = 0; i < mr; ++i)
      t = (br[0] + br[1]) / 2, br += 2, a[2 * i] += q0 * t, a[2 * i + 1] += q1 * t;
  }
  for (i = 0; i < mr; ++i) {
    double a0 = a[2 * i], a1 = a[2 * i + 1];
    for (k = 0; k < m; ++k)
      b[1] = b[0] = a0 + a1 * h[k], b += 2;
  }
  br = br_;
  for (j = 0; j < n; ++j, lbnp += 4 * m) {
    double zj = z[j];
    b = b_;
    for (i = 0; i < mr; ++i) {
      double t = a[2 * i] + a[2 * i + 1] * zj;
      double w0 = *br++ - t;
      double w1 = *br++ - t;
      if (w0 >= 0)
        for (k = 0; k < m; ++k)
          *b++ += w0 * lbnp[4 * k + 1] + w1 * lbnp[4 * k + 0], *b++ += w1 * lbnp[4 * k + 3] + w0 * lbnp[4 * k + 2];
      else if (w1 <= 0)
        for (k = 0; k < m; ++k)
          *b++ += w0 * lbnp[4 * k + 3] + w1 * lbnp[4 * k + 2], *b++ += w1 * lbnp[4 * k + 1] + w0 * lbnp[4 * k + 0];
      else
        for (k = 0; k < m; ++k)
          *b++ += w0 * lbnp[4 * k + 3] + w1 * lbnp[4 * k + 0], *b++ += w1 * lbnp[4 * k + 3] + w0 * lbnp[4 * k + 0];
    }
  }
}
```