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

```asm
<wavelet_adjoint3>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %r8d,%r13d
  push   %r12
  mov    %r9d,%r12d
  mov    $0x7,%r9d
  push   %rbp
  push   %rbx
  mov    %ecx,%ebx
  sub    $0x78,%rsp
  mov    %rdi,0x28(%rsp)
  mov    0xb8(%rsp),%ebp
  lea    0x50(%rsp),%rsi
  lea    0x60(%rsp),%rcx
  mov    %rdx,(%rsp)
  lea    0x48(%rsp),%rdi
  lea    0x58(%rsp),%rdx
  mov    0xc0(%rsp),%r14d
  movsbl 0xb0(%rsp),%r8d
  mov    %fs:0x28,%rax
  mov    %rax,0x68(%rsp)
  xor    %eax,%eax
  call   <getadjointfilter@plt>
  mov    %ebp,%eax
  cmp    %bpl,%r14b
  cmovg  %eax,%r14d
  lea    -0x1(%rbx),%eax
  mov    %eax,0x10(%rsp)
  lea    -0x1(%r13),%eax
  mov    %eax,0x14(%rsp)
  lea    -0x1(%r12),%eax
  movsbl %r14b,%esi
  mov    %eax,0x30(%rsp)
  movsbl %bpl,%eax
  mov    %eax,%ebp
  mov    %esi,0x38(%rsp)
  mov    %eax,0x3c(%rsp)
  sub    $0x1,%ebp
  js     L104
  lea    -0x1(%rsi),%eax
  mov    %ebp,%ebx
  xor    %r8d,%r8d
  mov    %eax,0x34(%rsp)
  cmp    %ebp,%esi
  jg     L87
  mov    %ebp,0x8(%rsp)
  mov    %r8,%rsi
  jmp    L103
  nopl   0x0(%rax,%rax,1)
L91:
  sub    $0x8,%rsp
  mov    %r13,%rdx
  xor    %eax,%eax
  sub    $0x1,%ebx
  push   $0x1
  push   %rcx
  mov    %r14,%rcx
  push   %r9
  mov    %r12,%r9
  push   %r8
  mov    %r15,%r8
  push   0x28(%rsp)
  push   %rsi
  mov    %r10,%rsi
  push   %rbp
  mov    %r10,0x60(%rsp)
  mov    %rdi,0x58(%rsp)
  call   *0x88(%rsp)
  add    $0x40,%rsp
  mov    0x20(%rsp),%r10
  mov    0x18(%rsp),%rdi
  cmp    0x34(%rsp),%ebx
  je     L35
L117:
  mov    (%rsp),%rsi
L103:
  mov    0x14(%rsp),%eax
  mov    %ebx,%ecx
  mov    0x30(%rsp),%r11d
  mov    0x10(%rsp),%edx
  sar    %cl,%eax
  sar    %cl,%r11d
  lea    0x1(%rax),%r9d
  add    $0x2,%eax
  sar    %cl,%edx
  lea    0x1(%r11),%ecx
  sar    %eax
  add    $0x2,%r11d
  mov    %r9d,%ebp
  mov    %ecx,%r14d
  sar    %r11d
  mov    %eax,%r10d
  lea    0x1(%rdx),%r8d
  sub    %eax,%ebp
  imul   %r11d,%r10d
  mov    %ebp,%edi
  sub    %r11d,%r14d
  add    $0x2,%edx
  imul   %r11d,%edi
  sar    %edx
  mov    %r8d,%r15d
  imul   %r14d,%eax
  sub    %edx,%r15d
  mov    %r10d,%r12d
  imul   %r14d,%ebp
  imul   %edx,%r12d
  mov    %edi,%r13d
  imul   %r15d,%r10d
  mov    %eax,%r11d
  imul   %edx,%r13d
  imul   %r15d,%edi
  movslq %r12d,%r12
  imul   %r15d,%eax
  mov    0x28(%rsp),%r15
  movslq %r10d,%r10
  imul   %edx,%r11d
  movslq %r13d,%r13
  imul   %ebp,%edx
  lea    (%r15,%r12,4),%rbp
  cmp    %ebx,0x8(%rsp)
  movslq %edi,%rdi
  lea    0x0(%rbp,%r10,4),%r12
  cmove  %r15,%rsi
  cltq
  cmp    %ebx,0x8(%rsp)
  lea    (%r12,%r13,4),%r15
  movslq %r11d,%r11
  cmovg  (%rsp),%rsi
  lea    (%r15,%rdi,4),%r14
  movslq %edx,%rdx
  lea    (%r14,%r11,4),%r13
  lea    0x0(%r13,%rax,4),%r10
  lea    (%r10,%rdx,4),%rdi
  test   %ebx,%ebx
  jne    L91
  sub    $0x8,%rsp
  mov    %r13,%rdx
  xor    %eax,%eax
  sub    $0x1,%ebx
  push   $0x1
  push   %rcx
  mov    %r14,%rcx
  push   %r9
  mov    %r12,%r9
  push   %r8
  mov    %r15,%r8
  push   0x28(%rsp)
  push   %rsi
  mov    %r10,%rsi
  push   %rbp
  mov    %r10,0x60(%rsp)
  mov    %rdi,0x58(%rsp)
  call   *0x98(%rsp)
  add    $0x40,%rsp
  mov    0x18(%rsp),%rdi
  mov    0x20(%rsp),%r10
  cmp    0x34(%rsp),%ebx
  jne    L117
L35:
  mov    %rbp,0x18(%rsp)
  mov    0x8(%rsp),%ebp
  mov    %r12,%r9
  mov    0x38(%rsp),%eax
  sub    0x3c(%rsp),%eax
  add    %eax,%ebp
  js     L11
  mov    %r9,%rbx
  mov    %r10,%rsi
  mov    %rdi,%r12
  jmp    L41
  nopl   0x0(%rax)
L105:
  sub    $0x8,%rsp
  mov    %rbx,%r9
  mov    %r15,%r8
  mov    %r12,%rdi
  push   $0x0
  sub    $0x1,%ebp
  push   %rcx
  mov    %r14,%rcx
  push   %rdx
  mov    %r13,%rdx
  push   %rax
  mov    0x28(%rsp),%rax
  push   %rax
  push   %rax
  xor    %eax,%eax
  push   0x50(%rsp)
  mov    %rsi,0x48(%rsp)
  call   *0x88(%rsp)
  mov    0x48(%rsp),%rsi
  add    $0x40,%rsp
L41:
  mov    0x30(%rsp),%edi
  mov    %ebp,%ecx
  mov    0x10(%rsp),%eax
  mov    0x14(%rsp),%edx
  sar    %cl,%edi
  sar    %cl,%eax
  sar    %cl,%edx
  mov    %edi,%ecx
  add    $0x1,%eax
  add    $0x1,%edx
  add    $0x1,%ecx
  test   %ebp,%ebp
  jne    L105
  sub    $0x8,%rsp
  mov    %r12,%rdi
  mov    %rbx,%r9
  mov    %r15,%r8
  push   $0x0
  push   %rcx
  mov    %r14,%rcx
  push   %rdx
  mov    %r13,%rdx
  push   %rax
  mov    0x28(%rsp),%rax
  push   %rax
  push   %rax
  xor    %eax,%eax
  push   0x50(%rsp)
  call   *0x98(%rsp)
  add    $0x40,%rsp
L11:
  mov    0x68(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L57
  add    $0x78,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L104:
  imul   %r13d,%ebx
  mov    0x28(%rsp),%rsi
  mov    (%rsp),%rdi
  imul   %r12d,%ebx
  lea    0x0(,%rbx,4),%edx
  call   <memcpy@plt>
  jmp    L11
  nopw   0x0(%rax,%rax,1)
L87:
  xor    %r12d,%r12d
  xor    %r9d,%r9d
  xor    %r13d,%r13d
  xor    %r15d,%r15d
  movq   $0x0,0x18(%rsp)
  xor    %r14d,%r14d
  mov    %r9,%rbx
  mov    %r12,%rsi
  jmp    L41
L57:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
int wavelet_adjoint3(float *reccovector, int colength, float *outvector, int Xlength, int Ylength, int Zlength, char Filterlength, char Levels, char Skip, char ifnotSilent) {
  float *vector = NULL;
  int scale, xlength, ylength, zlength, maxscale;
  int mxlength, mylength, mzlength;
  float *LLL = NULL, *LHH = NULL, *LLH = NULL, *LHL = NULL;
  float *HLL = NULL, *HHH = NULL, *HLH = NULL, *HHL = NULL;
  float *LLLv = NULL;
  int skipsize;
  int Size;
  int LLLsubsize;
  int LLHsubsize;
  int LHLsubsize;
  int LHHsubsize;
  int HLHsubsize;
  int HHLsubsize;
  int HHHsubsize;
  int HLLsubsize;
  int finalxlength, finalylength, finalzlength;
  int finalLLLsize;
  int ifnotallskip = 1;
  void (*adjbioD_3d)();
  void (*adjbioD_skip_3d)();
  void (*adjbioD_3d_char)();
  void (*adjbioD_skip_3d_char)();
  getadjointfilter(&adjbioD_3d, &adjbioD_skip_3d, &adjbioD_3d_char, &adjbioD_skip_3d_char, Filterlength, 7);
  Skip = (Skip < Levels ? Skip : Levels);
  Size = Xlength * Ylength * Zlength;
  xlength = 1 + ((Xlength - 1) >> Skip);
  ylength = 1 + ((Ylength - 1) >> Skip);
  zlength = 1 + ((Zlength - 1) >> Skip);
  skipsize = xlength * ylength * zlength;
  maxscale = Levels - 1;
  finalxlength = 1 + ((Xlength - 1) >> (maxscale + 1));
  finalylength = 1 + ((Ylength - 1) >> (maxscale + 1));
  finalzlength = 1 + ((Zlength - 1) >> (maxscale + 1));
  finalLLLsize = finalxlength * finalylength * finalzlength;
  if (Levels == 0)
    skipsize = Xlength * Ylength * Zlength;
  vector = outvector;
  if (maxscale < 0) {
    memcpy(outvector, reccovector, Size * sizeof(float));
    return 0;
  }
  scale = maxscale;
  LLLv = reccovector + skipsize;
  xlength = Xlength;
  ylength = Ylength;
  zlength = Zlength;
  LLLsubsize = zlength * ylength * xlength;
  while (scale >= Skip) {
    xlength = 1 + ((Xlength - 1) >> scale);
    ylength = 1 + ((Ylength - 1) >> scale);
    zlength = 1 + ((Zlength - 1) >> scale);
    mxlength = (xlength + 1) >> 1;
    mylength = (ylength + 1) >> 1;
    mzlength = (zlength + 1) >> 1;
    vector = outvector;
    LLLsubsize = mzlength * mylength * mxlength;
    LLHsubsize = mzlength * mylength * (xlength - mxlength);
    LHLsubsize = mzlength * (ylength - mylength) * mxlength;
    LHHsubsize = mzlength * (ylength - mylength) * (xlength - mxlength);
    HLLsubsize = (zlength - mzlength) * mylength * mxlength;
    HLHsubsize = (zlength - mzlength) * mylength * (xlength - mxlength);
    HHLsubsize = (zlength - mzlength) * (ylength - mylength) * mxlength;
    HHHsubsize = (zlength - mzlength) * (ylength - mylength) * (xlength - mxlength);
    if (scale == maxscale)
      LLL = reccovector;
    LLH = reccovector + LLLsubsize;
    LHL = LLH + LLHsubsize;
    LHH = LHL + LHLsubsize;
    HLL = LHH + LHHsubsize;
    HLH = HLL + HLLsubsize;
    HHL = HLH + HLHsubsize;
    HHH = HHL + HHLsubsize;
    ifnotallskip = (scale < Skip ? 0 : 1);
    if (scale < maxscale)
      LLL = outvector;
    zlength = 1 + ((Zlength - 1) >> scale);
    vector = outvector;
    if (scale == 0) {
      adjbioD_3d_char(HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, outvector, xlength, ylength, zlength, ifnotallskip);
    } else {
      adjbioD_3d(HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, vector, xlength, ylength, zlength, ifnotallskip);
    }
    LLL = vector;
    scale--;
  }
  LLL = outvector;
  while (scale >= 0) {
    ifnotallskip = 0;
    xlength = 1 + ((Xlength - 1) >> scale);
    ylength = 1 + ((Ylength - 1) >> scale);
    zlength = 1 + ((Zlength - 1) >> scale);
    if (scale == 0) {
      adjbioD_3d_char(HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, outvector, xlength, ylength, zlength, ifnotallskip);
    } else {
      adjbioD_3d(HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, vector, xlength, ylength, zlength, ifnotallskip);
    }
    LLL = vector;
    scale--;
  }
  return 0;
}
```