user: 请将以下的汇编代码转换成C语言代码:

```asm
<kern_r8_fma_o_f>:
  endbr64
  push   %r15
  mov    %rdi,%rax
  mov    %rsi,%r10
  mov    %rdx,%r11
  push   %r14
  mov    %rcx,%r14
  push   %r13
  mov    %r9,%r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xc8,%rsp
  mov    0x168(%rsp),%r15
  mov    %r8,-0x58(%rsp)
  mov    0x110(%rsp),%r12
  mov    0x118(%rsp),%rbp
  mov    0x120(%rsp),%rbx
  mov    0x128(%rsp),%rdx
  mov    0x130(%rsp),%rcx
  mov    0x138(%rsp),%rsi
  mov    0x140(%rsp),%rdi
  mov    0x148(%rsp),%r8
  imul   0x160(%rsp),%r15
  cmpq   $0x0,0x160(%rsp)
  lea    0x1(%r15,%r15,1),%r9
  mov    %r15,-0x48(%rsp)
  mov    0x150(%rsp),%r15
  movsd  (%r15,%r9,8),%xmm15
  jle    L4
  mov    0x160(%rsp),%r15
  xor    %r9d,%r9d
  shl    $0x4,%r15
  mov    %r15,-0x50(%rsp)
  nopl   (%rax)
L46:
  movsd  0x0(%rbp,%r9,1),%xmm10
  movsd  (%rsi,%r9,1),%xmm0
  movsd  0x8(%rsi,%r9,1),%xmm1
  movsd  (%rdx,%r9,1),%xmm9
  movapd %xmm10,%xmm5
  subsd  %xmm0,%xmm10
  movsd  0x8(%rdx,%r9,1),%xmm4
  movsd  (%rcx,%r9,1),%xmm2
  addsd  %xmm0,%xmm5
  movapd %xmm9,%xmm12
  movsd  (%rdi,%r9,1),%xmm8
  movsd  0x8(%rbp,%r9,1),%xmm0
  movapd %xmm4,%xmm7
  mov    -0x58(%rsp),%r15
  movapd %xmm0,%xmm6
  subsd  %xmm1,%xmm0
  addsd  %xmm1,%xmm6
  movsd  (%r8,%r9,1),%xmm1
  subsd  %xmm1,%xmm9
  addsd  %xmm1,%xmm12
  movsd  0x8(%r8,%r9,1),%xmm1
  addsd  %xmm1,%xmm7
  subsd  %xmm1,%xmm4
  movsd  (%r12,%r9,1),%xmm1
  movsd  %xmm12,-0x70(%rsp)
  movapd %xmm1,%xmm13
  subsd  %xmm2,%xmm1
  movsd  %xmm7,-0x78(%rsp)
  movsd  0x8(%r12,%r9,1),%xmm7
  addsd  %xmm2,%xmm13
  movsd  0x8(%rcx,%r9,1),%xmm2
  movapd %xmm7,%xmm3
  subsd  %xmm2,%xmm7
  addsd  %xmm2,%xmm3
  movsd  (%rbx,%r9,1),%xmm2
  movapd %xmm2,%xmm14
  subsd  %xmm8,%xmm2
  addsd  %xmm8,%xmm14
  movsd  0x8(%rbx,%r9,1),%xmm8
  movapd %xmm8,%xmm11
  movapd %xmm14,%xmm12
  movsd  0x8(%rdi,%r9,1),%xmm14
  addsd  %xmm14,%xmm11
  subsd  %xmm14,%xmm8
  movapd %xmm10,%xmm14
  addsd  %xmm0,%xmm14
  subsd  %xmm10,%xmm0
  movapd %xmm4,%xmm10
  xorpd  D8(%rip),%xmm4
  subsd  %xmm9,%xmm10
  subsd  %xmm9,%xmm4
  movapd %xmm13,%xmm9
  subsd  %xmm12,%xmm13
  addsd  %xmm12,%xmm9
  movsd  -0x70(%rsp),%xmm12
  movsd  %xmm10,-0x68(%rsp)
  movapd %xmm7,%xmm10
  subsd  %xmm2,%xmm10
  addsd  %xmm7,%xmm2
  movsd  -0x78(%rsp),%xmm7
  movsd  %xmm4,-0x60(%rsp)
  movapd %xmm3,%xmm4
  subsd  %xmm11,%xmm3
  addsd  %xmm11,%xmm4
  movapd %xmm1,%xmm11
  subsd  %xmm8,%xmm1
  addsd  %xmm8,%xmm11
  movapd %xmm12,%xmm8
  addsd  %xmm6,%xmm7
  addsd  %xmm5,%xmm8
  subsd  %xmm12,%xmm5
  movapd %xmm9,%xmm12
  subsd  -0x78(%rsp),%xmm6
  subsd  %xmm8,%xmm9
  addsd  %xmm8,%xmm12
  movsd  %xmm9,(%r15,%r9,1)
  movsd  -0x68(%rsp),%xmm9
  movsd  %xmm12,(%rax,%r9,1)
  movapd %xmm4,%xmm12
  subsd  %xmm7,%xmm4
  addsd  %xmm7,%xmm12
  movapd %xmm9,%xmm7
  addsd  %xmm14,%xmm7
  subsd  %xmm9,%xmm14
  movsd  %xmm4,0x8(%r15,%r9,1)
  movsd  %xmm12,0x8(%rax,%r9,1)
  movsd  -0x60(%rsp),%xmm12
  mulsd  %xmm15,%xmm7
  mulsd  %xmm15,%xmm14
  movapd %xmm7,%xmm4
  addsd  %xmm11,%xmm4
  subsd  %xmm7,%xmm11
  movsd  %xmm4,(%r10,%r9,1)
  movapd %xmm12,%xmm4
  addsd  %xmm0,%xmm4
  subsd  %xmm12,%xmm0
  mulsd  %xmm15,%xmm4
  mulsd  %xmm15,%xmm0
  movapd %xmm4,%xmm8
  addsd  %xmm10,%xmm8
  subsd  %xmm4,%xmm10
  movapd %xmm13,%xmm4
  addsd  %xmm6,%xmm4
  subsd  %xmm6,%xmm13
  movsd  %xmm8,0x8(%r10,%r9,1)
  mov    0x100(%rsp),%r15
  movsd  %xmm4,(%r11,%r9,1)
  movapd %xmm3,%xmm4
  addsd  %xmm5,%xmm3
  subsd  %xmm5,%xmm4
  movsd  %xmm13,(%r15,%r9,1)
  movsd  %xmm11,0x0(%r13,%r9,1)
  movsd  %xmm3,0x8(%r15,%r9,1)
  movapd %xmm0,%xmm3
  mov    0x108(%rsp),%r15
  addsd  %xmm1,%xmm3
  subsd  %xmm0,%xmm1
  movsd  %xmm10,0x8(%r13,%r9,1)
  movsd  %xmm4,0x8(%r11,%r9,1)
  movsd  %xmm3,(%r14,%r9,1)
  movapd %xmm2,%xmm3
  subsd  %xmm14,%xmm3
  addsd  %xmm2,%xmm14
  movsd  %xmm1,(%r15,%r9,1)
  movsd  %xmm3,0x8(%r14,%r9,1)
  movsd  %xmm14,0x8(%r15,%r9,1)
  add    $0x10,%r9
  cmp    %r9,-0x50(%rsp)
  jne    L46
L4:
  mov    0x160(%rsp),%r9
  shl    $0x7,%r9
  lea    (%rbx,%r9,1),%r15
  add    %r9,%r12
  add    %r9,%rbp
  add    %r9,%rdx
  add    %r9,%rcx
  add    %r9,%rsi
  add    %r9,%rdi
  add    %r9,%r8
  mov    0x160(%rsp),%rbx
  mov    %r9,0x70(%rsp)
  shl    $0x4,%rbx
  add    %rbx,%rax
  add    %rbx,%r10
  add    %rbx,%r11
  add    %rbx,%r14
  mov    %rax,0x78(%rsp)
  mov    -0x58(%rsp),%rax
  add    %rbx,%r13
  mov    %rbx,-0x40(%rsp)
  add    %rbx,%rax
  mov    %r10,0x80(%rsp)
  mov    %rax,0x98(%rsp)
  mov    0x100(%rsp),%rax
  mov    %r11,0x88(%rsp)
  add    %rbx,%rax
  mov    %r14,0x90(%rsp)
  mov    %rax,0xa8(%rsp)
  mov    0x108(%rsp),%rax
  mov    %r13,0xa0(%rsp)
  add    %rbx,%rax
  cmpq   $0x1,0x168(%rsp)
  mov    %rax,0xb0(%rsp)
  jle    L5
  mov    0x160(%rsp),%rax
  mov    %rbx,0x50(%rsp)
  mov    0x160(%rsp),%r14
  movq   $0x0,0x68(%rsp)
  movq   $0x1,0x40(%rsp)
  shl    $0x6,%rax
  mov    %rax,0xc0(%rsp)
  add    0x150(%rsp),%rax
  shl    $0x5,%r14
  mov    %rax,0x60(%rsp)
  mov    -0x48(%rsp),%rax
  add    0x160(%rsp),%rax
  mov    %r14,0xb8(%rsp)
  shl    $0x4,%rax
  mov    %r14,0x58(%rsp)
  mov    %rax,0x48(%rsp)
  nopl   0x0(%rax,%rax,1)
L18:
  mov    0x68(%rsp),%rax
  mov    0x90(%rsp),%rbx
  mov    0xb0(%rsp),%r14
  mov    0x88(%rsp),%r9
  add    %rax,%rbx
  mov    0xa0(%rsp),%r10
  mov    0x78(%rsp),%r13
  mov    %rbx,0x20(%rsp)
  mov    0xa8(%rsp),%rbx
  add    %rax,%r9
  add    %rax,%r14
  mov    0x80(%rsp),%r11
  add    %rax,%r10
  add    %rax,%r13
  mov    %r9,0x30(%rsp)
  add    %rax,%rbx
  mov    0x150(%rsp),%r9
  mov    %rbx,0x28(%rsp)
  mov    0x98(%rsp),%rbx
  add    %rax,%r11
  add    %rax,%rbx
  mov    0x60(%rsp),%rax
  movsd  (%rax),%xmm0
  movsd  0x8(%rax),%xmm6
  mov    0x58(%rsp),%rax
  movsd  %xmm0,-0x38(%rsp)
  movsd  (%r9,%rax,1),%xmm0
  movsd  %xmm0,-0x30(%rsp)
  movsd  0x8(%r9,%rax,1),%xmm0
  mov    0x158(%rsp),%r9
  movsd  %xmm0,-0x28(%rsp)
  movsd  (%r9,%rax,1),%xmm0
  movsd  %xmm0,-0x20(%rsp)
  movsd  0x8(%r9,%rax,1),%xmm0
  mov    0x50(%rsp),%rax
  mov    0x150(%rsp),%r9
  movsd  %xmm0,-0x18(%rsp)
  movsd  (%r9,%rax,1),%xmm0
  movsd  %xmm0,-0x10(%rsp)
  movsd  0x8(%r9,%rax,1),%xmm0
  movsd  %xmm0,-0x78(%rsp)
  mov    0x158(%rsp),%r9
  cmpq   $0x0,0x160(%rsp)
  movsd  (%r9,%rax,1),%xmm0
  movsd  %xmm0,-0x8(%rsp)
  movsd  0x8(%r9,%rax,1),%xmm0
  mov    0x48(%rsp),%rax
  mov    0x150(%rsp),%r9
  movsd  %xmm0,(%rsp)
  movsd  (%r9,%rax,1),%xmm0
  movsd  %xmm0,0x8(%rsp)
  movsd  0x8(%r9,%rax,1),%xmm0
  mov    0x158(%rsp),%r9
  movsd  %xmm0,-0x70(%rsp)
  movsd  (%r9,%rax,1),%xmm0
  movsd  %xmm0,0x10(%rsp)
  movsd  0x8(%r9,%rax,1),%xmm0
  movsd  %xmm0,0x18(%rsp)
  jle    L59
  xor    %eax,%eax
  movsd  %xmm6,0x38(%rsp)
  nopl   0x0(%rax,%rax,1)
L29:
  movsd  -0x30(%rsp),%xmm14
  movsd  (%r15,%rax,1),%xmm2
  movsd  0x8(%r15,%rax,1),%xmm5
  movsd  -0x38(%rsp),%xmm15
  movapd %xmm14,%xmm0
  movsd  (%rsi,%rax,1),%xmm1
  movsd  0x8(%rsi,%rax,1),%xmm3
  mulsd  %xmm2,%xmm0
  movsd  0x38(%rsp),%xmm6
  movsd  0x8(%rbp,%rax,1),%xmm9
  movsd  (%r8,%rax,1),%xmm7
  movsd  0x8(%r8,%rax,1),%xmm4
  movsd  0x0(%rbp,%rax,1),%xmm10
  movsd  (%r12,%rax,1),%xmm8
  movsd  0x8(%rdi,%rax,1),%xmm12
  movapd %xmm10,%xmm13
  addsd  %xmm5,%xmm0
  mulsd  %xmm14,%xmm5
  movsd  (%rcx,%rax,1),%xmm14
  subsd  %xmm2,%xmm5
  movsd  0x8(%rcx,%rax,1),%xmm2
  movsd  %xmm5,-0x68(%rsp)
  movapd %xmm15,%xmm5
  mulsd  %xmm1,%xmm5
  addsd  %xmm3,%xmm5
  mulsd  %xmm15,%xmm3
  mulsd  %xmm6,%xmm5
  subsd  %xmm1,%xmm3
  movapd %xmm9,%xmm1
  mulsd  %xmm6,%xmm3
  addsd  %xmm5,%xmm13
  subsd  %xmm5,%xmm10
  movsd  (%rdx,%rax,1),%xmm5
  movapd %xmm13,%xmm11
  movapd %xmm5,%xmm13
  addsd  %xmm3,%xmm1
  subsd  %xmm3,%xmm9
  movapd %xmm15,%xmm3
  mulsd  %xmm7,%xmm3
  addsd  %xmm4,%xmm3
  mulsd  %xmm15,%xmm4
  mulsd  %xmm6,%xmm3
  subsd  %xmm7,%xmm4
  movsd  0x8(%r12,%rax,1),%xmm7
  subsd  %xmm3,%xmm5
  mulsd  %xmm6,%xmm4
  addsd  %xmm3,%xmm13
  movsd  %xmm5,-0x50(%rsp)
  movsd  0x8(%rdx,%rax,1),%xmm5
  movsd  %xmm13,-0x58(%rsp)
  movapd %xmm5,%xmm3
  movapd %xmm5,%xmm13
  movapd %xmm15,%xmm5
  addsd  %xmm4,%xmm3
  subsd  %xmm4,%xmm13
  movapd %xmm15,%xmm4
  mulsd  %xmm14,%xmm4
  movapd %xmm8,%xmm15
  addsd  %xmm2,%xmm4
  mulsd  %xmm5,%xmm2
  mulsd  %xmm6,%xmm4
  subsd  %xmm14,%xmm2
  movapd %xmm7,%xmm14
  mulsd  %xmm6,%xmm2
  addsd  %xmm4,%xmm15
  subsd  %xmm4,%xmm8
  movsd  (%rdi,%rax,1),%xmm4
  movsd  -0x18(%rsp),%xmm6
  subsd  %xmm2,%xmm7
  addsd  %xmm2,%xmm14
  movsd  -0x20(%rsp),%xmm2
  mulsd  %xmm2,%xmm4
  movapd %xmm4,%xmm5
  movapd %xmm12,%xmm4
  addsd  %xmm12,%xmm5
  movsd  -0x10(%rsp),%xmm12
  mulsd  %xmm2,%xmm4
  subsd  (%rdi,%rax,1),%xmm4
  movapd %xmm12,%xmm2
  mulsd  %xmm11,%xmm2
  mulsd  %xmm6,%xmm4
  mulsd  %xmm6,%xmm5
  movsd  -0x28(%rsp),%xmm6
  addsd  %xmm1,%xmm2
  mulsd  %xmm12,%xmm1
  subsd  %xmm11,%xmm1
  movsd  0x8(%rsp),%xmm11
  movsd  %xmm1,-0x60(%rsp)
  movapd %xmm11,%xmm1
  mulsd  %xmm10,%xmm1
  addsd  %xmm9,%xmm1
  mulsd  %xmm11,%xmm9
  movsd  -0x68(%rsp),%xmm11
  addsd  %xmm4,%xmm11
  mulsd  %xmm6,%xmm11
  subsd  %xmm10,%xmm9
  movapd %xmm5,%xmm10
  addsd  %xmm0,%xmm10
  subsd  %xmm5,%xmm0
  movapd %xmm7,%xmm5
  movsd  %xmm9,-0x48(%rsp)
  mulsd  %xmm6,%xmm0
  movapd %xmm10,%xmm12
  mulsd  %xmm6,%xmm12
  movapd %xmm11,%xmm9
  addsd  %xmm14,%xmm9
  subsd  %xmm11,%xmm14
  movsd  -0x68(%rsp),%xmm11
  subsd  %xmm4,%xmm11
  subsd  %xmm0,%xmm5
  addsd  %xmm7,%xmm0
  movapd %xmm12,%xmm10
  movapd %xmm11,%xmm4
  addsd  %xmm15,%xmm10
  subsd  %xmm12,%xmm15
  mulsd  %xmm6,%xmm4
  movapd %xmm4,%xmm11
  addsd  %xmm8,%xmm11
  subsd  %xmm4,%xmm8
  movsd  %xmm8,-0x68(%rsp)
  movsd  -0x8(%rsp),%xmm4
  movsd  -0x58(%rsp),%xmm7
  movsd  (%rsp),%xmm12
  movsd  0x18(%rsp),%xmm6
  movapd %xmm4,%xmm8
  mov    0x30(%rsp),%r9
  mulsd  %xmm7,%xmm8
  addsd  %xmm3,%xmm8
  mulsd  %xmm4,%xmm3
  movsd  0x10(%rsp),%xmm4
  mulsd  %xmm12,%xmm8
  subsd  %xmm7,%xmm3
  movapd %xmm4,%xmm7
  mulsd  %xmm12,%xmm3
  movsd  -0x50(%rsp),%xmm12
  mulsd  %xmm12,%xmm7
  addsd  %xmm13,%xmm7
  mulsd  %xmm4,%xmm13
  mulsd  %xmm6,%xmm7
  movapd %xmm13,%xmm4
  movapd %xmm8,%xmm13
  subsd  %xmm12,%xmm4
  addsd  %xmm2,%xmm13
  subsd  %xmm8,%xmm2
  mulsd  %xmm6,%xmm4
  movsd  -0x78(%rsp),%xmm6
  mulsd  %xmm6,%xmm13
  movapd %xmm13,%xmm12
  addsd  %xmm10,%xmm12
  subsd  %xmm13,%xmm10
  movsd  -0x48(%rsp),%xmm13
  movsd  %xmm12,0x0(%r13,%rax,1)
  movsd  -0x60(%rsp),%xmm12
  movsd  %xmm10,(%rbx,%rax,1)
  addsd  %xmm3,%xmm12
  mulsd  %xmm6,%xmm12
  movapd %xmm12,%xmm6
  addsd  %xmm9,%xmm6
  subsd  %xmm12,%xmm9
  movsd  -0x70(%rsp),%xmm12
  movsd  %xmm6,0x8(%r13,%rax,1)
  movapd %xmm7,%xmm6
  addsd  %xmm1,%xmm6
  movsd  %xmm9,0x8(%rbx,%rax,1)
  movapd %xmm13,%xmm9
  subsd  %xmm7,%xmm1
  addsd  %xmm4,%xmm9
  movapd %xmm6,%xmm10
  mulsd  %xmm12,%xmm10
  mulsd  %xmm12,%xmm9
  movapd %xmm10,%xmm6
  addsd  %xmm11,%xmm6
  subsd  %xmm10,%xmm11
  movsd  -0x78(%rsp),%xmm10
  mulsd  %xmm10,%xmm2
  movsd  %xmm6,(%r11,%rax,1)
  movapd %xmm9,%xmm6
  addsd  %xmm5,%xmm6
  subsd  %xmm9,%xmm5
  movsd  %xmm11,(%r10,%rax,1)
  movsd  %xmm5,0x8(%r10,%rax,1)
  movsd  -0x60(%rsp),%xmm5
  movsd  %xmm6,0x8(%r11,%rax,1)
  subsd  %xmm3,%xmm5
  movapd %xmm5,%xmm3
  movsd  -0x70(%rsp),%xmm5
  mulsd  %xmm10,%xmm3
  mulsd  %xmm5,%xmm1
  movapd %xmm3,%xmm6
  addsd  %xmm15,%xmm6
  subsd  %xmm3,%xmm15
  movsd  %xmm6,(%r9,%rax,1)
  movapd %xmm14,%xmm6
  subsd  %xmm2,%xmm6
  addsd  %xmm14,%xmm2
  movsd  %xmm6,0x8(%r9,%rax,1)
  mov    0x28(%rsp),%r9
  movsd  %xmm2,0x8(%r9,%rax,1)
  movapd %xmm13,%xmm2
  subsd  %xmm4,%xmm2
  movsd  -0x68(%rsp),%xmm4
  movsd  %xmm15,(%r9,%rax,1)
  mov    0x20(%rsp),%r9
  movapd %xmm4,%xmm3
  mulsd  %xmm5,%xmm2
  addsd  %xmm2,%xmm3
  subsd  %xmm2,%xmm4
  movsd  %xmm3,(%r9,%rax,1)
  movapd %xmm0,%xmm3
  subsd  %xmm1,%xmm3
  addsd  %xmm0,%xmm1
  movsd  %xmm4,(%r14,%rax,1)
  movsd  %xmm3,0x8(%r9,%rax,1)
  movsd  %xmm1,0x8(%r14,%rax,1)
  add    $0x10,%rax
  cmp    %rax,-0x40(%rsp)
  jne    L29
L59:
  mov    0x70(%rsp),%rax
  mov    -0x40(%rsp),%rbx
  addq   $0x1,0x40(%rsp)
  mov    0xc0(%rsp),%r14
  add    %rax,%r12
  add    %rax,%rbp
  add    %rax,%r15
  add    %rax,%rdx
  add    %rax,%rcx
  add    %rax,%rsi
  add    %rax,%rdi
  add    %rax,%r8
  mov    0xb8(%rsp),%r10
  add    %rbx,0x68(%rsp)
  add    %r14,0x60(%rsp)
  mov    0x40(%rsp),%rax
  add    %r10,0x58(%rsp)
  add    %rbx,0x50(%rsp)
  add    %rbx,0x48(%rsp)
  cmp    %rax,0x168(%rsp)
  jne    L18
L5:
  add    $0xc8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D8
datatype: byte[16]

tool: {"D8": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
void kern_r8_fma_o_f(double *restrict o0, double *restrict o1, double *restrict o2, double *restrict o3, double *restrict o4, double *restrict o5, double *restrict o6, double *restrict o7, double *i0, double *i1, double *i2, double *i3, double *i4, double *i5, double *i6, double *i7, double *w, double *w2, long bsize, long blocks) {
  long i, j;
  double w5;
  w5 = w[2 * blocks * bsize + 1];
  for (j = 0; j < bsize; j++) {
    double sr5, si5;
    double sr7, si7;
    double qr0, qi0;
    double qr1, qi1;
    double qr2, qi2;
    double qr3, qi3;
    double qr4, qi4;
    double qr5, qi5;
    double qr6, qi6;
    double qr7, qi7;
    double tr0, ti0;
    double tr1, ti1;
    double tr2, ti2;
    double tr3, ti3;
    double tr4, ti4;
    double tr5, ti5;
    double tr6, ti6;
    double tr7, ti7;
    qr1 = i1[2 * j] + i5[2 * j];
    sr5 = i1[2 * j] - i5[2 * j];
    qi1 = i1[2 * j + 1] + i5[2 * j + 1];
    si5 = i1[2 * j + 1] - i5[2 * j + 1];
    qr3 = i3[2 * j] + i7[2 * j];
    sr7 = i3[2 * j] - i7[2 * j];
    qi3 = i3[2 * j + 1] + i7[2 * j + 1];
    si7 = i3[2 * j + 1] - i7[2 * j + 1];
    qr0 = i0[2 * j] + i4[2 * j];
    qr4 = i0[2 * j] - i4[2 * j];
    qi0 = i0[2 * j + 1] + i4[2 * j + 1];
    qi4 = i0[2 * j + 1] - i4[2 * j + 1];
    qr2 = i2[2 * j] + i6[2 * j];
    qr6 = i2[2 * j] - i6[2 * j];
    qi2 = i2[2 * j + 1] + i6[2 * j + 1];
    qi6 = i2[2 * j + 1] - i6[2 * j + 1];
    qr5 = sr5 + si5;
    qi5 = si5 - sr5;
    qr7 = -sr7 + si7;
    qi7 = -si7 - sr7;
    tr0 = qr0 + qr2;
    ti0 = qi0 + qi2;
    tr2 = qr0 - qr2;
    ti2 = qi0 - qi2;
    tr4 = qr4 + qi6;
    ti4 = qi4 - qr6;
    tr6 = qr4 - qi6;
    ti6 = qi4 + qr6;
    tr1 = qr1 + qr3;
    ti1 = qi1 + qi3;
    tr3 = qr1 - qr3;
    ti3 = qi1 - qi3;
    tr5 = qr5 + qr7;
    ti5 = qi5 + qi7;
    tr7 = qr5 - qr7;
    ti7 = qi5 - qi7;
    o0[2 * j] = tr0 + tr1;
    o0[2 * j + 1] = ti0 + ti1;
    o4[2 * j] = tr0 - tr1;
    o4[2 * j + 1] = ti0 - ti1;
    o1[2 * j] = tr4 + tr5 * w5;
    o1[2 * j + 1] = ti4 + ti5 * w5;
    o5[2 * j] = tr4 - tr5 * w5;
    o5[2 * j + 1] = ti4 - ti5 * w5;
    o2[2 * j] = tr2 + ti3;
    o2[2 * j + 1] = ti2 - tr3;
    o6[2 * j] = tr2 - ti3;
    o6[2 * j + 1] = ti2 + tr3;
    o3[2 * j] = tr6 + ti7 * w5;
    o3[2 * j + 1] = ti6 - tr7 * w5;
    o7[2 * j] = tr6 - ti7 * w5;
    o7[2 * j + 1] = ti6 + tr7 * w5;
  }
  i0 += 16 * bsize;
  i1 += 16 * bsize;
  i2 += 16 * bsize;
  i3 += 16 * bsize;
  i4 += 16 * bsize;
  i5 += 16 * bsize;
  i6 += 16 * bsize;
  i7 += 16 * bsize;
  o0 += 2 * bsize;
  o1 += 2 * bsize;
  o2 += 2 * bsize;
  o3 += 2 * bsize;
  o4 += 2 * bsize;
  o5 += 2 * bsize;
  o6 += 2 * bsize;
  o7 += 2 * bsize;
  for (i = 1; i < blocks; i++) {
    double wr1, wi1;
    double wr2, wi2;
    double wr3, wi3;
    double wr4, wi4;
    double wr5, wi5;
    double wr6, wi6;
    double wr7, wi7;
    wr4 = w[8 * i * bsize];
    wi4 = w[8 * i * bsize + 1];
    wr2 = w[4 * i * bsize];
    wi2 = w[4 * i * bsize + 1];
    wr6 = w2[4 * i * bsize];
    wi6 = w2[4 * i * bsize + 1];
    wr1 = w[2 * i * bsize];
    wi1 = w[2 * i * bsize + 1];
    wr3 = w2[2 * i * bsize];
    wi3 = w2[2 * i * bsize + 1];
    wr5 = w[2 * i * bsize + 2 * blocks * bsize];
    wi5 = w[2 * i * bsize + 2 * blocks * bsize + 1];
    wr7 = w2[2 * i * bsize + 2 * blocks * bsize];
    wi7 = w2[2 * i * bsize + 2 * blocks * bsize + 1];
    for (j = 0; j < bsize; j++) {
      double pr2, pi2;
      double pr4, pi4;
      double pr5, pi5;
      double pr6, pi6;
      double pr7, pi7;
      double sr1, si1;
      double sr3, si3;
      double sr5, si5;
      double sr7, si7;
      double qr0, qi0;
      double qr1, qi1;
      double qr2, qi2;
      double qr3, qi3;
      double qr4, qi4;
      double qr5, qi5;
      double qr6, qi6;
      double qr7, qi7;
      double tr0, ti0;
      double tr1, ti1;
      double tr2, ti2;
      double tr3, ti3;
      double tr4, ti4;
      double tr5, ti5;
      double tr6, ti6;
      double tr7, ti7;
      pr5 = i5[2 * j] * wr4 + i5[2 * j + 1];
      pi5 = i5[2 * j + 1] * wr4 - i5[2 * j];
      pr7 = i7[2 * j] * wr4 + i7[2 * j + 1];
      pi7 = i7[2 * j + 1] * wr4 - i7[2 * j];
      pr2 = i2[2 * j] * wr2 + i2[2 * j + 1];
      pi2 = i2[2 * j + 1] * wr2 - i2[2 * j];
      pr4 = i4[2 * j] * wr4 + i4[2 * j + 1];
      pi4 = i4[2 * j + 1] * wr4 - i4[2 * j];
      pr6 = i6[2 * j] * wr6 + i6[2 * j + 1];
      pi6 = i6[2 * j + 1] * wr6 - i6[2 * j];
      sr1 = i1[2 * j] + pr5 * wi4;
      sr5 = i1[2 * j] - pr5 * wi4;
      si1 = i1[2 * j + 1] + pi5 * wi4;
      si5 = i1[2 * j + 1] - pi5 * wi4;
      sr3 = i3[2 * j] + pr7 * wi4;
      sr7 = i3[2 * j] - pr7 * wi4;
      si3 = i3[2 * j + 1] + pi7 * wi4;
      si7 = i3[2 * j + 1] - pi7 * wi4;
      qr0 = i0[2 * j] + pr4 * wi4;
      qr4 = i0[2 * j] - pr4 * wi4;
      qi0 = i0[2 * j + 1] + pi4 * wi4;
      qi4 = i0[2 * j + 1] - pi4 * wi4;
      qr2 = pr2 + pr6 * wi6;
      qr6 = pr2 - pr6 * wi6;
      qi2 = pi2 + pi6 * wi6;
      qi6 = pi2 - pi6 * wi6;
      qr1 = sr1 * wr1 + si1;
      qi1 = si1 * wr1 - sr1;
      qr3 = sr3 * wr3 + si3;
      qi3 = si3 * wr3 - sr3;
      qr5 = sr5 * wr5 + si5;
      qi5 = si5 * wr5 - sr5;
      qr7 = sr7 * wr7 + si7;
      qi7 = si7 * wr7 - sr7;
      tr0 = qr0 + qr2 * wi2;
      ti0 = qi0 + qi2 * wi2;
      tr2 = qr0 - qr2 * wi2;
      ti2 = qi0 - qi2 * wi2;
      tr4 = qr4 + qi6 * wi2;
      ti4 = qi4 - qr6 * wi2;
      tr6 = qr4 - qi6 * wi2;
      ti6 = qi4 + qr6 * wi2;
      tr1 = qr1 + qr3 * wi3;
      ti1 = qi1 + qi3 * wi3;
      tr3 = qr1 - qr3 * wi3;
      ti3 = qi1 - qi3 * wi3;
      tr5 = qr5 + qr7 * wi7;
      ti5 = qi5 + qi7 * wi7;
      tr7 = qr5 - qr7 * wi7;
      ti7 = qi5 - qi7 * wi7;
      o0[2 * j] = tr0 + tr1 * wi1;
      o0[2 * j + 1] = ti0 + ti1 * wi1;
      o4[2 * j] = tr0 - tr1 * wi1;
      o4[2 * j + 1] = ti0 - ti1 * wi1;
      o1[2 * j] = tr4 + tr5 * wi5;
      o1[2 * j + 1] = ti4 + ti5 * wi5;
      o5[2 * j] = tr4 - tr5 * wi5;
      o5[2 * j + 1] = ti4 - ti5 * wi5;
      o2[2 * j] = tr2 + ti3 * wi1;
      o2[2 * j + 1] = ti2 - tr3 * wi1;
      o6[2 * j] = tr2 - ti3 * wi1;
      o6[2 * j + 1] = ti2 + tr3 * wi1;
      o3[2 * j] = tr6 + ti7 * wi5;
      o3[2 * j + 1] = ti6 - tr7 * wi5;
      o7[2 * j] = tr6 - ti7 * wi5;
      o7[2 * j + 1] = ti6 + tr7 * wi5;
    }
    i0 += 16 * bsize;
    i1 += 16 * bsize;
    i2 += 16 * bsize;
    i3 += 16 * bsize;
    i4 += 16 * bsize;
    i5 += 16 * bsize;
    i6 += 16 * bsize;
    i7 += 16 * bsize;
    o0 += 2 * bsize;
    o1 += 2 * bsize;
    o2 += 2 * bsize;
    o3 += 2 * bsize;
    o4 += 2 * bsize;
    o5 += 2 * bsize;
    o6 += 2 * bsize;
    o7 += 2 * bsize;
  }
}
```