user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<test_spg_get_ir_reciprocal_mesh>:
  endbr64
  push   %rbp
  pxor   %xmm0,%xmm0
  mov    %rsp,%rbp
  push   %r13
  push   %r12
  sub    $0x130,%rsp
  movq   D165(%rip),%xmm1
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  movaps %xmm0,-0xf0(%rbp)
  mov    D7(%rip),%rax
  movaps %xmm0,-0xd0(%rbp)
  movaps %xmm0,-0xb0(%rbp)
  movapd D95(%rip),%xmm0
  mov    %rax,-0xc0(%rbp)
  movaps %xmm0,-0xa0(%rbp)
  movapd D51(%rip),%xmm0
  movaps %xmm1,-0x100(%rbp)
  movaps %xmm0,-0x90(%rbp)
  movapd D155(%rip),%xmm0
  movaps %xmm1,-0xe0(%rbp)
  movaps %xmm0,-0x80(%rbp)
  movapd D120(%rip),%xmm0
  movaps %xmm0,-0x70(%rbp)
  movq   D10(%rip),%xmm0
  movaps %xmm0,-0x60(%rbp)
  movapd D32(%rip),%xmm0
  movaps %xmm0,-0x50(%rbp)
  movapd D87(%rip),%xmm0
  movaps %xmm0,-0x40(%rbp)
  movapd D9(%rip),%xmm0
  movaps %xmm0,-0x30(%rbp)
  movdqa D47(%rip),%xmm0
  movaps %xmm0,-0x120(%rbp)
  mov    D88(%rip),%rax
  movl   $0x28,-0x138(%rbp)
  mov    %rax,-0x110(%rbp)
  mov    D111(%rip),%rax
  movl   $0x1,-0x128(%rbp)
  mov    %rax,-0x140(%rbp)
  mov    D47(%rip),%rax
  mov    %rax,-0x130(%rbp)
  lea    -0xbb000(%rsp),%rax
  cmp    %rax,%rsp
  je     L9
L11:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rax,%rsp
  jne    L11
L9:
  sub    $0x800,%rsp
  orq    $0x0,0x7f8(%rsp)
  lea    -0x3e000(%rsp),%rax
  mov    %rsp,%r12
  cmp    %rax,%rsp
  je     L14
L34:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rax,%rsp
  jne    L34
L14:
  sub    $0x800,%rsp
  orq    $0x0,0x7f8(%rsp)
  lea    D125(%rip),%rdi
  call   <puts@plt>
  lea    -0x120(%rbp),%rax
  mov    %rsp,%r13
  push   $0x0
  push   $0x6
  mov    %r13,%rsi
  mov    %r12,%rdi
  lea    -0x130(%rbp),%rcx
  push   %rax
  lea    -0xb0(%rbp),%rax
  lea    -0x140(%rbp),%rdx
  mov    $0x1,%r8d
  push   %rax
  lea    -0x100(%rbp),%r9
  call   <spg_get_ir_reciprocal_mesh@plt>
  lea    D115(%rip),%rdi
  mov    %r13,%rsp
  mov    %eax,%r12d
  call   <puts@plt>
  xor    %eax,%eax
  mov    %r12d,%esi
  lea    D17(%rip),%rdi
  call   <printf@plt>
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L30
  lea    -0x10(%rbp),%rsp
  pop    %r12
  pop    %r13
  pop    %rbp
  ret
L30:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: label: D165
type: byte[8]
name: disassemble_data
arguments: label: D7
type: byte[8]
name: disassemble_data
arguments: label: D95
type: byte[16]
name: disassemble_data
arguments: label: D51
type: f64[2]
name: disassemble_data
arguments: label: D155
type: f64[2]
name: disassemble_data
arguments: label: D120
type: byte[16]
name: disassemble_data
arguments: label: D10
type: f64
name: disassemble_data
arguments: label: D32
type: f64[2]
name: disassemble_data
arguments: label: D87
type: f64[2]
name: disassemble_data
arguments: label: D9
type: f64[2]
name: disassemble_data
arguments: label: D47
type: byte[16]
name: disassemble_data
arguments: label: D88
type: byte[8]
name: disassemble_data
arguments: label: D111
type: byte[8]
name: disassemble_data
arguments: label: D125
type: string
name: disassemble_data
arguments: label: D115
type: string
name: disassemble_data
arguments: label: D17
type: string

tool: {"D165": [0x00,0x00,0x00,0x00,0x00,0x00,0x10,0x40]}

tool: {"D7": [0x00,0x00,0x00,0x00,0x00,0x00,0x08,0x40]}

tool: {"D95": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

tool: {"D51": [0.5,0.5]}

tool: {"D155": [0.3,0.3]}

tool: {"D120": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x66,0x66,0x66,0x66,0x66,0x66,0xe6,0x3f]}

tool: D10: 0.7

tool: {"D32": [0.2,0.8]}

tool: {"D87": [0.5,0.8]}

tool: {"D9": [0.2,0.5]}

tool: {"D47": [0x01,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x02,0x00,0x00,0x00]}

tool: {"D88": [0x02,0x00,0x00,0x00,0x02,0x00,0x00,0x00]}

tool: {"D111": [0x28,0x00,0x00,0x00,0x28,0x00,0x00,0x00]}

tool: D125: *** Example of spg_get_ir_reciprocal_mesh of Rutile structure ***:

tool: D115: Number of irreducible k-points of Rutile with

tool: D17: 40x40x40 Monkhorst-Pack mesh is %d (4200).


assistant: ```c
void test_spg_get_ir_reciprocal_mesh(void) {
  double lattice[3][3] = {{4, 0, 0}, {0, 4, 0}, {0, 0, 3}};
  double position[][3] = {
      {0, 0, 0}, {0.5, 0.5, 0.5}, {0.3, 0.3, 0}, {0.7, 0.7, 0}, {0.2, 0.8, 0.5}, {0.8, 0.2, 0.5},
  };
  int types[] = {1, 1, 2, 2, 2, 2};
  int num_atom = 6;
  int m = 40;
  int mesh[] = {m, m, m};
  int is_shift[] = {1, 1, 1};
  int grid_address[m * m * m][3];
  int grid_mapping_table[m * m * m];
  printf("*** Example of spg_get_ir_reciprocal_mesh of Rutile structure ***:\n");
  int num_ir = spg_get_ir_reciprocal_mesh(grid_address, grid_mapping_table, mesh, is_shift, 1, lattice, position, types, num_atom, 1e-5);
  printf("Number of irreducible k-points of Rutile with\n");
  printf("40x40x40 Monkhorst-Pack mesh is %d (4200).\n", num_ir);
}
```