user: Please show me the C implementation that matches the assembly code listed below:

```asm
<test_spg_standardize_cell_corundum>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x3d0,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D34(%rip),%xmm0
  movsd  %xmm0,-0x3b0(%rbp)
  movsd  D1(%rip),%xmm0
  movsd  %xmm0,-0x3a8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x3a0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x398(%rbp)
  movsd  D29(%rip),%xmm0
  movsd  %xmm0,-0x390(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x388(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x380(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x378(%rbp)
  movsd  D11(%rip),%xmm0
  movsd  %xmm0,-0x370(%rbp)
  lea    -0x2e0(%rbp),%rax
  lea    D42(%rip),%rdx
  mov    $0x5a,%ecx
  mov    %rax,%rdi
  mov    %rdx,%rsi
  rep movsq %ds:(%rsi),%es:(%rdi)
  movl   $0x1e,-0x3bc(%rbp)
  movsd  D73(%rip),%xmm0
  movsd  %xmm0,-0x3b8(%rbp)
  movl   $0x0,-0x3c8(%rbp)
  jmp    L42
L6:
  mov    -0x3c8(%rbp),%eax
  cltq
  movl   $0x1,-0x360(%rbp,%rax,4)
  addl   $0x1,-0x3c8(%rbp)
L42:
  cmpl   $0xb,-0x3c8(%rbp)
  jle    L6
  movl   $0xc,-0x3c8(%rbp)
  jmp    L39
L89:
  mov    -0x3c8(%rbp),%eax
  cltq
  movl   $0x2,-0x360(%rbp,%rax,4)
  addl   $0x1,-0x3c8(%rbp)
L39:
  cmpl   $0x1d,-0x3c8(%rbp)
  jle    L89
  lea    D36(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  lea    D6(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  movl   $0x0,-0x3c4(%rbp)
  jmp    L50
L10:
  movl   $0x0,-0x3c0(%rbp)
  jmp    L28
L41:
  mov    -0x3c0(%rbp),%r9d
  mov    -0x3c4(%rbp),%r8d
  mov    -0x3b8(%rbp),%rdi
  mov    -0x3bc(%rbp),%ecx
  lea    -0x360(%rbp),%rdx
  lea    -0x2e0(%rbp),%rsi
  lea    -0x3b0(%rbp),%rax
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  call   <sub_spg_standardize_cell@plt>
  lea    D6(%rip),%rax
  mov    %rax,%rdi
  call   <printf@plt>
  addl   $0x1,-0x3c0(%rbp)
L28:
  cmpl   $0x1,-0x3c0(%rbp)
  jle    L41
  addl   $0x1,-0x3c4(%rbp)
L50:
  cmpl   $0x1,-0x3c4(%rbp)
  jle    L10
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L55
  call   <__stack_chk_fail@plt>
L55:
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_tag: D34
data_type: f64
name: parse_data
arguments: data_tag: D1
data_type: byte[8]
name: parse_data
arguments: data_tag: D29
data_type: f64
name: parse_data
arguments: data_tag: D11
data_type: f64
name: parse_data
arguments: data_tag: D42
data_type: f64
name: parse_data
arguments: data_tag: D73
data_type: byte[8]
name: parse_data
arguments: data_tag: D36
data_type: string
name: parse_data
arguments: data_tag: D6
data_type: string

tool: D34: 4.8076344022756095

tool: {"D1": [0xcf,0x97,0x43,0x83,0x04,0x3b,0x03,0xc0]}

tool: D29: 4.163533524478696

tool: D11: 13.117269919812754

tool: D42: 0

tool: {"D73": [0xf1,0x68,0xe3,0x88,0xb5,0xf8,0xe4,0x3e]}

tool: D36: *** Example of spg_standardize_cell (Corundum) ***:


tool: D6: ------------------------------------------------------


assistant: ```c
void test_spg_standardize_cell_corundum(void) {
  double lattice[3][3] = {{4.8076344022756095, -2.4038172011378047, 0}, {0, 4.1635335244786962, 0}, {0, 0, 13.1172699198127543}};
  double position[][3] = {
      {0.0000000000000000, 0.0000000000000000, 0.3521850942289043}, {0.6666666666666643, 0.3333333333333357, 0.6855184275622400}, {0.3333333333333357, 0.6666666666666643, 0.0188517608955686}, {0.0000000000000000, 0.0000000000000000, 0.6478149057711028}, {0.6666666666666643, 0.3333333333333357, 0.9811482391044314}, {0.3333333333333357, 0.6666666666666643, 0.3144815724377600}, {0.0000000000000000, 0.0000000000000000, 0.1478149057710957}, {0.6666666666666643, 0.3333333333333357, 0.4811482391044314}, {0.3333333333333357, 0.6666666666666643, 0.8144815724377600}, {0.0000000000000000, 0.0000000000000000, 0.8521850942288972}, {0.6666666666666643, 0.3333333333333357, 0.1855184275622400}, {0.3333333333333357, 0.6666666666666643, 0.5188517608955686}, {0.3061673906454899, 0.0000000000000000, 0.2500000000000000}, {0.9728340573121541, 0.3333333333333357, 0.5833333333333357}, {0.6395007239788255, 0.6666666666666643, 0.9166666666666643}, {0.6938326093545102, 0.0000000000000000, 0.7500000000000000}, {0.3604992760211744, 0.3333333333333357, 0.0833333333333357}, {0.0271659426878458, 0.6666666666666643, 0.4166666666666643}, {0.0000000000000000, 0.3061673906454899, 0.2500000000000000}, {0.6666666666666643, 0.6395007239788255, 0.5833333333333357}, {0.3333333333333357, 0.9728340573121541, 0.9166666666666643}, {0.0000000000000000, 0.6938326093545102, 0.7500000000000000}, {0.6666666666666643, 0.0271659426878458, 0.0833333333333357}, {0.3333333333333357, 0.3604992760211744, 0.4166666666666643}, {0.6938326093545102, 0.6938326093545102, 0.2500000000000000}, {0.3604992760211744, 0.0271659426878458, 0.5833333333333357}, {0.0271659426878458, 0.3604992760211744, 0.9166666666666643}, {0.3061673906454899, 0.3061673906454899, 0.7500000000000000}, {0.9728340573121541, 0.6395007239788255, 0.0833333333333357}, {0.6395007239788255, 0.9728340573121541, 0.4166666666666643},
  };
  int types[30];
  int i, j, k, num_atom = 30;
  double symprec = 1e-5;
  for (i = 0; i < 12; i++) {
    types[i] = 1;
  }
  for (i = 12; i < 30; i++) {
    types[i] = 2;
  }
  printf("*** Example of spg_standardize_cell (Corundum) ***:\n");
  printf("------------------------------------------------------\n");
  for (j = 0; j < 2; j++) {
    for (k = 0; k < 2; k++) {
      sub_spg_standardize_cell(lattice, position, types, num_atom, symprec, j, k);
      printf("------------------------------------------------------\n");
    }
  }
}
```