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

```asm
<dip_vacancy>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x428,%rsp
  mov    %edi,%r13d
  mov    %esi,%r14d
  mov    %rdx,%rbx
  mov    %rcx,%rbp
  mov    %r8d,%r15d
  movss  %xmm0,0x4(%rsp)
  mov    0x460(%rsp),%r12d
  mov    %fs:0x28,%rax
  mov    %rax,0x418(%rsp)
  xor    %eax,%eax
  test   %r12d,%r12d
  jle    L146
  lea    -0x1(%r12),%ecx
  mov    %r12d,%edx
  pxor   %xmm0,%xmm0
L144:
  addss  (%r9,%rax,4),%xmm0
  movl   $0x0,0x10(%rsp,%rax,4)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L144
  pxor   %xmm2,%xmm2
  cvtsi2ss %r12d,%xmm2
  divss  %xmm2,%xmm0
  movss  %xmm0,0x8(%rsp)
  mov    %r9,%rax
  lea    0x4(%r9,%rcx,4),%rdx
  pxor   %xmm0,%xmm0
L153:
  movss  (%rax),%xmm1
  subss  0x8(%rsp),%xmm1
  mulss  %xmm1,%xmm1
  addss  %xmm1,%xmm0
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L153
L17:
  divss  %xmm2,%xmm0
  call   <sqrt@plt>
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  addss  %xmm0,%xmm0
  movss  0x8(%rsp),%xmm7
  subss  %xmm0,%xmm7
  movss  %xmm7,0xc(%rsp)
  test   %r15d,%r15d
  jle    L50
  lea    -0x1(%r15),%eax
  lea    0x4(%rbx,%rax,4),%r15
  pxor   %xmm2,%xmm2
  movss  %xmm2,0x8(%rsp)
  jmp    L118
L164:
  add    $0x4,%rbx
  add    $0x4,%rbp
  cmp    %r15,%rbx
  je     L150
L118:
  movss  (%rbx),%xmm0
  divss  D24(%rip),%xmm0
  movss  0x4(%rsp),%xmm1
  call   <fmod@plt>
  imul   %r12d,%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  divss  0x4(%rsp),%xmm0
  cvttss2si %xmm0,%eax
  movslq %eax,%rdx
  movss  D19(%rip),%xmm0
  addss  0x10(%rsp,%rdx,4),%xmm0
  movss  %xmm0,0x10(%rsp,%rdx,4)
  cmp    %eax,%r13d
  jg     L164
  cmp    %eax,%r14d
  jle    L164
  movss  0x0(%rbp),%xmm0
  comiss 0xc(%rsp),%xmm0
  jbe    L164
  movss  0x8(%rsp),%xmm6
  addss  D19(%rip),%xmm6
  movss  %xmm6,0x8(%rsp)
  jmp    L164
L50:
  pxor   %xmm7,%xmm7
  movss  %xmm7,0x8(%rsp)
L150:
  test   %r12d,%r12d
  jle    L139
  lea    0x10(%rsp),%rax
  lea    -0x1(%r12),%edx
  lea    0x14(%rsp,%rdx,4),%rdx
  pxor   %xmm0,%xmm0
L137:
  movss  (%rax),%xmm1
  maxss  %xmm0,%xmm1
  movaps %xmm1,%xmm0
  add    $0x4,%rax
  cmp    %rax,%rdx
  jne    L137
  comiss D23(%rip),%xmm1
  ja     L9
  movss  D19(%rip),%xmm0
L115:
  mov    0x418(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L73
  add    $0x428,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L9:
  sub    %r13d,%r14d
  add    $0x1,%r14d
  pxor   %xmm1,%xmm1
  cvtsi2ss %r14d,%xmm1
  mulss  %xmm0,%xmm1
  movss  0x8(%rsp),%xmm0
  divss  %xmm1,%xmm0
  jmp    L115
L139:
  movss  D19(%rip),%xmm0
  jmp    L115
L146:
  pxor   %xmm2,%xmm2
  cvtsi2ss %r12d,%xmm2
  pxor   %xmm0,%xmm0
  movaps %xmm0,%xmm7
  divss  %xmm2,%xmm7
  movss  %xmm7,0x8(%rsp)
  jmp    L17
L73:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_tag: D24
datatype: byte[4]
name: parse_data
arguments: data_tag: D19
datatype: f32
name: parse_data
arguments: data_tag: D23
datatype: f32

tool: {"D24": [0x00,0xc0,0xa8,0x47]}

tool: D19: 1

tool: D23: 0

assistant: ```c
float dip_vacancy(int start_index, int end_index, float timestamp[], float series[], int series_length, float period_days, float curve[], int curve_length) {
  int i, index;
  float days, max_samples = 0;
  float density = 0, curve_average_mag = 0;
  float curve_variance = 0, min_curve_mag;
  float den[256];
  for (i = 0; i < curve_length; i++) {
    curve_average_mag += curve[i];
    den[i] = 0;
  }
  curve_average_mag /= curve_length;
  for (i = 0; i < curve_length; i++) {
    curve_variance += (curve[i] - curve_average_mag) * (curve[i] - curve_average_mag);
  }
  curve_variance = (float)sqrt(curve_variance / curve_length);
  min_curve_mag = curve_average_mag - (curve_variance * 2.0f);
  for (i = 0; i < series_length; i++) {
    days = timestamp[i] / (60.0f * 60.0f * 24.0f);
    index = (int)(fmod(days, period_days) * curve_length / period_days);
    den[index]++;
    if ((index >= start_index) && (index < end_index)) {
      if (series[i] > min_curve_mag) {
        density++;
      }
    }
  }
  for (i = 0; i < curve_length; i++) {
    if (den[i] > max_samples) {
      max_samples = den[i];
    }
  }
  if (max_samples > 0) {
    return density / (max_samples * ((end_index - start_index) + 1));
  }
  return 1.0;
}
```