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

```asm
<Inverse_SO3_Naive>:
  endbr64
  lea    (%rdi,%rdi,1),%eax
  push   %r15
  mov    %r9,%r15
  push   %r14
  mov    %eax,%r10d
  push   %r13
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  mov    %rdx,%rbp
  mov    %eax,%edx
  push   %rbx
  imul   %eax,%edx
  mov    %edi,%ebx
  sub    $0xc8,%rsp
  mov    %edx,0xa0(%rsp)
  imul   %eax,%edx
  mov    0x100(%rsp),%r14
  cltq
  mov    %rcx,0x90(%rsp)
  mov    %r8,0x98(%rsp)
  mov    %edx,0xa4(%rsp)
  movslq %edx,%rdx
  mov    %rdx,0xb8(%rsp)
  shl    $0x3,%rdx
  add    %rdx,%r15
  mov    %rdx,0xa8(%rsp)
  lea    (%r15,%rdx,1),%rsi
  mov    %r9,0x28(%rsp)
  mov    %rsi,0x40(%rsp)
  add    %rdx,%rsi
  lea    0x0(,%rax,8),%rdx
  add    %rdx,%r14
  mov    %rsi,0x48(%rsp)
  mov    0x100(%rsp),%rsi
  mov    %r14,%rax
  mov    %rdx,0x8(%rsp)
  add    %rdx,%rax
  mov    %r10d,0x5c(%rsp)
  mov    %rax,0x10(%rsp)
  mov    %rax,%rcx
  mov    %r10d,%eax
  imul   %edi,%eax
  add    %rdx,%rcx
  mov    %r10d,%edi
  mov    %r14,0x20(%rsp)
  mov    %rcx,%r13
  mov    %rcx,0x18(%rsp)
  add    %rdx,%r13
  cltq
  lea    0x0(%r13,%rax,8),%rax
  mov    %rax,(%rsp)
  call   <SinEvalPts@plt>
  mov    %r14,%rsi
  mov    0x5c(%rsp),%r14d
  mov    %r14d,%edi
  call   <CosEvalPts@plt>
  mov    0x10(%rsp),%rsi
  mov    %r14d,%edi
  call   <SinEvalPts2@plt>
  mov    0x18(%rsp),%rsi
  mov    %r14d,%edi
  mov    %r14d,0x5c(%rsp)
  call   <CosEvalPts2@plt>
  mov    %ebx,%edi
  call   <totalCoeffs_so3@plt>
  mov    0x28(%rsp),%rdi
  mov    %r12,%rsi
  lea    0x0(,%rax,8),%r14d
  mov    %r14d,%edx
  call   <memcpy@plt>
  mov    %r14d,%edx
  mov    %rbp,%rsi
  mov    %r15,%rdi
  call   <memcpy@plt>
  test   %ebx,%ebx
  jle    L248
  mov    $0x1,%eax
  mov    0x8(%rsp),%rdx
  mov    0x5c(%rsp),%r10d
  xor    %r12d,%r12d
  sub    %ebx,%eax
  mov    %eax,0x70(%rsp)
  movslq %ebx,%rax
  lea    -0x1(%r10),%esi
  imul   %rdx,%rax
  lea    0x8(,%rsi,8),%rcx
  mov    %rsi,0xb0(%rsp)
  mov    %rcx,0x68(%rsp)
  mov    %rax,0x78(%rsp)
  add    %rcx,%rax
  mov    %rax,0x88(%rsp)
  lea    -0x2(%rbx),%eax
  add    $0x1,%rax
  imul   %rdx,%rax
  mov    %rax,0x80(%rsp)
  nopl   0x0(%rax)
L64:
  mov    0x48(%rsp),%rax
  mov    0x40(%rsp),%r14
  xor    %ebp,%ebp
  mov    %rax,0x30(%rsp)
  nopl   0x0(%rax)
L311:
  sub    $0x8,%rsp
  mov    %ebp,%esi
  mov    %ebx,%edx
  mov    %r12d,%edi
  push   0x8(%rsp)
  push   %r13
  push   0x30(%rsp)
  mov    0x30(%rsp),%r9
  mov    0x40(%rsp),%r8
  mov    0x120(%rsp),%rcx
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r14,%r9
  mov    %ebp,%esi
  push   0x8(%rsp)
  mov    0x38(%rsp),%rcx
  mov    %r13,%r8
  mov    %ebx,%edx
  mov    %r12d,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %rsi
  push   0x8(%rsp)
  mov    %r13,%r8
  mov    0x40(%rsp),%r9
  mov    %r15,%rcx
  mov    %ebp,%esi
  mov    %ebx,%edx
  mov    %r12d,%edi
  call   <wigNaiveSynthesis@plt>
  mov    %ebp,%esi
  mov    %ebx,%edx
  mov    %r12d,%edi
  call   <howMany_so3@plt>
  mov    0x38(%rsp),%rcx
  mov    %ebp,%esi
  mov    %r12d,%edi
  cltq
  lea    (%rcx,%rax,8),%rdx
  mov    %rdx,0x38(%rsp)
  mov    %ebx,%edx
  call   <howMany_so3@plt>
  mov    0x18(%rsp),%rcx
  add    %rcx,0x40(%rsp)
  cltq
  pop    %rdi
  pop    %r8
  lea    (%r15,%rax,8),%r15
  mov    %ebp,%eax
  add    $0x1,%ebp
  add    %rcx,%r14
  cmp    %ebp,%ebx
  jne    L311
  mov    0x48(%rsp),%rdx
  mov    0x78(%rsp),%rcx
  mov    %eax,0x50(%rsp)
  mov    0x40(%rsp),%rax
  mov    %rdx,%r14
  mov    %rax,%rdi
  add    %rcx,%r14
  add    %rcx,%rdi
  mov    0x5c(%rsp),%ecx
  test   %ecx,%ecx
  jle    L293
  mov    0x88(%rsp),%rcx
  add    %rcx,%rax
  add    %rcx,%rdx
  mov    %rax,0x40(%rsp)
  mov    %rdx,0x48(%rsp)
  cmp    %rax,%r14
  jae    L211
  cmp    %rdx,%rdi
  jb     L220
L211:
  mov    0x68(%rsp),%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    0x68(%rsp),%rdx
  xor    %esi,%esi
  mov    %r14,%rdi
  call   <memset@plt>
L233:
  mov    0x70(%rsp),%r14d
  test   %r14d,%r14d
  je     L248
  mov    0x48(%rsp),%rax
  mov    %ebx,0x60(%rsp)
  mov    %r14d,%ebx
  mov    0x28(%rsp),%r14
  mov    %rax,0x38(%rsp)
  mov    0x40(%rsp),%rax
  mov    %rax,0x30(%rsp)
  nopl   0x0(%rax)
L93:
  sub    $0x8,%rsp
  mov    %ebp,%edx
  mov    %ebx,%esi
  mov    %r12d,%edi
  push   0x8(%rsp)
  push   %r13
  push   0x30(%rsp)
  mov    0x40(%rsp),%r8
  mov    0x120(%rsp),%rcx
  mov    0x30(%rsp),%r9
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r14,%rcx
  mov    %r13,%r8
  push   0x8(%rsp)
  mov    0x40(%rsp),%r9
  mov    %ebp,%edx
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %rax
  push   0x8(%rsp)
  mov    %r15,%rcx
  mov    0x48(%rsp),%r9
  mov    %r13,%r8
  mov    %ebp,%edx
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <wigNaiveSynthesis@plt>
  mov    %ebp,%edx
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <howMany_so3@plt>
  mov    %ebp,%edx
  mov    %ebx,%esi
  mov    %r12d,%edi
  cltq
  lea    (%r14,%rax,8),%r14
  call   <howMany_so3@plt>
  mov    0x18(%rsp),%rcx
  add    %rcx,0x40(%rsp)
  cltq
  add    %rcx,0x48(%rsp)
  lea    (%r15,%rax,8),%r15
  pop    %rax
  pop    %rdx
  add    $0x1,%ebx
  jne    L93
  mov    %r14,0x28(%rsp)
  mov    0x60(%rsp),%ebx
  lea    0x1(%r12),%eax
  mov    0x80(%rsp),%rdx
  add    %rdx,0x40(%rsp)
  add    %rdx,0x48(%rsp)
  cmp    0x50(%rsp),%r12d
  je     L248
  mov    %eax,%r12d
  jmp    L64
  nopl   0x0(%rax)
L248:
  mov    0xa0(%rsp),%eax
  test   %eax,%eax
  je     L109
  movslq %eax,%rdx
  mov    0x40(%rsp),%rcx
  mov    0x48(%rsp),%rsi
  lea    0x0(,%rdx,8),%r12
  lea    (%rcx,%r12,1),%rax
  lea    (%rsi,%r12,1),%rdi
  mov    %rax,0x38(%rsp)
  mov    %rdi,0x50(%rsp)
  cmp    %rsi,%rax
  jbe    L310
  cmp    %rcx,%rdi
  ja     L127
L310:
  mov    0x40(%rsp),%rdi
  mov    %r12,%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    0x48(%rsp),%rdi
  mov    %r12,%rdx
  xor    %esi,%esi
  call   <memset@plt>
L104:
  mov    $0x1,%eax
  sub    %ebx,%eax
  mov    %eax,0x68(%rsp)
  mov    %eax,%ebp
  jns    L18
  mov    0x8(%rsp),%rcx
  lea    -0x1(%rbx),%eax
  mov    0x28(%rsp),%r12
  lea    0x1(%rax),%rdx
  imul   %rcx,%rdx
  imul   %rcx,%rax
  mov    %rdx,0x78(%rsp)
  mov    0x5c(%rsp),%edx
  mov    %rax,0x70(%rsp)
  sub    $0x1,%edx
  mov    %rdx,0x80(%rsp)
  lea    0x8(,%rdx,8),%rdx
  mov    %rdx,0x60(%rsp)
  nopl   0x0(%rax)
L267:
  mov    0x50(%rsp),%rax
  xor    %r14d,%r14d
  mov    %rax,0x30(%rsp)
  mov    0x38(%rsp),%rax
  mov    %rax,0x28(%rsp)
  test   %ebx,%ebx
  jle    L209
  nopw   0x0(%rax,%rax,1)
L155:
  sub    $0x8,%rsp
  mov    %r14d,%esi
  mov    %ebx,%edx
  mov    %ebp,%edi
  push   0x8(%rsp)
  push   %r13
  push   0x30(%rsp)
  mov    0x40(%rsp),%r8
  mov    0x120(%rsp),%rcx
  mov    0x30(%rsp),%r9
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r12,%rcx
  mov    %r14d,%esi
  push   0x8(%rsp)
  mov    0x38(%rsp),%r9
  mov    %r13,%r8
  mov    %ebx,%edx
  mov    %ebp,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %r10
  push   0x8(%rsp)
  mov    %r15,%rcx
  mov    0x40(%rsp),%r9
  mov    %r13,%r8
  mov    %r14d,%esi
  mov    %ebx,%edx
  mov    %ebp,%edi
  call   <wigNaiveSynthesis@plt>
  mov    %r14d,%esi
  mov    %ebx,%edx
  mov    %ebp,%edi
  call   <howMany_so3@plt>
  mov    %r14d,%esi
  mov    %ebx,%edx
  mov    %ebp,%edi
  cltq
  add    $0x1,%r14d
  lea    (%r12,%rax,8),%r12
  call   <howMany_so3@plt>
  mov    0x18(%rsp),%rsi
  add    %rsi,0x38(%rsp)
  cltq
  add    %rsi,0x40(%rsp)
  pop    %r11
  lea    (%r15,%rax,8),%r15
  pop    %rax
  cmp    %r14d,%ebx
  jne    L155
  mov    0x78(%rsp),%rdx
  add    %rdx,0x38(%rsp)
  add    %rdx,0x50(%rsp)
L209:
  mov    0x5c(%rsp),%r9d
  test   %r9d,%r9d
  jle    L87
  mov    0x38(%rsp),%rax
  mov    0x50(%rsp),%rsi
  mov    0x60(%rsp),%rdi
  mov    %rax,%rdx
  mov    %rsi,%rcx
  add    %rdi,%rdx
  add    %rdi,%rcx
  mov    %rdx,0x40(%rsp)
  mov    %rcx,0x48(%rsp)
  cmp    %rdx,%rsi
  jae    L202
  cmp    %rcx,%rax
  jb     L73
L202:
  mov    0x60(%rsp),%r14
  mov    0x38(%rsp),%rdi
  xor    %esi,%esi
  mov    %r14,%rdx
  call   <memset@plt>
  mov    0x50(%rsp),%rdi
  mov    %r14,%rdx
  xor    %esi,%esi
  call   <memset@plt>
L128:
  mov    0x48(%rsp),%rax
  mov    0x68(%rsp),%r14d
  mov    %rax,0x30(%rsp)
  mov    0x40(%rsp),%rax
  mov    %rax,0x28(%rsp)
  nopw   0x0(%rax,%rax,1)
L320:
  sub    $0x8,%rsp
  mov    %ebx,%edx
  mov    %r14d,%esi
  mov    %ebp,%edi
  push   0x8(%rsp)
  push   %r13
  push   0x30(%rsp)
  mov    0x40(%rsp),%r8
  mov    0x120(%rsp),%rcx
  mov    0x30(%rsp),%r9
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r12,%rcx
  mov    %r13,%r8
  push   0x8(%rsp)
  mov    0x38(%rsp),%r9
  mov    %ebx,%edx
  mov    %r14d,%esi
  mov    %ebp,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %rsi
  push   0x8(%rsp)
  mov    %r13,%r8
  mov    0x40(%rsp),%r9
  mov    %r15,%rcx
  mov    %ebx,%edx
  mov    %r14d,%esi
  mov    %ebp,%edi
  call   <wigNaiveSynthesis@plt>
  mov    %ebx,%edx
  mov    %r14d,%esi
  mov    %ebp,%edi
  call   <howMany_so3@plt>
  mov    %ebp,%edi
  mov    %ebx,%edx
  mov    %r14d,%esi
  cltq
  lea    (%r12,%rax,8),%r12
  call   <howMany_so3@plt>
  mov    0x18(%rsp),%rdi
  add    %rdi,0x38(%rsp)
  cltq
  add    %rdi,0x40(%rsp)
  pop    %rdi
  lea    (%r15,%rax,8),%r15
  pop    %r8
  add    $0x1,%r14d
  jne    L320
  mov    0x70(%rsp),%rdx
  mov    0x40(%rsp),%rax
  add    %rdx,%rax
  mov    %rax,0x38(%rsp)
  mov    0x48(%rsp),%rax
  add    %rdx,%rax
  mov    %rax,0x50(%rsp)
  add    $0x1,%ebp
  jne    L267
  mov    %r12,0x28(%rsp)
L18:
  mov    %ebx,%edi
  call   <totalCoeffs_so3@plt>
  mov    0x28(%rsp),%rbp
  mov    %ebx,%edi
  cltq
  shl    $0x3,%rax
  sub    %rax,%rbp
  call   <totalCoeffs_so3@plt>
  mov    0x38(%rsp),%r14
  mov    0x50(%rsp),%r13
  mov    %rbp,%rsi
  cltq
  mov    0xa0(%rsp),%ecx
  shl    $0x3,%rax
  sub    %rax,%r15
  mov    0xa8(%rsp),%rax
  mov    %r15,%r12
  mov    0x5c(%rsp),%r15d
  sub    %rax,%r14
  sub    %rax,%r13
  mov    %r14,%rdi
  mov    %r15d,%edx
  call   <transpose_so3@plt>
  mov    %r13,%rdi
  mov    %r15d,%edx
  mov    %r12,%rsi
  mov    0xa0(%rsp),%ecx
  call   <transpose_so3@plt>
  sub    $0x8,%rsp
  mov    %r15d,%r9d
  mov    %r13,%rcx
  push   0x8(%rsp)
  mov    %r14,%rdx
  mov    %r12,%rsi
  mov    %rbp,%rdi
  mov    0xb0(%rsp),%r8d
  call   <grid_fourier_so3@plt>
  mov    %r14,%rdi
  mov    %r15d,%edx
  mov    %rbp,%rsi
  mov    0xb0(%rsp),%ecx
  call   <transpose_so3@plt>
  mov    %r13,%rdi
  mov    %r15d,%edx
  mov    %r12,%rsi
  mov    0xb0(%rsp),%r14d
  mov    %r14d,%ecx
  call   <transpose_so3@plt>
  pop    %rax
  mov    %r14d,%r8d
  push   0x8(%rsp)
  mov    0xa8(%rsp),%r13
  mov    %r15d,%r9d
  mov    %r12,%rsi
  mov    %rbp,%rdi
  mov    0xa0(%rsp),%r14
  mov    %r13,%rcx
  mov    %r14,%rdx
  call   <grid_fourier_so3@plt>
  pxor   %xmm1,%xmm1
  pop    %rdx
  pop    %rcx
  cvtsi2sd %r15d,%xmm1
  mov    0xa4(%rsp),%edx
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  divsd  D8(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  test   %edx,%edx
  jle    L255
  lea    0xf(%r14),%rax
  sub    %r13,%rax
  cmp    $0x1e,%rax
  jbe    L326
  shr    %edx
  unpcklpd %xmm0,%xmm0
  xor    %eax,%eax
  mov    %r14,%rcx
  shl    $0x4,%rdx
  mov    %r13,%rsi
  xchg   %ax,%ax
L172:
  movupd (%rcx,%rax,1),%xmm1
  mulpd  %xmm0,%xmm1
  movups %xmm1,(%rcx,%rax,1)
  movupd (%rsi,%rax,1),%xmm1
  mulpd  %xmm0,%xmm1
  movups %xmm1,(%rsi,%rax,1)
  add    $0x10,%rax
  cmp    %rdx,%rax
  jne    L172
L255:
  add    $0xc8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L87:
  mov    0x50(%rsp),%rax
  mov    %rax,0x48(%rsp)
  mov    0x38(%rsp),%rax
  mov    %rax,0x40(%rsp)
  jmp    L128
  nopl   0x0(%rax)
L293:
  mov    %r14,0x48(%rsp)
  mov    %rdi,0x40(%rsp)
  jmp    L233
L73:
  mov    0x80(%rsp),%rdi
  mov    %rax,%rcx
  xor    %eax,%eax
  nopl   0x0(%rax)
L304:
  mov    %rax,%rdx
  movq   $0x0,(%rcx,%rax,8)
  movq   $0x0,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %rdi,%rdx
  jne    L304
  jmp    L128
L220:
  mov    0xb0(%rsp),%rcx
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L153:
  mov    %rax,%rdx
  movq   $0x0,(%rdi,%rax,8)
  movq   $0x0,(%r14,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rcx
  jne    L153
  jmp    L233
L127:
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L35:
  movq   $0x0,(%rcx,%rax,8)
  movq   $0x0,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L35
  jmp    L104
L326:
  mov    0x90(%rsp),%rdx
  mov    0x98(%rsp),%rcx
  xor    %eax,%eax
  mov    0xb8(%rsp),%rsi
  nopl   0x0(%rax,%rax,1)
L70:
  movsd  (%rdx,%rax,8),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,(%rdx,%rax,8)
  movsd  (%rcx,%rax,8),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rsi
  jne    L70
  jmp    L255
L109:
  mov    0x48(%rsp),%rax
  mov    %rax,0x50(%rsp)
  mov    0x40(%rsp),%rax
  mov    %rax,0x38(%rsp)
  jmp    L104
```

