user: 这段汇编代码用C语言表示是什么样的?

```asm
<linetrend>:
  endbr64
  push   %r13
  mov    %edx,%r13d
  sub    $0x40,%rsp
  test   %edx,%edx
  jle    L98
  pxor   %xmm5,%xmm5
  lea    0x1(%rdx),%edx
  mov    $0x1,%eax
  movaps %xmm5,%xmm1
  movaps %xmm5,%xmm4
  movaps %xmm5,%xmm3
  movaps %xmm5,%xmm2
  nopl   (%rax)
L12:
  movss  (%rdi,%rax,4),%xmm0
  addss  %xmm0,%xmm2
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm4
  movss  (%rsi,%rax,4),%xmm0
  add    $0x1,%rax
  addss  %xmm0,%xmm3
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  cmp    %rax,%rdx
  jne    L12
L52:
  pxor   %xmm8,%xmm8
  movaps %xmm2,%xmm7
  movaps %xmm3,%xmm6
  cvtsi2ss %r13d,%xmm8
  movaps %xmm4,%xmm11
  divss  %xmm8,%xmm7
  divss  %xmm8,%xmm6
  movaps %xmm7,%xmm0
  mulss  %xmm7,%xmm0
  divss  %xmm8,%xmm11
  divss  %xmm8,%xmm1
  subss  %xmm0,%xmm11
  movaps %xmm6,%xmm0
  mulss  %xmm6,%xmm0
  comiss %xmm11,%xmm5
  subss  %xmm0,%xmm1
  jae    L113
  comiss %xmm1,%xmm5
  jb     L50
L113:
  movaps %xmm5,%xmm10
  movaps %xmm5,%xmm6
L1:
  movss  %xmm5,0x4(%rcx)
  movss  %xmm10,0x8(%rcx)
  movss  %xmm6,0xc(%rcx)
  add    $0x40,%rsp
  pop    %r13
  ret
  nopl   0x0(%rax,%rax,1)
L50:
  ucomiss %xmm11,%xmm5
  ja     L36
  sqrtss %xmm11,%xmm11
L18:
  ucomiss %xmm1,%xmm5
  ja     L55
  sqrtss %xmm1,%xmm1
L30:
  test   %r13d,%r13d
  jle    L114
  lea    0x1(%r13),%edx
  mov    $0x1,%eax
  movaps %xmm5,%xmm10
  nopw   0x0(%rax,%rax,1)
L94:
  movss  (%rdi,%rax,4),%xmm9
  movss  (%rsi,%rax,4),%xmm12
  add    $0x1,%rax
  movaps %xmm9,%xmm0
  movaps %xmm12,%xmm13
  subss  %xmm7,%xmm0
  subss  %xmm6,%xmm13
  mulss  %xmm12,%xmm9
  mulss  %xmm13,%xmm0
  addss  %xmm9,%xmm10
  addss  %xmm0,%xmm5
  cmp    %rax,%rdx
  jne    L94
L14:
  mulss  %xmm2,%xmm3
  sub    $0x1,%r13d
  pxor   %xmm0,%xmm0
  mulss  %xmm2,%xmm2
  cvtsi2ss %r13d,%xmm0
  mulss  %xmm1,%xmm11
  divss  %xmm8,%xmm3
  divss  %xmm8,%xmm2
  subss  %xmm3,%xmm10
  divss  %xmm0,%xmm5
  subss  %xmm2,%xmm4
  divss  %xmm4,%xmm10
  divss  %xmm11,%xmm5
  mulss  %xmm10,%xmm7
  subss  %xmm7,%xmm6
  jmp    L1
  nopw   0x0(%rax,%rax,1)
L98:
  pxor   %xmm5,%xmm5
  movaps %xmm5,%xmm1
  movaps %xmm5,%xmm4
  movaps %xmm5,%xmm3
  movaps %xmm5,%xmm2
  jmp    L52
L114:
  movaps %xmm5,%xmm10
  jmp    L14
L36:
  movaps %xmm11,%xmm0
  mov    %rcx,0x30(%rsp)
  mov    %rsi,0x28(%rsp)
  mov    %rdi,0x20(%rsp)
  movss  %xmm5,0x3c(%rsp)
  movss  %xmm3,0x38(%rsp)
  movss  %xmm4,0x1c(%rsp)
  movss  %xmm2,0x18(%rsp)
  movss  %xmm1,0x14(%rsp)
  movss  %xmm6,0x10(%rsp)
  movss  %xmm7,0xc(%rsp)
  movss  %xmm8,0x8(%rsp)
  call   <sqrtf@plt>
  movss  0x3c(%rsp),%xmm5
  mov    0x30(%rsp),%rcx
  mov    0x28(%rsp),%rsi
  mov    0x20(%rsp),%rdi
  movaps %xmm0,%xmm11
  movss  0x38(%rsp),%xmm3
  movss  0x1c(%rsp),%xmm4
  movss  0x18(%rsp),%xmm2
  movss  0x14(%rsp),%xmm1
  movss  0x10(%rsp),%xmm6
  movss  0xc(%rsp),%xmm7
  movss  0x8(%rsp),%xmm8
  jmp    L18
L55:
  movaps %xmm1,%xmm0
  mov    %rcx,0x30(%rsp)
  mov    %rsi,0x28(%rsp)
  mov    %rdi,0x20(%rsp)
  movss  %xmm5,0x3c(%rsp)
  movss  %xmm3,0x38(%rsp)
  movss  %xmm4,0x1c(%rsp)
  movss  %xmm2,0x18(%rsp)
  movss  %xmm11,0x14(%rsp)
  movss  %xmm6,0x10(%rsp)
  movss  %xmm7,0xc(%rsp)
  movss  %xmm8,0x8(%rsp)
  call   <sqrtf@plt>
  movss  0x3c(%rsp),%xmm5
  mov    0x30(%rsp),%rcx
  mov    0x28(%rsp),%rsi
  mov    0x20(%rsp),%rdi
  movaps %xmm0,%xmm1
  movss  0x38(%rsp),%xmm3
  movss  0x1c(%rsp),%xmm4
  movss  0x18(%rsp),%xmm2
  movss  0x14(%rsp),%xmm11
  movss  0x10(%rsp),%xmm6
  movss  0xc(%rsp),%xmm7
  movss  0x8(%rsp),%xmm8
  jmp    L30
```

