user: 请将以下的汇编代码转换成C语言代码：

```asm
<least_square>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  mov    %rdx,%r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x68,%rsp
  mov    %rcx,-0x88(%rbp)
  mov    %edi,-0x4c(%rbp)
  mov    %rsp,%rcx
  mov    %esi,-0x50(%rbp)
  mov    %r8,-0x90(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  lea    0x1(%rsi),%eax
  mov    %eax,-0x70(%rbp)
  imul   %eax,%eax
  cltq
  lea    0xf(,%rax,8),%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rdx
  sub    %rdx,%rcx
  and    $0xfffffffffffffff0,%rax
  cmp    %rcx,%rsp
  je     L74
L93:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L93
L74:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  jne    L18
L34:
  movslq -0x70(%rbp),%rax
  mov    %rsp,%rcx
  mov    %rsp,-0x78(%rbp)
  lea    0xf(,%rax,8),%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rdx
  sub    %rdx,%rcx
  and    $0xfffffffffffffff0,%rax
  cmp    %rcx,%rsp
  je     L189
L11:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L11
L189:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  jne    L131
L168:
  movslq -0x4c(%rbp),%r13
  mov    %rsp,%rcx
  mov    %rsp,-0x80(%rbp)
  lea    0xf(,%r13,8),%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rdx
  sub    %rdx,%rcx
  and    $0xfffffffffffffff0,%rax
  cmp    %rcx,%rsp
  je     L79
L109:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rcx,%rsp
  jne    L109
L79:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  je     L118
  orq    $0x0,-0x8(%rsp,%rax,1)
L118:
  mov    -0x50(%rbp),%ecx
  mov    %rsp,%r14
  test   %ecx,%ecx
  js     L8
  movl   $0x0,-0x60(%rbp)
  movq   $0x0,-0x68(%rbp)
  nopw   0x0(%rax,%rax,1)
L185:
  mov    -0x68(%rbp),%rax
  mov    -0x78(%rbp),%rsi
  movq   $0x0,-0x48(%rbp)
  mov    %eax,-0x5c(%rbp)
  mov    %eax,-0x6c(%rbp)
  movslq -0x60(%rbp),%rax
  lea    (%rsi,%rax,8),%rax
  mov    %rax,-0x58(%rbp)
  nopw   0x0(%rax,%rax,1)
L47:
  mov    -0x4c(%rbp),%edx
  mov    -0x5c(%rbp),%ebx
  xor    %r12d,%r12d
  add    -0x48(%rbp),%ebx
  test   %edx,%edx
  jle    L175
  nopl   0x0(%rax,%rax,1)
L141:
  movsd  (%r15,%r12,8),%xmm0
  mov    %ebx,%edi
  call   <pow@plt>
  movsd  %xmm0,(%r14,%r12,8)
  add    $0x1,%r12
  cmp    %r13,%r12
  jne    L141
L175:
  mov    -0x4c(%rbp),%esi
  mov    %r14,%rdi
  call   <get_sum@plt>
  mov    -0x48(%rbp),%rax
  mov    -0x58(%rbp),%rcx
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  mov    %rax,-0x48(%rbp)
  cmp    %eax,-0x50(%rbp)
  jge    L47
  mov    -0x4c(%rbp),%eax
  mov    -0x6c(%rbp),%r12d
  xor    %ebx,%ebx
  test   %eax,%eax
  jle    L178
  xchg   %ax,%ax
L176:
  movsd  (%r15,%rbx,8),%xmm0
  mov    %r12d,%edi
  call   <pow@plt>
  movsd  %xmm0,(%r14,%rbx,8)
  add    $0x1,%rbx
  cmp    %rbx,%r13
  jne    L176
L178:
  mov    -0x4c(%rbp),%edx
  mov    -0x88(%rbp),%rdi
  mov    %r14,%rsi
  call   <get_crossed_sum@plt>
  mov    -0x68(%rbp),%rax
  mov    -0x80(%rbp),%rsi
  mov    -0x70(%rbp),%ecx
  add    %ecx,-0x60(%rbp)
  movsd  %xmm0,(%rsi,%rax,8)
  add    $0x1,%rax
  mov    %rax,-0x68(%rbp)
  cmp    %eax,-0x50(%rbp)
  jge    L185
L8:
  mov    -0x90(%rbp),%rcx
  mov    -0x80(%rbp),%rdx
  mov    -0x78(%rbp),%rsi
  mov    -0x70(%rbp),%edi
  call   <gauss_jordan_iteration@plt>
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L38
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L18:
  orq    $0x0,-0x8(%rsp,%rax,1)
  jmp    L34
L131:
  orq    $0x0,-0x8(%rsp,%rax,1)
  jmp    L168
L38:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void least_square(uint32_t plot_num, uint32_t order, double x[], double y[], double a[]) {
  double c[(order + 1) * (order + 1)];
  double d[order + 1];
  double temp[plot_num];
  uint32_t i = 0;
  uint32_t j = 0;
  uint32_t k = 0;
  for (k = 0; k <= order; k++) {
    for (j = 0; j <= order; j++) {
      for (i = 0; i < plot_num; i++)
        temp[i] = pow(x[i], k + j);
      c[k * (order + 1) + j] = get_sum(temp, plot_num);
    }
    for (i = 0; i < plot_num; i++)
      temp[i] = pow(x[i], k);
    d[k] = get_crossed_sum(y, temp, plot_num);
  }
  gauss_jordan_iteration(order + 1, c, d, a);
}
```