user: 下面的汇编指令对应的C源代码是什么:

```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,0x90(%rsp)
  imul   %eax,%edx
  mov    0x100(%rsp),%r14
  cltq
  mov    %rcx,0xa0(%rsp)
  mov    0x100(%rsp),%rsi
  mov    %r8,0xa8(%rsp)
  mov    %edx,0x94(%rsp)
  movslq %edx,%rdx
  mov    %rdx,0xb0(%rsp)
  shl    $0x3,%rdx
  add    %rdx,%r15
  mov    %rdx,0x98(%rsp)
  lea    (%r15,%rdx,1),%rcx
  mov    %r9,0x30(%rsp)
  mov    %rcx,0x48(%rsp)
  add    %rdx,%rcx
  lea    0x0(,%rax,8),%rdx
  add    %rdx,%r14
  mov    %rcx,0x50(%rsp)
  mov    %r14,%rax
  mov    %r10d,0x60(%rsp)
  add    %rdx,%rax
  mov    %rdx,0x10(%rsp)
  mov    %rax,0x18(%rsp)
  mov    %rax,%rcx
  mov    %r10d,%eax
  imul   %edi,%eax
  add    %rdx,%rcx
  mov    %r10d,%edi
  mov    %r14,0x28(%rsp)
  mov    %rcx,%r13
  mov    %rcx,0x20(%rsp)
  add    %rdx,%r13
  cltq
  lea    0x0(%r13,%rax,8),%rax
  mov    %rax,0x8(%rsp)
  call   <SinEvalPts@plt>
  mov    %r14,%rsi
  mov    0x60(%rsp),%r14d
  mov    %r14d,%edi
  call   <CosEvalPts@plt>
  mov    0x18(%rsp),%rsi
  mov    %r14d,%edi
  call   <SinEvalPts2@plt>
  mov    0x20(%rsp),%rsi
  mov    %r14d,%edi
  mov    %r14d,0x60(%rsp)
  call   <CosEvalPts2@plt>
  mov    %ebx,%edi
  call   <totalCoeffs_so3@plt>
  mov    0x30(%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>
  mov    $0x1,%eax
  sub    %ebx,%eax
  mov    %eax,0x64(%rsp)
  test   %ebx,%ebx
  jle    L282
  mov    0x10(%rsp),%rdx
  mov    0x60(%rsp),%r10d
  movslq %ebx,%rax
  xor    %r12d,%r12d
  imul   %rdx,%rax
  lea    -0x1(%r10),%ecx
  lea    0x8(,%rcx,8),%rsi
  mov    %rcx,0xb8(%rsp)
  mov    %rsi,0x70(%rsp)
  mov    %rax,0x78(%rsp)
  add    %rsi,%rax
  mov    %rax,0x88(%rsp)
  lea    -0x2(%rbx),%eax
  add    $0x1,%rax
  imul   %rdx,%rax
  mov    %rax,0x80(%rsp)
  nopl   (%rax)
L23:
  mov    0x50(%rsp),%rax
  mov    0x48(%rsp),%r14
  xor    %ebp,%ebp
  mov    %rax,0x38(%rsp)
  nopl   0x0(%rax)
L246:
  sub    $0x8,%rsp
  mov    %ebp,%esi
  mov    %ebx,%edx
  mov    %r12d,%edi
  push   0x10(%rsp)
  push   %r13
  push   0x38(%rsp)
  mov    0x38(%rsp),%r9
  mov    0x48(%rsp),%r8
  mov    0x120(%rsp),%rcx
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r14,%r9
  mov    %ebp,%esi
  push   0x10(%rsp)
  mov    0x40(%rsp),%rcx
  mov    %r13,%r8
  mov    %ebx,%edx
  mov    %r12d,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %r9
  push   0x10(%rsp)
  mov    %r15,%rcx
  mov    0x48(%rsp),%r9
  mov    %r13,%r8
  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    0x40(%rsp),%rcx
  mov    %ebp,%esi
  mov    %r12d,%edi
  cltq
  lea    (%rcx,%rax,8),%rdx
  mov    %rdx,0x40(%rsp)
  mov    %ebx,%edx
  call   <howMany_so3@plt>
  mov    0x20(%rsp),%rcx
  add    %rcx,0x48(%rsp)
  cltq
  pop    %r10
  pop    %r11
  lea    (%r15,%rax,8),%r15
  mov    %ebp,%eax
  add    $0x1,%ebp
  add    %rcx,%r14
  cmp    %ebp,%ebx
  jne    L246
  mov    0x50(%rsp),%rdx
  mov    0x78(%rsp),%rcx
  mov    %eax,0x58(%rsp)
  mov    0x48(%rsp),%rax
  mov    0x60(%rsp),%r8d
  mov    %rdx,%r14
  mov    %rax,%rdi
  add    %rcx,%r14
  add    %rcx,%rdi
  test   %r8d,%r8d
  jle    L97
  mov    0x88(%rsp),%rcx
  add    %rcx,%rax
  add    %rcx,%rdx
  mov    %rax,0x48(%rsp)
  mov    %rdx,0x50(%rsp)
  cmp    %rax,%r14
  jae    L268
  cmp    %rdx,%rdi
  jb     L277
L268:
  mov    0x70(%rsp),%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    0x70(%rsp),%rdx
  xor    %esi,%esi
  mov    %r14,%rdi
  call   <memset@plt>
L61:
  mov    0x64(%rsp),%r14d
  test   %r14d,%r14d
  je     L282
  mov    0x50(%rsp),%rax
  mov    %ebx,0x68(%rsp)
  mov    %r14d,%ebx
  mov    0x30(%rsp),%r14
  mov    %rax,0x40(%rsp)
  mov    0x48(%rsp),%rax
  mov    %rax,0x38(%rsp)
  nopl   0x0(%rax,%rax,1)
L298:
  sub    $0x8,%rsp
  mov    %ebp,%edx
  mov    %ebx,%esi
  mov    %r12d,%edi
  push   0x10(%rsp)
  push   %r13
  push   0x38(%rsp)
  mov    0x48(%rsp),%r8
  mov    0x120(%rsp),%rcx
  mov    0x38(%rsp),%r9
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r14,%rcx
  mov    %r13,%r8
  push   0x10(%rsp)
  mov    0x48(%rsp),%r9
  mov    %ebp,%edx
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %rcx
  push   0x10(%rsp)
  mov    %r15,%rcx
  mov    0x50(%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    %ebx,%esi
  mov    %r12d,%edi
  mov    %ebp,%edx
  cltq
  lea    (%r14,%rax,8),%r14
  call   <howMany_so3@plt>
  mov    0x20(%rsp),%rcx
  add    %rcx,0x48(%rsp)
  cltq
  add    %rcx,0x50(%rsp)
  pop    %rsi
  lea    (%r15,%rax,8),%r15
  pop    %rdi
  add    $0x1,%ebx
  jne    L298
  mov    %r14,0x30(%rsp)
  mov    0x68(%rsp),%ebx
  lea    0x1(%r12),%eax
  mov    0x80(%rsp),%rdx
  add    %rdx,0x48(%rsp)
  add    %rdx,0x50(%rsp)
  cmp    0x58(%rsp),%r12d
  je     L282
  mov    %eax,%r12d
  jmp    L23
  nopl   0x0(%rax)
L282:
  mov    0x90(%rsp),%eax
  test   %eax,%eax
  je     L76
  movslq %eax,%rdx
  mov    0x48(%rsp),%rsi
  mov    0x50(%rsp),%rcx
  lea    0x0(,%rdx,8),%r12
  lea    (%rsi,%r12,1),%rax
  lea    (%rcx,%r12,1),%rdi
  mov    %rax,0x40(%rsp)
  mov    %rdi,0x58(%rsp)
  cmp    %rcx,%rax
  jbe    L191
  cmp    %rdi,%rsi
  jb     L234
L191:
  mov    0x48(%rsp),%rdi
  mov    %r12,%rdx
  xor    %esi,%esi
  call   <memset@plt>
  mov    0x50(%rsp),%rdi
  mov    %r12,%rdx
  xor    %esi,%esi
  call   <memset@plt>
L235:
  mov    0x64(%rsp),%ebp
  test   %ebp,%ebp
  jns    L91
  mov    0x10(%rsp),%rsi
  lea    -0x1(%rbx),%eax
  mov    0x30(%rsp),%r12
  lea    0x1(%rax),%rdx
  imul   %rsi,%rdx
  imul   %rsi,%rax
  mov    %rdx,0x78(%rsp)
  mov    0x60(%rsp),%edx
  mov    %rax,0x70(%rsp)
  sub    $0x1,%edx
  mov    %rdx,0x80(%rsp)
  lea    0x8(,%rdx,8),%rdx
  mov    %rdx,0x68(%rsp)
  nopl   (%rax)
L270:
  mov    0x58(%rsp),%rax
  xor    %r14d,%r14d
  mov    %rax,0x38(%rsp)
  mov    0x40(%rsp),%rax
  mov    %rax,0x30(%rsp)
  test   %ebx,%ebx
  jle    L221
  nopw   0x0(%rax,%rax,1)
L49:
  sub    $0x8,%rsp
  mov    %r14d,%esi
  mov    %ebx,%edx
  mov    %ebp,%edi
  push   0x10(%rsp)
  push   %r13
  push   0x38(%rsp)
  mov    0x48(%rsp),%r8
  mov    0x120(%rsp),%rcx
  mov    0x38(%rsp),%r9
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r12,%rcx
  mov    %r14d,%esi
  push   0x10(%rsp)
  mov    0x40(%rsp),%r9
  mov    %r13,%r8
  mov    %ebx,%edx
  mov    %ebp,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %r11
  push   0x10(%rsp)
  mov    %r15,%rcx
  mov    0x48(%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    0x20(%rsp),%rsi
  add    %rsi,0x40(%rsp)
  cltq
  add    %rsi,0x48(%rsp)
  lea    (%r15,%rax,8),%r15
  pop    %rax
  pop    %rdx
  cmp    %r14d,%ebx
  jne    L49
  mov    0x78(%rsp),%rdx
  add    %rdx,0x40(%rsp)
  add    %rdx,0x58(%rsp)
L221:
  mov    0x60(%rsp),%r10d
  test   %r10d,%r10d
  jle    L171
  mov    0x40(%rsp),%rax
  mov    0x58(%rsp),%rsi
  mov    0x68(%rsp),%rdi
  mov    %rax,%rdx
  mov    %rsi,%rcx
  add    %rdi,%rdx
  add    %rdi,%rcx
  mov    %rdx,0x48(%rsp)
  mov    %rcx,0x50(%rsp)
  cmp    %rdx,%rsi
  jae    L31
  cmp    %rcx,%rax
  jb     L81
L31:
  mov    0x68(%rsp),%r14
  mov    0x40(%rsp),%rdi
  xor    %esi,%esi
  mov    %r14,%rdx
  call   <memset@plt>
  mov    0x58(%rsp),%rdi
  mov    %r14,%rdx
  xor    %esi,%esi
  call   <memset@plt>
L291:
  mov    0x50(%rsp),%rax
  mov    0x64(%rsp),%r14d
  mov    %rax,0x38(%rsp)
  mov    0x48(%rsp),%rax
  mov    %rax,0x30(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L204:
  sub    $0x8,%rsp
  mov    %ebx,%edx
  mov    %r14d,%esi
  mov    %ebp,%edi
  push   0x10(%rsp)
  push   %r13
  push   0x38(%rsp)
  mov    0x48(%rsp),%r8
  mov    0x120(%rsp),%rcx
  mov    0x38(%rsp),%r9
  call   <genWigTrans_L2@plt>
  add    $0x18,%rsp
  mov    %r12,%rcx
  mov    %r13,%r8
  push   0x10(%rsp)
  mov    0x40(%rsp),%r9
  mov    %ebx,%edx
  mov    %r14d,%esi
  mov    %ebp,%edi
  call   <wigNaiveSynthesis@plt>
  pop    %rdi
  push   0x10(%rsp)
  mov    %r13,%r8
  mov    0x48(%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    0x20(%rsp),%rdi
  add    %rdi,0x40(%rsp)
  cltq
  add    %rdi,0x48(%rsp)
  pop    %r8
  lea    (%r15,%rax,8),%r15
  pop    %r9
  add    $0x1,%r14d
  jne    L204
  mov    0x70(%rsp),%rdx
  mov    0x48(%rsp),%rax
  add    %rdx,%rax
  mov    %rax,0x40(%rsp)
  mov    0x50(%rsp),%rax
  add    %rdx,%rax
  mov    %rax,0x58(%rsp)
  add    $0x1,%ebp
  jne    L270
  mov    %r12,0x30(%rsp)
L91:
  mov    %ebx,%edi
  call   <totalCoeffs_so3@plt>
  mov    0x30(%rsp),%rbp
  mov    %ebx,%edi
  cltq
  shl    $0x3,%rax
  sub    %rax,%rbp
  call   <totalCoeffs_so3@plt>
  mov    0x40(%rsp),%r14
  mov    0x58(%rsp),%r13
  mov    %rbp,%rsi
  cltq
  mov    0x90(%rsp),%ecx
  shl    $0x3,%rax
  sub    %rax,%r15
  mov    0x98(%rsp),%rax
  mov    %r15,%r12
  mov    0x60(%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    0x90(%rsp),%ecx
  call   <transpose_so3@plt>
  sub    $0x8,%rsp
  mov    %r15d,%r9d
  mov    %r13,%rcx
  push   0x10(%rsp)
  mov    %r14,%rdx
  mov    %r12,%rsi
  mov    %rbp,%rdi
  mov    0xa0(%rsp),%r8d
  call   <grid_fourier_so3@plt>
  mov    %r14,%rdi
  mov    %r15d,%edx
  mov    %rbp,%rsi
  mov    0xa0(%rsp),%ecx
  call   <transpose_so3@plt>
  mov    %r13,%rdi
  mov    %r15d,%edx
  mov    %r12,%rsi
  mov    0xa0(%rsp),%r14d
  mov    %r14d,%ecx
  call   <transpose_so3@plt>
  pop    %rax
  mov    %r14d,%r8d
  push   0x10(%rsp)
  mov    0xb8(%rsp),%r13
  mov    %r12,%rsi
  mov    %r15d,%r9d
  mov    %rbp,%rdi
  mov    0xb0(%rsp),%r14
  mov    %r13,%rcx
  mov    %r14,%rdx
  call   <grid_fourier_so3@plt>
  pxor   %xmm1,%xmm1
  pop    %rdx
  pop    %rcx
  cvtsi2sd %r15d,%xmm1
  mov    0x94(%rsp),%esi
  pxor   %xmm0,%xmm0
  xor    %eax,%eax
  cvtsi2sd %ebx,%xmm0
  divsd  D5(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  test   %esi,%esi
  jle    L230
  mov    0xb0(%rsp),%rsi
  mov    %r14,%rdx
  mov    %r13,%rcx
  nopl   0x0(%rax)
L43:
  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    L43
L230:
  add    $0xc8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L171:
  mov    0x58(%rsp),%rax
  mov    %rax,0x50(%rsp)
  mov    0x40(%rsp),%rax
  mov    %rax,0x48(%rsp)
  jmp    L291
  nopl   0x0(%rax)
L97:
  mov    %r14,0x50(%rsp)
  mov    %rdi,0x48(%rsp)
  jmp    L61
L81:
  mov    0x40(%rsp),%rcx
  mov    0x58(%rsp),%rsi
  xor    %eax,%eax
  mov    0x80(%rsp),%rdi
  nopl   0x0(%rax,%rax,1)
L148:
  mov    %rax,%rdx
  movq   $0x0,(%rcx,%rax,8)
  movq   $0x0,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %rdi,%rdx
  jne    L148
  jmp    L291
L277:
  mov    0xb8(%rsp),%rcx
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L257:
  mov    %rax,%rdx
  movq   $0x0,(%rdi,%rax,8)
  movq   $0x0,(%r14,%rax,8)
  add    $0x1,%rax
  cmp    %rcx,%rdx
  jne    L257
  jmp    L61
L234:
  mov    0x48(%rsp),%rcx
  mov    0x50(%rsp),%rsi
  xor    %eax,%eax
  nopl   (%rax)
L307:
  movq   $0x0,(%rcx,%rax,8)
  movq   $0x0,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L307
  jmp    L235
L76:
  mov    0x50(%rsp),%rax
  mov    %rax,0x58(%rsp)
  mov    0x48(%rsp),%rax
  mov    %rax,0x40(%rsp)
  jmp    L235
```

tool_calls: name: disassemble_data
arguments: label: D5
type: f64

tool: D5: 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;
  }
}
```