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

```asm
<Forward_SO3_Naive>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x98,%rsp
  mov    %edi,%r15d
  mov    %rsi,%r14
  mov    %rdx,0x30(%rsp)
  mov    %rcx,0x78(%rsp)
  mov    %r8,0x80(%rsp)
  mov    %r9,%r12
  lea    (%rdi,%rdi,1),%ebx
  mov    %ebx,%eax
  imul   %ebx,%eax
  mov    %eax,0x88(%rsp)
  imul   %ebx,%eax
  cltq
  shl    $0x3,%rax
  lea    (%r9,%rax,1),%rbp
  lea    0x0(%rbp,%rax,1),%rdi
  mov    %rax,%r10
  mov    %rdi,0x40(%rsp)
  add    %rdi,%r10
  mov    %r10,0x48(%rsp)
  movslq %ebx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdi
  add    0xd0(%rsp),%rdi
  mov    %rdi,0x28(%rsp)
  mov    %rdi,%r11
  add    %rax,%r11
  mov    %r11,0x10(%rsp)
  mov    %r11,%r10
  add    %rax,%r10
  mov    %r10,%r11
  mov    %r10,0x18(%rsp)
  add    %rax,%r11
  mov    %r11,%r9
  mov    %r11,0x20(%rsp)
  mov    %rax,0x8(%rsp)
  add    %rax,%r9
  mov    %r9,(%rsp)
  mov    %r15d,%r13d
  imul   %ebx,%r13d
  mov    %r13d,0x8c(%rsp)
  movslq %r13d,%rax
  lea    (%r9,%rax,8),%r13
  mov    %rdi,%rsi
  mov    %ebx,%edi
  call   <SinEvalPts@plt>
  mov    0x10(%rsp),%rsi
  mov    %ebx,%edi
  call   <CosEvalPts@plt>
  mov    0x18(%rsp),%rsi
  mov    %ebx,%edi
  call   <SinEvalPts2@plt>
  mov    0x20(%rsp),%rsi
  mov    %ebx,%edi
  call   <CosEvalPts2@plt>
  mov    0xd0(%rsp),%rsi
  mov    %r15d,%edi
  call   <makeweights2@plt>
  sub    $0x8,%rsp
  push   %r13
  mov    %ebx,%r9d
  mov    0x98(%rsp),%r8d
  mov    %rbp,%rcx
  mov    %r12,%rdx
  mov    0x40(%rsp),%rsi
  mov    %r14,%rdi
  call   <grid_invfourier_so3@plt>
  mov    %ebx,%ecx
  mov    0x98(%rsp),%r14d
  mov    %r14d,%edx
  mov    0x50(%rsp),%rsi
  mov    %r12,%rdi
  call   <transpose_so3@plt>
  mov    %ebx,%ecx
  mov    %r14d,%edx
  mov    0x58(%rsp),%rsi
  mov    %rbp,%rdi
  call   <transpose_so3@plt>
  mov    %r13,(%rsp)
  mov    %ebx,%r9d
  mov    %r14d,%r8d
  mov    %rbp,%rcx
  mov    %r12,%rdx
  mov    0x58(%rsp),%rsi
  mov    0x50(%rsp),%rdi
  call   <grid_invfourier_so3@plt>
  mov    %ebx,%ecx
  mov    %r14d,%edx
  mov    0x50(%rsp),%rsi
  mov    %r12,%rdi
  call   <transpose_so3@plt>
  mov    %ebx,%ecx
  mov    %r14d,%edx
  mov    0x58(%rsp),%rsi
  mov    %rbp,%rdi
  call   <transpose_so3@plt>
  add    $0x10,%rsp
  test   %r15d,%r15d
  jle    L107
  mov    %r15d,%eax
  mov    0x8(%rsp),%rbx
  imul   %rbx,%rax
  mov    $0x1,%edx
  sub    %r15d,%edx
  mov    %edx,0x60(%rsp)
  lea    -0x2(%r15),%edx
  add    $0x1,%rdx
  imul   %rbx,%rdx
  mov    %rdx,0x70(%rsp)
  mov    0x80(%rsp),%r12
  mov    0x78(%rsp),%r14
  mov    $0x0,%ebp
  add    %rbx,%rax
  mov    %rax,0x68(%rsp)
  mov    %r15d,0x58(%rsp)
L120:
  mov    0x48(%rsp),%r15
  mov    0x40(%rsp),%rax
  mov    %rax,0x30(%rsp)
  mov    $0x0,%ebx
  mov    %r15,0x38(%rsp)
  mov    0x58(%rsp),%r15d
L64:
  sub    $0x8,%rsp
  push   %r13
  push   0x10(%rsp)
  push   0x38(%rsp)
  mov    0x38(%rsp),%r9
  mov    0x30(%rsp),%r8
  mov    0x48(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %ebp,%edi
  call   <genWig_L2@plt>
  add    $0x20,%rsp
  push   %r13
  push   %r14
  mov    0xe0(%rsp),%r9
  mov    0x10(%rsp),%r8
  mov    0x40(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %ebp,%edi
  call   <wigNaiveAnalysis@plt>
  push   %r13
  push   %r12
  mov    0xf0(%rsp),%r9
  mov    0x20(%rsp),%r8
  mov    0x58(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %ebp,%edi
  call   <wigNaiveAnalysis@plt>
  mov    0x28(%rsp),%rdi
  add    %rdi,0x50(%rsp)
  add    %rdi,0x58(%rsp)
  add    $0x20,%rsp
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %ebp,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r14,%rax,8),%r14
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %ebp,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r12,%rax,8),%r12
  mov    %ebx,%eax
  add    $0x1,%ebx
  cmp    %ebx,%r15d
  jne    L64
  mov    %eax,0x54(%rsp)
  mov    0x68(%rsp),%rdx
  add    %rdx,0x40(%rsp)
  mov    0x40(%rsp),%rax
  add    %rdx,0x48(%rsp)
  mov    0x48(%rsp),%rcx
  mov    0x60(%rsp),%r15d
  test   %r15d,%r15d
  jns    L31
  mov    %rcx,0x38(%rsp)
  mov    %rax,0x30(%rsp)
L152:
  sub    $0x8,%rsp
  push   %r13
  push   0x10(%rsp)
  push   0x38(%rsp)
  mov    0x38(%rsp),%r9
  mov    0x30(%rsp),%r8
  mov    0x48(%rsp),%rcx
  mov    %ebx,%edx
  mov    %r15d,%esi
  mov    %ebp,%edi
  call   <genWig_L2@plt>
  add    $0x20,%rsp
  push   %r13
  push   %r14
  mov    0xe0(%rsp),%r9
  mov    0x10(%rsp),%r8
  mov    0x40(%rsp),%rcx
  mov    %ebx,%edx
  mov    %r15d,%esi
  mov    %ebp,%edi
  call   <wigNaiveAnalysis@plt>
  push   %r13
  push   %r12
  mov    0xf0(%rsp),%r9
  mov    0x20(%rsp),%r8
  mov    0x58(%rsp),%rcx
  mov    %ebx,%edx
  mov    %r15d,%esi
  mov    %ebp,%edi
  call   <wigNaiveAnalysis@plt>
  mov    0x28(%rsp),%rdx
  add    %rdx,0x50(%rsp)
  add    %rdx,0x58(%rsp)
  add    $0x20,%rsp
  mov    %ebx,%edx
  mov    %r15d,%esi
  mov    %ebp,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r14,%rax,8),%r14
  mov    %ebx,%edx
  mov    %r15d,%esi
  mov    %ebp,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r12,%rax,8),%r12
  add    $0x1,%r15d
  jne    L152
  mov    0x70(%rsp),%rcx
  add    %rcx,0x40(%rsp)
  add    %rcx,0x48(%rsp)
L31:
  lea    0x1(%rbp),%eax
  cmp    0x54(%rsp),%ebp
  je     L82
  mov    %eax,%ebp
  jmp    L120
L107:
  mov    0x80(%rsp),%r12
  mov    0x78(%rsp),%r14
  jmp    L145
L82:
  mov    0x58(%rsp),%r15d
L145:
  movslq 0x88(%rsp),%rax
  shl    $0x3,%rax
  mov    0x40(%rsp),%rcx
  add    %rax,%rcx
  mov    %rcx,0x40(%rsp)
  add    0x48(%rsp),%rax
  mov    %rax,0x48(%rsp)
  mov    $0x1,%eax
  sub    %r15d,%eax
  mov    %eax,%ebx
  mov    %eax,0x54(%rsp)
  jns    L13
  lea    -0x1(%r15),%eax
  lea    0x1(%rax),%rdx
  mov    0x8(%rsp),%rcx
  imul   %rcx,%rdx
  mov    %rdx,0x60(%rsp)
  imul   %rcx,%rax
  mov    %rax,0x58(%rsp)
  jmp    L75
L147:
  sub    $0x8,%rsp
  push   %r13
  push   0x10(%rsp)
  push   0x38(%rsp)
  mov    0x38(%rsp),%r9
  mov    0x30(%rsp),%r8
  mov    0x48(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <genWig_L2@plt>
  add    $0x20,%rsp
  push   %r13
  push   %r14
  mov    0xe0(%rsp),%r9
  mov    0x10(%rsp),%r8
  mov    0x40(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <wigNaiveAnalysis@plt>
  push   %r13
  push   %r12
  mov    0xf0(%rsp),%r9
  mov    0x20(%rsp),%r8
  mov    0x58(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <wigNaiveAnalysis@plt>
  mov    0x28(%rsp),%rsi
  add    %rsi,0x50(%rsp)
  add    %rsi,0x58(%rsp)
  add    $0x20,%rsp
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r14,%rax,8),%r14
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r12,%rax,8),%r12
  add    $0x1,%ebp
  cmp    %ebp,%r15d
  jne    L147
  mov    0x60(%rsp),%rcx
  add    %rcx,0x40(%rsp)
  add    %rcx,0x48(%rsp)
L149:
  mov    0x40(%rsp),%rax
  mov    0x8(%rsp),%rsi
  add    %rsi,%rax
  mov    %rax,0x40(%rsp)
  mov    0x48(%rsp),%rcx
  add    %rsi,%rcx
  mov    %rcx,0x48(%rsp)
  mov    %rcx,0x38(%rsp)
  mov    %rax,0x30(%rsp)
  mov    0x54(%rsp),%ebp
L74:
  sub    $0x8,%rsp
  push   %r13
  push   0x10(%rsp)
  push   0x38(%rsp)
  mov    0x38(%rsp),%r9
  mov    0x30(%rsp),%r8
  mov    0x48(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <genWig_L2@plt>
  add    $0x20,%rsp
  push   %r13
  push   %r14
  mov    0xe0(%rsp),%r9
  mov    0x10(%rsp),%r8
  mov    0x40(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <wigNaiveAnalysis@plt>
  push   %r13
  push   %r12
  mov    0xf0(%rsp),%r9
  mov    0x20(%rsp),%r8
  mov    0x58(%rsp),%rcx
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <wigNaiveAnalysis@plt>
  mov    0x28(%rsp),%rdi
  add    %rdi,0x50(%rsp)
  add    %rdi,0x58(%rsp)
  add    $0x20,%rsp
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r14,%rax,8),%r14
  mov    %r15d,%edx
  mov    %ebp,%esi
  mov    %ebx,%edi
  call   <howMany_so3@plt>
  cltq
  lea    (%r12,%rax,8),%r12
  add    $0x1,%ebp
  jne    L74
  mov    0x40(%rsp),%rax
  mov    0x58(%rsp),%rcx
  add    %rcx,%rax
  mov    %rax,0x40(%rsp)
  mov    0x48(%rsp),%rax
  add    %rcx,%rax
  mov    %rax,0x48(%rsp)
  add    $0x1,%ebx
  je     L13
L75:
  mov    0x48(%rsp),%rax
  mov    %rax,0x38(%rsp)
  mov    0x40(%rsp),%rax
  mov    %rax,0x30(%rsp)
  mov    $0x0,%ebp
  test   %r15d,%r15d
  jg     L147
  jmp    L149
L13:
  mov    %r15d,%edi
  call   <totalCoeffs_so3@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sdl 0x8c(%rsp),%xmm1
  movsd  D1(%rip),%xmm0
  divsd  %xmm1,%xmm0
  test   %eax,%eax
  jle    L33
  mov    %eax,%edx
  mov    $0x0,%eax
  mov    0x78(%rsp),%rcx
  mov    0x80(%rsp),%rsi
L21:
  movapd %xmm0,%xmm1
  mulsd  (%rcx,%rax,8),%xmm1
  movsd  %xmm1,(%rcx,%rax,8)
  movapd %xmm0,%xmm1
  mulsd  (%rsi,%rax,8),%xmm1
  movsd  %xmm1,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L21
L33:
  add    $0x98,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D1
type_guessed: f64

tool: D1: 3.141592653589793

assistant: ```c
void Forward_SO3_Naive(int bw, double *rdata, double *idata, double *rcoeffs, double *icoeffs, double *workspace1, double *workspace2) {
  int j, n;
  int m1, m2;
  double *t1r, *t1i, *t2r, *t2i;
  double *sinPts, *cosPts, *sinPts2, *cosPts2;
  double *wigners, *scratch;
  double *weights;
  double *rcoeffsPtr, *icoeffsPtr;
  double scalar;
  int tmpInt;
  n = 2 * bw;
  t1r = workspace1;
  t1i = workspace1 + (n * n * n);
  t2r = t1i + (n * n * n);
  t2i = t2r + (n * n * n);
  rcoeffsPtr = rcoeffs;
  icoeffsPtr = icoeffs;
  weights = workspace2;
  sinPts = weights + n;
  cosPts = sinPts + n;
  sinPts2 = cosPts + n;
  cosPts2 = sinPts2 + n;
  wigners = cosPts2 + n;
  scratch = wigners + (bw * n);
  SinEvalPts(n, sinPts);
  CosEvalPts(n, cosPts);
  SinEvalPts2(n, sinPts2);
  CosEvalPts2(n, cosPts2);
  makeweights2(bw, weights);
  grid_invfourier_so3(rdata, idata, t1r, t1i, n * n, n, scratch);
  transpose_so3(t1r, t2r, n * n, n);
  transpose_so3(t1i, t2i, n * n, n);
  grid_invfourier_so3(t2r, t2i, t1r, t1i, n * n, n, scratch);
  transpose_so3(t1r, t2r, n * n, n);
  transpose_so3(t1i, t2i, n * n, n);
  for (m1 = 0; m1 < bw; m1++) {
    for (m2 = 0; m2 < bw; m2++) {
      genWig_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2r, wigners, weights, rcoeffs, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2i, wigners, weights, icoeffs, scratch);
      t2r += n;
      t2i += n;
      rcoeffs += howMany_so3(m1, m2, bw);
      icoeffs += howMany_so3(m1, m2, bw);
    }
    t2r += n;
    t2i += n;
    for (m2 = -(bw - 1); m2 < 0; m2++) {
      genWig_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2r, wigners, weights, rcoeffs, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2i, wigners, weights, icoeffs, scratch);
      t2r += n;
      t2i += n;
      rcoeffs += howMany_so3(m1, m2, bw);
      icoeffs += howMany_so3(m1, m2, bw);
    }
  }
  t2r += n * n;
  t2i += n * n;
  for (m1 = -(bw - 1); m1 < 0; m1++) {
    for (m2 = 0; m2 < bw; m2++) {
      genWig_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2r, wigners, weights, rcoeffs, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2i, wigners, weights, icoeffs, scratch);
      t2r += n;
      t2i += n;
      rcoeffs += howMany_so3(m1, m2, bw);
      icoeffs += howMany_so3(m1, m2, bw);
    }
    t2r += n;
    t2i += n;
    for (m2 = -(bw - 1); m2 < 0; m2++) {
      genWig_L2(m1, m2, bw, sinPts, cosPts, sinPts2, cosPts2, wigners, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2r, wigners, weights, rcoeffs, scratch);
      wigNaiveAnalysis(m1, m2, bw, t2i, wigners, weights, icoeffs, scratch);
      t2r += n;
      t2i += n;
      rcoeffs += howMany_so3(m1, m2, bw);
      icoeffs += howMany_so3(m1, m2, bw);
    }
  }
  tmpInt = totalCoeffs_so3(bw);
  scalar = 3.14159265358979323846 / ((double)(bw * n));
  for (j = 0; j < tmpInt; j++) {
    rcoeffsPtr[j] *= scalar;
    icoeffsPtr[j] *= scalar;
  }
}
```