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

```asm
<kern_r4_u1_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,-0xa8(%rbp)
  jmp    L8
L9:
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  -0x98(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  subsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  mulsd  -0x90(%rbp),%xmm0
  movsd  -0x98(%rbp),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x88(%rbp),%xmm1
  movsd  -0x98(%rbp),%xmm0
  mulsd  -0x90(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  mov    %rax,%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0xa0(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  lea    0x8(%rax),%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x98(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  mov    %rax,%rdx
  mov    -0xe0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x98(%rbp),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%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    $0x6,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x90(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  lea    0x8(%rax),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  mov    %rax,%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x88(%rbp),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%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    $0x6,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  lea    0x8(%rax),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x78(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  mov    %rax,%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x78(%rbp),%xmm1
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%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    $0x6,%rax
  mov    %rax,%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x60(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  mov    %rax,%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  subsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  lea    0x8(%rax),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x58(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x6,%rax
  lea    0x8(%rax),%rdx
  mov    -0xd8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  subsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x70(%rbp),%xmm0
  addsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x70(%rbp),%xmm0
  subsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x68(%rbp),%xmm0
  addsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x68(%rbp),%xmm0
  subsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x40(%rbp),%xmm0
  addsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xb8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x30(%rbp),%xmm0
  addsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x40(%rbp),%xmm0
  subsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xc8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x30(%rbp),%xmm0
  subsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x38(%rbp),%xmm0
  subsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xc0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x38(%rbp),%xmm0
  addsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xa8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xd0(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x28(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0xa8(%rbp)
L8:
  mov    -0xa8(%rbp),%rax
  cmp    0x30(%rbp),%rax
  jl     L9
  nop
  nop
  leave
  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;
  }
}
```