user: What is the c source code of the assembly code below?

```asm
<enorm>:
  endbr64
  movsd  D29(%rip),%xmm10
  pxor   %xmm0,%xmm0
  cvtsi2sd %edi,%xmm0
  divsd  %xmm0,%xmm10
  test   %edi,%edi
  jle    L82
  pxor   %xmm6,%xmm6
  lea    -0x1(%rdi),%eax
  sub    $0x18,%rsp
  movq   D2(%rip),%xmm7
  movsd  D40(%rip),%xmm4
  movapd %xmm6,%xmm3
  movapd %xmm6,%xmm2
  movapd %xmm6,%xmm5
  lea    0x8(%rsi,%rax,8),%rax
  movapd %xmm6,%xmm9
  movapd %xmm6,%xmm8
  movsd  D5(%rip),%xmm11
  jmp    L176
  nopl   0x0(%rax,%rax,1)
L174:
  comisd %xmm10,%xmm12
  jae    L227
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm5
L83:
  add    $0x8,%rsi
  cmp    %rsi,%rax
  je     L30
L176:
  movsd  (%rsi),%xmm1
  movapd %xmm1,%xmm12
  andpd  %xmm7,%xmm12
  comisd %xmm12,%xmm4
  jb     L174
L227:
  comisd %xmm4,%xmm12
  jbe    L3
  comisd %xmm2,%xmm12
  jbe    L138
  divsd  %xmm1,%xmm2
  add    $0x8,%rsi
  mulsd  %xmm2,%xmm2
  mulsd  %xmm2,%xmm8
  movapd %xmm12,%xmm2
  addsd  %xmm11,%xmm8
  cmp    %rsi,%rax
  jne    L176
  nopl   0x0(%rax,%rax,1)
L30:
  ucomisd %xmm6,%xmm8
  jp     L103
  jne    L103
  ucomisd %xmm6,%xmm5
  jp     L135
  je     L156
L135:
  mulsd  %xmm3,%xmm9
  comisd %xmm3,%xmm5
  jb     L37
  movapd %xmm3,%xmm0
  divsd  %xmm5,%xmm0
  mulsd  %xmm9,%xmm0
  addsd  D5(%rip),%xmm0
  mulsd  %xmm5,%xmm0
  ucomisd %xmm0,%xmm6
  ja     L9
L186:
  sqrtsd %xmm0,%xmm0
  add    $0x18,%rsp
  ret
  nopl   0x0(%rax,%rax,1)
L3:
  comisd %xmm3,%xmm12
  jbe    L140
  divsd  %xmm1,%xmm3
  mulsd  %xmm3,%xmm3
  mulsd  %xmm3,%xmm9
  movapd %xmm12,%xmm3
  addsd  %xmm11,%xmm9
  jmp    L83
  nopl   0x0(%rax,%rax,1)
L138:
  divsd  %xmm2,%xmm1
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm8
  jmp    L83
  nopw   0x0(%rax,%rax,1)
L140:
  ucomisd %xmm6,%xmm1
  jp     L236
  je     L83
L236:
  divsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm9
  jmp    L83
  xchg   %ax,%ax
L103:
  divsd  %xmm2,%xmm5
  divsd  %xmm2,%xmm5
  addsd  %xmm8,%xmm5
  ucomisd %xmm5,%xmm6
  ja     L190
  sqrtsd %xmm5,%xmm5
L17:
  movapd %xmm2,%xmm0
  add    $0x18,%rsp
  mulsd  %xmm5,%xmm0
  ret
  nopl   0x0(%rax,%rax,1)
L37:
  divsd  %xmm3,%xmm5
  addsd  %xmm9,%xmm5
  mulsd  %xmm3,%xmm5
  ucomisd %xmm5,%xmm6
  movapd %xmm5,%xmm0
  jbe    L186
L9:
  add    $0x18,%rsp
  jmp    1050 <sqrt@plt>
  nopl   0x0(%rax)
L156:
  ucomisd %xmm9,%xmm6
  ja     L226
  movapd %xmm9,%xmm0
  sqrtsd %xmm0,%xmm0
L238:
  mulsd  %xmm3,%xmm0
  add    $0x18,%rsp
  ret
  nopl   0x0(%rax)
L82:
  pxor   %xmm3,%xmm3
  movapd %xmm3,%xmm0
  sqrtsd %xmm0,%xmm0
  mulsd  %xmm3,%xmm0
  ret
L226:
  movapd %xmm9,%xmm0
  movsd  %xmm3,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm3
  jmp    L238
L190:
  movapd %xmm5,%xmm0
  movsd  %xmm2,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm2
  movapd %xmm0,%xmm5
  jmp    L17
```

tool_calls: name: disassemble_data
arguments: data_label: D29
datatype: byte[8]
name: disassemble_data
arguments: data_label: D2
datatype: qword
name: disassemble_data
arguments: data_label: D40
datatype: byte[8]
name: disassemble_data
arguments: data_label: D5
datatype: f64

tool: {"D29": [0xcc,0xcc,0xcc,0xcc,0xcc,0xcc,0xec,0x5f]}

tool: D2: 9223372036854775807

tool: {"D40": [0xcd,0xcc,0xcc,0xcc,0xcc,0xcc,0xfc,0x1f]}

tool: D5: 1

assistant: ```c
double enorm(int n, const double *x) {
  double ret_val, d1;
  int i;
  double s1, s2, s3, xabs, x1max, x3max, agiant, floatn;
  s1 = 0.;
  s2 = 0.;
  s3 = 0.;
  x1max = 0.;
  x3max = 0.;
  floatn = (double)(n);
  agiant = (1.3407807929942596e+154 * 0.9) / floatn;
  for (i = 0; i < n; ++i) {
    xabs = fabs(x[i]);
    if (xabs <= (1.4916681462400413e-154 * 0.9) || xabs >= agiant) {
      if (xabs > (1.4916681462400413e-154 * 0.9)) {
        if (xabs > x1max) {
          d1 = x1max / xabs;
          s1 = 1. + s1 * (d1 * d1);
          x1max = xabs;
        } else {
          d1 = xabs / x1max;
          s1 += d1 * d1;
        }
      } else {
        if (xabs > x3max) {
          d1 = x3max / xabs;
          s3 = 1. + s3 * (d1 * d1);
          x3max = xabs;
        } else {
          if (xabs != 0.) {
            d1 = xabs / x3max;
            s3 += d1 * d1;
          }
        }
      }
    } else {
      s2 += xabs * xabs;
    }
  }
  if (s1 != 0.) {
    ret_val = x1max * sqrt(s1 + (s2 / x1max) / x1max);
  } else {
    if (s2 != 0.) {
      if (s2 >= x3max) {
        ret_val = sqrt(s2 * (1. + (x3max / s2) * (x3max * s3)));
      } else {
        ret_val = sqrt(x3max * ((s2 / x3max) + (x3max * s3)));
      }
    } else {
      ret_val = x3max * sqrt(s3);
    }
  }
  return ret_val;
}
```