tool_calls: name: read_data
arguments: label: D8
type: f64

tool: D8: 3.141592653589793

assistant: ```c
void Inverse_SO3_Naive(int bw, double *rcoeffs, double *icoeffs, double *rdata, double *idata, double *workspace1, double *workspace2) {
  int j, n;
  int m1, m2;
  double *t1r, *t1i, *t2r, *t2i;
  double *sinPts, *cosPts, *sinPts2, *cosPts2;
  double *wignersTrans, *scratch;
  double *rdataPtr, *idataPtr;
  double *rcoeffsPtr, *icoeffsPtr;
  double dn;
  int tmpInt;
  n = 2 * bw;
  t1r = workspace1;
  t1i = workspace1 + (n * n * n);
  t2r = t1i + (n * n * n);
  t2i = t2r + (n * n * n);
  rdataPtr = rdata;
  idataPtr = idata;
  rcoeffsPtr = rcoeffs;
  icoeffsPtr = icoeffs;
  sinPts = workspace2;
  cosPts = sinPts + n;
  sinPts2 = cosPts + n;
  cosPts2 = sinPts2 + n;
  wignersTrans = cosPts2 + n;
  scratch = wignersTrans + (bw * n);
  SinEvalPts(n, sinPts);
  CosEvalPts(n, cosPts);
  SinEvalPts2(n, sinPts2);
  CosEvalPts2(n, cosPts2);
  tmpInt = totalCoeffs_so3(bw);
  memcpy(t1r, rcoeffs, sizeof(double) * tmpInt);
  memcpy(t1i, icoeffs, sizeof(double) * tmpInt);
  for (m1 = 0; m1 < bw; m1++) {
    for (m2 = 0; m2 < bw; m2++) {
      genWigTrans_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wignersTrans, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1r, wignersTrans, t2r, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1i, wignersTrans, t2i, scratch);
      t1r += howMany_so3(m1, m2, bw);
      t1i += howMany_so3(m1, m2, bw);
      t2r += n;
      t2i += n;
    }
    for (j = 0; j < n; j++) {
      *t2r++ = 0.0;
      *t2i++ = 0.0;
    }
    for (m2 = -(bw - 1); m2 < 0; m2++) {
      genWigTrans_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wignersTrans, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1r, wignersTrans, t2r, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1i, wignersTrans, t2i, scratch);
      t1r += howMany_so3(m1, m2, bw);
      t1i += howMany_so3(m1, m2, bw);
      t2r += n;
      t2i += n;
    }
  }
  for (j = 0; j < n * n; j++) {
    *t2r++ = 0.0;
    *t2i++ = 0.0;
  }
  for (m1 = -(bw - 1); m1 < 0; m1++) {
    for (m2 = 0; m2 < bw; m2++) {
      genWigTrans_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wignersTrans, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1r, wignersTrans, t2r, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1i, wignersTrans, t2i, scratch);
      t1r += howMany_so3(m1, m2, bw);
      t1i += howMany_so3(m1, m2, bw);
      t2r += n;
      t2i += n;
    }
    for (j = 0; j < n; j++) {
      *t2r++ = 0.0;
      *t2i++ = 0.0;
    }
    for (m2 = -(bw - 1); m2 < 0; m2++) {
      genWigTrans_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wignersTrans, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1r, wignersTrans, t2r, scratch);
      wigNaiveSynthesis(m1, m2, bw, t1i, wignersTrans, t2i, scratch);
      t1r += howMany_so3(m1, m2, bw);
      t1i += howMany_so3(m1, m2, bw);
      t2r += n;
      t2i += n;
    }
  }
  t1r -= totalCoeffs_so3(bw);
  t1i -= totalCoeffs_so3(bw);
  t2r -= (n * n * n);
  t2i -= (n * n * n);
  transpose_so3(t2r, t1r, n, n * n);
  transpose_so3(t2i, t1i, n, n * n);
  grid_fourier_so3(t1r, t1i, t2r, t2i, n * n, n, scratch);
  transpose_so3(t2r, t1r, n, n * n);
  transpose_so3(t2i, t1i, n, n * n);
  grid_fourier_so3(t1r, t1i, rdata, idata, n * n, n, scratch);
  dn = (double)n;
  dn *= ((double)bw) / 3.14159265358979323846;
  for (j = 0; j < n * n * n; j++) {
    rdata[j] *= dn;
    idata[j] *= dn;
  }
}
```