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

```asm
<cpu_cal_align_score>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %ecx,%r13d
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    0x90(%rsp),%eax
  mov    0xb0(%rsp),%ebx
  mov    %rdx,0x18(%rsp)
  mov    %rdi,0x40(%rsp)
  mov    0xa8(%rsp),%r8d
  cltd
  mov    %rsi,0x48(%rsp)
  idiv   %ebx
  mov    %r9d,0x24(%rsp)
  mov    %edx,0x20(%rsp)
  mov    %eax,0xc(%rsp)
  test   %edx,%edx
  je     L21
  add    $0x1,%eax
  mov    %eax,0xc(%rsp)
L4:
  mov    0xa0(%rsp),%eax
  cmp    %eax,0x98(%rsp)
  jge    L31
  lea    0x1(%r13),%eax
  movl   $0x0,0x34(%rsp)
  mov    %eax,0x3c(%rsp)
  imul   0x98(%rsp),%eax
  mov    %eax,0x38(%rsp)
  mov    %r8d,%eax
  imul   %ebx,%eax
  lea    (%rax,%rax,4),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,0x28(%rsp)
  mov    0xc(%rsp),%eax
  sub    $0x1,%eax
  mov    %eax,0x30(%rsp)
  nopl   0x0(%rax,%rax,1)
L44:
  mov    0xc(%rsp),%eax
  test   %eax,%eax
  jle    L9
  mov    0x34(%rsp),%r14d
  movslq 0x38(%rsp),%rbp
  xor    %r15d,%r15d
  add    0x40(%rsp),%rbp
  mov    0x48(%rsp),%r12
  mov    %rbp,0x10(%rsp)
  mov    %ebx,%ebp
  mov    %r14d,%ebx
  mov    %r8d,%r14d
  nopl   0x0(%rax)
L6:
  mov    0x20(%rsp),%r9d
  cmp    0x30(%rsp),%r15d
  mov    %r12,%rsi
  mov    %r14d,%r8d
  mov    %r13d,%edx
  cmovne %ebp,%r9d
  sub    $0x8,%rsp
  add    $0x1,%r15d
  push   0xc0(%rsp)
  push   0x28(%rsp)
  push   %rbx
  add    %ebp,%ebx
  mov    0x44(%rsp),%ecx
  mov    0x30(%rsp),%rdi
  call   <align_cpu@plt>
  add    0x48(%rsp),%r12
  add    $0x20,%rsp
  cmp    %r15d,0xc(%rsp)
  jne    L6
  mov    %r14d,%r8d
  mov    %ebp,%ebx
L9:
  addl   $0x1,0x98(%rsp)
  mov    0x90(%rsp),%ecx
  mov    0x3c(%rsp),%esi
  add    %ecx,0x34(%rsp)
  add    %esi,0x38(%rsp)
  mov    0x98(%rsp),%eax
  cmp    %eax,0xa0(%rsp)
  jne    L44
L31:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L21:
  mov    %ebx,0x20(%rsp)
  jmp    L4
```

assistant: ```c
void cpu_cal_align_score(char *content, uint64_t *preprocess_reads, int8_t *align_results, int ref_len, int ref_count, int read_len, int read_count, int ref_start, int ref_end, int word_num, int chunk_read_num, uint64_t *dvdh_bit_mem) {
  int i = 0, j = 0;
  char *ref_itr;
  uint64_t *read_itr;
  int64_t align_index_start;
  int loop_num = read_count / 1;
  int loop_num_temp = loop_num;
  int index = 0;
  int rest_read_num;
  int tmp_read_num;
  if (loop_num % chunk_read_num == 0) {
    loop_num = loop_num / chunk_read_num;
    rest_read_num = chunk_read_num;
  } else {
    rest_read_num = loop_num % chunk_read_num;
    loop_num = loop_num / chunk_read_num + 1;
  }
#pragma omp parallel for default(none) num_threads(cpu_threads) private(i, j, read_itr, ref_itr, align_index_start, tmp_read_num) shared(ref_len, read_len, ref_start, ref_end, content, align_results, preprocess_reads, loop_num, loop_num_temp, word_num, dvdh_bit_mem, chunk_read_num, rest_read_num) schedule(guided) collapse(2)
  for (i = ref_start; i < ref_end; i++) {
    for (j = 0; j < loop_num; j++) {
      if (j == loop_num - 1) {
        tmp_read_num = rest_read_num;
      } else {
        tmp_read_num = chunk_read_num;
      }
      ref_itr = &content[i * (ref_len + 1)];
      align_index_start = (i - ref_start) * loop_num_temp;
      read_itr = &preprocess_reads[5 * 1 * j * word_num * chunk_read_num];
      align_cpu(ref_itr, read_itr, ref_len, read_len, word_num, tmp_read_num, align_index_start + j * chunk_read_num, align_results, dvdh_bit_mem);
    }
  }
}
```