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

```asm
<solv>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %rsi,%r12
  mov    $0x8,%esi
  push   %rbp
  mov    %rdi,%rbp
  mov    %edx,%edi
  push   %rbx
  mov    %edx,%ebx
  sub    $0x68,%rsp
  call   <calloc@plt>
  mov    %rax,0x50(%rsp)
  test   %ebx,%ebx
  jle    L115
  movslq %ebx,%rdi
  mov    %rax,%r13
  mov    %rbp,0x30(%rsp)
  pxor   %xmm3,%xmm3
  shl    $0x3,%rdi
  lea    0x8(%r13),%rcx
  movl   $0x0,0x2c(%rsp)
  movq   D15(%rip),%xmm2
  lea    0x8(%rdi),%rsi
  mov    %r12,0x38(%rsp)
  mov    %rbp,%r15
  movapd %xmm3,%xmm4
  mov    %rsi,0x20(%rsp)
  mov    %ebx,%esi
  mov    %rcx,%r12
  movsd  D24(%rip),%xmm5
  mov    %rsi,0x10(%rsp)
  shl    $0x3,%rsi
  movsd  D26(%rip),%xmm6
  mov    $0x1,%r14d
  mov    %rsi,0x40(%rsp)
  add    %rax,%rsi
  lea    0x0(%rbp,%rdi,1),%rax
  mov    %rax,0x58(%rsp)
  sub    %r13,%rax
  add    $0x8,%rax
  mov    %rsi,0x18(%rsp)
  mov    %rbp,%rsi
  mov    %rax,0x48(%rsp)
L170:
  movsd  (%r15),%xmm7
  lea    -0x1(%r14),%r8d
  mov    %r14d,%ebp
  movapd %xmm7,%xmm1
  andpd  %xmm2,%xmm1
  cmp    %r14d,%ebx
  jle    L60
  mov    %r15,%rdx
  mov    %r8d,%ecx
  mov    %r14d,%eax
  nop
L179:
  add    %rdi,%rdx
  movsd  (%rdx),%xmm0
  andpd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  maxsd  %xmm1,%xmm0
  cmova  %eax,%ecx
  add    $0x1,%eax
  movapd %xmm0,%xmm1
  cmp    %eax,%ebx
  jne    L179
L233:
  maxsd  %xmm1,%xmm4
  movapd %xmm4,%xmm0
  mulsd  %xmm5,%xmm0
  comisd %xmm1,%xmm0
  ja     L214
  cmp    %r8d,%ecx
  jne    L146
  movapd %xmm6,%xmm1
  divsd  %xmm7,%xmm1
  cmp    0x10(%rsp),%r14
  je     L98
L159:
  mov    %r15,%rdx
  mov    %ebp,%eax
  nopl   0x0(%rax,%rax,1)
L210:
  add    %rdi,%rdx
  mov    %eax,%r11d
  add    $0x1,%eax
  movsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx)
  cmp    %eax,%ebx
  jne    L210
  lea    0x8(%rsi),%rcx
  mov    0x18(%rsp),%r9
  add    0x20(%rsp),%r15
  mov    %r13,%r8
  mov    %rcx,0x8(%rsp)
  mov    %rcx,%rdx
L22:
  movsd  (%rdx),%xmm0
  mov    %r8,%rcx
  add    $0x8,%r8
  add    %rdi,%rdx
  movsd  %xmm0,(%rcx)
  cmp    %r8,%r9
  jne    L22
  cmp    $0x1,%eax
  je     L216
  mov    0x48(%rsp),%rax
  sub    0x8(%rsp),%rax
  mov    %r12,%r10
  mov    $0x1,%r9d
  lea    (%rsi,%rax,1),%rcx
  nopl   (%rax)
L42:
  cmp    %r9d,%ebp
  mov    %r9d,%eax
  movapd %xmm3,%xmm1
  cmovle %ebp,%eax
  sub    $0x1,%eax
  lea    (%r12,%rax,8),%rsi
  mov    %r13,%rax
  nop
L130:
  mov    %rax,%rdx
  movsd  (%rax,%rcx,1),%xmm0
  add    $0x8,%rax
  mulsd  (%rdx),%xmm0
  addsd  %xmm0,%xmm1
  cmp    %rax,%rsi
  jne    L130
  movsd  (%r10),%xmm0
  lea    0x1(%r9),%eax
  add    $0x8,%r10
  add    %rdi,%rcx
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r10)
  cmp    %r9d,%r11d
  je     L216
  mov    %eax,%r9d
  jmp    L42
  nopl   (%rax)
L216:
  mov    0x8(%rsp),%rdx
  mov    %r13,%rax
  nopl   0x0(%rax,%rax,1)
L165:
  mov    %rax,%rcx
  add    $0x8,%rax
  movsd  (%rcx),%xmm0
  movsd  %xmm0,(%rdx)
  add    %rdi,%rdx
  cmp    %r8,%rax
  jne    L165
  add    %ebx,0x2c(%rsp)
  mov    0x8(%rsp),%rsi
  add    $0x1,%r14
  jmp    L170
  nopl   0x0(%rax)
L146:
  mov    0x38(%rsp),%rdx
  movslq %ecx,%rax
  imul   %ebx,%ecx
  mov    0x30(%rsp),%r10
  lea    (%rdx,%rax,8),%rax
  movsd  -0x8(%rdx,%r14,8),%xmm0
  movsd  (%rax),%xmm1
  movsd  %xmm1,-0x8(%rdx,%r14,8)
  movslq %ecx,%rdx
  mov    0x40(%rsp),%rcx
  movsd  %xmm0,(%rax)
  movslq 0x2c(%rsp),%rax
  lea    (%r10,%rdx,8),%rdx
  lea    (%r10,%rax,8),%rax
  add    %rax,%rcx
  nopw   0x0(%rax,%rax,1)
L220:
  movsd  (%rax),%xmm0
  movsd  (%rdx),%xmm1
  add    $0x8,%rax
  add    $0x8,%rdx
  movsd  %xmm1,-0x8(%rax)
  movsd  %xmm0,-0x8(%rdx)
  cmp    %rax,%rcx
  jne    L220
  movsd  (%r15),%xmm7
  movapd %xmm6,%xmm1
  divsd  %xmm7,%xmm1
  cmp    0x10(%rsp),%r14
  jne    L159
L98:
  mov    0x38(%rsp),%r12
  mov    0x30(%rsp),%rbp
  movsd  (%r12),%xmm2
  lea    0x8(%r12),%rsi
  cmp    $0x1,%ebx
  je     L67
  mov    0x58(%rsp),%rcx
  mov    $0x1,%edx
  nopw   0x0(%rax,%rax,1)
L8:
  xor    %eax,%eax
  movapd %xmm3,%xmm1
  cs nopw 0x0(%rax,%rax,1)
L221:
  movsd  (%rcx,%rax,8),%xmm0
  mulsd  (%r12,%rax,8),%xmm0
  add    $0x1,%rax
  addsd  %xmm0,%xmm1
  cmp    %rax,%rdx
  jne    L221
  movsd  -0x8(%rsi,%rdx,8),%xmm2
  add    %rdi,%rcx
  subsd  %xmm1,%xmm2
  movsd  %xmm2,-0x8(%rsi,%rdx,8)
  add    $0x1,%rdx
  cmp    %edx,%ebx
  jg     L8
  lea    -0x2(%rbx),%eax
  xor    %edx,%edx
  shl    $0x3,%rax
  cmp    $0x1,%ebx
  cmovle %rdx,%rax
  mov    %ebx,%edx
  imul   %ebx,%edx
  add    %rax,%rsi
  lea    -0x1(%rbx),%eax
  movslq %edx,%rdx
  lea    -0x8(%rbp,%rdx,8),%rbp
L53:
  mov    $0xfffffffffffffff8,%r8
  lea    -0x1(%rbx),%r9d
  sub    %rdi,%r8
  lea    0x1(%rax),%edi
  cltq
  lea    (%r12,%rax,8),%rdx
  nopl   (%rax)
L15:
  cmp    %edi,%ebx
  jle    L166
  mov    %r9d,%eax
  movapd %xmm3,%xmm1
  sub    %edi,%eax
  lea    0x10(,%rax,8),%rcx
  mov    $0x8,%eax
  nopw   0x0(%rax,%rax,1)
L135:
  movsd  0x0(%rbp,%rax,1),%xmm0
  mulsd  (%rdx,%rax,1),%xmm0
  add    $0x8,%rax
  addsd  %xmm0,%xmm1
  cmp    %rax,%rcx
  jne    L135
  subsd  %xmm1,%xmm2
L166:
  movsd  %xmm2,(%rsi)
  lea    -0x1(%rdi),%eax
  sub    $0x8,%rsi
  divsd  0x0(%rbp),%xmm2
  sub    $0x8,%rdx
  add    %r8,%rbp
  sub    $0x2,%edi
  movsd  %xmm2,0x8(%rsi)
  js     L115
  movsd  (%rsi),%xmm2
  mov    %eax,%edi
  jmp    L15
L60:
  mov    %r8d,%ecx
  jmp    L233
L115:
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  add    $0x68,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L214:
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  add    $0x68,%rsp
  mov    $0xffffffff,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L67:
  mov    %r12,%rsi
  xor    %eax,%eax
  jmp    L53
```

