user: Please show me the C implementation that matches the assembly code listed below:

```asm
<__jaciter_kernel_parallel>:
  endbr64
  push   %r15
  mov    %r9,%r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rcx,%rbx
  sub    $0x38,%rsp
  mov    %rdi,(%rsp)
  mov    0x88(%rsp),%r12d
  mov    %rsi,0x8(%rsp)
  mov    0x80(%rsp),%r14d
  mov    %rdx,0x10(%rsp)
  mov    %r12d,%edi
  mov    0x70(%rsp),%rbp
  mov    %r8,0x18(%rsp)
  mov    0x78(%rsp),%r13
  mov    %fs:0x28,%rax
  mov    %rax,0x28(%rsp)
  xor    %eax,%eax
  mov    D45(%rip),%rax
  mov    %rax,0x20(%rsp)
  call   <omp_set_num_threads@plt>
  mov    0x10(%rsp),%rdx
  mov    %r12d,%r9d
  mov    %rbp,%r8
  mov    D45(%rip),%rax
  mov    %r15,%rcx
  mov    %r14d,%esi
  mov    %r14d,%edi
  pxor   %xmm1,%xmm1
  movq   %rax,%xmm0
  call   <dgemv_seq@plt>
  test   %r14d,%r14d
  jle    L47
  lea    0xf(%rbp),%rax
  mov    0x18(%rsp),%r10
  sub    %rbx,%rax
  cmp    $0x1e,%rax
  seta   %dl
  cmp    $0x1,%r14d
  setne  %al
  test   %al,%dl
  je     L30
  lea    0x8(%r10),%rax
  cmp    %rax,%rbp
  je     L30
  mov    %r14d,%edx
  movapd D54(%rip),%xmm1
  xor    %eax,%eax
  shr    %edx
  shl    $0x4,%rdx
  nopl   (%rax)
L73:
  movupd (%r10,%rax,1),%xmm0
  movupd 0x0(%rbp,%rax,1),%xmm2
  xorpd  %xmm1,%xmm0
  mulpd  %xmm2,%xmm0
  movups %xmm0,0x0(%rbp,%rax,1)
  movupd (%rbx,%rax,1),%xmm3
  addpd  %xmm3,%xmm0
  movups %xmm0,0x0(%rbp,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%rdx
  jne    L73
  mov    %r14d,%eax
  and    $0xfffffffe,%eax
  test   $0x1,%r14b
  je     L47
  lea    0x0(%rbp,%rax,8),%rdx
  movsd  (%r10,%rax,8),%xmm0
  xorpd  D0(%rip),%xmm0
  mulsd  (%rdx),%xmm0
  movsd  %xmm0,(%rdx)
  addsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx)
L47:
  mov    %r15,%rdx
  mov    %rbp,%rsi
  mov    %r14d,%edi
  call   <dcopy_seq@plt>
  mov    (%rsp),%rdx
  mov    %r12d,%r9d
  mov    %r13,%r8
  mov    D45(%rip),%rax
  mov    %r15,%rcx
  mov    %r14d,%esi
  mov    %r14d,%edi
  pxor   %xmm1,%xmm1
  movq   %rax,%xmm0
  call   <dgemv_seq@plt>
  mov    0x8(%rsp),%rsi
  mov    %r12d,%ecx
  mov    %r13,%rdx
  movsd  D20(%rip),%xmm0
  mov    %r14d,%edi
  call   <daxpy_seq@plt>
  lea    0x20(%rsp),%rdx
  mov    %r12d,%ecx
  mov    %r13,%rsi
  mov    %r14d,%edi
  call   <dnrm2_seq@plt>
  movsd  0x20(%rsp),%xmm0
  mov    0x28(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L37
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L30:
  movq   D0(%rip),%xmm1
  movslq %r14d,%rdx
  xor    %eax,%eax
  nopl   (%rax)
L54:
  movsd  (%r10,%rax,8),%xmm0
  xorpd  %xmm1,%xmm0
  mulsd  0x0(%rbp,%rax,8),%xmm0
  movsd  %xmm0,0x0(%rbp,%rax,8)
  addsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,0x0(%rbp,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L54
  jmp    L47
L37:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D45
type_infered: f64
name: read_data
arguments: data_label: D54
type_infered: f32[4]
name: read_data
arguments: data_label: D0
type_infered: f32[4]
name: read_data
arguments: data_label: D20
type_infered: f32[2]

tool: D45: 1

tool: {"D54": [0.0,-0.0,0.0,-0.0]}

tool: {"D0": [0.0,-0.0,0.0,0.0]}

tool: {"D20": [0.0,-1.875]}

assistant: ```c
double __jaciter_kernel_parallel(double *A, double *b, double *R, double *C, double *Dinv, double *xk, double *xkp1, double *xconv, int n, int p) {
  double e = 1.0;
  int i;
  omp_set_num_threads(p);
  dgemv_seq(n, n, 1.0, R, xk, 0.0, xkp1, p);
#pragma omp parallel for schedule(static) private(i) shared(xkp1, Dinv, C, n)
  for (i = 0; i < n; i++) {
    xkp1[i] *= -Dinv[i];
    xkp1[i] += C[i];
  }
  dcopy_seq(n, xkp1, xk);
  dgemv_seq(n, n, 1.0, A, xk, 0.0, xconv, p);
  daxpy_seq(n, -1.0, b, xconv, p);
  dnrm2_seq(n, xconv, &e, p);
  return e;
}
```