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

```asm
<fwt53>:
  endbr64
  push   %r15
  mov    %ecx,%r15d
  push   %r14
  mov    %rsi,%r14
  push   %r13
  mov    %ecx,%r13d
  imul   %edx,%r13d
  push   %r12
  mov    %edx,%r12d
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  shl    $0x3,%r13d
  mov    %r13d,%edi
  sub    $0x28,%rsp
  mov    %ecx,0xc(%rsp)
  call   <malloc@plt>
  mov    %r13d,%edx
  mov    %r14,%rsi
  mov    %rax,%rdi
  mov    %rax,%rbp
  call   <memcpy@plt>
  test   %r15d,%r15d
  jle    L58
  lea    -0x3(%r12),%r15d
  movslq %r12d,%r9
  movsd  D21(%rip),%xmm6
  movsd  D59(%rip),%xmm5
  mov    %r15d,%ecx
  mov    %rbx,0x18(%rsp)
  lea    -0x1(%r12),%esi
  xor    %r10d,%r10d
  shr    $0x2,%ecx
  movapd D59(%rip),%xmm3
  movsd  D9(%rip),%xmm4
  lea    0x0(,%r9,8),%r13
  lea    -0x1(%rcx),%edi
  lea    0x8(%rbp),%rax
  xor    %r14d,%r14d
  mov    %r15d,%ebx
  lea    0x1(,%rcx,4),%r11d
  lea    0x2(,%rcx,4),%ecx
  add    $0x1,%rdi
  movapd D9(%rip),%xmm2
  mov    %ecx,0x14(%rsp)
  lea    0x10(%rbp),%rdx
  shl    $0x5,%rdi
  mov    %r11d,0x10(%rsp)
  xor    %r11d,%r11d
  xchg   %ax,%ax
L95:
  cmp    $0x1,%esi
  jle    L30
  cmp    $0x3,%ebx
  jbe    L109
  lea    -0x8(%rax),%r8
  xor    %ecx,%ecx
  nopl   0x0(%rax,%rax,1)
L51:
  movupd (%rdx,%rcx,1),%xmm0
  movupd (%r8,%rcx,1),%xmm1
  movhpd 0x10(%rdx,%rcx,1),%xmm0
  movhpd 0x10(%r8,%rcx,1),%xmm1
  addpd  %xmm1,%xmm0
  movupd 0x10(%r8,%rcx,1),%xmm1
  movlpd 0x8(%r8,%rcx,1),%xmm1
  mulpd  %xmm3,%xmm0
  addpd  %xmm1,%xmm0
  movlpd %xmm0,(%rax,%rcx,1)
  add    $0x20,%rcx
  movhpd %xmm0,-0x10(%rcx,%rax,1)
  cmp    %rcx,%rdi
  jne    L51
  mov    0x10(%rsp),%ecx
L161:
  movslq %ecx,%r8
  add    $0x2,%ecx
  add    %r10,%r8
  shl    $0x3,%r8
  lea    0x0(%rbp,%r8,1),%r15
  movsd  0x8(%rbp,%r8,1),%xmm0
  movsd  -0x8(%r15),%xmm1
  addsd  %xmm0,%xmm1
  mulsd  %xmm5,%xmm1
  addsd  (%r15),%xmm1
  movsd  %xmm1,(%r15)
  cmp    %ecx,%esi
  jle    L30
  add    %r11d,%ecx
  movslq %ecx,%rcx
  shl    $0x3,%rcx
  addsd  0x8(%rbp,%rcx,1),%xmm0
  lea    0x0(%rbp,%rcx,1),%r8
  mulsd  %xmm5,%xmm0
  addsd  (%r8),%xmm0
  movsd  %xmm0,(%r8)
L30:
  movsd  -0x10(%rax,%r9,8),%xmm0
  subsd  -0x18(%rax,%r9,8),%xmm0
  movsd  %xmm0,-0x10(%rax,%r9,8)
  cmp    $0x2,%r12d
  jle    L59
  cmp    $0x3,%ebx
  jbe    L133
  lea    0x10(%rax),%r8
  xor    %ecx,%ecx
L141:
  movupd (%r8,%rcx,1),%xmm0
  movupd (%rax,%rcx,1),%xmm1
  movhpd 0x10(%r8,%rcx,1),%xmm0
  movhpd 0x10(%rax,%rcx,1),%xmm1
  addpd  %xmm1,%xmm0
  movupd 0x10(%rax,%rcx,1),%xmm1
  movlpd 0x8(%rax,%rcx,1),%xmm1
  mulpd  %xmm2,%xmm0
  addpd  %xmm1,%xmm0
  movlpd %xmm0,(%rdx,%rcx,1)
  add    $0x20,%rcx
  movhpd %xmm0,-0x10(%rcx,%rdx,1)
  cmp    %rcx,%rdi
  jne    L141
  mov    0x14(%rsp),%ecx
L124:
  movslq %ecx,%r8
  add    $0x2,%ecx
  add    %r10,%r8
  shl    $0x3,%r8
  lea    0x0(%rbp,%r8,1),%r15
  movsd  0x8(%rbp,%r8,1),%xmm0
  movsd  -0x8(%r15),%xmm1
  addsd  %xmm0,%xmm1
  mulsd  %xmm4,%xmm1
  addsd  (%r15),%xmm1
  movsd  %xmm1,(%r15)
  cmp    %r12d,%ecx
  jge    L59
  add    %r11d,%ecx
  movslq %ecx,%rcx
  shl    $0x3,%rcx
  addsd  0x8(%rbp,%rcx,1),%xmm0
  lea    0x0(%rbp,%rcx,1),%r8
  mulsd  %xmm4,%xmm0
  addsd  (%r8),%xmm0
  movsd  %xmm0,(%r8)
L59:
  movsd  (%rax),%xmm0
  lea    0x1(%r14),%ecx
  add    %r12d,%r11d
  add    %r9,%r10
  add    %r13,%rdx
  mulsd  %xmm6,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  add    %r13,%rax
  cmp    %ecx,0xc(%rsp)
  je     L5
  mov    %ecx,%r14d
  jmp    L95
  nopl   0x0(%rax)
L58:
  add    $0x28,%rsp
  mov    %rbp,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <free@plt>
  xchg   %ax,%ax
L5:
  mov    %r12d,%r8d
  mov    0x18(%rsp),%rbx
  sar    %r8d
  test   %r12d,%r12d
  jle    L58
  movsd  D45(%rip),%xmm2
  movsd  D44(%rip),%xmm1
  mov    %rbp,%r9
  xor    %edi,%edi
  nopl   (%rax)
L138:
  xor    %edx,%edx
  jmp    L43
  nopl   0x0(%rax)
L126:
  mulsd  %xmm1,%xmm0
  imul   %ecx,%eax
  add    %edi,%eax
  cltq
  movsd  %xmm0,(%rbx,%rax,8)
  lea    0x1(%rdx),%rax
  cmp    %rdx,%rsi
  je     L172
L108:
  mov    %rax,%rdx
L43:
  mov    %edx,%eax
  movsd  (%r9,%rdx,8),%xmm0
  sar    %eax
  test   $0x1,%dl
  je     L126
  mulsd  %xmm2,%xmm0
  add    %r8d,%eax
  imul   %ecx,%eax
  add    %edi,%eax
  cltq
  movsd  %xmm0,(%rbx,%rax,8)
  lea    0x1(%rdx),%rax
  cmp    %rdx,%rsi
  jne    L108
L172:
  lea    0x1(%rdi),%eax
  add    %r13,%r9
  cmp    %edi,%r14d
  je     L58
  mov    %eax,%edi
  jmp    L138
L133:
  mov    $0x2,%ecx
  jmp    L124
L109:
  mov    $0x1,%ecx
  jmp    L161
```

