user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<genWigAllTrans>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  mov    %edi,-0x44(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %rdx,-0x58(%rbp)
  mov    -0x44(%rbp),%eax
  add    %eax,%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x58(%rbp),%rax
  mov    %rax,-0x28(%rbp)
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x20(%rbp)
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x28(%rbp),%rdx
  mov    -0x2c(%rbp),%eax
  mov    %rdx,%rsi
  mov    %eax,%edi
  call   <SinEvalPts@plt>
  mov    -0x20(%rbp),%rdx
  mov    -0x2c(%rbp),%eax
  mov    %rdx,%rsi
  mov    %eax,%edi
  call   <CosEvalPts@plt>
  mov    -0x18(%rbp),%rdx
  mov    -0x2c(%rbp),%eax
  mov    %rdx,%rsi
  mov    %eax,%edi
  call   <SinEvalPts2@plt>
  mov    -0x10(%rbp),%rdx
  mov    -0x2c(%rbp),%eax
  mov    %rdx,%rsi
  mov    %eax,%edi
  call   <CosEvalPts2@plt>
  mov    -0x18(%rbp),%rsi
  mov    -0x20(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    -0x44(%rbp),%eax
  sub    $0x8,%rsp
  push   -0x8(%rbp)
  push   -0x50(%rbp)
  push   -0x10(%rbp)
  mov    %rsi,%r9
  mov    %rcx,%r8
  mov    %rdx,%rcx
  mov    %eax,%edx
  mov    $0x0,%esi
  mov    $0x0,%edi
  call   <genWigTrans_L2@plt>
  add    $0x20,%rsp
  mov    -0x44(%rbp),%eax
  imul   -0x2c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x50(%rbp)
  movl   $0x1,-0x34(%rbp)
  jmp    L34
L64:
  mov    -0x18(%rbp),%r8
  mov    -0x20(%rbp),%rdi
  mov    -0x28(%rbp),%rcx
  mov    -0x44(%rbp),%edx
  mov    -0x34(%rbp),%esi
  mov    -0x34(%rbp),%eax
  sub    $0x8,%rsp
  push   -0x8(%rbp)
  push   -0x50(%rbp)
  push   -0x10(%rbp)
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <genWigTrans_L2@plt>
  add    $0x20,%rsp
  mov    -0x44(%rbp),%eax
  sub    -0x34(%rbp),%eax
  imul   -0x2c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x50(%rbp)
  addl   $0x1,-0x34(%rbp)
L34:
  mov    -0x34(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L64
  movl   $0x1,-0x34(%rbp)
  jmp    L12
L55:
  mov    -0x18(%rbp),%rdi
  mov    -0x20(%rbp),%rsi
  mov    -0x28(%rbp),%rcx
  mov    -0x44(%rbp),%edx
  mov    -0x34(%rbp),%eax
  sub    $0x8,%rsp
  push   -0x8(%rbp)
  push   -0x50(%rbp)
  push   -0x10(%rbp)
  mov    %rdi,%r9
  mov    %rsi,%r8
  mov    $0x0,%esi
  mov    %eax,%edi
  call   <genWigTrans_L2@plt>
  add    $0x20,%rsp
  mov    -0x44(%rbp),%eax
  sub    -0x34(%rbp),%eax
  imul   -0x2c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x50(%rbp)
  addl   $0x1,-0x34(%rbp)
L12:
  mov    -0x34(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L55
  movl   $0x1,-0x34(%rbp)
  jmp    L28
L58:
  mov    -0x34(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x30(%rbp)
  jmp    L56
L6:
  mov    -0x18(%rbp),%r8
  mov    -0x20(%rbp),%rdi
  mov    -0x28(%rbp),%rcx
  mov    -0x44(%rbp),%edx
  mov    -0x30(%rbp),%esi
  mov    -0x34(%rbp),%eax
  sub    $0x8,%rsp
  push   -0x8(%rbp)
  push   -0x50(%rbp)
  push   -0x10(%rbp)
  mov    %r8,%r9
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <genWigTrans_L2@plt>
  add    $0x20,%rsp
  mov    -0x44(%rbp),%eax
  sub    -0x30(%rbp),%eax
  imul   -0x2c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,-0x50(%rbp)
  addl   $0x1,-0x30(%rbp)
L56:
  mov    -0x30(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L6
  addl   $0x1,-0x34(%rbp)
L28:
  mov    -0x34(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L58
  nop
  nop
  leave
  ret
```

assistant: ```c
void genWigAllTrans(int bw, double *wigners, double *workspace) {
  int m1, m2, n;
  double *sinPts, *cosPts;
  double *sinPts2, *cosPts2;
  double *scratch;
  n = 2 * bw;
  sinPts = workspace;
  cosPts = sinPts + n;
  sinPts2 = cosPts + n;
  cosPts2 = sinPts2 + n;
  scratch = cosPts2 + n;
  SinEvalPts(n, sinPts);
  CosEvalPts(n, cosPts);
  SinEvalPts2(n, sinPts2);
  CosEvalPts2(n, cosPts2);
  genWigTrans_L2(0, 0, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
  wigners += bw * n;
  for (m1 = 1; m1 < bw; m1++) {
    genWigTrans_L2(m1, m1, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
    wigners += n * (bw - m1);
  }
  for (m1 = 1; m1 < bw; m1++) {
    genWigTrans_L2(m1, 0, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
    wigners += n * (bw - m1);
  }
  for (m1 = 1; m1 < bw; m1++)
    for (m2 = m1 + 1; m2 < bw; m2++) {
      genWigTrans_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
      wigners += n * (bw - m2);
    }
}
```