user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<hinit>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,%r12d
  movsd  %xmm0,0x30(%rsp)
  mov    %rsi,%rbx
  mov    %rdx,0x38(%rsp)
  movsd  %xmm1,0x20(%rsp)
  mov    %rcx,%rbp
  mov    %r8,%r13
  mov    %r9,%r14
  movsd  %xmm2,0x28(%rsp)
  mov    0x90(%rsp),%r15
  mov    0x88(%rsp),%rax
  movsd  (%rax),%xmm5
  movsd  %xmm5,0x10(%rsp)
  movsd  (%r15),%xmm6
  movsd  %xmm6,0x18(%rsp)
  cmpl   $0x0,0x98(%rsp)
  je     L93
  test   %edi,%edi
  je     L183
  mov    %edi,%edx
  mov    $0x0,%eax
  pxor   %xmm3,%xmm3
  movsd  %xmm3,(%rsp)
  movq   D62(%rip),%xmm4
L245:
  movsd  (%rbx,%rax,8),%xmm0
  movapd %xmm0,%xmm1
  andpd  %xmm4,%xmm1
  mulsd  (%r15,%rax,8),%xmm1
  mov    0x88(%rsp),%rcx
  addsd  (%rcx,%rax,8),%xmm1
  movsd  0x0(%rbp,%rax,8),%xmm2
  divsd  %xmm1,%xmm2
  mulsd  %xmm2,%xmm2
  addsd  (%rsp),%xmm2
  movsd  %xmm2,(%rsp)
  divsd  %xmm1,%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm3
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L245
  jmp    L36
L93:
  test   %edi,%edi
  je     L183
  mov    %edi,%edx
  mov    $0x0,%eax
  pxor   %xmm3,%xmm3
  movsd  %xmm3,(%rsp)
  movq   D62(%rip),%xmm4
L232:
  movsd  (%rbx,%rax,8),%xmm0
  movapd %xmm0,%xmm1
  andpd  %xmm4,%xmm1
  mulsd  0x18(%rsp),%xmm1
  addsd  0x10(%rsp),%xmm1
  movsd  0x0(%rbp,%rax,8),%xmm2
  divsd  %xmm1,%xmm2
  mulsd  %xmm2,%xmm2
  addsd  (%rsp),%xmm2
  movsd  %xmm2,(%rsp)
  divsd  %xmm1,%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm3
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L232
L36:
  movsd  D120(%rip),%xmm0
  comisd (%rsp),%xmm0
  jae    L58
  movsd  D159(%rip),%xmm0
  movsd  D120(%rip),%xmm1
  comisd %xmm3,%xmm1
  jb     L40
L239:
  movsd  0x28(%rsp),%xmm1
  mov    $0x2,%eax
  call   <min_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  0x20(%rsp),%xmm1
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm6,%xmm6
  cvtsi2sd %eax,%xmm6
  movsd  %xmm6,0x8(%rsp)
  test   %r12d,%r12d
  je     L161
L208:
  mov    $0x0,%eax
