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

```asm
<model_>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x78,%rsp
  mov    %rcx,%rdx
  mov    %rcx,-0x78(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x90(%rbp)
  mov    0x20(%rbp),%rcx
  mov    0x28(%rbp),%rbx
  mov    0x30(%rbp),%r15
  mov    0x38(%rbp),%r13
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  movslq (%rdx),%rax
  lea    0xf(,%rax,8),%rax
  mov    %rax,%rsi
  and    $0xfffffffffffffff0,%rsi
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rdx
  sub    %rax,%rdx
L17:
  cmp    %rdx,%rsp
  je     L65
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L17
L65:
  mov    %rsi,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L16
  orq    $0x0,-0x8(%rsp,%rax,1)
L16:
  mov    %rsp,-0x88(%rbp)
  movsd  (%rcx),%xmm7
  movsd  %xmm7,-0x98(%rbp)
  movsd  0x8(%rcx),%xmm6
  movsd  %xmm6,-0x48(%rbp)
  movsd  0x10(%rcx),%xmm1
  movsd  D24(%rip),%xmm0
  divsd  %xmm7,%xmm0
  movsd  %xmm0,(%rbx)
  movsd  D16(%rip),%xmm0
  call   <pow@plt>
  movsd  %xmm0,-0x80(%rbp)
  movsd  %xmm0,0x8(%rbx)
  mov    -0x78(%rbp),%rax
  mov    (%rax),%r12d
  mov    $0x1,%ebx
  pxor   %xmm6,%xmm6
  movsd  %xmm6,-0x50(%rbp)
  test   %r12d,%r12d
  jg     L20
L70:
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L51
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L12:
  pxor   %xmm7,%xmm7
  movsd  %xmm7,-0x50(%rbp)
L4:
  add    $0x1,%rbx
L20:
  mov    -0x88(%rbp),%rax
  movsd  -0x50(%rbp),%xmm3
  movsd  %xmm3,-0x8(%rax,%rbx,8)
  mov    %rbx,%r14
  shl    $0x4,%r14
  movsd  -0x8(%r15,%r14,1),%xmm2
  movsd  %xmm2,-0x58(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movapd %xmm2,%xmm0
  call   <pow@plt>
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  -0x80(%rbp),%xmm0
  call   <pow@plt>
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  -0x58(%rbp),%xmm0
  call   <pow@plt>
  movsd  -0x68(%rbp),%xmm1
  addsd  %xmm0,%xmm1
  movsd  -0x60(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r13,%rbx,8)
  lea    -0x1(%r12,%rbx,1),%eax
  cltq
  movsd  D19(%rip),%xmm4
  movsd  %xmm4,0x0(%r13,%rax,8)
  mov    -0x78(%rbp),%rax
  mov    (%rax),%r12d
  cmp    %ebx,%r12d
  jle    L70
  lea    (%rbx,%rbx,1),%eax
  test   %ebx,%ebx
  je     L12
  mov    -0x90(%rbp),%rdx
  mov    (%rdx),%esi
  mov    %esi,-0x58(%rbp)
  sub    %esi,-0x6c(%rbp)
  mov    -0x6c(%rbp),%edi
  sub    %esi,%eax
  cltq
  movsd  (%r15,%r14,1),%xmm0
  subsd  (%r15,%rax,8),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  movslq %edi,%rax
  movsd  (%r15,%rax,8),%xmm1
  movsd  -0x50(%rbp),%xmm5
  subsd  %xmm5,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  %xmm5,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  jmp    L4
L51:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: tag: D24
data_type: byte[8]
name: retrieve_data
arguments: tag: D16
data_type: f64
name: retrieve_data
arguments: tag: D19
data_type: f64

tool: {"D24": [0xef,0x39,0xfa,0xfe,0x42,0x2e,0xe6,0x3f]}

tool: D16: 10

tool: D19: 1

assistant: ```c
void model_(int *CurSet, int *NoOfParams, int *NoOfDerived, int *TotalDataValues, int *MaxNoOfDataValues, int *NoOfDataCols, int *NoOfAdirissaCols, int *NoOfModelVectors, double Params[], double Derived[], double Adirissa[], double Signal[]) {
  int CurEntry, time_j, cp_j, time_previous, cp_previous, emax_i, e0_i;
  double Ce[*TotalDataValues], ke0, thalf_effect, nHill, q, ec50, e0, emax, dt;
  ke0 = Params[0];
  nHill = Params[1];
  q = Params[2];
  e0 = Params[3];
  emax = Params[4];
  thalf_effect = log(2) / ke0;
  Derived[0] = thalf_effect;
  ec50 = pow(10.0, q);
  Derived[1] = ec50;
  double pct_effect(double conc) { return (pow(conc, nHill) / (pow(ec50, nHill) + pow(conc, nHill))); }
  for (CurEntry = 0; CurEntry < *TotalDataValues; CurEntry++) {
    time_j = 2 * CurEntry;
    cp_j = 2 * CurEntry + 1;
    if (CurEntry == 0) {
      dt = Adirissa[time_j];
      Ce[CurEntry] = 0;
    } else {
      time_previous = time_j - *NoOfAdirissaCols;
      cp_previous = cp_previous - *NoOfAdirissaCols;
      dt = Adirissa[time_j] - Adirissa[time_previous];
      Ce[CurEntry] = Ce[CurEntry - 1] + dt * ke0 * (Adirissa[cp_previous] - Ce[CurEntry - 1]);
    }
    emax_i = CurEntry;
    e0_i = CurEntry + 1 * (*TotalDataValues);
    Signal[emax_i] = pct_effect(Adirissa[cp_j]);
    Signal[e0_i] = 1.0;
  }
  return;
}
```