user: What would the C source look like for the assembly code provided here:

```asm
<inverse_dit_butterfly>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x28,%rsp
  mov    %rdi,-0x88(%rbp)
  mov    %rsi,-0x90(%rbp)
  mov    %rdx,-0x98(%rbp)
  mov    %rcx,-0xa0(%rbp)
  mov    -0x90(%rbp),%rax
  sar    %rax
  mov    %rax,-0x70(%rbp)
  mov    -0x90(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x50(%rbp)
  movq   $0x2,-0x68(%rbp)
  jmp    L43
L51:
  mov    -0x88(%rbp),%rax
  mov    %rax,-0x60(%rbp)
  mov    -0x68(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x58(%rbp)
  jmp    L52
L30:
  mov    -0x58(%rbp),%rax
  mov    %rax,-0x48(%rbp)
  movq   $0x0,-0x78(%rbp)
  jmp    L28
L29:
  mov    -0x78(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0x78(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x60(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x60(%rbp),%rax
  movsd  0x8(%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0x58(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x58(%rbp),%rax
  movsd  0x8(%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x20(%rbp),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  movsd  -0x18(%rbp),%xmm1
  mulsd  -0x38(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x38(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x30(%rbp),%xmm0
  addsd  -0x10(%rbp),%xmm0
  mov    -0x60(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x60(%rbp),%rax
  add    $0x8,%rax
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x30(%rbp),%xmm0
  subsd  -0x10(%rbp),%xmm0
  mov    -0x58(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x58(%rbp),%rax
  add    $0x8,%rax
  movsd  -0x28(%rbp),%xmm0
  subsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x70(%rbp),%rax
  add    %rax,-0x78(%rbp)
  addq   $0x10,-0x60(%rbp)
  addq   $0x10,-0x58(%rbp)
L28:
  mov    -0x60(%rbp),%rax
  cmp    -0x48(%rbp),%rax
  jb     L29
  mov    -0x58(%rbp),%rax
  mov    %rax,-0x60(%rbp)
  mov    -0x68(%rbp),%rax
  shl    $0x3,%rax
  add    %rax,-0x58(%rbp)
L52:
  mov    -0x58(%rbp),%rax
  cmp    -0x50(%rbp),%rax
  jb     L30
  shlq   -0x68(%rbp)
  sarq   -0x70(%rbp)
L43:
  cmpq   $0x0,-0x70(%rbp)
  jg     L51
  nop
  nop
  leave
  ret
```

assistant: ```c
void inverse_dit_butterfly(double *data, long size, double *cos_ang, double *sin_ang) {
  long angle, astep, dl;
  double xr, yr, xi, yi, wr, wi, dr, di;
  double *l1, *l2, *end, *ol2;
  astep = size >> 1;
  end = data + size + size;
  for (dl = 2; astep > 0; dl += dl, astep >>= 1) {
    l1 = data;
    l2 = data + dl;
    for (; l2 < end; l1 = l2, l2 = l2 + dl) {
      ol2 = l2;
      for (angle = 0; l1 < ol2; l1 += 2, l2 += 2) {
        wr = cos_ang[angle];
        wi = sin_ang[angle];
        xr = *l1;
        xi = *(l1 + 1);
        yr = *l2;
        yi = *(l2 + 1);
        dr = yr * wr - yi * wi;
        di = yr * wi + yi * wr;
        *(l1) = xr + dr;
        *(l1 + 1) = xi + di;
        *(l2) = xr - dr;
        *(l2 + 1) = xi - di;
        angle += astep;
      }
    }
  }
}
```