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

```asm
<matched_filter_precise>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %rcx,0x38(%rsp)
  mov    0xc0(%rsp),%ecx
  mov    %rsi,0x30(%rsp)
  mov    %rdx,0x60(%rsp)
  mov    %r8,0x40(%rsp)
  test   %ecx,%ecx
  jle    L34
  mov    0xc8(%rsp),%ecx
  pxor   %xmm2,%xmm2
  mov    %rdi,%rbx
  mov    0xb8(%rsp),%eax
  sub    0xb0(%rsp),%eax
  cvtsi2ss %r9d,%xmm2
  mov    %r9d,%ebp
  movq   $0x0,0x20(%rsp)
  imul   0xd0(%rsp),%ecx
  mov    %eax,0x68(%rsp)
  mov    0xb0(%rsp),%eax
  movl   $0x0,(%rsp)
  movl   $0x0,0x48(%rsp)
  imul   %ecx,%eax
  mov    %ecx,0x4c(%rsp)
  movss  %xmm2,0x6c(%rsp)
  cltq
  shl    $0x2,%rax
  mov    %rax,0x50(%rsp)
  movslq %ecx,%rax
  shl    $0x2,%rax
  mov    %rax,0x58(%rsp)
  lea    -0x1(%rcx),%eax
  lea    0x4(%rdx,%rax,4),%rax
  mov    %rax,0x28(%rsp)
  movslq %r9d,%rax
  shl    $0x2,%rax
  mov    %rax,0x18(%rsp)
  nopl   0x0(%rax)
L8:
  mov    0x4c(%rsp),%eax
  test   %eax,%eax
  jle    L29
  mov    0x60(%rsp),%r12
  mov    0x28(%rsp),%rcx
  xor    %r13d,%r13d
  xor    %edi,%edi
  add    0x20(%rsp),%r12
  mov    %r12,%rdx
  nopl   0x0(%rax,%rax,1)
L50:
  mov    (%rdx),%eax
  cmp    %eax,%edi
  cmovg  %eax,%edi
  cmp    %eax,%r13d
  cmovl  %eax,%r13d
  add    $0x4,%rdx
  cmp    %rcx,%rdx
  jne    L50
L66:
  mov    0x20(%rsp),%rcx
  mov    0x40(%rsp),%rax
  add    %rcx,%rax
  mov    %rax,0x10(%rsp)
  mov    0x30(%rsp),%rax
  add    %rcx,%rax
  mov    %rax,0x8(%rsp)
  call   <abs@plt>
  divss  0x6c(%rsp),%xmm0
  call   <ceilf@plt>
  mov    0x68(%rsp),%ecx
  mov    0x38(%rsp),%rdi
  mov    %eax,%r15d
  imul   %ebp,%r15d
  sub    %r13d,%ecx
  mov    %ecx,0x4(%rsp)
  movslq %r15d,%rax
  lea    (%rdi,%rax,4),%r13
  cmp    %r15d,%ecx
  jl     L7
  mov    %r12,%rax
  mov    %r15d,%r14d
  mov    %ebp,%r12d
  mov    %r13,%r15
  mov    %rax,%rbp
  nopl   0x0(%rax)
L32:
  mov    %r14d,%eax
  mov    %r15,%rcx
  mov    %rbx,%rdi
  add    %r12d,%r14d
  cltd
  mov    0xe0(%rsp),%rsi
  idiv   %r12d
  mov    %rbp,%rdx
  add    (%rsp),%eax
  cltq
  lea    (%rsi,%rax,4),%r13
  mov    0xe8(%rsp),%eax
  push   %rax
  mov    0xd8(%rsp),%eax
  push   %rax
  mov    0xd8(%rsp),%eax
  push   %rax
  mov    0xd0(%rsp),%eax
  push   %rax
  mov    0xd0(%rsp),%r9d
  mov    0x30(%rsp),%r8
  mov    0x28(%rsp),%rsi
  call   <network_corr_precise@plt>
  add    0x38(%rsp),%r15
  add    $0x20,%rsp
  movss  %xmm0,0x0(%r13)
  cmp    %r14d,0x4(%rsp)
  jge    L32
  mov    %r12d,%ebp
L7:
  mov    0x58(%rsp),%rsi
  addl   $0x1,0x48(%rsp)
  mov    0xd8(%rsp),%edi
  add    %rsi,0x28(%rsp)
  add    %edi,(%rsp)
  mov    0x48(%rsp),%eax
  add    %rsi,0x20(%rsp)
  add    0x50(%rsp),%rbx
  cmp    %eax,0xc0(%rsp)
  jne    L8
L34:
  add    $0x78,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L29:
  mov    0x60(%rsp),%rax
  mov    0x20(%rsp),%rsi
  xor    %r13d,%r13d
  xor    %edi,%edi
  lea    (%rax,%rsi,1),%r12
  jmp    L66
```

assistant: ```c
void matched_filter_precise(float *templates, float *sum_square_templates, int *moveouts, float *data, float *weights, int step, int n_samples_template, int n_samples_data, int n_templates, int n_stations, int n_components, int n_corr, float *cc_sum, int normalize) {
  int t, ch, i;
  int start_i, stop_i, cc_i;
  int min_moveout, max_moveout;
  int network_offset, station_offset, cc_sum_offset;
  int *moveouts_t = NULL;
  float *templates_t = NULL, *sum_square_templates_t = NULL, *weights_t = NULL;
  for (t = 0; t < n_templates; t++) {
    network_offset = t * n_stations * n_components;
    station_offset = t * n_stations;
    cc_sum_offset = t * n_corr;
    min_moveout = 0;
    max_moveout = 0;
    for (ch = 0; ch < (n_stations * n_components); ch++) {
      if (moveouts[network_offset + ch] < min_moveout)
        min_moveout = moveouts[network_offset + ch];
      if (moveouts[network_offset + ch] > max_moveout)
        max_moveout = moveouts[network_offset + ch];
    }
    templates_t = templates + network_offset * n_samples_template;
    moveouts_t = moveouts + network_offset;
    weights_t = weights + network_offset;
    sum_square_templates_t = sum_square_templates + network_offset;
    start_i = (int)(ceilf(abs(min_moveout) / (float)step)) * step;
    stop_i = 1 + (n_samples_data - n_samples_template - max_moveout);
#pragma omp parallel for private(i, cc_i)
    for (i = start_i; i < stop_i; i += step) {
      cc_i = i / step;
      cc_sum[cc_sum_offset + cc_i] = network_corr_precise(templates_t, sum_square_templates_t, moveouts_t, data + i, weights_t, n_samples_template, n_samples_data, n_stations, n_components, normalize);
    }
  }
}
```