user: Please show me the C implementation that matches the assembly code listed below:

```asm
<anigauss>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xe0,%rsp
  mov    %rdi,-0x98(%rbp)
  mov    %rsi,-0xa0(%rbp)
  mov    %edx,-0xa4(%rbp)
  mov    %ecx,-0xa8(%rbp)
  movsd  %xmm0,-0xb0(%rbp)
  movsd  %xmm1,-0xb8(%rbp)
  movsd  %xmm2,-0xc0(%rbp)
  mov    %r8d,-0xc4(%rbp)
  mov    %r9d,-0xc8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  -0xb8(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0xb0(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0xc0(%rbp),%xmm1
  movsd  D9(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D0(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movsd  %xmm0,-0xd0(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  mulsd  -0xd0(%rbp),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  %xmm0,-0xd8(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  mulsd  -0xd8(%rbp),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  addsd  -0xd0(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movsd  %xmm0,-0xd0(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0xd0(%rbp),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x88(%rbp),%xmm0
  subsd  -0x80(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movsd  %xmm0,-0xd0(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  mulsd  -0xd0(%rbp),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  %xmm0,-0xd8(%rbp)
  mov    -0x78(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  mulsd  -0xd8(%rbp),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  addsd  -0xd0(%rbp),%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x68(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x60(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x58(%rbp)
  movsd  -0x68(%rbp),%xmm0
  divsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x60(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x48(%rbp)
  lea    -0x40(%rbp),%rdx
  mov    -0x58(%rbp),%rax
  mov    %rdx,%rdi
  movq   %rax,%xmm0
  call   <YvVfilterCoef@plt>
  lea    -0x40(%rbp),%rdi
  mov    -0xa8(%rbp),%ecx
  mov    -0xa4(%rbp),%edx
  mov    -0xa0(%rbp),%rsi
  mov    -0x98(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <f_iir_xline_filter@plt>
  lea    -0x40(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  mov    %rdx,%rdi
  movq   %rax,%xmm0
  call   <YvVfilterCoef@plt>
  pxor   %xmm0,%xmm0
  ucomisd -0x50(%rbp),%xmm0
  jp     L52
  pxor   %xmm0,%xmm0
  ucomisd -0x50(%rbp),%xmm0
  je     L12
L52:
  mov    -0x50(%rbp),%rdi
  lea    -0x40(%rbp),%r8
  mov    -0xa8(%rbp),%ecx
  mov    -0xa4(%rbp),%edx
  mov    -0xa0(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  call   <f_iir_tline_filter@plt>
  jmp    L15
L12:
  lea    -0x40(%rbp),%rdi
  mov    -0xa8(%rbp),%ecx
  mov    -0xa4(%rbp),%edx
  mov    -0xa0(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <f_iir_yline_filter@plt>
L15:
  movl   $0x0,-0x8c(%rbp)
  jmp    L50
L25:
  movsd  -0x78(%rbp),%xmm0
  movsd  D8(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rdi
  mov    -0xa8(%rbp),%ecx
  mov    -0xa4(%rbp),%edx
  mov    -0xa0(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  mov    $0x1,%r8d
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  call   <f_iir_derivative_filter@plt>
  addl   $0x1,-0x8c(%rbp)
L50:
  mov    -0x8c(%rbp),%eax
  cmp    -0xc4(%rbp),%eax
  jl     L25
  movl   $0x0,-0x8c(%rbp)
  jmp    L33
L73:
  mov    -0x78(%rbp),%rdi
  mov    -0xa8(%rbp),%ecx
  mov    -0xa4(%rbp),%edx
  mov    -0xa0(%rbp),%rsi
  mov    -0xa0(%rbp),%rax
  mov    $0x1,%r8d
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  call   <f_iir_derivative_filter@plt>
  addl   $0x1,-0x8c(%rbp)
L33:
  mov    -0x8c(%rbp),%eax
  cmp    -0xc8(%rbp),%eax
  jl     L73
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L59
  call   <__stack_chk_fail@plt>
L59:
  leave
  ret
```

tool_calls: name: parse_data
arguments: tag: D9
data_type: f64
name: parse_data
arguments: tag: D0
data_type: f64
name: parse_data
arguments: tag: D8
data_type: byte[8]

tool: D9: 3.141592653589793

tool: D0: 180

tool: {"D8": [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);
}
```