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

```asm
<hinit>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xb0,%rsp
  mov    %edi,-0x64(%rbp)
  movsd  %xmm0,-0x70(%rbp)
  mov    %rsi,-0x78(%rbp)
  mov    %rdx,-0x80(%rbp)
  movsd  %xmm1,-0x88(%rbp)
  mov    %rcx,-0x90(%rbp)
  mov    %r8,-0x98(%rbp)
  mov    %r9,-0xa0(%rbp)
  movsd  %xmm2,-0xa8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    0x18(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  mov    0x20(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  cmpl   $0x0,0x28(%rbp)
  jne    L177
  movl   $0x0,-0x54(%rbp)
  jmp    L96
L91:
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D11(%rip),%xmm1
  andpd  %xmm1,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  movsd  -0x28(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0x50(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0x48(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  addl   $0x1,-0x54(%rbp)
L96:
  mov    -0x54(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jb     L91
  jmp    L41
L177:
  movl   $0x0,-0x54(%rbp)
  jmp    L66
L37:
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D11(%rip),%xmm3
  andpd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0x50(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  divsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0x48(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  addl   $0x1,-0x54(%rbp)
L66:
  mov    -0x54(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jb     L37
L41:
  movsd  D132(%rip),%xmm0
  comisd -0x50(%rbp),%xmm0
  jae    L10
  movsd  D132(%rip),%xmm0
  comisd -0x48(%rbp),%xmm0
  jb     L112
L10:
  movsd  D147(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  jmp    L59
L112:
  movsd  -0x48(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movsd  D104(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
L59:
  movsd  -0xa8(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <min_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x88(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movl   $0x0,-0x54(%rbp)
  jmp    L107
L170:
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x40(%rbp),%xmm0
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x54(%rbp)
L107:
  mov    -0x54(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jb     L170
  movsd  -0x70(%rbp),%xmm0
  addsd  -0x40(%rbp),%xmm0
  movq   %xmm0,%rdi
  mov    -0x98(%rbp),%rcx
  mov    -0x80(%rbp),%rdx
  mov    -0xa0(%rbp),%rsi
  mov    -0x64(%rbp),%eax
  movq   %rdi,%xmm0
  mov    %eax,%edi
  mov    $0x1,%eax
  call   <safeOdeFunc@plt>
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  cmpl   $0x0,0x28(%rbp)
  jne    L40
  movl   $0x0,-0x54(%rbp)
  jmp    L122
L4:
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D11(%rip),%xmm1
  andpd  %xmm1,%xmm0
  mulsd  -0x20(%rbp),%xmm0
  movsd  -0x28(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  divsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0x30(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  addl   $0x1,-0x54(%rbp)
L122:
  mov    -0x54(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jb     L4
  jmp    L27
L40:
  movl   $0x0,-0x54(%rbp)
  jmp    L175
L155:
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D11(%rip),%xmm3
  andpd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x54(%rbp),%eax
  lea    0x0(,%rax,8),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  divsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  -0x30(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  addl   $0x1,-0x54(%rbp)
L175:
  mov    -0x54(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jb     L155
L27:
  mov    -0x30(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  divsd  -0x40(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rdx
  movsd  -0x30(%rbp),%xmm0
  movq   D11(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rdx,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <max_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  D70(%rip),%xmm0
  comisd -0x8(%rbp),%xmm0
  jb     L95
  movsd  -0x40(%rbp),%xmm0
  movq   D11(%rip),%xmm1
  andpd  %xmm0,%xmm1
  movsd  D119(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mov    D147(%rip),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <max_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  jmp    L65
L95:
  pxor   %xmm1,%xmm1
  cvtsi2sdl 0x10(%rbp),%xmm1
  movsd  D107(%rip),%xmm0
  movapd %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  movsd  D104(%rip),%xmm0
  divsd  -0x8(%rbp),%xmm0
  movq   %xmm0,%rax
  movapd %xmm2,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x38(%rbp)
L65:
  movsd  -0xa8(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <min_d@plt>
  mov    %eax,%edx
  movsd  -0x40(%rbp),%xmm1
  movsd  D158(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  mov    %edx,%edi
  movq   %rax,%xmm0
  mov    $0x1,%eax
  call   <min_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm0
  ucomisd -0x40(%rbp),%xmm0
  jnp    L45
  int3
L45:
  movsd  -0x88(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: label: D11
datatype: qword
name: disassemble_data
arguments: label: D132
datatype: byte[8]
name: disassemble_data
arguments: label: D147
datatype: byte[8]
name: disassemble_data
arguments: label: D104
datatype: f64
name: disassemble_data
arguments: label: D70
datatype: byte[8]
name: disassemble_data
arguments: label: D119
datatype: byte[8]
name: disassemble_data
arguments: label: D107
datatype: f64
name: disassemble_data
arguments: label: D158
datatype: f64

tool: D11: 9223372036854775807

tool: {"D132": [0xbb,0xbd,0xd7,0xd9,0xdf,0x7c,0xdb,0x3d]}

tool: {"D147": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xb0,0x3e]}

tool: D104: 0.01

tool: {"D70": [0x16,0x56,0xe7,0x9e,0xaf,0x03,0xd2,0x3c]}

tool: {"D119": [0xfc,0xa9,0xf1,0xd2,0x4d,0x62,0x50,0x3f]}

tool: D107: 1

tool: D158: 100

assistant: ```c
double hinit(unsigned n, double x, double *y, double *pars, double posneg, double *f0, double *f1, double *yy1, int iord, double hmax, double *atoler, double *rtoler, int itoler) {
  double dnf, dny, atoli, rtoli, sk, h, h1, der2, der12, sqr;
  unsigned i;
  dnf = 0.0;
  dny = 0.0;
  atoli = atoler[0];
  rtoli = rtoler[0];
  if (!itoler)
    for (i = 0; i < n; i++) {
      sk = atoli + rtoli * fabs(y[i]);
      sqr = f0[i] / sk;
      dnf += sqr * sqr;
      sqr = y[i] / sk;
      dny += sqr * sqr;
    }
  else
    for (i = 0; i < n; i++) {
      sk = atoler[i] + rtoler[i] * fabs(y[i]);
      sqr = f0[i] / sk;
      dnf += sqr * sqr;
      sqr = y[i] / sk;
      dny += sqr * sqr;
    }
  if ((dnf <= 1.0E-10) || (dny <= 1.0E-10))
    h = 1.0E-6;
  else
    h = sqrt(dny / dnf) * 0.01;
  h = min_d(h, hmax);
  h = sign(h, posneg);
  for (i = 0; i < n; i++)
    yy1[i] = y[i] + h * f0[i];
  safeOdeFunc(n, x + h, yy1, pars, f1);
  der2 = 0.0;
  if (!itoler)
    for (i = 0; i < n; i++) {
      sk = atoli + rtoli * fabs(y[i]);
      sqr = (f1[i] - f0[i]) / sk;
      der2 += sqr * sqr;
    }
  else
    for (i = 0; i < n; i++) {
      sk = atoler[i] + rtoler[i] * fabs(y[i]);
      sqr = (f1[i] - f0[i]) / sk;
      der2 += sqr * sqr;
    }
  der2 = sqrt(der2) / h;
  der12 = max_d(fabs(der2), sqrt(dnf));
  if (der12 <= 1.0E-15)
    h1 = max_d(1.0E-6, fabs(h) * 1.0E-3);
  else
    h1 = pow(0.01 / der12, 1.0 / (double)iord);
  h = min_d(100.0 * h, min_d(h1, hmax));
  if (isnan(h)) {
    __asm__("int $3");
  };
  return sign(h, posneg);
}
```