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

```asm
<kern_r4_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    L3
L35:
  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)
L3:
  mov    -0xe8(%rbp),%rax
  cmp    0x28(%rbp),%rax
  jl     L35
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x118(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x120(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,0x10(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,0x18(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0xf8(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x100(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x108(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x110(%rbp)
  movq   $0x1,-0xf0(%rbp)
  jmp    L58
L43:
  mov    -0xf0(%rbp),%rax
  imul   0x28(%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   0x28(%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   0x28(%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   0x28(%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   0x28(%rbp),%rax
  mov    %rax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  mov    -0xf0(%rbp),%rax
  imul   0x28(%rbp),%rax
  mov    %rax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  add    %rax,%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movq   $0x0,-0xe8(%rbp)
  jmp    L10
L30:
  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
  mulsd  -0xd8(%rbp),%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),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xd8(%rbp),%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
  mulsd  -0xc8(%rbp),%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),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xc8(%rbp),%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
  mulsd  -0xb8(%rbp),%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),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0xb8(%rbp),%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),%xmm0
  movsd  -0xa0(%rbp),%xmm1
  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
  subsd  -0xa0(%rbp),%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),%xmm0
  movsd  -0x98(%rbp),%xmm1
  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
  subsd  -0x98(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0xb0(%rbp),%xmm0
  addsd  -0x90(%rbp),%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0xb0(%rbp),%xmm0
  subsd  -0x90(%rbp),%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  addsd  -0x88(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  subsd  -0x88(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x80(%rbp),%xmm0
  addsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0xf8(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x70(%rbp),%xmm0
  addsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x80(%rbp),%xmm0
  subsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x108(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x70(%rbp),%xmm0
  subsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x78(%rbp),%xmm0
  subsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x100(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x68(%rbp),%xmm0
  addsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  mov    %rax,%rdx
  mov    -0x110(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x78(%rbp),%xmm0
  addsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xe8(%rbp),%rax
  shl    $0x4,%rax
  lea    0x8(%rax),%rdx
  mov    -0x110(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x68(%rbp),%xmm0
  subsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0xe8(%rbp)
L10:
  mov    -0xe8(%rbp),%rax
  cmp    0x28(%rbp),%rax
  jl     L30
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x118(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,-0x120(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,0x10(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x6,%rax
  add    %rax,0x18(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0xf8(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x100(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x108(%rbp)
  mov    0x28(%rbp),%rax
  shl    $0x4,%rax
  add    %rax,-0x110(%rbp)
  addq   $0x1,-0xf0(%rbp)
L58:
  mov    -0xf0(%rbp),%rax
  cmp    0x30(%rbp),%rax
  jl     L43
  nop
  nop
  leave
  ret
```

assistant: ```c
void kern_r4_o_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, 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, wi3;
    wr1 = w[2 * i * bsize];
    wi1 = w[2 * i * bsize + 1];
    wr2 = w[4 * i * bsize];
    wi2 = w[4 * i * bsize + 1];
    wr3 = w[6 * i * bsize];
    wi3 = w[6 * 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] * wi1;
      pi1 = i1[2 * j] * wi1 + i1[2 * j + 1] * wr1;
      pr2 = i2[2 * j] * wr2 - i2[2 * j + 1] * wi2;
      pi2 = i2[2 * j] * wi2 + i2[2 * j + 1] * wr2;
      pr3 = i3[2 * j] * wr3 - i3[2 * j + 1] * wi3;
      pi3 = i3[2 * j] * wi3 + i3[2 * j + 1] * wr3;
      qr0 = i0[2 * j] + pr2;
      qr2 = i0[2 * j] - pr2;
      qi0 = i0[2 * j + 1] + pi2;
      qi2 = i0[2 * j + 1] - pi2;
      qr1 = pr1 + pr3;
      qr3 = pr1 - pr3;
      qi1 = pi1 + pi3;
      qi3 = pi1 - pi3;
      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;
  }
}
```