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

```asm
<kern_r4_fma_n_b>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x68,%rsp
  mov    %rdi,-0xb8(%rbp)
  mov    %rsi,-0xc0(%rbp)
  mov    %rdx,-0xc8(%rbp)
  mov    %rcx,-0xd0(%rbp)
  mov    %r8,-0xd8(%rbp)
  mov    %r9,-0xe0(%rbp)
  movq   $0x0,-0xb0(%rbp)
  jmp    L38
L26:
  mov    -0xb0(%rbp),%rax
  imul   0x30(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  mov    -0xb0(%rbp),%rax
  imul   0x30(%rbp),%rax
  add    %rax,%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  mov    -0xb0(%rbp),%rax
  imul   0x30(%rbp),%rax
  shl    $0x5,%rax
  mov    %rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x90(%rbp)
  mov    -0xb0(%rbp),%rax
  imul   0x30(%rbp),%rax
  shl    $0x2,%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  mov    -0xb0(%rbp),%rax
  imul   0x30(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  mov    -0xb0(%rbp),%rax
  imul   0x30(%rbp),%rax
  add    %rax,%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movq   $0x0,-0xa8(%rbp)
  jmp    L21
L14:
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xa0(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xa0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x90(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x90(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x60(%rbp),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x60(%rbp),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x58(%rbp),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x78(%rbp),%xmm0
  movsd  -0x70(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x78(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x48(%rbp),%xmm0
  mulsd  -0x78(%rbp),%xmm0
  movsd  -0x68(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x78(%rbp),%xmm1
  movsd  -0x68(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x20(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x10(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x98(%rbp),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x40(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x98(%rbp),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x30(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x98(%rbp),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x38(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x18(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x8(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x18(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x98(%rbp),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0xa8(%rbp)
L21:
  mov    -0xa8(%rbp),%rax
  cmp    0x30(%rbp),%rax
  jl     L14
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0xd8(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0xe0(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,0x10(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,0x18(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0xb8(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0xc0(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0xc8(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0xd0(%rbp)
  addq   $0x1,-0xb0(%rbp)
L38:
  mov    -0xb0(%rbp),%rax
  cmp    0x38(%rbp),%rax
  jl     L26
  nop
  nop
  leave
  ret
```

assistant: ```c
void kern_r4_fma_n_b(double *restrict o0, double *restrict o1, double *restrict o2, double *restrict o3, double *i0, double *i1, double *i2, double *i3, double *w, double *w2, long bsize, long blocks) {
  long i, j;
  for (i = 0; i < blocks; i++) {
    double wr1, wi1;
    double wr2, wi2;
    double wr3, wi13;
    wr1 = w[2 * i * bsize];
    wi1 = w[2 * i * bsize + 1];
    wr2 = w[4 * i * bsize];
    wi2 = w[4 * i * bsize + 1];
    wr3 = w2[2 * i * bsize];
    wi13 = w2[2 * i * bsize + 1];
    for (j = 0; j < bsize; j++) {
      double pr1, pi1;
      double pr2, pi2;
      double pr3, pi3;
      double qr0, qi0;
      double qr1, qi1;
      double qr2, qi2;
      double qr3, qi3;
      pr1 = i1[2 * j] * wr1 - i1[2 * j + 1];
      pi1 = i1[2 * j] + i1[2 * j + 1] * wr1;
      pr2 = i2[2 * j] * wr2 - i2[2 * j + 1];
      pi2 = i2[2 * j] + i2[2 * j + 1] * wr2;
      pr3 = i3[2 * j] * wr3 - i3[2 * j + 1];
      pi3 = i3[2 * j] + i3[2 * j + 1] * wr3;
      qr0 = i0[2 * j] + pr2 * wi2;
      qr2 = i0[2 * j] - pr2 * wi2;
      qi0 = i0[2 * j + 1] + pi2 * wi2;
      qi2 = i0[2 * j + 1] - pi2 * wi2;
      qr1 = pr1 + pr3 * wi13;
      qr3 = pr1 - pr3 * wi13;
      qi1 = pi1 + pi3 * wi13;
      qi3 = pi1 - pi3 * wi13;
      o0[2 * j] = qr0 + qr1 * wi1;
      o0[2 * j + 1] = qi0 + qi1 * wi1;
      o2[2 * j] = qr0 - qr1 * wi1;
      o2[2 * j + 1] = qi0 - qi1 * wi1;
      o1[2 * j] = qr2 - qi3 * wi1;
      o1[2 * j + 1] = qi2 + qr3 * wi1;
      o3[2 * j] = qr2 + qi3 * wi1;
      o3[2 * j + 1] = qi2 - qr3 * wi1;
    }
    i0 += 8 * bsize;
    i1 += 8 * bsize;
    i2 += 8 * bsize;
    i3 += 8 * bsize;
    o0 += 2 * bsize;
    o1 += 2 * bsize;
    o2 += 2 * bsize;
    o3 += 2 * bsize;
  }
}
```