user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<glassofast>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  lea    -0x19a000(%rsp),%r11
L20:
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  cmp    %r11,%rsp
  jne    L20
  sub    $0x338,%rsp
  mov    0x19a370(%rsp),%rax
  mov    %rsi,0x78(%rsp)
  mov    %r8d,%r10d
  mov    0x19a378(%rsp),%r8
  mov    %rdx,0x80(%rsp)
  mov    %ecx,0x88(%rsp)
  mov    %rax,0x60(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x19a328(%rsp)
  xor    %eax,%eax
  test   %edi,%edi
  jle    L352
  lea    -0x1(%rdi),%eax
  mov    %r9d,%ecx
  pxor   %xmm3,%xmm3
  mov    %edi,%r13d
  lea    0x0(,%rax,8),%r12
  mov    %eax,0x8c(%rsp)
  movq   D50(%rip),%xmm5
  movapd %xmm3,%xmm2
  lea    0x8(%r12),%r9
  lea    0x8(%rsi,%r12,1),%rdi
  mov    %r9,0x90(%rsp)
L137:
  mov    (%rsi),%rax
  lea    (%rax,%r9,1),%rdx
  cs nopw 0x0(%rax,%rax,1)
L213:
  movsd  (%rax),%xmm1
  add    $0x8,%rax
  andpd  %xmm5,%xmm1
  addsd  %xmm1,%xmm2
  cmp    %rax,%rdx
  jne    L213
  add    $0x8,%rsi
  cmp    %rdi,%rsi
  jne    L137
  mov    %r13d,%ebp
  mov    0x78(%rsp),%rsi
  lea    0x0(,%rbp,8),%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
L91:
  mov    (%rsi,%rax,1),%rdx
  movsd  (%rdx,%rax,1),%xmm1
  add    $0x8,%rax
  andpd  %xmm5,%xmm1
  subsd  %xmm1,%xmm2
  cmp    %rax,0x58(%rsp)
  jne    L91
  ucomisd %xmm3,%xmm2
  jp     L5
  je     L96
L5:
  pxor   %xmm1,%xmm1
  cvtsi2sdl 0x8c(%rsp),%xmm1
  divsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,0x98(%rsp)
  movapd %xmm0,%xmm8
  pxor   %xmm0,%xmm0
  cvtsi2sd %r13d,%xmm0
  divsd  %xmm0,%xmm8
  movsd  D30(%rip),%xmm0
  maxsd  %xmm8,%xmm0
  movapd %xmm0,%xmm8
  test   %ecx,%ecx
  jne    L390
  xor    %edx,%edx
L335:
  mov    0x78(%rsp),%rax
  mov    (%r8,%rdx,8),%rsi
  mov    (%rax,%rdx,8),%rdi
  mov    0x60(%rsp),%rax
  mov    (%rax,%rdx,8),%rcx
  xor    %eax,%eax
L364:
  movsd  (%rdi,%rax,8),%xmm0
  movsd  %xmm0,(%rsi,%rax,8)
  movq   $0x0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rbp
  jne    L364
  add    $0x1,%rdx
  cmp    %rdx,%rbp
  jne    L335
  movq   D56(%rip),%xmm9
L92:
  mov    0x78(%rsp),%rdi
  mov    0x80(%rsp),%r9
  xor    %eax,%eax
  lea    0xa0(%rsp),%rcx
L410:
  mov    (%r9,%rax,1),%rdx
  mov    (%rdi,%rax,1),%rsi
  movsd  (%rsi,%rax,1),%xmm0
  addsd  (%rdx,%rax,1),%xmm0
  mov    (%r8,%rax,1),%rdx
  movsd  %xmm0,(%rcx,%rax,1)
  movsd  %xmm0,(%rdx,%rax,1)
  add    $0x8,%rax
  cmp    %rax,0x58(%rsp)
  jne    L410
  mov    0x88(%rsp),%edx
  test   %edx,%edx
  jle    L318
  test   %r13d,%r13d
  mov    $0x8,%eax
  cmovg  0x58(%rsp),%rax
  movl   $0x1,0x74(%rsp)
  lea    0xcd1e0(%rsp),%rbx
  lea    0xcd1e8(%rsp,%r12,1),%r9
  mov    %rax,0x68(%rsp)
L253:
  mov    0x80(%rsp),%r14
  mov    0x78(%rsp),%r15
  xor    %r12d,%r12d
  movsd  %xmm3,0x8(%rsp)
  nopl   0x0(%rax)
L138:
  mov    0x58(%rsp),%rdx
  xor    %esi,%esi
  mov    %rbx,%rdi
  mov    %rcx,0x50(%rsp)
  mov    %r8,0x38(%rsp)
  mov    %r10d,0x30(%rsp)
  mov    %r9,0x28(%rsp)
  movaps %xmm9,0x10(%rsp)
  movsd  %xmm8,0x20(%rsp)
  call   <memset@plt>
  mov    0x60(%rsp),%rax
  mov    0x28(%rsp),%r9
  lea    0x0(,%r12,8),%r11
  mov    0x30(%rsp),%r10d
  mov    0x38(%rsp),%r8
  xor    %esi,%esi
  pxor   %xmm3,%xmm3
  mov    (%rax,%r12,8),%rdi
  mov    D11(%rip),%rax
  mov    0x50(%rsp),%rcx
  movsd  0x20(%rsp),%xmm8
  movq   %rax,%xmm4
  mov    D70(%rip),%rax
  movq   D50(%rip),%xmm5
  movapd 0x10(%rsp),%xmm9
  movq   %rax,%xmm7
  jmp    L454
  nopl   0x0(%rax,%rax,1)
L61:
  add    $0x1,%rsi
  cmp    %rsi,%rbp
  je     L290
L454:
  movsd  (%rdi,%rsi,8),%xmm1
  ucomisd %xmm3,%xmm1
  jp     L8
  je     L61
L8:
  mov    (%r8,%rsi,8),%rdx
  mov    %rbx,%rax
  nopl   (%rax)
L362:
  movsd  (%rdx),%xmm0
  add    $0x8,%rax
  add    $0x8,%rdx
  mulsd  %xmm1,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%r9
  jne    L362
  add    $0x1,%rsi
  cmp    %rsi,%rbp
  jne    L454
L290:
  xor    %eax,%eax
  movapd %xmm3,%xmm6
  cs nopw 0x0(%rax,%rax,1)
L531:
  cmp    %eax,%r12d
  je     L64
  mov    (%r14),%rdx
  movsd  (%rdx,%rax,8),%xmm1
  comisd %xmm1,%xmm4
  jbe    L64
  movsd  (%rcx,%rax,8),%xmm12
  movsd  (%rdi,%rax,8),%xmm11
  movapd %xmm3,%xmm10
  mov    (%r15),%rdx
  movapd %xmm12,%xmm2
  mulsd  %xmm11,%xmm2
  movsd  (%rdx,%rax,8),%xmm0
  subsd  (%rbx,%rax,8),%xmm0
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm2
  andpd  %xmm5,%xmm2
  subsd  %xmm1,%xmm2
  comisd %xmm2,%xmm3
  jae    L7
  comisd %xmm3,%xmm0
  jb     L157
  divsd  %xmm12,%xmm2
  movapd %xmm2,%xmm10
L7:
  movapd %xmm10,%xmm1
  subsd  %xmm11,%xmm1
  ucomisd %xmm3,%xmm1
  jp     L380
  je     L64
L380:
  movapd %xmm1,%xmm2
  andpd  %xmm5,%xmm2
  test   %r10d,%r10d
  je     L345
  comisd %xmm7,%xmm2
  ja     L345
  nopw   0x0(%rax,%rax,1)
L64:
  add    $0x1,%rax
  cmp    %rax,%rbp
  jne    L531
  comisd %xmm6,%xmm8
  jbe    L290
  movsd  (%rcx,%r12,8),%xmm0
  mov    (%r8,%r12,8),%rdi
  xor    %eax,%eax
  movapd %xmm3,%xmm1
  movsd  %xmm0,(%rbx,%r12,8)
  xchg   %ax,%ax
L578:
  movsd  (%rbx,%rax,8),%xmm0
  subsd  (%rdi,%rax,8),%xmm0
  add    $0x1,%rax
  andpd  %xmm5,%xmm0
  addsd  %xmm0,%xmm1
  cmp    %eax,%r13d
  jg     L578
  maxsd  0x8(%rsp),%xmm1
  mov    0x68(%rsp),%rdx
  mov    %rbx,%rsi
  mov    %r11,0x38(%rsp)
  mov    %r9,0x30(%rsp)
  mov    %r10d,0x28(%rsp)
  mov    %r8,0x10(%rsp)
  mov    %rcx,0x20(%rsp)
  movaps %xmm9,0x40(%rsp)
  movsd  %xmm8,0x50(%rsp)
  movsd  %xmm1,0x8(%rsp)
  call   <memcpy@plt>
  mov    0x38(%rsp),%r11
  mov    0x30(%rsp),%r9
  xor    %eax,%eax
  mov    0x28(%rsp),%r10d
  mov    0x10(%rsp),%r8
  pxor   %xmm3,%xmm3
  mov    0x20(%rsp),%rcx
  movapd 0x40(%rsp),%xmm9
  movsd  0x50(%rsp),%xmm8
  nopw   0x0(%rax,%rax,1)
L587:
  mov    (%r8,%rax,8),%rdx
  movsd  (%rbx,%rax,8),%xmm0
  add    $0x1,%rax
  movsd  %xmm0,(%rdx,%r11,1)
  cmp    %eax,%r13d
  jg     L587
  add    $0x1,%r12
  add    $0x8,%r14
  add    $0x8,%r15
  cmp    %r12,%rbp
  jne    L138
  movsd  0x98(%rsp),%xmm6
  comisd 0x8(%rsp),%xmm6
  jae    L263
  addl   $0x1,0x74(%rsp)
  mov    0x74(%rsp),%eax
  cmp    %eax,0x88(%rsp)
  jge    L253
L263:
  movsd  D1(%rip),%xmm4
  xor    %edi,%edi
L577:
  mov    0x60(%rsp),%rax
  mov    (%r8,%rdi,1),%rsi
  movapd %xmm3,%xmm1
  mov    (%rax,%rdi,1),%rdx
  xor    %eax,%eax
  nopl   0x0(%rax)
L127:
  movsd  (%rdx,%rax,8),%xmm0
  mulsd  (%rsi,%rax,8),%xmm0
  add    $0x1,%rax
  addsd  %xmm0,%xmm1
  cmp    %rax,%rbp
  jne    L127
  movsd  (%rcx,%rdi,1),%xmm0
  movapd %xmm4,%xmm2
  mov    0x90(%rsp),%rsi
  mov    %rdx,%rax
  subsd  %xmm1,%xmm0
  add    %rdx,%rsi
  divsd  %xmm0,%xmm2
  movapd %xmm2,%xmm1
  xorpd  %xmm9,%xmm1
  nopw   0x0(%rax,%rax,1)
L501:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rsi
  jne    L501
  movsd  %xmm2,(%rdx,%rdi,1)
  add    $0x8,%rdi
  cmp    %rdi,0x58(%rsp)
  jne    L577
  mov    0x8c(%rsp),%eax
  test   %eax,%eax
  je     L555
  lea    -0x2(%r13),%r8d
  mov    0x60(%rsp),%r9
  xor    %esi,%esi
  mov    $0x1,%edi
  movsd  D107(%rip),%xmm1
  add    $0x2,%r8
L367:
  cmp    %edi,%r13d
  jle    L314
  mov    (%r9,%rsi,1),%rcx
  mov    %rdi,%rax
  nopl   (%rax)
L257:
  mov    (%r9,%rax,8),%rdx
  movsd  (%rcx,%rax,8),%xmm0
  add    %rsi,%rdx
  addsd  (%rdx),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  movsd  %xmm0,(%rdx)
  cmp    %eax,%r13d
  jg     L257
L314:
  add    $0x1,%rdi
  add    $0x8,%rsi
  cmp    %rdi,%r8
  jne    L367
L555:
  mov    0x19a328(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L310
  mov    0x74(%rsp),%eax
  add    $0x19a338,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L390:
  movq   D56(%rip),%xmm9
  mov    0x60(%rsp),%rdi
  xor    %ecx,%ecx
L596:
  mov    (%rdi,%rcx,1),%rax
  mov    0x90(%rsp),%rsi
  lea    (%rax,%rcx,1),%rdx
  add    %rax,%rsi
L537:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  xorpd  %xmm9,%xmm0
  divsd  (%rdx),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rsi
  jne    L537
  movq   $0x0,(%rdx)
  add    $0x8,%rcx
  cmp    %rcx,0x58(%rsp)
  jne    L596
  jmp    L92
  nopl   (%rax)
L157:
  xorpd  %xmm9,%xmm2
  movapd %xmm2,%xmm10
  divsd  %xmm12,%xmm10
  jmp    L7
  nopl   0x0(%rax)
L345:
  mov    (%r8,%rax,8),%rsi
  movsd  %xmm10,(%rdi,%rax,8)
  xor    %edx,%edx
  nopl   0x0(%rax)
L382:
  movsd  (%rsi,%rdx,8),%xmm0
  mulsd  %xmm1,%xmm0
  addsd  (%rbx,%rdx,8),%xmm0
  movsd  %xmm0,(%rbx,%rdx,8)
  add    $0x1,%rdx
  cmp    %edx,%r13d
  jg     L382
  maxsd  %xmm6,%xmm2
  movapd %xmm2,%xmm6
  jmp    L64
L96:
  mov    0x58(%rsp),%r14
  xor    %ebx,%ebx
  mov    %r8,%r15
L371:
  mov    0x60(%rsp),%rax
  mov    (%r15,%rbx,8),%r13
  mov    (%rax,%rbx,8),%rdi
  lea    0x0(%r13,%r14,1),%rax
  cmp    %rax,%rdi
  jae    L395
  lea    (%rdi,%r14,1),%rax
  cmp    %rax,%r13
  jb     L441
L395:
  mov    %r14,%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    %r14,%rdx
  xor    %esi,%esi
  mov    %r13,%rdi
  call   <memset@plt>
L80:
  add    $0x1,%rbx
  cmp    %rbx,%rbp
  jne    L371
  mov    0x80(%rsp),%rsi
  mov    %r15,%r8
  xor    %eax,%eax
L522:
  mov    (%r8,%rax,1),%rdx
  mov    (%rsi,%rax,1),%rcx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  addsd  (%rcx,%rax,1),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,(%rdx)
  cmp    %rax,0x58(%rsp)
  jne    L522
  mov    0x60(%rsp),%rax
  lea    0x8(%rax,%r12,1),%rbp
  mov    %rax,%rbx
  mov    %r8,%r12
L26:
  mov    (%rbx),%rdi
  mov    0x58(%rsp),%rdx
  xor    %esi,%esi
  add    $0x8,%rbx
  call   <memset@plt>
  cmp    %rbp,%rbx
  jne    L26
  movsd  D1(%rip),%xmm2
  mov    0x60(%rsp),%rcx
  mov    %r12,%r8
  xor    %eax,%eax
  jmp    L73
L341:
  movapd %xmm2,%xmm0
  divsd  %xmm1,%xmm0
L481:
  mov    (%rcx,%rax,1),%rdx
  movsd  %xmm0,(%rdx,%rax,1)
  add    $0x8,%rax
  cmp    %rax,0x58(%rsp)
  je     L352
L73:
  mov    (%r8,%rax,1),%rdx
  movsd  (%rdx,%rax,1),%xmm1
  comisd D68(%rip),%xmm1
  ja     L341
  mov    D47(%rip),%rsi
  movq   %rsi,%xmm0
  jmp    L481
L352:
  movl   $0x0,0x74(%rsp)
  jmp    L555
L441:
  xor    %eax,%eax
L331:
  movq   $0x0,(%rdi,%rax,8)
  movq   $0x0,0x0(%r13,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rbp
  jne    L331
  jmp    L80
L318:
  movl   $0x1,0x74(%rsp)
  jmp    L263
L310:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D50
type: qword
name: extract_symbolic_data
arguments: label: D30
type: byte[8]
name: extract_symbolic_data
arguments: label: D56
type: f32[2]
name: extract_symbolic_data
arguments: label: D11
type: f64
name: extract_symbolic_data
arguments: label: D70
type: byte[8]
name: extract_symbolic_data
arguments: label: D1
type: f64
name: extract_symbolic_data
arguments: label: D107
type: f64
name: extract_symbolic_data
arguments: label: D68
type: byte[16]
name: extract_symbolic_data
arguments: label: D47
type: byte[8]

tool: D50: 9223372036854775807

tool: {"D30": [0xb1,0xde,0x64,0xc8,0xda,0xd0,0xe3,0x3c]}

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

tool: D11: 1000000000

tool: {"D70": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xb0,0x3e]}

tool: D1: 1

tool: D107: 0.5

tool: {"D68": [0xb1,0xde,0x64,0xc8,0xda,0xd0,0xd3,0x3c,0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

tool: {"D47": [0xe9,0xa2,0xa3,0xc5,0x81,0xd6,0x09,0x43]}

assistant: ```c
int glassofast(const int n, double **S, double **L, const double thr, const int maxit, int approxflg, int warm, double **X, double **W) {
  int i, j, ii, iter, jj;
  double a, b, c, delta, dlx, dw, shr, sum, thrlasso, tmp, wd[5000 * 21], wxj[5000 * 21];
  for (shr = ii = 0; ii < n; ii++)
    for (jj = 0; jj < n; jj++)
      shr += fabs(S[ii][jj]);
  for (i = 0; i < n; i++)
    shr -= fabs(S[i][i]);
  if (shr == 0.0) {
    for (ii = 0; ii < n; ii++)
      for (jj = 0; jj < n; jj++)
        W[ii][jj] = X[ii][jj] = 0.0;
    for (i = 0; i < n; i++)
      W[i][i] = W[i][i] + L[i][i];
    for (ii = 0; ii < n; ii++)
      for (jj = 0; jj < n; jj++)
        X[ii][jj] = 0.0;
    for (i = 0; i < n; i++)
      X[i][i] = 1.0 / ((W[i][i]) > ((1.1e-15)) ? (W[i][i]) : ((1.1e-15)));
    return 0;
  }
  shr *= thr / (n - 1);
  thrlasso = shr / n;
  if (thrlasso < 2 * (1.1e-15))
    thrlasso = 2 * (1.1e-15);
  if (!warm) {
    for (ii = 0; ii < n; ii++)
      for (jj = 0; jj < n; jj++) {
        W[ii][jj] = S[ii][jj];
        X[ii][jj] = 0.0;
      }
  } else {
    for (i = 0; i < n; i++) {
      for (ii = 0; ii < n; ii++)
        X[i][ii] = -X[i][ii] / X[i][i];
      X[i][i] = 0.0;
    }
  }
  for (i = 0; i < n; i++) {
    wd[i] = S[i][i] + L[i][i];
    W[i][i] = wd[i];
  }
  for (iter = 1; iter <= maxit; iter++) {
    dw = 0.0;
#pragma omp parallel for default(shared) private(i, j, ii, wxj, a, b, c, dlx, delta, sum)
    for (j = 0; j < n; j++) {
      for (ii = 0; ii < n; ii++)
        wxj[ii] = 0.0;
      for (i = 0; i < n; i++)
        if (X[j][i] != 0.0)
          for (ii = 0; ii < n; ii++)
            wxj[ii] += W[i][ii] * X[j][i];
      for (;;) {
        dlx = 0.0;
        for (i = 0; i < n; i++) {
          if (i != j && L[j][i] < (1e9)) {
            a = S[j][i] - wxj[i] + wd[i] * X[j][i];
            b = fabs(a) - L[j][i];
            if (b <= 0.0)
              c = 0.0;
            else if (a >= 0.0)
              c = b / wd[i];
            else
              c = -b / wd[i];
            delta = c - X[j][i];
            if (delta != 0.0 && (!approxflg || fabs(delta) > 1e-6)) {
              X[j][i] = c;
              for (ii = 0; ii < n; ii++)
                wxj[ii] += W[i][ii] * delta;
              if (fabs(delta) > dlx)
                dlx = fabs(delta);
            }
          }
        }
        if (dlx < thrlasso)
          break;
      }
      wxj[j] = wd[j];
      for (sum = ii = 0; ii < n; ii++)
        sum += fabs(wxj[ii] - W[j][ii]);
#pragma omp critical
      if (sum > dw)
        dw = sum;
      for (ii = 0; ii < n; ii++)
        W[j][ii] = wxj[ii];
      for (ii = 0; ii < n; ii++)
        W[ii][j] = wxj[ii];
    }
    if (dw <= shr)
      break;
  }
  for (i = 0; i < n; i++) {
    for (sum = ii = 0; ii < n; ii++)
      sum += X[i][ii] * W[i][ii];
    tmp = 1.0 / (wd[i] - sum);
    for (ii = 0; ii < n; ii++)
      X[i][ii] = -tmp * X[i][ii];
    X[i][i] = tmp;
  }
  for (i = 0; i < n - 1; i++) {
    for (ii = i + 1; ii < n; ii++) {
      X[i][ii] = 0.5 * (X[i][ii] + X[ii][i]);
      X[ii][i] = X[i][ii];
    }
  }
  return iter;
}
```