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

```asm
<kern_r4_fma_o_b>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xa8,%rsp
  mov    %rdi,-0xf8(%rbp)
  mov    %rsi,-0x100(%rbp)
  mov    %rdx,-0x108(%rbp)
  mov    %rcx,-0x110(%rbp)
  mov    %r8,-0x118(%rbp)
  mov    %r9,-0x120(%rbp)
  movq   $0x0,-0xe8(%rbp)
  jmp    L21
L39:
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x40(%rbp),%xmm0
  addsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x30(%rbp),%xmm0
  addsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x40(%rbp),%xmm0
  subsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x30(%rbp),%xmm0
  subsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x38(%rbp),%xmm0
  subsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x110(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x38(%rbp),%xmm0
  addsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x110(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0xe8(%rbp)
L21:
  mov    -0xe8(%rbp),%rax
  cmp    0x30(%rbp),%rax
  jl     L39
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x118(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x120(%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,-0xf8(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x100(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x108(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x110(%rbp)
  movq   $0x1,-0xf0(%rbp)
  jmp    L5
L36:
  mov    -0xf0(%rbp),%rax
  imul   0x30(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  mov    -0xf0(%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,-0xd8(%rbp)
  mov    -0xf0(%rbp),%rax
  imul   0x30(%rbp),%rax
  shl    $0x5,%rax
  mov    %rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  mov    -0xf0(%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,-0xc8(%rbp)
  mov    -0xf0(%rbp),%rax
  imul   0x30(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  mov    -0xf0(%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,-0xb8(%rbp)
  movq   $0x0,-0xe8(%rbp)
  jmp    L56
L22:
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xe0(%rbp),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x120(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xe0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xd0(%rbp),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xd0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xc0(%rbp),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xc0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0xa0(%rbp),%xmm0
  mulsd  -0xc8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0xa0(%rbp),%xmm1
  mulsd  -0xc8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x98(%rbp),%xmm0
  mulsd  -0xc8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x118(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x98(%rbp),%xmm1
  mulsd  -0xc8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x90(%rbp),%xmm0
  mulsd  -0xb8(%rbp),%xmm0
  movsd  -0xb0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x90(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xb8(%rbp),%xmm1
  movsd  -0xb0(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0x88(%rbp),%xmm0
  mulsd  -0xb8(%rbp),%xmm0
  movsd  -0xa8(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x88(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xb8(%rbp),%xmm1
  movsd  -0xa8(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x60(%rbp),%xmm0
  mulsd  -0xd8(%rbp),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0xd8(%rbp),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x60(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xd8(%rbp),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x80(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x50(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xd8(%rbp),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x70(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xd8(%rbp),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x78(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x58(%rbp),%xmm0
  mulsd  -0xd8(%rbp),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x68(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x48(%rbp),%xmm0
  mulsd  -0xd8(%rbp),%xmm0
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x110(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x78(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xd8(%rbp),%xmm1
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x110(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x68(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0xe8(%rbp)
L56:
  mov    -0xe8(%rbp),%rax
  cmp    0x30(%rbp),%rax
  jl     L22
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x118(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x120(%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,-0xf8(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x100(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x108(%rbp)
  mov    0x30(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x110(%rbp)
  addq   $0x1,-0xf0(%rbp)
L5:
  mov    -0xf0(%rbp),%rax
  cmp    0x38(%rbp),%rax
  jl     L36
  nop
  nop
  leave
  ret
```

assistant: ```c
void kern_r4_fma_o_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 (j = 0; j < bsize; j++) {
    double qr0, qi0;
    double qr1, qi1;
    double qr2, qi2;
    double qr3, qi3;
    qr0 = i0[2 * j] + i2[2 * j];
    qr2 = i0[2 * j] - i2[2 * j];
    qi0 = i0[2 * j + 1] + i2[2 * j + 1];
    qi2 = i0[2 * j + 1] - i2[2 * j + 1];
    qr1 = i1[2 * j] + i3[2 * j];
    qr3 = i1[2 * j] - i3[2 * j];
    qi1 = i1[2 * j + 1] + i3[2 * j + 1];
    qi3 = i1[2 * j + 1] - i3[2 * j + 1];
    o0[2 * j] = qr0 + qr1;
    o0[2 * j + 1] = qi0 + qi1;
    o2[2 * j] = qr0 - qr1;
    o2[2 * j + 1] = qi0 - qi1;
    o1[2 * j] = qr2 - qi3;
    o1[2 * j + 1] = qi2 + qr3;
    o3[2 * j] = qr2 + qi3;
    o3[2 * j + 1] = qi2 - qr3;
  }
  i0 += 8 * bsize;
  i1 += 8 * bsize;
  i2 += 8 * bsize;
  i3 += 8 * bsize;
  o0 += 2 * bsize;
  o1 += 2 * bsize;
  o2 += 2 * bsize;
  o3 += 2 * bsize;
  for (i = 1; 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;
  }
}
```