user: What is the c source code of the assembly code below?

```asm
<dwt_decode_1>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rdi,%rbx
  mov    %esi,0x1c(%rsp)
  mov    %edx,%r14d
  mov    %esi,%ebp
  shr    $0x1f,%ebp
  add    %esi,%ebp
  sar    %ebp
  lea    -0x1(%rsi),%eax
  mov    %eax,%r12d
  shr    $0x1f,%r12d
  add    %eax,%r12d
  sar    %r12d
  lea    0x1(%r12),%r13d
  mov    $0x0,%eax
  call   <dwt_interleave@plt>
  test   %r13d,%r13d
  jle    L173
  mov    %r14d,%esi
  lea    -0x1(%rbp,%rbp,1),%eax
  imul   %r14d,%eax
  cltq
  lea    (%rbx,%rax,4),%r15
  movslq %r14d,%rax
  lea    (%rbx,%rax,4),%rax
  mov    %rax,0x10(%rsp)
  lea    (%r14,%r14,1),%r11d
  movslq %r11d,%rax
  shl    $0x2,%rax
  mov    %rax,0x8(%rsp)
  mov    %rbx,%rcx
  mov    %r11d,%eax
  neg    %eax
  mov    %eax,0x18(%rsp)
  mov    $0x0,%eax
  jmp    L24
L95:
  mov    0x10(%rsp),%r10
  mov    (%r10),%r9d
  jns    L182
  lea    0x2(%r9,%r9,1),%edi
  sar    $0x2,%edi
  sub    %edi,%edx
  mov    %edx,(%rcx)
  jmp    L100
L69:
  mov    0x18(%rsp),%r10d
  lea    (%rsi,%r10,1),%r9d
  movslq %r9d,%r9
  mov    (%rbx,%r9,4),%r9d
L182:
  cmp    %edi,%ebp
  jle    L122
  movslq %esi,%r10
  mov    (%rbx,%r10,4),%r10d
L179:
  lea    0x2(%r9,%r10,1),%r9d
  sar    $0x2,%r9d
  sub    %r9d,%edx
  mov    %edx,(%r8)
  cmp    %edi,%r12d
  jle    L173
L100:
  add    $0x1,%eax
  add    0x8(%rsp),%rcx
  add    %r11d,%esi
L24:
  mov    %rcx,%r8
  mov    (%rcx),%edx
  mov    %eax,%edi
  test   %eax,%eax
  jle    L95
  cmp    %eax,%ebp
  jge    L69
  mov    (%r15),%r9d
L122:
  mov    (%r15),%r10d
  jmp    L179
L173:
  cmpl   $0x1,0x1c(%rsp)
  jle    L79
  lea    (%r14,%r14,1),%r11d
  mov    %r12d,%eax
  imul   %r11d,%eax
  cltq
  lea    (%rbx,%rax,4),%r15
  movslq %r11d,%rax
  shl    $0x2,%rax
  mov    %rax,0x8(%rsp)
  movslq %r14d,%r14
  lea    (%rbx,%r14,4),%rcx
  mov    %r11d,%r14d
  neg    %r14d
  mov    %r11d,%edi
  mov    $0x1,%r8d
  mov    $0x0,%edx
  jmp    L35
L55:
  mov    (%rbx),%esi
  cmp    $0xffffffff,%edx
  jge    L126
  add    %esi,%esi
  sar    %esi
  add    %r10d,%esi
  mov    %esi,(%rcx)
  jmp    L87
L108:
  lea    (%rdi,%r14,1),%esi
  movslq %esi,%rsi
  mov    (%rbx,%rsi,4),%esi
  jmp    L126
L184:
  movslq %edi,%rax
  mov    (%rbx,%rax,4),%eax
L63:
  add    %esi,%eax
  sar    %eax
  add    %r10d,%eax
  mov    %eax,(%r9)
  cmp    %r8d,%ebp
  jle    L79
L87:
  add    $0x1,%edx
  add    0x8(%rsp),%rcx
  add    $0x1,%r8d
  add    %r11d,%edi
L35:
  mov    %rcx,%r9
  mov    (%rcx),%r10d
  mov    %edx,%eax
  test   %edx,%edx
  js     L55
  cmp    %edx,%r13d
  jg     L108
  mov    (%r15),%esi
L126:
  cmp    %eax,%r12d
  jg     L184
  mov    (%r15),%eax
  jmp    L63
L79:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void dwt_decode_1(int *a, int n, int x) {
  int dn, sn, i;
  dn = n / 2;
  sn = 1 + (n - 1) / 2;
  dwt_interleave(a, n, x);
  for (i = 0; i < sn; i++)
    a[x * (i)*2] -= (((i - 1) < 0 ? a[x * (1 + (0) * 2)] : ((i - 1) >= dn ? a[x * (1 + (dn - 1) * 2)] : a[x * (1 + (i - 1) * 2)])) + ((i) < 0 ? a[x * (1 + (0) * 2)] : ((i) >= dn ? a[x * (1 + (dn - 1) * 2)] : a[x * (1 + (i)*2)])) + 2) >> 2;
  for (i = 0; i < dn; i++)
    a[x * (1 + (i)*2)] += (((i) < 0 ? a[x * (0) * 2] : ((i) >= sn ? a[x * (sn - 1) * 2] : a[x * (i)*2])) + ((i + 1) < 0 ? a[x * (0) * 2] : ((i + 1) >= sn ? a[x * (sn - 1) * 2] : a[x * (i + 1) * 2]))) >> 1;
}
```