assistant: ```c
void linetrend(float *x, float *y, int n, float *res) {
  float mean_x, mean_y, std_x, std_y;
  float sum_x, sum_y, sum_x2, sum_y2;
  float cov, sum_xy;
  int i;
  sum_x = 0;
  sum_y = 0;
  sum_x2 = 0;
  sum_y2 = 0;
  for (i = 1; i <= n; i++) {
    sum_x += x[i];
    sum_x2 += (x[i] * x[i]);
    sum_y += y[i];
    sum_y2 += (y[i] * y[i]);
  }
  mean_x = sum_x / n;
  mean_y = sum_y / n;
  std_x = sum_x2 / n - mean_x * mean_x;
  std_y = sum_y2 / n - mean_y * mean_y;
  if ((std_x <= 0) | (std_y <= 0)) {
    res[1] = 0;
    res[2] = 0;
    res[3] = 0;
  } else {
    std_x = sqrt(std_x);
    std_y = sqrt(std_y);
    cov = 0;
    sum_xy = 0;
    for (i = 1; i <= n; i++) {
      cov += (x[i] - mean_x) * (y[i] - mean_y);
      sum_xy += (x[i] * y[i]);
    }
    cov = cov / (n - 1);
    res[1] = cov / (std_x * std_y);
    res[2] = (sum_xy - sum_x * sum_y / n) / (sum_x2 - sum_x * sum_x / n);
    res[3] = mean_y - (res[2]) * mean_x;
  }
}
```