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

```asm
<sfa_special>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  add    $0xffffffffffffff80,%rsp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %rcx,-0x50(%rbp)
  mov    %r8,-0x58(%rbp)
  mov    %r9,-0x60(%rbp)
  movsd  %xmm0,-0x68(%rbp)
  movsd  %xmm1,-0x80(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0x70(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0x78(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    0x10(%rbp),%eax
  cltq
  shl    $0x2,%eax
  mov    %eax,%edi
  call   <malloc@plt>
  mov    %rax,-0x18(%rbp)
  movsd  D31(%rip),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movl   $0xffffffff,-0x28(%rbp)
  movl   $0xffffffff,-0x24(%rbp)
  mov    -0x40(%rbp),%rax
  movsd  D31(%rip),%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x1,-0x2c(%rbp)
  jmp    L30
L35:
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm0
  addsd  %xmm0,%xmm0
  mov    -0x50(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x60(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x78(%rbp),%rax
  movsd  %xmm0,(%rax)
  movl   $0x1,-0x30(%rbp)
  jmp    L1
L107:
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x30(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x30(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x30(%rbp)
L1:
  mov    0x10(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0x30(%rbp)
  jl     L107
  mov    0x10(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    0x10(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x10(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    0x10(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  mov    0x10(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    0x10(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x28(%rbp),%eax
  mov    %eax,-0x24(%rbp)
  mov    -0x68(%rbp),%rdi
  mov    -0x18(%rbp),%r8
  mov    -0x78(%rbp),%rcx
  mov    -0x50(%rbp),%rdx
  mov    -0x58(%rbp),%rsi
  mov    -0x38(%rbp),%rax
  mov    0x10(%rbp),%r9d
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  call   <supportSet@plt>
  mov    %eax,-0x28(%rbp)
  movl   $0x0,-0x30(%rbp)
  jmp    L68
L112:
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x30(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x30(%rbp)
L68:
  mov    -0x30(%rbp),%eax
  cmp    0x10(%rbp),%eax
  jl     L112
  mov    -0x70(%rbp),%rdx
  mov    -0x50(%rbp),%rsi
  lea    -0x20(%rbp),%rax
  mov    0x10(%rbp),%ecx
  mov    %rax,%rdi
  call   <Thomas@plt>
  movl   $0x0,-0x30(%rbp)
  jmp    L108
L131:
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x68(%rbp),%xmm0
  jbe    L44
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x68(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L87
L44:
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x68(%rbp),%xmm0
  movq   D25(%rip),%xmm2
  xorpd  %xmm2,%xmm0
  comisd %xmm1,%xmm0
  jbe    L87
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x68(%rbp),%xmm0
  movq   D25(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L87:
  addl   $0x1,-0x30(%rbp)
L108:
  mov    -0x30(%rbp),%eax
  cmp    0x10(%rbp),%eax
  jl     L131
  mov    -0x2c(%rbp),%eax
  cltd
  idivl  0x30(%rbp)
  mov    %edx,%eax
  test   %eax,%eax
  jne    L19
  mov    -0x40(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x68(%rbp),%rdi
  mov    -0x60(%rbp),%r8
  mov    -0x70(%rbp),%rcx
  mov    -0x78(%rbp),%rdx
  mov    -0x50(%rbp),%rsi
  mov    -0x40(%rbp),%rax
  mov    0x10(%rbp),%r9d
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  call   <dualityGap@plt>
  mov    -0x40(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x80(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jae    L81
  mov    -0x40(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D19(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L19
  mov    -0x28(%rbp),%eax
  cmp    -0x24(%rbp),%eax
  jne    L19
  mov    -0x40(%rbp),%rax
  movsd  (%rax),%xmm0
  ucomisd -0x10(%rbp),%xmm0
  setnp  %al
  mov    $0x0,%edx
  ucomisd -0x10(%rbp),%xmm0
  cmovne %edx,%eax
  movzbl %al,%eax
  mov    %eax,%edi
  call   <fabs@plt>
  test   %rax,%rax
  jne    L47
L19:
  addl   $0x1,-0x2c(%rbp)
L30:
  mov    -0x2c(%rbp),%eax
  cmp    0x18(%rbp),%eax
  jle    L35
  jmp    L56
L81:
  nop
  jmp    L56
L47:
  nop
L56:
  mov    -0x18(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    -0x48(%rbp),%rax
  mov    -0x28(%rbp),%edx
  mov    %edx,(%rax)
  mov    -0x2c(%rbp),%eax
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L82
  call   <__stack_chk_fail@plt>
L82:
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D31
datatype: byte[8]
name: extract_symbolic_data
arguments: data_label: D25
datatype: byte[8]
name: extract_symbolic_data
arguments: data_label: D19
datatype: f64

tool: {"D31": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: {"D25": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

tool: D19: 1

assistant: ```c
int sfa_special(double *x, double *gap, int *activeS, double *z, double *v, double *Av, double lambda, int nn, int maxStep, double *s, double *g, double tol, int tau) {
  int i, iterStep;
  double temp;
  int *S = (int *)malloc(sizeof(int) * nn);
  double gapp = -1;
  int numS = -1, numSp = -1;
  *gap = -1;
  for (iterStep = 1; iterStep <= maxStep; iterStep++) {
    g[0] = z[0] + z[0] - z[1] - Av[0];
    for (i = 1; i < nn - 1; i++) {
      g[i] = -z[i - 1] + z[i] + z[i] - z[i + 1] - Av[i];
    }
    g[nn - 1] = -z[nn - 2] + z[nn - 1] + z[nn - 1] - Av[nn - 1];
    numSp = numS;
    numS = supportSet(x, v, z, g, S, lambda, nn);
    for (i = 0; i < nn; i++)
      s[i] = Av[i] - x[i + 1] + x[i];
    Thomas(&temp, z, s, nn);
    for (i = 0; i < nn; i++) {
      if (z[i] > lambda)
        z[i] = lambda;
      else if (z[i] < -lambda)
        z[i] = -lambda;
    }
    if (iterStep % tau == 0) {
      gapp = *gap;
      dualityGap(gap, z, g, s, Av, lambda, nn);
      if (*gap <= tol) {
        break;
      }
      if ((*gap < 1) && (numS == numSp) && fabs(gapp == *gap)) {
        break;
      }
    }
  }
  free(S);
  *activeS = numS;
  return (iterStep);
}
```