tool_calls: name: unpack_data
arguments: label: D21
datatype: f64
name: unpack_data
arguments: label: D59
datatype: byte[16]
name: unpack_data
arguments: label: D9
datatype: f64[2]
name: unpack_data
arguments: label: D45
datatype: byte[8]
name: unpack_data
arguments: label: D44
datatype: byte[8]

tool: D21: 0.5

tool: {"D59": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0xbf,0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0xbf]}

tool: {"D9": [0.25,0.25]}

tool: {"D45": [0xcd,0x3b,0x7f,0x66,0x9e,0xa0,0xe6,0x3f]}

tool: {"D44": [0xcd,0x3b,0x7f,0x66,0x9e,0xa0,0xf6,0x3f]}

assistant: ```c
void fwt53(double *wav, const double *x, const int nx, const int ny) {
  const double a0 = -0.5;
  const double a1 = 0.25;
  const double s0 = sqrt(2.0);
  const double s1 = sqrt(2.0) * 0.5;
  register int i, j, off;
  double *tempx = malloc(nx * ny * sizeof(double));
  memcpy(tempx, x, nx * ny * sizeof(double));
  for (j = 0; j < ny; ++j) {
    off = nx * j;
    for (i = 1; i < nx - 1; i += 2) {
      tempx[off + i] += a0 * (tempx[off + i - 1] + tempx[off + i + 1]);
    }
    tempx[off + nx - 1] += 2 * a0 * tempx[off + nx - 2];
    for (i = 2; i < nx; i += 2) {
      tempx[off + i] += a1 * (tempx[off + i - 1] + tempx[off + i + 1]);
    }
    tempx[off] += 2 * a1 * tempx[off + 1];
  }
  for (j = 0; j < ny; ++j) {
    off = nx * j;
    for (i = 0; i < nx; ++i) {
      if (i % 2 == 0)
        wav[(i / 2) * ny + j] = s0 * tempx[off + i];
      else
        wav[(nx / 2 + i / 2) * ny + j] = s1 * tempx[off + i];
    }
  }
  free(tempx);
}
```