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

```asm
<fe_SigmaB2_at>:
  endbr64
  push   %r15
  mov    %rdi,%r15
  lea    D80(%rip),%rdi
  push   %r14
  push   %r13
  mov    $0x1,%r13d
  push   %r12
  push   %rbp
  mov    $0x2,%ebp
  push   %rbx
  mov    %edx,%ebx
  sub    $0x78,%rsp
  mov    %rsi,0x40(%rsp)
  mov    %ecx,%esi
  sub    $0x1,%ecx
  shl    %cl,%ebp
  shl    %cl,%r13d
  movsd  %xmm0,0x8(%rsp)
  lea    -0x1(%rbp),%eax
  mov    %edx,0x68(%rsp)
  mov    %eax,0x6c(%rsp)
  xor    %eax,%eax
  call   <printf@plt>
  mov    %r13d,%esi
  lea    D84(%rip),%rdi
  xor    %eax,%eax
  call   <printf@plt>
  mov    %ebp,%esi
  lea    D59(%rip),%rdi
  xor    %eax,%eax
  call   <printf@plt>
  lea    D39(%rip),%rdi
  xor    %eax,%eax
  mov    %ebx,%esi
  call   <printf@plt>
  movsd  0x8(%rsp),%xmm0
  mov    $0x1,%eax
  lea    D58(%rip),%rdi
  call   <printf@plt>
  test   %ebx,%ebx
  jle    L5
  mov    0x6c(%rsp),%eax
  pxor   %xmm7,%xmm7
  xor    %r14d,%r14d
  mov    %r13d,0x64(%rsp)
  cvtsi2sd %r13d,%xmm7
  imul   %ebp,%eax
  mov    %eax,0x54(%rsp)
  lea    -0x1(%r13),%eax
  mov    %eax,0x60(%rsp)
  movslq 0x68(%rsp),%rax
  movsd  %xmm7,0x58(%rsp)
  mov    %rax,0x48(%rsp)
  mov    %r15,%rax
  mov    %r14,%r15
  mov    %rax,%r14
  jmp    L66
  nopw   0x0(%rax,%rax,1)
L75:
  lea    -0x1(%r13),%esi
  lea    (%rbx,%rbx,1),%r12d
  movapd %xmm4,%xmm0
  add    $0x1,%r15
  imul   %ebp,%esi
  lea    (%rsi,%r13,1),%eax
  add    %r13d,%r13d
  add    %esi,%ebx
  mov    %esi,0x10(%rsp)
  add    %eax,%eax
  add    %ebx,%ebx
  cltq
  movslq %ebx,%rbx
  lea    -0x10(%r14,%rax,8),%rax
  addsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  lea    -0x1(%r12),%eax
  imul   %ebp,%eax
  lea    (%r12,%rax,1),%edi
  mov    %eax,0x20(%rsp)
  movslq %edi,%rdi
  lea    -0x8(%r14,%rdi,8),%rdi
  addsd  (%rdi),%xmm2
  movsd  %xmm2,(%rdi)
  lea    (%rax,%r13,1),%edi
  movslq %edi,%rdi
  lea    -0x10(%r14,%rdi,8),%rdi
  addsd  (%rdi),%xmm3
  movsd  %xmm3,(%rdi)
  movsd  %xmm3,-0x8(%r14,%rbx,8)
  call   <hat@plt>
  mov    0x20(%rsp),%eax
  mov    0x54(%rsp),%ebx
  mulsd  0x8(%rsp),%xmm0
  add    %ebp,%eax
  lea    (%r12,%rbx,1),%edx
  cltq
  movslq %edx,%rdx
  lea    -0x8(%r14,%rax,8),%rax
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  movsd  %xmm0,-0x8(%r14,%rdx,8)
  movsd  0x18(%rsp),%xmm0
  call   <hat@plt>
  mov    0x10(%rsp),%esi
  mulsd  0x8(%rsp),%xmm0
  lea    0x0(%r13,%rbx,1),%ecx
  movslq %ecx,%rcx
  lea    0x0(%rbp,%rsi,2),%eax
  cltq
  lea    -0x8(%r14,%rax,8),%rax
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  movsd  %xmm0,-0x10(%r14,%rcx,8)
  cmp    %r15,0x48(%rsp)
  je     L80
L66:
  mov    0x40(%rsp),%rax
  movsd  0x58(%rsp),%xmm5
  mulsd  (%rax,%r15,8),%xmm5
  movapd %xmm5,%xmm0
  movsd  %xmm5,0x10(%rsp)
  call   <ceil@plt>
  mov    0x64(%rsp),%esi
  mov    %eax,%edi
  call   <mini@plt>
  mov    $0x1,%esi
  mov    %eax,%edi
  call   <maxi@plt>
  movsd  0x10(%rsp),%xmm0
  subsd  D3(%rip),%xmm0
  mov    %eax,%r13d
  call   <ceil@plt>
  mov    0x60(%rsp),%esi
  mov    %eax,%edi
  call   <mini@plt>
  mov    $0x1,%esi
  mov    %eax,%edi
  call   <maxi@plt>
  pxor   %xmm1,%xmm1
  movsd  0x10(%rsp),%xmm0
  movsd  D12(%rip),%xmm6
  cvtsi2sd %r13d,%xmm1
  mov    %eax,%ebx
  subsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm6
  movapd %xmm6,%xmm0
  movsd  %xmm6,0x18(%rsp)
  call   <hat@plt>
  call   <sq@plt>
  movsd  0x8(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm4
  mulsd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  subsd  %xmm0,%xmm4
  addsd  D3(%rip),%xmm4
  movsd  %xmm1,0x38(%rsp)
  movapd %xmm4,%xmm0
  movsd  %xmm4,0x28(%rsp)
  call   <hat@plt>
  call   <sq@plt>
  movsd  0x8(%rsp),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  0x18(%rsp),%xmm0
  movsd  %xmm2,0x30(%rsp)
  call   <hat@plt>
  movsd  0x28(%rsp),%xmm4
  movsd  %xmm0,0x20(%rsp)
  movapd %xmm4,%xmm0
  call   <hat@plt>
  movsd  0x20(%rsp),%xmm3
  cmp    $0x9,%r15
  movsd  0x28(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm2
  movsd  0x38(%rsp),%xmm1
  mulsd  %xmm0,%xmm3
  mulsd  0x8(%rsp),%xmm3
  ja     L75
  movsd  0x10(%rsp),%xmm0
  mov    %ebx,%edx
  mov    %r13d,%esi
  mov    $0x4,%eax
  lea    D60(%rip),%rdi
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm3,0x30(%rsp)
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm1,0x20(%rsp)
  call   <printf@plt>
  movsd  0x38(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm3
  movsd  0x28(%rsp),%xmm2
  movsd  0x20(%rsp),%xmm1
  jmp    L75
  nopl   0x0(%rax)
L80:
  mov    %r14,%r15
L5:
  pxor   %xmm0,%xmm0
  mov    %ebp,%eax
  mov    0x6c(%rsp),%r14d
  xor    %ebx,%ebx
  cvtsi2sdl 0x68(%rsp),%xmm0
  mulsd  0x8(%rsp),%xmm0
  lea    D63(%rip),%r13
  imul   %ebp,%eax
  movslq %eax,%r12
  movsd  %xmm0,-0x8(%r15,%r12,8)
  test   %eax,%eax
  jne    L14
  jmp    L77
  xchg   %ax,%ax
L8:
  add    $0x1,%rbx
  cmp    %r12,%rbx
  je     L77
L14:
  movsd  (%r15,%rbx,8),%xmm0
  mov    %r13,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  mov    %ebx,%eax
  cltd
  idiv   %ebp
  cmp    %edx,%r14d
  jne    L8
  mov    $0xa,%edi
  add    $0x1,%rbx
  call   <putchar@plt>
  cmp    %r12,%rbx
  jne    L14
L77:
  add    $0x78,%rsp
  mov    $0xa,%edi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10b0 <putchar@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D80
type: string
name: extract_symbolic_data
arguments: label: D84
type: string
name: extract_symbolic_data
arguments: label: D59
type: string
name: extract_symbolic_data
arguments: label: D39
type: string
name: extract_symbolic_data
arguments: label: D58
type: string
name: extract_symbolic_data
arguments: label: D3
type: f64
name: extract_symbolic_data
arguments: label: D12
type: byte[8]
name: extract_symbolic_data
arguments: label: D60
type: string
name: extract_symbolic_data
arguments: label: D63
type: string

tool: D80: L: %i


tool: D84: n: %i


tool: D59: d: %i


tool: D39: N: %i


tool: D58: dt: %f


tool: D3: 0.5

tool: {"D12": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

tool: D60: %5.3f(%i, %i), %5.6f %5.6f %5.6f 


tool: D63:  %10.3f

assistant: ```c
void fe_SigmaB2_at(double *S, double *y, int N, double dt, int L) {
  int n = 1 << (L - 1);
  int K = 1;
  int d = 2 * n - 1 + K;
  int t;
  int i, j;
  double yn;
  printf("L: %i\n", L);
  printf("n: %i\n", n);
  printf("d: %i\n", d);
  printf("N: %i\n", N);
  printf("dt: %f\n", dt);
  double a, b, c;
  for (t = 0; t < N; t++) {
    yn = y[t] * (double)n;
    i = maxi(mini(((int)ceil(yn)), n), 1);
    j = maxi(mini(((int)ceil(yn - 0.5)), n - 1), 1);
    a = sq(hat(yn - i + 1)) * dt;
    b = sq(hat(yn - j + 0.5)) * dt;
    c = hat(yn - i + 1) * hat(yn - j + 0.5) * dt;
    if (t < 10)
      printf("%5.3f(%i, %i), %5.6f %5.6f %5.6f \n", yn, i, j, a, b, c);
    S[(2 * i - 2) * d + 2 * i - 2] += a;
    S[(2 * j - 1) * d + 2 * j - 1] += b;
    S[(2 * i - 2) * d + 2 * j - 1] = S[(2 * j - 1) * d + 2 * i - 2] += c;
    S[(d - 1) * d + 2 * j - 1] = S[(2 * j - 1) * d + d - 1] += hat(yn - j + 0.5) * dt;
    S[(d - 1) * d + 2 * i - 2] = S[(2 * i - 2) * d + d - 1] += hat(yn - i + 1) * dt;
  }
  S[d * d - 1] = N * dt;
  for (i = 0; i < d * d; i++) {
    printf(" %10.3f", S[i]);
    if (i % d == d - 1)
      printf("\n");
  }
  printf("\n");
}
```