user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<main>:
  endbr64
  push   %r15
  mov    $0x10,%r8d
  mov    $0xe,%ecx
  mov    $0x2,%esi
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  mov    D66(%rip),%rax
  lea    0x20(%rsp),%r13
  mov    %r13,%rdx
  mov    %rax,0x18(%rsp)
  lea    0x18(%rsp),%rax
  mov    %rax,%rdi
  mov    %rax,0x8(%rsp)
  call   <create_omega@plt>
  mov    $0x4,%esi
  mov    $0x4,%edi
  call   <alloc_array_z@plt>
  mov    $0x4,%esi
  mov    $0x4,%edi
  mov    %rax,%rbp
  call   <alloc_array_z@plt>
  mov    $0x4,%edx
  mov    $0x4,%esi
  mov    $0x1,%edi
  call   <test_mat@plt>
  mov    $0x4,%esi
  mov    $0x4,%edi
  mov    %rax,%r14
  call   <alloc_array_z@plt>
  mov    $0x4,%esi
  mov    $0x4,%edi
  mov    %rax,%r12
  call   <alloc_array_z@plt>
  movl   $0x3e8,0x4(%rsp)
  mov    %rax,%rbx
  nopl   0x0(%rax)
L15:
  movsd  0x0(%rbp),%xmm0
  mov    $0x4,%edx
  mov    $0x4,%esi
  mov    %rbx,%rdi
  movsd  %xmm0,(%rbx)
  movsd  0x8(%rbp),%xmm0
  movsd  %xmm0,0x8(%rbx)
  movsd  0x10(%rbp),%xmm0
  movsd  %xmm0,0x10(%rbx)
  movsd  0x18(%rbp),%xmm0
  movsd  %xmm0,0x18(%rbx)
  movsd  0x20(%rbp),%xmm0
  movsd  %xmm0,0x20(%rbx)
  movsd  0x28(%rbp),%xmm0
  movsd  %xmm0,0x28(%rbx)
  movsd  0x30(%rbp),%xmm0
  movsd  %xmm0,0x30(%rbx)
  movsd  0x38(%rbp),%xmm0
  movsd  %xmm0,0x38(%rbx)
  movsd  0x40(%rbp),%xmm0
  movsd  %xmm0,0x40(%rbx)
  movsd  0x48(%rbp),%xmm0
  movsd  %xmm0,0x48(%rbx)
  movsd  0x50(%rbp),%xmm0
  movsd  %xmm0,0x50(%rbx)
  movsd  0x58(%rbp),%xmm0
  movsd  %xmm0,0x58(%rbx)
  movsd  0x60(%rbp),%xmm0
  movsd  %xmm0,0x60(%rbx)
  movsd  0x68(%rbp),%xmm0
  movsd  %xmm0,0x68(%rbx)
  movsd  0x70(%rbp),%xmm0
  movsd  %xmm0,0x70(%rbx)
  movsd  0x78(%rbp),%xmm0
  movsd  %xmm0,0x78(%rbx)
  mov    D43(%rip),%rax
  movq   %rax,%xmm0
  call   <shrink@plt>
  mov    %r12,%r8
  mov    $0xe,%ecx
  mov    %r13,%rdx
  mov    %rax,%rsi
  mov    %r14,%rdi
  mov    %rax,%r15
  call   <Proj_sub@plt>
  sub    $0x8,%rsp
  mov    %r12,%rcx
  mov    %rbp,%rdi
  mov    D68(%rip),%rax
  push   %rbp
  mov    $0x4,%edx
  mov    $0x4,%r9d
  push   $0xe
  mov    $0x4,%r8d
  mov    $0x4,%esi
  push   %r13
  movq   %rax,%xmm1
  mov    D26(%rip),%rax
  movq   %rax,%xmm0
  call   <ma@plt>
  add    $0x20,%rsp
  subl   $0x1,0x4(%rsp)
  jne    L15
  mov    0x8(%rsp),%rdx
  mov    $0x2,%ecx
  mov    %r15,%rsi
  mov    %r14,%rdi
  call   <RMSE2@plt>
  lea    D65(%rip),%rdi
  call   <printf@plt>
  movsd  (%r14),%xmm0
  lea    D30(%rip),%rdi
  call   <printf@plt>
  movsd  (%r15),%xmm0
  lea    D23(%rip),%rdi
  call   <printf@plt>
  mov    %rbp,%rdi
  call   <free_array@plt>
  mov    %r14,%rdi
  call   <free_array@plt>
  mov    %rbx,%rdi
  call   <free_array@plt>
  mov    %r15,%rdi
  call   <free_array@plt>
  mov    %r12,%rdi
  call   <free_array@plt>
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L14
  add    $0x68,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L14:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_label: D66
type: byte[8]
name: retrieve_data
arguments: data_label: D43
type: byte[8]
name: retrieve_data
arguments: data_label: D68
type: f64
name: retrieve_data
arguments: data_label: D26
type: f64
name: retrieve_data
arguments: data_label: D65
type: string
name: retrieve_data
arguments: data_label: D30
type: string
name: retrieve_data
arguments: data_label: D23
type: string

tool: {"D66": [0x00,0x00,0x00,0x00,0x05,0x00,0x00,0x00]}

tool: {"D43": [0x00,0x00,0x00,0x00,0x00,0x00,0x18,0x40]}

tool: D68: 1.4

tool: D26: 1

tool: D65: Error = %f



tool: D30: Known value = %f


tool: D23: Approximation = %f


assistant: ```c
int main(void) {
  int nrows = 4;
  int ncols = 4;
  int rank = 1;
  int numIter = 1000;
  int i, j;
  double tau = 1.5 * nrows;
  double delta = 1.4;
  int ku = 2;
  int omega_c[] = {0, 5};
  int kn = (nrows * ncols) - ku;
  int omega[kn];
  create_omega(omega_c, ku, omega, kn, nrows * ncols);
  double *Y = alloc_array_z(nrows, ncols);
  double *Z = alloc_array_z(nrows, ncols);
  double *M = test_mat(rank, nrows, ncols);
  double *dummyMatrix = alloc_array_z(nrows, ncols);
  double *dummy2 = alloc_array_z(nrows, ncols);
  for (i = 0; i < numIter; i++) {
    for (j = 0; j < ncols * nrows; j++) {
      dummy2[j] = Y[j];
    }
    Z = shrink(dummy2, tau, nrows, ncols);
    Proj_sub(M, Z, omega, kn, dummyMatrix);
    ma(Y, ncols, nrows, 1.0, dummyMatrix, ncols, nrows, delta, omega, kn, Y);
  }
  printf("Error = %f\n\n", RMSE2(M, Z, omega_c, ku));
  printf("Known value = %f\n", M[0]);
  printf("Approximation = %f\n", Z[0]);
  free_array(Y);
  free_array(M);
  free_array(dummy2);
  free_array(Z);
  free_array(dummyMatrix);
  return 0;
}
```