tool_calls: name: disassemble_data
arguments: tag: D15
type: qword
name: disassemble_data
arguments: tag: D24
type: byte[8]
name: disassemble_data
arguments: tag: D26
type: f64

tool: D15: 9223372036854775807

tool: {"D24": [0x16,0x56,0xe7,0x9e,0xaf,0x03,0xd2,0x3c]}

tool: D26: 1

assistant: ```c
int solv(double *a, double *b, int n) {
  int i, j, k, lc;
  double *ps, *p, *q, *pa, *pd;
  double *q0, s, t, tq = 0., zr = 1.e-15;
  q0 = (double *)calloc(n, sizeof(double));
  for (j = 0, pa = a, pd = a; j < n; ++j, ++pa, pd += n + 1) {
    if (j) {
      for (i = 0, q = q0, p = pa; i < n; ++i, p += n)
        *q++ = *p;
      for (i = 1; i < n; ++i) {
        lc = i < j ? i : j;
        for (k = 0, p = pa + i * n - j, q = q0, t = 0.; k < lc; ++k)
          t += *p++ * *q++;
        q0[i] -= t;
      }
      for (i = 0, q = q0, p = pa; i < n; ++i, p += n)
        *p = *q++;
    }
    s = fabs(*pd);
    lc = j;
    for (k = j + 1, ps = pd; k < n; ++k) {
      if ((t = fabs(*(ps += n))) > s) {
        s = t;
        lc = k;
      }
    }
    tq = tq > s ? tq : s;
    if (s < zr * tq) {
      free(q0);
      return -1;
    }
    if (lc != j) {
      t = b[j];
      b[j] = b[lc];
      b[lc] = t;
      for (k = 0, p = a + n * j, q = a + n * lc; k < n; ++k) {
        t = *p;
        *p++ = *q;
        *q++ = t;
      }
    }
    for (k = j + 1, ps = pd, t = 1. / *pd; k < n; ++k)
      *(ps += n) *= t;
  }
  for (j = 1, ps = b + 1; j < n; ++j) {
    for (k = 0, p = a + n * j, q = b, t = 0.; k < j; ++k)
      t += *p++ * *q++;
    *ps++ -= t;
  }
  for (j = n - 1, --ps, pd = a + n * n - 1; j >= 0; --j, pd -= n + 1) {
    for (k = j + 1, p = pd, q = b + j, t = 0.; k < n; ++k)
      t += *++p * *++q;
    *ps -= t;
    *ps-- /= *pd;
  }
  free(q0);
  return 0;
}
```