L149:
  movsd  0x8(%rsp),%xmm0
  mulsd  0x0(%rbp,%rax,8),%xmm0
  addsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,(%r14,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%r12d
  ja     L149
  movsd  0x30(%rsp),%xmm0
  addsd  0x8(%rsp),%xmm0
  mov    %r13,%rcx
  mov    0x38(%rsp),%rdx
  mov    %r14,%rsi
  mov    %r12d,%edi
  mov    $0x1,%eax
  call   <safeOdeFunc@plt>
  cmpl   $0x0,0x98(%rsp)
  je     L151
  test   %r12d,%r12d
  je     L49
  mov    %r12d,%r12d
  mov    $0x0,%eax
  pxor   %xmm0,%xmm0
  movq   D62(%rip),%xmm3
L15:
  movsd  0x0(%r13,%rax,8),%xmm1
  subsd  0x0(%rbp,%rax,8),%xmm1
  movsd  (%rbx,%rax,8),%xmm2
  andpd  %xmm3,%xmm2
  mulsd  (%r15,%rax,8),%xmm2
  mov    0x88(%rsp),%rsi
  addsd  (%rsi,%rax,8),%xmm2
  divsd  %xmm2,%xmm1
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm0
  add    $0x1,%rax
  cmp    %rax,%r12
  jne    L15
  jmp    L203
L40:
  divsd  (%rsp),%xmm3
  movapd %xmm3,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm3,%xmm1
  ja     L126
  sqrtsd %xmm0,%xmm0
L228:
  mulsd  D77(%rip),%xmm0
  jmp    L239
L126:
  call   <sqrt@plt>
  jmp    L228
L151:
  mov    $0x0,%eax
  pxor   %xmm0,%xmm0
  movq   D62(%rip),%xmm3
L138:
  movsd  0x0(%r13,%rax,8),%xmm1
  subsd  0x0(%rbp,%rax,8),%xmm1
  movsd  (%rbx,%rax,8),%xmm2
  andpd  %xmm3,%xmm2
  mulsd  0x18(%rsp),%xmm2
  addsd  0x10(%rsp),%xmm2
  divsd  %xmm2,%xmm1
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm0
  add    $0x1,%rax
  cmp    %eax,%r12d
  ja     L138
L203:
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L4
L198:
  sqrtsd %xmm0,%xmm0
L127:
  divsd  0x8(%rsp),%xmm0
  movq   %xmm0,%rbx
  pxor   %xmm0,%xmm0
  movsd  (%rsp),%xmm4
  ucomisd %xmm4,%xmm0
  ja     L240
  sqrtsd %xmm4,%xmm4
  movapd %xmm4,%xmm1
L199:
  movq   %rbx,%xmm0
  andpd  D62(%rip),%xmm0
  mov    $0x2,%eax
  call   <max_d@plt>
  test   %eax,%eax
  jle    L73
  pxor   %xmm2,%xmm2
  cvtsi2sdl 0x80(%rsp),%xmm2
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D77(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  D55(%rip),%xmm1
  divsd  %xmm2,%xmm1
  call   <pow@plt>
L155:
  movsd  0x28(%rsp),%xmm1
  mov    $0x2,%eax
  call   <min_d@plt>
  mov    %eax,%edi
  movsd  0x8(%rsp),%xmm0
  mulsd  D24(%rip),%xmm0
  mov    $0x1,%eax
  call   <min_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  0x20(%rsp),%xmm1
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L4:
  call   <sqrt@plt>
  jmp    L127
L240:
  movsd  (%rsp),%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  jmp    L199
L73:
  movsd  0x8(%rsp),%xmm1
  andpd  D62(%rip),%xmm1
  mulsd  D114(%rip),%xmm1
  movsd  D159(%rip),%xmm0
  mov    $0x2,%eax
  call   <max_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  jmp    L155
L58:
  movsd  0x28(%rsp),%xmm1
  movsd  D159(%rip),%xmm0
  mov    $0x2,%eax
  call   <min_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  0x20(%rsp),%xmm1
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm7,%xmm7
  cvtsi2sd %eax,%xmm7
  movsd  %xmm7,0x8(%rsp)
  jmp    L208
L183:
  movsd  0x28(%rsp),%xmm1
  movsd  D159(%rip),%xmm0
  mov    $0x2,%eax
  call   <min_d@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  0x20(%rsp),%xmm1
  mov    $0x2,%eax
  call   <sign@plt>
  pxor   %xmm4,%xmm4
  cvtsi2sd %eax,%xmm4
  movsd  %xmm4,0x8(%rsp)
  pxor   %xmm7,%xmm7
  movsd  %xmm7,(%rsp)
L161:
  movsd  0x30(%rsp),%xmm0
  addsd  0x8(%rsp),%xmm0
  mov    %r13,%rcx
  mov    0x38(%rsp),%rdx
  mov    %r14,%rsi
  mov    $0x0,%edi
  mov    $0x1,%eax
  call   <safeOdeFunc@plt>
L49:
  pxor   %xmm0,%xmm0
  jmp    L198
```

tool_calls: name: read_data
arguments: data_label: D62
type: qword[2]
name: read_data
arguments: data_label: D120
type: byte[8]
name: read_data
arguments: data_label: D159
type: byte[8]
name: read_data
arguments: data_label: D77
type: f64
name: read_data
arguments: data_label: D55
type: f64
name: read_data
arguments: data_label: D24
type: f64
name: read_data
arguments: data_label: D114
type: byte[8]

tool: {"D62": [0x7fffffffffffffff,0x0000000000000000]}

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

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

tool: D77: 0.01

tool: D55: 1

tool: D24: 100

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

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);
}
```