user: Demonstrate the C code that generates the assembly output below:

```asm
<anigauss>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,0x28(%rsp)
  mov    %rsi,%rbx
  mov    %edx,%ebp
  mov    %ecx,%r12d
  mov    %r8d,%r14d
  mov    %r9d,%r13d
  mov    %fs:0x28,%rax
  mov    %rax,0x78(%rsp)
  xor    %eax,%eax
  mulsd  %xmm1,%xmm1
  movsd  %xmm1,0x10(%rsp)
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,0x18(%rsp)
  mulsd  D23(%rip),%xmm2
  divsd  D27(%rip),%xmm2
  movsd  %xmm2,0x8(%rsp)
  lea    0x38(%rsp),%rdi
  lea    0x30(%rsp),%rsi
  movapd %xmm2,%xmm0
  call   <sincos@plt>
  movsd  0x30(%rsp),%xmm1
  movsd  0x38(%rsp),%xmm3
  movapd %xmm1,%xmm6
  mulsd  %xmm1,%xmm6
  movapd %xmm6,%xmm0
  movapd %xmm3,%xmm2
  mulsd  %xmm3,%xmm2
  mulsd  %xmm3,%xmm1
  movsd  0x10(%rsp),%xmm4
  movapd %xmm4,%xmm3
  movsd  0x18(%rsp),%xmm5
  subsd  %xmm5,%xmm3
  movapd %xmm1,%xmm6
  mulsd  %xmm3,%xmm6
  movsd  %xmm6,0x10(%rsp)
  movapd %xmm5,%xmm3
  mulsd  %xmm0,%xmm3
  movapd %xmm4,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm3
  movsd  %xmm3,0x18(%rsp)
  mulsd  %xmm4,%xmm0
  mulsd  %xmm5,%xmm2
  addsd  %xmm2,%xmm0
  movapd %xmm6,%xmm1
  mulsd  %xmm6,%xmm1
  divsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L84
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x20(%rsp)
L115:
  movsd  0x10(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm7
  divsd  %xmm7,%xmm4
  movsd  %xmm4,0x10(%rsp)
  pxor   %xmm0,%xmm0
  ucomisd %xmm7,%xmm0
  ja     L110
  sqrtsd %xmm7,%xmm7
  movsd  %xmm7,0x18(%rsp)
L89:
  lea    0x40(%rsp),%r15
  mov    %r15,%rdi
  movsd  0x20(%rsp),%xmm0
  call   <YvVfilterCoef@plt>
  mov    %r15,%r8
  mov    %r12d,%ecx
  mov    %ebp,%edx
  mov    %rbx,%rsi
  mov    0x28(%rsp),%rdi
  call   <f_iir_xline_filter@plt>
  mov    %r15,%rdi
  movsd  0x18(%rsp),%xmm0
  call   <YvVfilterCoef@plt>
  movsd  0x10(%rsp),%xmm5
  ucomisd D8(%rip),%xmm5
  jp     L8
  jne    L8
  lea    0x40(%rsp),%r8
  mov    %r12d,%ecx
  mov    %ebp,%edx
  mov    %rbx,%rsi
  mov    %rbx,%rdi
  call   <f_iir_yline_filter@plt>
  jmp    L7
L84:
  call   <sqrt@plt>
  movsd  %xmm0,0x20(%rsp)
  jmp    L115
L110:
  movsd  0x18(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  %xmm0,0x18(%rsp)
  jmp    L89
L8:
  movsd  0x10(%rsp),%xmm0
  lea    0x40(%rsp),%r8
  mov    %r12d,%ecx
  mov    %ebp,%edx
  mov    %rbx,%rsi
  mov    %rbx,%rdi
  call   <f_iir_tline_filter@plt>
L7:
  mov    $0x0,%r15d
  test   %r14d,%r14d
  jle    L40
  movsd  0x8(%rsp),%xmm4
  subsd  D2(%rip),%xmm4
  movsd  %xmm4,0x10(%rsp)
L98:
  mov    $0x1,%r8d
  movsd  0x10(%rsp),%xmm0
  mov    %r12d,%ecx
  mov    %ebp,%edx
  mov    %rbx,%rsi
  mov    %rbx,%rdi
  call   <f_iir_derivative_filter@plt>
  add    $0x1,%r15d
  cmp    %r15d,%r14d
  jne    L98
L40:
  test   %r13d,%r13d
  jle    L42
  mov    $0x0,%r14d
L112:
  mov    $0x1,%r8d
  movsd  0x8(%rsp),%xmm0
  mov    %r12d,%ecx
  mov    %ebp,%edx
  mov    %rbx,%rsi
  mov    %rbx,%rdi
  call   <f_iir_derivative_filter@plt>
  add    $0x1,%r14d
  cmp    %r14d,%r13d
  jne    L112
L42:
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L100
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L100:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: tag: D23
type: f64
name: unpack_data
arguments: tag: D27
type: f64
name: unpack_data
arguments: tag: D8
type: f64
name: unpack_data
arguments: tag: D2
type: byte[8]

tool: D23: 3.141592653589793

tool: D27: 180

tool: D8: 0

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