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

```asm
<vradfreq>:
  endbr64
  movsd  D7(%rip),%xmm6
  movapd %xmm0,%xmm5
  mov    %edi,%r11d
  divsd  %xmm6,%xmm5
  test   %edi,%edi
  jle    L55
  mov    %edx,%eax
  push   %r15
  movslq %esi,%rdx
  movslq %r11d,%r10
  push   %r14
  mov    %esi,%edi
  sub    $0x1,%r10
  lea    0x0(,%rdx,8),%rsi
  push   %r13
  push   %r12
  movslq %eax,%r12
  push   %rbp
  push   %rbx
  mov    %rcx,%rbx
  lea    0x0(,%r12,8),%rcx
  cmp    $0x1,%eax
  jne    L121
  cmp    $0x1,%edi
  jne    L121
  mov    %rcx,%rax
  xor    %r14d,%r14d
  mov    %rsi,%rdi
  imul   %r10,%rax
  mov    %r14,%rbp
  mov    %r14,%r15
  mov    %r14,%r13
  imul   %r10,%rdi
  test   %rcx,%rcx
  cmovs  %rax,%rbp
  add    %r8,%rbp
  test   %rcx,%rcx
  cmovns %rax,%r15
  add    %r8,%r15
  test   %rsi,%rsi
  cmovs  %rdi,%r13
  add    %rbx,%r13
  cmp    %r13,%r15
  setb   %r13b
  test   %rsi,%rsi
  cmovs  %r14,%rdi
  add    %rbx,%rdi
  cmp    %rdi,%rbp
  seta   %dil
  or     %edi,%r13d
  lea    -0x1(%r11),%edi
  cmp    $0x2d,%edi
  seta   %dil
  test   %dil,%r13b
  je     L29
  add    $0x4,%rax
  mov    $0x4,%edi
  test   %rcx,%rcx
  cmovs  %rdi,%rax
  lea    (%r9,%r10,4),%rdi
  add    %r8,%rax
  cmp    %rax,%r9
  seta   %al
  cmp    %rdi,%rbp
  seta   %dil
  or     %dil,%al
  je     L29
  mov    %r11d,%r14d
  mov    %rdx,%rbp
  mov    %r12,%r13
  shl    $0x4,%rdx
  shr    $0x2,%r14d
  movapd %xmm5,%xmm4
  lea    (%r12,%r12,2),%r12
  shl    $0x5,%rbp
  sub    $0x1,%r14d
  unpcklpd %xmm4,%xmm4
  shl    $0x5,%r13
  mov    %rbx,%r10
  shl    $0x4,%r14
  mov    %r8,%rax
  mov    %r9,%rdi
  add    %rbx,%rdx
  movapd D7(%rip),%xmm3
  lea    0x10(%r9,%r14,1),%r14
  shl    $0x3,%r12
  pxor   %xmm7,%xmm7
  nopl   0x0(%rax)
L49:
  movsd  (%r10),%xmm0
  movsd  (%rdx),%xmm2
  movapd %xmm3,%xmm1
  add    $0x10,%rdi
  movhpd (%r10,%rsi,1),%xmm0
  movhpd (%rdx,%rsi,1),%xmm2
  add    %rbp,%r10
  add    %rbp,%rdx
  subpd  %xmm0,%xmm1
  movapd %xmm3,%xmm0
  subpd  %xmm2,%xmm0
  mulpd  %xmm4,%xmm1
  mulpd  %xmm4,%xmm0
  movlpd %xmm1,(%rax)
  movhpd %xmm1,(%rax,%rcx,1)
  movlpd %xmm0,(%rax,%rcx,2)
  movhpd %xmm0,(%rax,%r12,1)
  add    %r13,%rax
  movups %xmm7,-0x10(%rdi)
  cmp    %rdi,%r14
  jne    L49
L120:
  mov    %r11d,%edi
  mov    %rsi,%rdx
  and    $0xfffffffc,%edi
  mov    %edi,%eax
  imul   %rax,%rdx
  lea    (%r9,%rax,4),%r9
  imul   %rcx,%rax
  add    %rbx,%rdx
  add    %r8,%rax
  cmp    %edi,%r11d
  je     L9
  movapd %xmm6,%xmm0
  subsd  (%rdx),%xmm0
  lea    0x1(%rdi),%r8d
  add    %rsi,%rdx
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,(%rax)
  add    %rcx,%rax
  movl   $0x0,(%r9)
  cmp    %r8d,%r11d
  jle    L9
  movapd %xmm6,%xmm0
  subsd  (%rdx),%xmm0
  add    $0x2,%edi
  add    %rsi,%rdx
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,(%rax)
  add    %rcx,%rax
  movl   $0x0,0x4(%r9)
  cmp    %edi,%r11d
  jle    L9
  subsd  (%rdx),%xmm6
  mulsd  %xmm5,%xmm6
  movsd  %xmm6,(%rax)
  movl   $0x0,0x8(%r9)
L9:
  pop    %rbx
  xor    %eax,%eax
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L121:
  mov    %rcx,%rax
  xor    %r14d,%r14d
  mov    %rsi,%rdi
  imul   %r10,%rax
  mov    %r14,%rbp
  mov    %r14,%r15
  mov    %r14,%r13
  imul   %r10,%rdi
  test   %rcx,%rcx
  cmovs  %rax,%rbp
  add    %r8,%rbp
  test   %rcx,%rcx
  cmovns %rax,%r15
  add    %r8,%r15
  test   %rsi,%rsi
  cmovs  %rdi,%r13
  add    %rbx,%r13
  cmp    %r13,%r15
  setb   %r13b
  test   %rsi,%rsi
  cmovs  %r14,%rdi
  add    %rbx,%rdi
  cmp    %rdi,%rbp
  seta   %dil
  or     %edi,%r13d
  lea    -0x1(%r11),%edi
  cmp    $0x2d,%edi
  seta   %dil
  test   %dil,%r13b
  je     L129
  add    $0x4,%rax
  mov    $0x4,%edi
  test   %rcx,%rcx
  cmovs  %rdi,%rax
  lea    (%r9,%r10,4),%rdi
  add    %r8,%rax
  cmp    %rax,%r9
  seta   %al
  cmp    %rdi,%rbp
  seta   %dil
  or     %dil,%al
  je     L129
  mov    %r11d,%r14d
  mov    %rdx,%rbp
  mov    %r12,%r13
  shl    $0x4,%rdx
  shr    $0x2,%r14d
  movapd %xmm5,%xmm4
  lea    (%r12,%r12,2),%r12
  shl    $0x5,%rbp
  sub    $0x1,%r14d
  unpcklpd %xmm4,%xmm4
  shl    $0x5,%r13
  mov    %rbx,%r10
  shl    $0x4,%r14
  mov    %r8,%rax
  mov    %r9,%rdi
  add    %rbx,%rdx
  movapd D7(%rip),%xmm3
  lea    0x10(%r9,%r14,1),%r14
  shl    $0x3,%r12
  pxor   %xmm7,%xmm7
  nopw   0x0(%rax,%rax,1)
L126:
  movsd  (%r10),%xmm0
  movsd  (%rdx),%xmm2
  movapd %xmm3,%xmm1
  add    $0x10,%rdi
  movhpd (%r10,%rsi,1),%xmm0
  movhpd (%rdx,%rsi,1),%xmm2
  add    %rbp,%r10
  add    %rbp,%rdx
  subpd  %xmm0,%xmm1
  movapd %xmm3,%xmm0
  subpd  %xmm2,%xmm0
  mulpd  %xmm4,%xmm1
  mulpd  %xmm4,%xmm0
  movlpd %xmm1,(%rax)
  movhpd %xmm1,(%rax,%rcx,1)
  movlpd %xmm0,(%rax,%rcx,2)
  movhpd %xmm0,(%rax,%r12,1)
  add    %r13,%rax
  movups %xmm7,-0x10(%rdi)
  cmp    %r14,%rdi
  jne    L126
  jmp    L120
  nopl   0x0(%rax,%rax,1)
L29:
  mov    %r11d,%r11d
  lea    (%r9,%r11,4),%rax
  nop
L96:
  movapd %xmm6,%xmm0
  subsd  (%rbx),%xmm0
  add    $0x4,%r9
  add    %rsi,%rbx
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,(%r8)
  add    %rcx,%r8
  movl   $0x0,-0x4(%r9)
  cmp    %rax,%r9
  jne    L96
  pop    %rbx
  xor    %eax,%eax
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L55:
  xor    %eax,%eax
  ret
L129:
  mov    %r11d,%r11d
  lea    (%r9,%r11,4),%rax
  nopw   0x0(%rax,%rax,1)
L151:
  movapd %xmm6,%xmm0
  subsd  (%rbx),%xmm0
  add    $0x4,%r9
  add    %rsi,%rbx
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,(%r8)
  add    %rcx,%r8
  movl   $0x0,-0x4(%r9)
  cmp    %rax,%r9
  jne    L151
  jmp    L9
```

tool_calls: name: read_data
arguments: data_label: D7
datatype: f64[2]

tool: {"D7": [299792458.0,299792458.0]}

assistant: ```c
int vradfreq(restfrq, nvrad, svrad, sfreq, vrad, freq, stat)
double restfrq;
int nvrad, svrad, sfreq;
const double vrad[];
double freq[];
int stat[];
{
  double r;
  register int ivrad, *statp;
  register const double *vradp;
  register double *freqp;
  r = restfrq / 2.99792458e8;
  vradp = vrad;
  freqp = freq;
  statp = stat;
  for (ivrad = 0; ivrad < nvrad; ivrad++) {
    *freqp = r * (2.99792458e8 - *vradp);
    *(statp++) = 0;
    vradp += svrad;
    freqp += sfreq;
  }
  return 0;
}
```