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

```asm
<main>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  add    $0xffffffffffffff80,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D108(%rip),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  D101(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  D149(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  D122(%rip),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  D18(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  D71(%rip),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  D116(%rip),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    D129(%rip),%rax
  movq   %rax,%xmm0
  lea    D83(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  lea    -0x78(%rbp),%rax
  mov    %rax,%rsi
  lea    D133(%rip),%rax
  mov    %rax,%rdi
  call   <scanf@plt>
  mov    %eax,-0x80(%rbp)
  cmpl   $0x1,-0x80(%rbp)
  je     L40
  lea    D48(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  mov    $0x1,%eax
  jmp    L48
L40:
  movsd  -0x20(%rbp),%xmm2
  movsd  -0x40(%rbp),%xmm3
  movsd  -0x78(%rbp),%xmm1
  movsd  D13(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  lea    -0x50(%rbp),%r8
  lea    -0x58(%rbp),%rcx
  lea    -0x60(%rbp),%rdx
  lea    -0x68(%rbp),%rsi
  lea    -0x70(%rbp),%rdi
  mov    $0x1,%r9d
  movapd %xmm3,%xmm1
  movq   %rax,%xmm0
  call   <exactP@plt>
  mov    %eax,-0x7c(%rbp)
  cmpl   $0x0,-0x7c(%rbp)
  je     L18
  lea    D45(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  mov    -0x7c(%rbp),%eax
  mov    %eax,%edi
  call   <error_message_testP@plt>
  mov    $0x1,%eax
  jmp    L48
L18:
  lea    D158(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  movsd  -0x50(%rbp),%xmm0
  movsd  D74(%rip),%xmm1
  movapd %xmm0,%xmm5
  divsd  %xmm1,%xmm5
  movsd  -0x58(%rbp),%xmm4
  movsd  -0x60(%rbp),%xmm3
  movsd  -0x68(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  -0x48(%rbp),%xmm2
  movsd  -0x70(%rbp),%xmm1
  movsd  D115(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D74(%rip),%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x70(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  lea    D14(%rip),%rax
  mov    %rax,%rdi
  mov    $0x6,%eax
  call   <printf@plt>
  mov    $0x0,%eax
L48:
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L7
  call   <__stack_chk_fail@plt>
L7:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D108
datatype: f64
name: retrieve_data
arguments: tag: D101
datatype: byte[8]
name: retrieve_data
arguments: tag: D149
datatype: byte[8]
name: retrieve_data
arguments: tag: D122
datatype: byte[8]
name: retrieve_data
arguments: tag: D18
datatype: byte[8]
name: retrieve_data
arguments: tag: D71
datatype: byte[8]
name: retrieve_data
arguments: tag: D116
datatype: byte[8]
name: retrieve_data
arguments: tag: D129
datatype: byte[8]
name: retrieve_data
arguments: tag: D83
datatype: string
name: retrieve_data
arguments: tag: D133
datatype: string
name: retrieve_data
arguments: tag: D48
datatype: string
name: retrieve_data
arguments: tag: D13
datatype: f64
name: retrieve_data
arguments: tag: D45
datatype: string
name: retrieve_data
arguments: tag: D158
datatype: string
name: retrieve_data
arguments: tag: D74
datatype: f64
name: retrieve_data
arguments: tag: D115
datatype: byte[8]
name: retrieve_data
arguments: tag: D14
datatype: string

tool: D108: 31556926

tool: {"D101": [0x11,0xea,0x2d,0x81,0x99,0x97,0x71,0x3d]}

tool: {"D149": [0x95,0xd6,0x26,0xe8,0x0b,0x2e,0x11,0x3e]}

tool: {"D122": [0x48,0xaf,0xbc,0x9a,0xf2,0xd7,0x7a,0x3e]}

tool: {"D18": [0x16,0x56,0xe7,0x9e,0xaf,0x03,0xd2,0x3c]}

tool: {"D71": [0x9b,0x2b,0xa1,0x86,0x9b,0x84,0x06,0x3d]}

tool: {"D116": [0x95,0x64,0x79,0xe1,0x7f,0xfd,0xa5,0x3d]}

tool: {"D129": [0x00,0x00,0x00,0x00,0x00,0x80,0x36,0x40]}

tool: D83: Enter  r  (in km; 0 <= r <= TESTP_L = %.1f):   

tool: D133: %lf

tool: D48: ... input error; exiting

tool: D13: 1000

tool: D45: 

simpleP ENDING because of ERROR from exactP():

tool: D158: Results from Test P:

tool: D74: 100000

tool: {"D115": [0xcd,0xcc,0xcc,0xcc,0x8c,0x6f,0xc1,0x40]}

tool: D14:     h = %.4f (m)  Po = %.7f (bar)  |vb| = %.5f (m/year)
    W_c = %.8f (m)  W = %.8f (m)  P = %.7f (bar)


assistant: ```c
int main() {
  double r, h, magvb, Wcrit, W, P;
  int scanret, ierr;
  const double secpera = 31556926.0;
  double EPS_ABS[] = {1.0e-12, 1.0e-9, 1.0e-7};
  double EPS_REL[] = {1.0e-15, 1.0e-14, 1.0e-11};
  printf("Enter  r  (in km; 0 <= r <= TESTP_L = %.1f):   ", 0.9 * 25000.0 / 1000.0);
  scanret = scanf("%lf", &r);
  if (scanret != 1) {
    printf("... input error; exiting\n");
    return 1;
  }
  ierr = exactP(r * 1000.0, &h, &magvb, &Wcrit, &W, &P, EPS_ABS[0], EPS_REL[0], 1);
  if (ierr) {
    printf("\n\nsimpleP ENDING because of ERROR from exactP():\n");
    error_message_testP(ierr);
    return 1;
  }
  printf("Results from Test P:\n");
  printf("    h = %.4f (m)  Po = %.7f (bar)  |vb| = %.5f (m/year)\n"
         "    W_c = %.8f (m)  W = %.8f (m)  P = %.7f (bar)\n",
         h, 910.0 * 9.81 * h / 1.0e5, magvb * secpera, Wcrit, W, P / 1.0e5);
  return 0;
}
```