user: 请从给定的汇编指令推导出C源代码：

```asm
<main>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  lea    D65(%rip),%rbp
  push   %rbx
  sub    $0x118,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x108(%rsp)
  xor    %eax,%eax
  lea    0x40(%rsp),%r14
  lea    0x44(%rsp),%r15
  movl   $0x45c7511d,0x58(%rsp)
  movl   $0x1,0x5c(%rsp)
  lea    0x3c(%rsp),%r13
  lea    0x58(%rsp),%r12
  call   <setpar@plt>
  mov    %rbp,%rsi
  mov    %r14,%rdx
  lea    0x70(%rsp),%rbx
  lea    D39(%rip),%rdi
  call   <mstpar@plt>
  mov    %rbp,%rsi
  mov    %r15,%rdx
  lea    D132(%rip),%rdi
  call   <mstpar@plt>
  mov    %rbp,%rsi
  mov    %r13,%rdx
  lea    D73(%rip),%rdi
  call   <mstpar@plt>
  mov    %rbp,%rsi
  lea    0x34(%rsp),%rdx
  lea    D57(%rip),%rdi
  call   <mstpar@plt>
  mov    %rbp,%rsi
  lea    0x38(%rsp),%rdx
  lea    D14(%rip),%rdi
  call   <mstpar@plt>
  mov    %rbp,%rsi
  mov    %r12,%rdx
  xor    %eax,%eax
  lea    D67(%rip),%rdi
  lea    0xc0(%rsp),%rbp
  call   <getpar@plt>
  lea    0x5c(%rsp),%rdx
  lea    D86(%rip),%rsi
  xor    %eax,%eax
  lea    D51(%rip),%rdi
  call   <getpar@plt>
  xor    %eax,%eax
  call   <endpar@plt>
  mov    %r15,%r8
  mov    %r14,%rcx
  mov    %r13,%rdx
  mov    %rbp,%rsi
  mov    %rbx,%rdi
  call   <gen_matrices@plt>
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  movsd  D106(%rip),%xmm2
  cvtss2sd 0x38(%rsp),%xmm0
  cvtss2sd 0x58(%rsp),%xmm1
  mov    0x5c(%rsp),%eax
  mulsd  %xmm2,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x60(%rsp)
  pxor   %xmm0,%xmm0
  cvtss2sd 0x34(%rsp),%xmm0
  mulsd  %xmm2,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x68(%rsp)
  test   %eax,%eax
  je     L16
  mov    %rbx,0x10(%rsp)
  lea    D113(%rip),%rax
  mov    %r12,%rbx
  lea    0x54(%rsp),%r14
  mov    %rax,0x8(%rsp)
  lea    0x50(%rsp),%r13
  lea    D124(%rip),%r15
  jmp    L2
  nopl   0x0(%rax)
L64:
  lea    0x60(%rsp),%r9
  lea    0x4c(%rsp),%rcx
  mov    %rbx,%r8
  mov    %r14,%rsi
  lea    0x48(%rsp),%r12
  push   $0x0
  mov    %r13,%rdi
  push   %rbp
  mov    %r12,%rdx
  push   0x20(%rsp)
  lea    0x80(%rsp),%r10
  push   %r10
  mov    %r10,0x48(%rsp)
  mov    %r9,0x40(%rsp)
  mov    %rcx,0x38(%rsp)
  call   <gcproj@plt>
  movss  0x6c(%rsp),%xmm1
  add    $0x20,%rsp
  mov    %r15,%rdi
  movss  0x48(%rsp),%xmm0
  call   <printf@plt>
  push   $0x1
  mov    %r13,%rdi
  mov    %rbx,%r8
  push   %rbp
  mov    %r12,%rdx
  mov    %r14,%rsi
  push   0x20(%rsp)
  mov    0x40(%rsp),%r10
  push   %r10
  mov    0x40(%rsp),%r9
  mov    0x38(%rsp),%rcx
  call   <gcproj@plt>
  movss  0x74(%rsp),%xmm1
  add    $0x20,%rsp
  mov    %r15,%rdi
  movss  0x50(%rsp),%xmm0
  call   <printf@plt>
L2:
  mov    0x8(%rsp),%rdi
  mov    %r14,%rdx
  mov    %r13,%rsi
  call   <scanf@plt>
  cmp    $0x2,%eax
  je     L64
L41:
  mov    0x108(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L39
  add    $0x118,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L16:
  lea    D113(%rip),%rax
  lea    0x4c(%rsp),%r14
  mov    %rax,0x8(%rsp)
  lea    0x48(%rsp),%r13
  lea    D124(%rip),%r15
  jmp    L1
  nopw   0x0(%rax,%rax,1)
L5:
  lea    0x50(%rsp),%r10
  lea    0x60(%rsp),%r9
  mov    %r12,%r8
  mov    %r14,%rcx
  lea    0x54(%rsp),%rsi
  push   $0x1
  mov    %r10,%rdi
  mov    %r13,%rdx
  push   %rbp
  push   %rbx
  lea    0x80(%rsp),%r11
  push   %r11
  mov    %r10,0x30(%rsp)
  mov    %r11,0x48(%rsp)
  mov    %r9,0x40(%rsp)
  mov    %rsi,0x38(%rsp)
  call   <gcproj@plt>
  movss  0x74(%rsp),%xmm1
  add    $0x20,%rsp
  mov    %r15,%rdi
  movss  0x50(%rsp),%xmm0
  call   <printf@plt>
  push   $0x0
  mov    %r12,%r8
  mov    %r14,%rcx
  push   %rbp
  mov    %r13,%rdx
  push   %rbx
  mov    0x40(%rsp),%r11
  push   %r11
  mov    0x30(%rsp),%rdi
  mov    0x40(%rsp),%r9
  mov    0x38(%rsp),%rsi
  call   <gcproj@plt>
  movss  0x6c(%rsp),%xmm1
  add    $0x20,%rsp
  mov    %r15,%rdi
  movss  0x48(%rsp),%xmm0
  call   <printf@plt>
L1:
  mov    0x8(%rsp),%rdi
  mov    %r14,%rdx
  mov    %r13,%rsi
  call   <scanf@plt>
  cmp    $0x2,%eax
  je     L5
  jmp    L41
L39:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_label: D65
type_guessed: string
name: retrieve_data
arguments: data_label: D39
type_guessed: string
name: retrieve_data
arguments: data_label: D132
type_guessed: string
name: retrieve_data
arguments: data_label: D73
type_guessed: string
name: retrieve_data
arguments: data_label: D57
type_guessed: string
name: retrieve_data
arguments: data_label: D14
type_guessed: string
name: retrieve_data
arguments: data_label: D67
type_guessed: string
name: retrieve_data
arguments: data_label: D86
type_guessed: string
name: retrieve_data
arguments: data_label: D51
type_guessed: string
name: retrieve_data
arguments: data_label: D106
type_guessed: f64
name: retrieve_data
arguments: data_label: D113
type_guessed: string
name: retrieve_data
arguments: data_label: D124
type_guessed: string

tool: D65: f

tool: D39: ref_lon

tool: D132: ref_lat

tool: D73: alpha

tool: D57: xlen

tool: D14: ylen

tool: D67: ref_rad

tool: D86: d

tool: D51: input_xy

tool: D106: 0.5

tool: D113: %f %f

tool: D124: %.5f %.5f


assistant: ```c
main(int ac, char **av) {
  FILE *fpw, *fpr, *fopfile();
  float xlen, ylen, alpha, ref_lon, ref_lat, rlon, rlat;
  float dx, dy, xf, yf;
  double g0, b0;
  double amat[9], ainv[9];
  float ref_rad = 6378.139;
  int input_xy = 1;
  setpar(ac, av);
  mstpar("ref_lon", "f", &ref_lon);
  mstpar("ref_lat", "f", &ref_lat);
  mstpar("alpha", "f", &alpha);
  mstpar("xlen", "f", &xlen);
  mstpar("ylen", "f", &ylen);
  getpar("ref_rad", "f", &ref_rad);
  getpar("input_xy", "d", &input_xy);
  endpar();
  gen_matrices(amat, ainv, &alpha, &ref_lon, &ref_lat);
  g0 = (double)(0.5 * ylen) / (double)(ref_rad);
  b0 = (double)(0.5 * xlen) / (double)(ref_rad);
  if (input_xy) {
    while (scanf("%f %f", &xf, &yf) == 2) {
      gcproj(&xf, &yf, &rlon, &rlat, &ref_rad, &g0, &b0, amat, ainv, 0);
      printf("%.5f %.5f\n", rlon, rlat);
      gcproj(&xf, &yf, &rlon, &rlat, &ref_rad, &g0, &b0, amat, ainv, 1);
      printf("%.5f %.5f\n", xf, yf);
    }
  } else {
    while (scanf("%f %f", &rlon, &rlat) == 2) {
      gcproj(&xf, &yf, &rlon, &rlat, &ref_rad, &g0, &b0, amat, ainv, 1);
      printf("%.5f %.5f\n", xf, yf);
      gcproj(&xf, &yf, &rlon, &rlat, &ref_rad, &g0, &b0, amat, ainv, 0);
      printf("%.5f %.5f\n", rlon, rlat);
    }
  }
}
```