user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<iceform>:
  endbr64
  push   %rbp
  movsd  D27(%rip),%xmm7
  mov    %rdx,%r11
  mov    %rdi,%r10
  push   %rbx
  mov    0x30(%rsp),%rbx
  movq   $0x0,(%rdi)
  subsd  %xmm1,%xmm7
  mov    0x18(%rsp),%rdx
  movq   $0x0,(%rbx)
  mov    0x20(%rsp),%rbp
  test   %ecx,%ecx
  jle    L100
  movslq %ecx,%rdi
  movapd %xmm0,%xmm5
  xor    %eax,%eax
  sub    $0x1,%ecx
  movsd  D31(%rip),%xmm8
  pxor   %xmm6,%xmm6
  movsd  D49(%rip),%xmm9
  jmp    L94
  nopl   0x0(%rax)
L1:
  mulsd  %xmm2,%xmm0
  mulsd  %xmm8,%xmm0
  mulsd  %xmm11,%xmm0
  mulsd  %xmm7,%xmm0
  cmp    %eax,%ecx
  je     L21
  addsd  0x8(%rdx,%rax,8),%xmm1
  mulsd  %xmm1,%xmm0
  mulsd  %xmm9,%xmm0
L130:
  addsd  %xmm0,%xmm6
  movsd  %xmm5,(%rsi,%rax,8)
L111:
  add    $0x1,%rax
  cmp    %rdi,%rax
  je     L31
L94:
  movsd  (%rsi,%rax,8),%xmm10
  comisd %xmm10,%xmm5
  jbe    L111
  movapd %xmm5,%xmm0
  movsd  (%r9,%rax,8),%xmm11
  movsd  (%rdx,%rax,8),%xmm1
  subsd  %xmm10,%xmm0
  test   %rax,%rax
  jne    L1
  mulsd  %xmm3,%xmm0
  addsd  0x8(%rdx),%xmm1
  movapd %xmm0,%xmm10
  movapd %xmm1,%xmm0
  mulsd  %xmm8,%xmm10
  mulsd  %xmm11,%xmm10
  mulsd  %xmm7,%xmm10
  mulsd  %xmm10,%xmm0
  mulsd  %xmm9,%xmm0
  jmp    L130
  cs nopw 0x0(%rax,%rax,1)
L31:
  movapd %xmm6,%xmm0
  divsd  D12(%rip),%xmm0
L41:
  comisd %xmm0,%xmm4
  movsd  %xmm0,(%rbx)
  jbe    L19
  imul   $0xe10,%r8d,%r8d
  pxor   %xmm1,%xmm1
  movapd %xmm6,%xmm0
  movsd  D62(%rip),%xmm2
  divsd  D15(%rip),%xmm6
  cvtsi2sd %r8d,%xmm1
  mulsd  (%rdx),%xmm1
  mulsd  %xmm7,%xmm1
  divsd  %xmm1,%xmm0
  movapd %xmm6,%xmm1
  divsd  %xmm2,%xmm1
L64:
  movsd  %xmm0,(%r10)
  movsd  %xmm2,0x0(%rbp)
  movsd  %xmm1,(%r11)
  movsd  (%rdx),%xmm0
  mulsd  %xmm7,%xmm0
  comisd %xmm0,%xmm1
  jbe    L50
  divsd  %xmm0,%xmm6
  movsd  %xmm6,0x0(%rbp)
  mulsd  (%rdx),%xmm7
  movsd  %xmm7,(%r11)
L50:
  pop    %rbx
  pop    %rbp
  ret
L21:
  mulsd  %xmm1,%xmm0
  jmp    L130
  nopl   0x0(%rax)
L19:
  pxor   %xmm3,%xmm3
  comisd %xmm3,%xmm4
  jbe    L58
  imul   $0xe10,%r8d,%r8d
  pxor   %xmm1,%xmm1
  movsd  %xmm4,(%rbx)
  movsd  D62(%rip),%xmm2
  mulsd  D31(%rip),%xmm4
  divsd  D76(%rip),%xmm4
  movapd %xmm4,%xmm0
  movapd %xmm4,%xmm6
  cvtsi2sd %r8d,%xmm1
  mulsd  (%rdx),%xmm1
  mulsd  %xmm7,%xmm1
  divsd  %xmm1,%xmm0
  movapd %xmm4,%xmm1
  divsd  %xmm2,%xmm1
  jmp    L64
  nopw   0x0(%rax,%rax,1)
L58:
  pxor   %xmm0,%xmm0
  movq   $0x0,(%rbx)
  movsd  D62(%rip),%xmm2
  movapd %xmm0,%xmm1
  movapd %xmm0,%xmm6
  jmp    L64
L100:
  pxor   %xmm0,%xmm0
  movapd %xmm0,%xmm6
  jmp    L41
```

tool_calls: name: parse_data
arguments: tag: D27
type: f64
name: parse_data
arguments: tag: D31
type: f64
name: parse_data
arguments: tag: D49
type: f32[2]
name: parse_data
arguments: tag: D12
type: byte[8]
name: parse_data
arguments: tag: D62
type: f64
name: parse_data
arguments: tag: D15
type: byte[8]
name: parse_data
arguments: tag: D76
type: f64

tool: D27: 1

tool: D31: 999.842594

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

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

tool: D62: 0.1

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

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