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

```asm
<anigauss>:
  endbr64
  push   %r15
  mulsd  %xmm1,%xmm1
  movapd %xmm0,%xmm5
  mov    %rsi,%r15
  push   %r14
  mulsd  %xmm0,%xmm5
  push   %r13
  mov    %r8d,%r13d
  push   %r12
  mov    %r9d,%r12d
  push   %rbp
  mov    %ecx,%ebp
  push   %rbx
  mov    %edx,%ebx
  sub    $0x78,%rsp
  mov    %rdi,0x8(%rsp)
  lea    0x20(%rsp),%rsi
  lea    0x28(%rsp),%rdi
  mulsd  D26(%rip),%xmm2
  divsd  D22(%rip),%xmm2
  movapd %xmm2,%xmm0
  mov    %fs:0x28,%rax
  mov    %rax,0x68(%rsp)
  xor    %eax,%eax
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm2,(%rsp)
  call   <sincos@plt>
  movsd  0x20(%rsp),%xmm2
  pxor   %xmm7,%xmm7
  movsd  0x28(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm5
  movapd %xmm2,%xmm3
  movapd %xmm4,%xmm0
  mulsd  %xmm2,%xmm3
  mulsd  %xmm4,%xmm2
  mulsd  %xmm4,%xmm0
  movapd %xmm1,%xmm4
  subsd  %xmm5,%xmm4
  mulsd  %xmm4,%xmm2
  movapd %xmm3,%xmm4
  mulsd  %xmm5,%xmm4
  movapd %xmm0,%xmm6
  mulsd  %xmm1,%xmm6
  mulsd  %xmm1,%xmm3
  mulsd  %xmm5,%xmm0
  movapd %xmm2,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm6,%xmm4
  addsd  %xmm0,%xmm3
  divsd  %xmm4,%xmm1
  subsd  %xmm1,%xmm3
  ucomisd %xmm3,%xmm7
  ja     L25
  sqrtsd %xmm3,%xmm3
  movapd %xmm3,%xmm1
L10:
  divsd  %xmm4,%xmm2
  ucomisd %xmm4,%xmm7
  ja     L57
  sqrtsd %xmm4,%xmm4
L22:
  lea    0x30(%rsp),%r14
  movapd %xmm1,%xmm0
  movsd  %xmm2,0x10(%rsp)
  mov    %r14,%rdi
  movsd  %xmm4,0x18(%rsp)
  call   <YvVfilterCoef@plt>
  mov    0x8(%rsp),%rdi
  mov    %r14,%r8
  mov    %ebp,%ecx
  mov    %ebx,%edx
  mov    %r15,%rsi
  call   <f_iir_xline_filter@plt>
  movsd  0x18(%rsp),%xmm4
  mov    %r14,%rdi
  movapd %xmm4,%xmm0
  call   <YvVfilterCoef@plt>
  movsd  0x10(%rsp),%xmm2
  pxor   %xmm6,%xmm6
  ucomisd %xmm6,%xmm2
  jp     L107
  jne    L107
  mov    %r14,%r8
  mov    %ebp,%ecx
  mov    %ebx,%edx
  mov    %r15,%rsi
  mov    %r15,%rdi
  call   <f_iir_yline_filter@plt>
L85:
  test   %r13d,%r13d
  jle    L1
  movsd  (%rsp),%xmm7
  subsd  D9(%rip),%xmm7
  xor    %r14d,%r14d
  movsd  %xmm7,0x8(%rsp)
  xchg   %ax,%ax
L54:
  movsd  0x8(%rsp),%xmm0
  mov    $0x1,%r8d
  mov    %ebp,%ecx
  mov    %ebx,%edx
  mov    %r15,%rsi
  mov    %r15,%rdi
  add    $0x1,%r14d
  call   <f_iir_derivative_filter@plt>
  cmp    %r14d,%r13d
  jne    L54
L1:
  xor    %r13d,%r13d
  test   %r12d,%r12d
  jle    L42
  nopl   0x0(%rax)
L21:
  movsd  (%rsp),%xmm0
  mov    $0x1,%r8d
  mov    %ebp,%ecx
  mov    %ebx,%edx
  mov    %r15,%rsi
  mov    %r15,%rdi
  add    $0x1,%r13d
  call   <f_iir_derivative_filter@plt>
  cmp    %r13d,%r12d
  jne    L21
L42:
  mov    0x68(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L20
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  cs nopw 0x0(%rax,%rax,1)
L107:
  movapd %xmm2,%xmm0
  mov    %r14,%r8
  mov    %ebp,%ecx
  mov    %ebx,%edx
  mov    %r15,%rsi
  mov    %r15,%rdi
  call   <f_iir_tline_filter@plt>
  jmp    L85
L25:
  movapd %xmm3,%xmm0
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm2,0x10(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm2
  pxor   %xmm7,%xmm7
  movapd %xmm0,%xmm1
  jmp    L10
L57:
  movapd %xmm4,%xmm0
  movsd  %xmm2,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  movapd %xmm0,%xmm4
  jmp    L22
L20:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_tag: D26
type_infered: f64
name: unpack_data
arguments: data_tag: D22
type_infered: f64
name: unpack_data
arguments: data_tag: D9
type_infered: byte[8]

tool: D26: 3.141592653589793

tool: D22: 180

tool: {"D9": [0x18,0x2d,0x44,0x54,0xfb,0x21,0xf9,0x3f]}

assistant: ```c
void anigauss(double *input, double *output, int sizex, int sizey, double sigmav, double sigmau, double phi, int orderv, int orderu) {
  double filter[7];
  double sigmax, sigmay, tanp;
  double su2, sv2;
  double phirad;
  double a11, a21, a22;
  int i;
  su2 = sigmau * sigmau;
  sv2 = sigmav * sigmav;
  phirad = phi * 3.14159265358979323846 / 180.;
  a11 = cos(phirad) * cos(phirad) * su2 + sin(phirad) * sin(phirad) * sv2;
  a21 = cos(phirad) * sin(phirad) * (su2 - sv2);
  a22 = cos(phirad) * cos(phirad) * sv2 + sin(phirad) * sin(phirad) * su2;
  sigmax = sqrt(a11 - a21 * a21 / a22);
  tanp = a21 / a22;
  sigmay = sqrt(a22);
  YvVfilterCoef(sigmax, filter);
  f_iir_xline_filter(input, output, sizex, sizey, filter);
  YvVfilterCoef(sigmay, filter);
  if (tanp != 0.0) {
    f_iir_tline_filter(output, output, sizex, sizey, filter, tanp);
  } else {
    f_iir_yline_filter(output, output, sizex, sizey, filter);
  }
  for (i = 0; i < orderv; i++)
    f_iir_derivative_filter(output, output, sizex, sizey, phirad - 3.14159265358979323846 / 2., 1);
  for (i = 0; i < orderu; i++)
    f_iir_derivative_filter(output, output, sizex, sizey, phirad, 1);
}
```