user: 下面的汇编指令对应的C源代码是什么:

```asm
<iceform>:
  endbr64
  push   %rbp
  push   %rbx
  mov    %rdi,%r10
  mov    %rdx,%r11
  mov    0x18(%rsp),%rdx
  mov    0x20(%rsp),%rbp
  mov    0x30(%rsp),%rbx
  movq   $0x0,(%rdi)
  movq   $0x0,(%rbx)
  test   %ecx,%ecx
  jle    L71
  movapd %xmm0,%xmm12
  mov    %ecx,%edi
  mov    $0x0,%eax
  pxor   %xmm5,%xmm5
  sub    $0x1,%ecx
  movsd  D137(%rip),%xmm8
  movsd  D83(%rip),%xmm7
  movsd  D124(%rip),%xmm6
  movapd %xmm8,%xmm11
  movapd %xmm7,%xmm10
  movapd %xmm6,%xmm9
  jmp    L95
L37:
  movapd %xmm12,%xmm0
  subsd  %xmm13,%xmm0
  mulsd  %xmm3,%xmm0
  mulsd  %xmm11,%xmm0
  mulsd  (%r9,%rax,8),%xmm0
  movapd %xmm10,%xmm13
  subsd  %xmm1,%xmm13
  mulsd  %xmm13,%xmm0
  movsd  (%rdx,%rax,8),%xmm13
  addsd  0x8(%rdx),%xmm13
  mulsd  %xmm13,%xmm0
  mulsd  %xmm9,%xmm0
  jmp    L57
L62:
  movapd %xmm12,%xmm0
  subsd  %xmm13,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  %xmm8,%xmm0
  mulsd  (%r9,%rax,8),%xmm0
  movapd %xmm7,%xmm13
  subsd  %xmm1,%xmm13
  mulsd  %xmm13,%xmm0
  movsd  (%rdx,%rax,8),%xmm13
  addsd  0x8(%rdx,%rax,8),%xmm13
  mulsd  %xmm13,%xmm0
  mulsd  %xmm6,%xmm0
L57:
  movsd  %xmm12,(%rsi,%rax,8)
  addsd  %xmm0,%xmm5
L114:
  add    $0x1,%rax
  cmp    %rax,%rdi
  je     L77
L95:
  movsd  (%rsi,%rax,8),%xmm13
  comisd %xmm13,%xmm12
  jbe    L114
  test   %eax,%eax
  je     L37
  cmp    %eax,%ecx
  jne    L62
  movapd %xmm12,%xmm0
  subsd  %xmm13,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  D137(%rip),%xmm0
  mulsd  (%r9,%rax,8),%xmm0
  movsd  D83(%rip),%xmm13
  subsd  %xmm1,%xmm13
  mulsd  %xmm13,%xmm0
  mulsd  (%rdx,%rax,8),%xmm0
  jmp    L57
L71:
  pxor   %xmm5,%xmm5
L77:
  movapd %xmm5,%xmm0
  divsd  D110(%rip),%xmm0
  movsd  %xmm0,(%rbx)
  comisd %xmm0,%xmm4
  jbe    L123
  imul   $0xe10,%r8d,%r8d
  pxor   %xmm0,%xmm0
  cvtsi2sd %r8d,%xmm0
  mulsd  (%rdx),%xmm0
  movsd  D83(%rip),%xmm2
  subsd  %xmm1,%xmm2
  mulsd  %xmm2,%xmm0
  movapd %xmm5,%xmm2
  divsd  %xmm0,%xmm2
  divsd  D139(%rip),%xmm5
L40:
  movsd  %xmm2,(%r10)
  movsd  D47(%rip),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movapd %xmm5,%xmm2
  divsd  %xmm0,%xmm2
  movsd  %xmm2,(%r11)
  movsd  D83(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  (%rdx),%xmm1
  comisd %xmm1,%xmm2
  jbe    L26
  divsd  %xmm1,%xmm5
  movsd  %xmm5,0x0(%rbp)
  mulsd  (%rdx),%xmm0
  movsd  %xmm0,(%r11)
L26:
  pop    %rbx
  pop    %rbp
  ret
L123:
  comisd D87(%rip),%xmm4
  jbe    L31
  movsd  %xmm4,(%rbx)
  mulsd  D137(%rip),%xmm4
  divsd  D77(%rip),%xmm4
  movapd %xmm4,%xmm5
  imul   $0xe10,%r8d,%r8d
  pxor   %xmm0,%xmm0
  cvtsi2sd %r8d,%xmm0
  mulsd  (%rdx),%xmm0
  movsd  D83(%rip),%xmm2
  subsd  %xmm1,%xmm2
  mulsd  %xmm2,%xmm0
  movapd %xmm4,%xmm2
  divsd  %xmm0,%xmm2
  jmp    L40
L31:
  movq   $0x0,(%rbx)
  pxor   %xmm2,%xmm2
  movapd %xmm2,%xmm5
  jmp    L40
```

tool_calls: name: retrieve_data
arguments: label: D137
type: f64
name: retrieve_data
arguments: label: D83
type: f64
name: retrieve_data
arguments: label: D124
type: f32[2]
name: retrieve_data
arguments: label: D110
type: byte[8]
name: retrieve_data
arguments: label: D139
type: byte[8]
name: retrieve_data
arguments: label: D47
type: f64
name: retrieve_data
arguments: label: D87
type: f32[4]
name: retrieve_data
arguments: label: D77
type: f64

tool: D137: 999.842594

tool: D83: 1

tool: {"D124": [0.0,1.75]}

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

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

tool: D47: 0.1

tool: {"D87": [0.0,0.0,-1.6012746345487803e-18,4.4764084815979]}

tool: D77: 917

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