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

```asm
<__jaciter_kernel_parallel>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %rdx,-0x38(%rbp)
  mov    %rcx,-0x40(%rbp)
  mov    %r8,-0x48(%rbp)
  mov    %r9,-0x50(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x58(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0x60(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D24(%rip),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    0x28(%rbp),%eax
  mov    %eax,%edi
  call   <omp_set_num_threads@plt>
  mov    0x28(%rbp),%r9d
  mov    -0x58(%rbp),%r8
  mov    -0x50(%rbp),%rcx
  mov    -0x38(%rbp),%rdx
  mov    D24(%rip),%rdi
  mov    0x20(%rbp),%esi
  mov    0x20(%rbp),%eax
  pxor   %xmm1,%xmm1
  movq   %rdi,%xmm0
  mov    %eax,%edi
  call   <dgemv_seq@plt>
  movl   $0x0,-0x14(%rbp)
  jmp    L2
L0:
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D21(%rip),%xmm2
  xorpd  %xmm2,%xmm0
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x40(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x14(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x14(%rbp)
L2:
  mov    -0x14(%rbp),%eax
  cmp    0x20(%rbp),%eax
  jl     L0
  mov    -0x50(%rbp),%rdx
  mov    -0x58(%rbp),%rcx
  mov    0x20(%rbp),%eax
  mov    %rcx,%rsi
  mov    %eax,%edi
  call   <dcopy_seq@plt>
  mov    0x28(%rbp),%r9d
  mov    -0x60(%rbp),%r8
  mov    -0x50(%rbp),%rcx
  mov    -0x28(%rbp),%rdx
  mov    D24(%rip),%rdi
  mov    0x20(%rbp),%esi
  mov    0x20(%rbp),%eax
  pxor   %xmm1,%xmm1
  movq   %rdi,%xmm0
  mov    %eax,%edi
  call   <dgemv_seq@plt>
  mov    0x28(%rbp),%ecx
  mov    -0x60(%rbp),%rdx
  mov    -0x30(%rbp),%rsi
  mov    D23(%rip),%rdi
  mov    0x20(%rbp),%eax
  movq   %rdi,%xmm0
  mov    %eax,%edi
  call   <daxpy_seq@plt>
  mov    0x28(%rbp),%ecx
  lea    -0x10(%rbp),%rdx
  mov    -0x60(%rbp),%rsi
  mov    0x20(%rbp),%eax
  mov    %eax,%edi
  call   <dnrm2_seq@plt>
  movsd  -0x10(%rbp),%xmm0
  movq   %xmm0,%rax
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L4
  call   <__stack_chk_fail@plt>
L4:
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: unpack_data
arguments: label: D24
data_type: f64
name: unpack_data
arguments: label: D21
data_type: f32[2]
name: unpack_data
arguments: label: D23
data_type: f32[2]

tool: D24: 1

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

tool: {"D23": [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;
}
```