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

```asm
<bitset>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  mov    %edx,%ebp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x858,%rsp
  mov    (%rsi),%r10d
  mov    %rcx,0x20(%rsp)
  mov    %r8d,0x14(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x848(%rsp)
  xor    %eax,%eax
  cmp    $0x1,%edx
  jle    L203
  lea    -0x2(%rdx),%eax
  lea    -0x1(%rdx),%ecx
  cmp    $0x2,%eax
  jbe    L106
  mov    %ecx,%edx
  movd   %r10d,%xmm5
  lea    0x4(%rsi),%rax
  shr    $0x2,%edx
  pshufd $0x0,%xmm5,%xmm1
  sub    $0x1,%edx
  movdqa %xmm1,%xmm3
  shl    $0x4,%rdx
  lea    0x14(%rsi,%rdx,1),%rdx
  nopl   0x0(%rax,%rax,1)
L173:
  movdqu (%rax),%xmm0
  movdqa %xmm1,%xmm2
  add    $0x10,%rax
  pcmpgtd %xmm0,%xmm2
  pand   %xmm2,%xmm1
  pandn  %xmm0,%xmm2
  por    %xmm2,%xmm1
  movdqa %xmm3,%xmm2
  pcmpgtd %xmm0,%xmm2
  pand   %xmm2,%xmm0
  pandn  %xmm3,%xmm2
  movdqa %xmm2,%xmm3
  por    %xmm0,%xmm3
  cmp    %rdx,%rax
  jne    L173
  movdqa %xmm3,%xmm2
  mov    %ecx,%edx
  psrldq $0x8,%xmm2
  and    $0xfffffffc,%edx
  movdqa %xmm2,%xmm0
  lea    0x1(%rdx),%eax
  pcmpgtd %xmm3,%xmm0
  pand   %xmm0,%xmm3
  pandn  %xmm2,%xmm0
  por    %xmm3,%xmm0
  movdqa %xmm0,%xmm3
  psrldq $0x4,%xmm3
  movdqa %xmm3,%xmm2
  pcmpgtd %xmm0,%xmm2
  pand   %xmm2,%xmm0
  pandn  %xmm3,%xmm2
  por    %xmm0,%xmm2
  movd   %xmm2,%r10d
  movdqa %xmm1,%xmm2
  psrldq $0x8,%xmm2
  movdqa %xmm2,%xmm0
  pcmpgtd %xmm1,%xmm0
  pand   %xmm0,%xmm2
  pandn  %xmm1,%xmm0
  por    %xmm2,%xmm0
  movdqa %xmm0,%xmm2
  psrldq $0x4,%xmm2
  movdqa %xmm2,%xmm1
  pcmpgtd %xmm0,%xmm1
  pand   %xmm1,%xmm2
  pandn  %xmm0,%xmm1
  por    %xmm2,%xmm1
  movd   %xmm1,0x28(%rsp)
  cmp    %edx,%ecx
  je     L60
L108:
  movslq %eax,%rdx
  mov    0x28(%rsp),%edi
  lea    0x0(,%rdx,4),%rcx
  mov    (%rbx,%rdx,4),%edx
  cmp    %edx,%edi
  cmovl  %edx,%edi
  cmp    %edx,%r10d
  cmovg  %edx,%r10d
  lea    0x1(%rax),%edx
  mov    %edi,0x28(%rsp)
  cmp    %edx,%ebp
  jle    L60
  mov    0x4(%rbx,%rcx,1),%edx
  cmp    %edx,%edi
  cmovl  %edx,%edi
  cmp    %edx,%r10d
  cmovg  %edx,%r10d
  add    $0x2,%eax
  mov    %edi,0x28(%rsp)
  cmp    %eax,%ebp
  jle    L60
  mov    0x8(%rbx,%rcx,1),%eax
  cmp    %eax,%edi
  cmovl  %eax,%edi
  cmp    %eax,%r10d
  cmovg  %eax,%r10d
  mov    %edi,0x28(%rsp)
L60:
  mov    0x28(%rsp),%edi
  cmp    %r10d,%edi
  jge    L35
  mov    %r10d,%eax
  mov    %edi,%r10d
  mov    %eax,0x28(%rsp)
L35:
  xor    %eax,%eax
  test   %r10d,%r10d
  mov    0x28(%rsp),%edi
  cmovns %r10d,%eax
  cmp    %r9d,%edi
  mov    %eax,%r15d
  lea    -0x1(%r9),%eax
  cmovl  %edi,%eax
  mov    %eax,0x28(%rsp)
  cmp    %r15d,%eax
  jl     L216
  lea    0x40(%rsp),%rax
  mov    0x14(%rsp),%r14d
  mov    %rax,0x18(%rsp)
  lea    0x44(%rsp),%rax
  imul   %r15d,%r14d
  mov    %rax,0x38(%rsp)
  jmp    L116
  nopl   (%rax)
L57:
  add    $0x1,%r15d
  add    0x14(%rsp),%r14d
  cmp    %r15d,0x28(%rsp)
  jl     L216
L116:
  mov    0x18(%rsp),%r8
  xor    %eax,%eax
  mov    %r15d,%ecx
  mov    %ebp,%edx
  mov    %rbx,%rsi
  mov    %r12,%rdi
  call   <xxsect@plt>
  cmp    $0x1,%eax
  jle    L57
  mov    0x38(%rsp),%rdi
  sub    $0x2,%eax
  mov    0x18(%rsp),%r13
  lea    (%rdi,%rax,4),%rax
  mov    %rax,0x8(%rsp)
  mov    0x14(%rsp),%eax
  sub    $0x1,%eax
  mov    %eax,0x2c(%rsp)
  movslq %r14d,%rax
  mov    %rax,0x30(%rsp)
  mov    %r14d,%eax
  mov    %r13,%r14
  mov    %r15d,%r13d
  mov    %eax,%r15d
  jmp    L244
  nop
L58:
  add    $0x4,%r14
  cmp    0x8(%rsp),%r14
  je     L187
L244:
  mov    0x4(%r14),%eax
  add    (%r14),%eax
  mov    %r13d,%r8d
  mov    %ebp,%edx
  mov    %eax,%ecx
  mov    %rbx,%rsi
  mov    %r12,%rdi
  shr    $0x1f,%ecx
  add    %eax,%ecx
  xor    %eax,%eax
  sar    %ecx
  call   <inside@plt>
  test   %eax,%eax
  je     L58
  mov    (%r14),%edx
  xor    %eax,%eax
  mov    0x4(%r14),%edi
  test   %edx,%edx
  cmovs  %eax,%edx
  cmp    0x14(%rsp),%edi
  cmovge 0x2c(%rsp),%edi
  cmp    %edx,%edi
  jl     L58
  mov    %edi,%eax
  sub    %edx,%eax
  lea    0x1(%rax),%r10d
  cmp    $0x2,%eax
  jbe    L55
  mov    0x20(%rsp),%rsi
  movslq %edx,%rax
  add    0x30(%rsp),%rax
  movaps D5(%rip),%xmm4
  lea    (%rsi,%rax,4),%rcx
  mov    %r10d,%eax
  shr    $0x2,%eax
  lea    -0x1(%rax),%esi
  lea    0x10(%rcx),%rax
  shl    $0x4,%rsi
  add    %rax,%rsi
  jmp    L19
  nopl   0x0(%rax,%rax,1)
L213:
  add    $0x10,%rax
L19:
  movups %xmm4,(%rcx)
  mov    %rax,%rcx
  cmp    %rax,%rsi
  jne    L213
  mov    %r10d,%eax
  and    $0xfffffffc,%eax
  add    %eax,%edx
  cmp    %r10d,%eax
  je     L58
L55:
  mov    0x20(%rsp),%rsi
  lea    (%r15,%rdx,1),%eax
  cltq
  movl   $0x3f800000,(%rsi,%rax,4)
  lea    0x1(%rdx),%eax
  cmp    %eax,%edi
  jl     L58
  add    %r15d,%eax
  add    $0x2,%edx
  cltq
  movl   $0x3f800000,(%rsi,%rax,4)
  cmp    %edi,%edx
  jg     L58
  add    %r15d,%edx
  add    $0x4,%r14
  movslq %edx,%rdx
  movl   $0x3f800000,(%rsi,%rdx,4)
  cmp    0x8(%rsp),%r14
  jne    L244
L187:
  mov    %r15d,%r14d
  mov    %r13d,%r15d
  add    0x14(%rsp),%r14d
  add    $0x1,%r15d
  cmp    %r15d,0x28(%rsp)
  jge    L116
L216:
  mov    0x848(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L129
  add    $0x858,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L203:
  mov    %r10d,0x28(%rsp)
  jmp    L35
L106:
  mov    %r10d,0x28(%rsp)
  mov    $0x1,%eax
  jmp    L108
L129:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D5
data_type: f32[4]

tool: {"D5": [1.0,1.0,1.0,1.0]}

assistant: ```c
void bitset(xverts, yverts, nverts, mask, npix, nrast) float mask[];
int xverts[], yverts[], nverts;
int npix, nrast;
{
  int xsects[512];
  int ifrast, ilrast;
  int ifpix, ilpix;
  int nxsect = 0;
  int i, j, k, n, ix, itemp;
  ifrast = ilrast = yverts[0];
  for (i = 1; i < nverts; i++) {
    ilrast = ((yverts[i]) > (ilrast) ? (yverts[i]) : (ilrast));
    ifrast = ((yverts[i]) < (ifrast) ? (yverts[i]) : (ifrast));
  }
  if (ilrast < ifrast)
    (itemp = (ifrast), (ifrast) = (ilrast), (ilrast) = itemp);
  ifrast = ((0) > (ifrast) ? (0) : (ifrast));
  ilrast = (((nrast - 1)) < (ilrast) ? ((nrast - 1)) : (ilrast));
  for (j = ifrast; j <= ilrast; j++)
    if ((nxsect = xxsect(xverts, yverts, nverts, j, xsects)) > 1)
      for (i = 1; i < nxsect; i++) {
        n = (j * npix);
        ix = (xsects[i - 1] + xsects[i]) / 2;
        if (inside(xverts, yverts, nverts, ix, j)) {
          ifpix = ((0) > (xsects[i - 1]) ? (0) : (xsects[i - 1]));
          ilpix = (((npix - 1)) < (xsects[i]) ? ((npix - 1)) : (xsects[i]));
          for (k = ifpix; k <= ilpix; k++)
            mask[n + k] = 1.0;
        }
      }
}
```