user: 请从给定的汇编指令推导出C源代码:

```asm
<iceform>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  movsd  %xmm0,-0x38(%rbp)
  movsd  %xmm1,-0x40(%rbp)
  mov    %rdx,-0x48(%rbp)
  mov    %ecx,-0x4c(%rbp)
  mov    %r8d,-0x50(%rbp)
  movsd  %xmm2,-0x58(%rbp)
  movsd  %xmm3,-0x60(%rbp)
  mov    %r9,-0x68(%rbp)
  movsd  %xmm4,-0x70(%rbp)
  mov    -0x28(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    0x28(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  movl   $0x0,-0x1c(%rbp)
  jmp    L40
L16:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x38(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L52
  cmpl   $0x0,-0x1c(%rbp)
  jne    L107
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x38(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movsd  D163(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D49(%rip),%xmm0
  subsd  -0x40(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D205(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  jmp    L93
L107:
  mov    -0x4c(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0x1c(%rbp)
  jne    L48
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x38(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x58(%rbp),%xmm1
  movsd  D163(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D49(%rip),%xmm0
  subsd  -0x40(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  jmp    L93
L48:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x38(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x58(%rbp),%xmm1
  movsd  D163(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D49(%rip),%xmm0
  subsd  -0x40(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  mov    -0x1c(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D205(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
L93:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x18(%rbp),%xmm0
  addsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
L52:
  addl   $0x1,-0x1c(%rbp)
L40:
  mov    -0x1c(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jl     L16
  movsd  -0x18(%rbp),%xmm0
  movsd  D184(%rip),%xmm1
  divsd  %xmm1,%xmm0
  mov    0x28(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    0x28(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x70(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L73
  mov    -0x50(%rbp),%eax
  imul   $0xe10,%eax,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mov    0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D49(%rip),%xmm0
  subsd  -0x40(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x18(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  mov    -0x28(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  -0x18(%rbp),%xmm0
  movsd  D62(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  jmp    L60
L73:
  movsd  -0x70(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L106
  mov    0x28(%rbp),%rax
  movsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    0x28(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D163(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D141(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    0x28(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D163(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D141(%rip),%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x50(%rbp),%eax
  imul   $0xe10,%eax,%eax
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  mov    0x10(%rbp),%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm1,%xmm2
  movsd  D49(%rip),%xmm1
  subsd  -0x40(%rbp),%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x28(%rbp),%rax
  movsd  %xmm0,(%rax)
  jmp    L60
L106:
  mov    0x28(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x28(%rbp),%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
L60:
  mov    0x18(%rbp),%rax
  movsd  D36(%rip),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x8(%rbp),%xmm0
  movsd  D36(%rip),%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x48(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x48(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  D49(%rip),%xmm1
  movapd %xmm1,%xmm2
  subsd  -0x40(%rbp),%xmm2
  mov    0x10(%rbp),%rax
  movsd  (%rax),%xmm1
  mulsd  %xmm2,%xmm1
  comisd %xmm1,%xmm0
  ja     L71
  jmp    L114
L71:
  movsd  D49(%rip),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x40(%rbp),%xmm1
  mov    0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x8(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  mov    0x18(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  D49(%rip),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x40(%rbp),%xmm1
  mov    0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x48(%rbp),%rax
  movsd  %xmm0,(%rax)
L114:
  nop
  pop    %rbp
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D163
type_infered: f64
name: retrieve_data
arguments: tag: D49
type_infered: f64
name: retrieve_data
arguments: tag: D205
type_infered: f64
name: retrieve_data
arguments: tag: D184
type_infered: byte[8]
name: retrieve_data
arguments: tag: D62
type_infered: byte[8]
name: retrieve_data
arguments: tag: D141
type_infered: f64
name: retrieve_data
arguments: tag: D36
type_infered: f64

tool: D163: 999.842594

tool: D49: 1

tool: D205: 2

tool: {"D184": [0x24,0x28,0x9e,0x71,0x0e,0xe3,0xb3,0x41]}

tool: {"D62": [0x00,0x00,0x00,0xd4,0x3b,0x3d,0xb2,0x41]}

tool: D141: 917

tool: D36: 0.1

assistant: ```c
void iceform(double *qfusion, double *T, double Tcutoff, double fracprv, double *areaadd, int numnod, int dt, double dz, double surfdz, double *cp, double *surface, double *new_ice_height, double *water_density, double *new_ice_water_eq, double lvolume) {
  double sum, extra;
  int j;
  double xfrac, di;
  *qfusion = 0.0;
  sum = 0.0;
  *new_ice_water_eq = 0.0;
  for (j = 0; j < numnod; j++) {
    if (T[j] < Tcutoff) {
      if (j == 0)
        extra = ((Tcutoff - T[j]) * surfdz * 999.842594 * cp[j] * (1.0 - fracprv) * (surface[j] + surface[j + 1]) / 2.);
      else if (j == numnod - 1)
        extra = ((Tcutoff - T[j]) * dz * 999.842594 * cp[j] * (1.0 - fracprv) * surface[j]);
      else
        extra = ((Tcutoff - T[j]) * dz * 999.842594 * cp[j] * (1.0 - fracprv) * (surface[j] + surface[j + 1]) / 2.);
      T[j] = Tcutoff;
      sum += extra;
    }
  }
  *new_ice_water_eq = sum / (999.842594 * 3.337e5);
  if (lvolume > *new_ice_water_eq) {
    *qfusion = (sum / (dt * 3600 * surface[0] * (1.0 - fracprv)));
    di = sum / (3.337e5 * 917.);
  } else if (lvolume > 0.0) {
    *new_ice_water_eq = lvolume;
    di = *new_ice_water_eq * 999.842594 / 917.;
    *qfusion = (*new_ice_water_eq * 999.842594 / 917.) / (dt * 3600 * surface[0] * (1.0 - fracprv));
  } else {
    *new_ice_water_eq = 0.0;
    di = 0.0;
    *qfusion = 0.0;
  }
  *new_ice_height = 0.10;
  *areaadd = di / (0.10);
  if (*areaadd > (1.0 - fracprv) * surface[0]) {
    *new_ice_height = di / ((1. - fracprv) * surface[0]);
    *areaadd = (1. - fracprv) * surface[0];
  }
}
```