user: 请从给定的汇编指令推导出C源代码:

```asm
<kern_r4_u1_b>:
  endbr64
  push   %rbp
  push   %rbx
  mov    0x38(%rsp),%rbp
  mov    0x18(%rsp),%r10
  mov    0x20(%rsp),%r11
  mov    0x28(%rsp),%rbx
  test   %rbp,%rbp
  jle    L12
  shl    $0x4,%rbp
  xor    %eax,%eax
  nopl   0x0(%rax)
L17:
  movsd  (%rbx,%rax,1),%xmm1
  movsd  0x8(%rbx,%rax,1),%xmm2
  movsd  0x8(%r9,%rax,4),%xmm4
  movsd  (%r9,%rax,4),%xmm5
  movapd %xmm1,%xmm8
  movapd %xmm2,%xmm0
  movapd %xmm1,%xmm6
  mulsd  %xmm2,%xmm0
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm7
  mulsd  %xmm1,%xmm8
  movapd %xmm4,%xmm9
  mulsd  %xmm2,%xmm9
  mulsd  %xmm1,%xmm4
  subsd  %xmm0,%xmm8
  movapd %xmm1,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  %xmm8,%xmm6
  addsd  %xmm0,%xmm0
  mulsd  %xmm0,%xmm3
  mulsd  %xmm0,%xmm7
  subsd  %xmm3,%xmm6
  movapd %xmm2,%xmm3
  mulsd  %xmm8,%xmm3
  mulsd  %xmm5,%xmm2
  addsd  %xmm3,%xmm7
  movapd %xmm5,%xmm3
  mulsd  %xmm1,%xmm3
  addsd  %xmm4,%xmm2
  movsd  0x8(%r10,%rax,4),%xmm1
  movsd  (%r10,%rax,4),%xmm4
  movapd %xmm1,%xmm5
  mulsd  %xmm8,%xmm1
  mulsd  %xmm0,%xmm5
  subsd  %xmm9,%xmm3
  movapd %xmm4,%xmm9
  mulsd  %xmm0,%xmm4
  movsd  0x8(%r11,%rax,4),%xmm0
  mulsd  %xmm8,%xmm9
  movsd  (%r8,%rax,4),%xmm8
  addsd  %xmm1,%xmm4
  movsd  (%r11,%rax,4),%xmm1
  subsd  %xmm5,%xmm9
  movapd %xmm0,%xmm5
  mulsd  %xmm7,%xmm5
  movapd %xmm1,%xmm10
  mulsd  %xmm6,%xmm10
  mulsd  %xmm6,%xmm0
  movapd %xmm2,%xmm6
  mulsd  %xmm7,%xmm1
  movapd %xmm3,%xmm7
  subsd  %xmm5,%xmm10
  movapd %xmm8,%xmm5
  subsd  %xmm9,%xmm8
  addsd  %xmm9,%xmm5
  addsd  %xmm0,%xmm1
  movsd  0x8(%r8,%rax,4),%xmm0
  addsd  %xmm10,%xmm7
  subsd  %xmm10,%xmm3
  movapd %xmm0,%xmm9
  subsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm6
  subsd  %xmm1,%xmm2
  movapd %xmm5,%xmm1
  addsd  %xmm4,%xmm9
  addsd  %xmm7,%xmm1
  subsd  %xmm7,%xmm5
  movsd  %xmm1,(%rdi,%rax,1)
  movapd %xmm9,%xmm1
  subsd  %xmm6,%xmm9
  addsd  %xmm6,%xmm1
  movsd  %xmm5,(%rdx,%rax,1)
  movsd  %xmm9,0x8(%rdx,%rax,1)
  movsd  %xmm1,0x8(%rdi,%rax,1)
  movapd %xmm8,%xmm1
  subsd  %xmm2,%xmm1
  addsd  %xmm8,%xmm2
  movsd  %xmm1,(%rsi,%rax,1)
  movapd %xmm0,%xmm1
  subsd  %xmm3,%xmm0
  addsd  %xmm3,%xmm1
  movsd  %xmm2,(%rcx,%rax,1)
  movsd  %xmm0,0x8(%rcx,%rax,1)
  movsd  %xmm1,0x8(%rsi,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%rbp
  jne    L17
L12:
  pop    %rbx
  pop    %rbp
  ret
```

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