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

```asm
<linear_GRP_solver_LAG>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x140,%rsp
  mov    %rdi,-0xe8(%rbp)
  mov    %rsi,-0xf0(%rbp)
  movsd  %xmm0,-0xf8(%rbp)
  movsd  %xmm1,-0x100(%rbp)
  movsd  %xmm2,-0x108(%rbp)
  movsd  %xmm3,-0x110(%rbp)
  movsd  %xmm4,-0x118(%rbp)
  movsd  %xmm5,-0x120(%rbp)
  movsd  %xmm6,-0x128(%rbp)
  movsd  %xmm7,-0x130(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  0x30(%rbp),%xmm0
  movsd  D69(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  0x30(%rbp),%xmm2
  movsd  D69(%rip),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  0x30(%rbp),%xmm0
  mulsd  0x10(%rbp),%xmm0
  divsd  -0xf8(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x80(%rbp)
  movsd  0x30(%rbp),%xmm0
  mulsd  0x18(%rbp),%xmm0
  divsd  -0x100(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x78(%rbp)
  movsd  -0xf8(%rbp),%xmm0
  mulsd  -0x80(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x100(%rbp),%xmm0
  mulsd  -0x78(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  0x38(%rbp),%xmm5
  lea    -0x10(%rbp),%rdx
  movsd  -0x78(%rbp),%xmm4
  movsd  -0x80(%rbp),%xmm3
  movsd  0x18(%rbp),%xmm2
  movsd  -0x120(%rbp),%xmm1
  movsd  -0x118(%rbp),%xmm0
  mov    0x30(%rbp),%rdi
  lea    -0xd0(%rbp),%rsi
  lea    -0xd8(%rbp),%rax
  mov    $0x1f4,%ecx
  movapd %xmm5,%xmm7
  movapd %xmm4,%xmm6
  movapd %xmm3,%xmm5
  movapd %xmm2,%xmm4
  movsd  0x10(%rbp),%xmm3
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  call   <Riemann_solver_exact@plt>
  movsd  -0xd0(%rbp),%xmm0
  comisd 0x10(%rbp),%xmm0
  jbe    L59
  movsd  -0x88(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  0x10(%rbp),%xmm1
  movsd  -0xd0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0xf8(%rbp),%xmm0
  movsd  -0xd0(%rbp),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  addsd  0x10(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  jmp    L12
L59:
  movsd  D69(%rip),%xmm0
  movapd %xmm0,%xmm1
  divsd  0x30(%rbp),%xmm1
  movsd  -0xd0(%rbp),%xmm0
  divsd  0x10(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <pow@plt>
  movsd  -0xf8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
L12:
  movsd  -0xd0(%rbp),%xmm0
  comisd 0x18(%rbp),%xmm0
  jbe    L46
  movsd  -0x88(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  0x18(%rbp),%xmm1
  movsd  -0xd0(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  -0xd0(%rbp),%xmm1
  mulsd  -0x88(%rbp),%xmm1
  addsd  0x18(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  jmp    L23
L46:
  movsd  D69(%rip),%xmm0
  movapd %xmm0,%xmm1
  divsd  0x30(%rbp),%xmm1
  movsd  -0xd0(%rbp),%xmm0
  divsd  0x18(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <pow@plt>
  movsd  -0x100(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xc0(%rbp)
L23:
  mov    -0xf0(%rbp),%rax
  add    $0x8,%rax
  movsd  -0xd8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xf0(%rbp),%rax
  add    $0x10,%rax
  movsd  -0xd0(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xf0(%rbp),%rax
  movsd  -0xc8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xf0(%rbp),%rax
  add    $0x18,%rax
  movsd  -0xc0(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0xd0(%rbp),%xmm0
  mulsd  0x30(%rbp),%xmm0
  divsd  -0xc8(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x60(%rbp)
  movsd  -0xd0(%rbp),%xmm0
  mulsd  0x30(%rbp),%xmm0
  divsd  -0xc0(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x58(%rbp)
  movsd  -0xc0(%rbp),%xmm0
  mulsd  -0x58(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0xc8(%rbp),%xmm0
  mulsd  -0x60(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x48(%rbp),%xmm0
  divsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x118(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x120(%rbp),%xmm1
  movsd  -0x118(%rbp),%xmm0
  subsd  -0x120(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  0x10(%rbp),%xmm0
  movapd %xmm0,%xmm2
  subsd  0x18(%rbp),%xmm2
  movsd  0x10(%rbp),%xmm0
  subsd  0x18(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x38(%rbp)
  movsd  0x38(%rbp),%xmm1
  movsd  D331(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  comisd -0x38(%rbp),%xmm0
  jbe    L3
  movsd  D69(%rip),%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movsd  D69(%rip),%xmm0
  divsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  movsd  D82(%rip),%xmm0
  divsd  -0x48(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x70(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  mulsd  -0x128(%rbp),%xmm0
  addsd  0x20(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  D82(%rip),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  D69(%rip),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  D82(%rip),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x68(%rbp),%xmm1
  movsd  -0x68(%rbp),%xmm0
  mulsd  -0x130(%rbp),%xmm0
  subsd  0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  jmp    L63
L3:
  movsd  D69(%rip),%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movsd  D69(%rip),%xmm0
  divsd  -0x48(%rbp),%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  movsd  0x20(%rbp),%xmm0
  divsd  -0x70(%rbp),%xmm0
  addsd  -0x128(%rbp),%xmm0
  movsd  %xmm0,-0x138(%rbp)
  movsd  D69(%rip),%xmm0
  divsd  -0x70(%rbp),%xmm0
  movsd  0x30(%rbp),%xmm2
  movsd  D216(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  movsd  D69(%rip),%xmm3
  movapd %xmm1,%xmm2
  subsd  %xmm3,%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  -0x80(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  -0x108(%rbp),%xmm0
  subsd  0x20(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,-0x140(%rbp)
  movsd  0x30(%rbp),%xmm1
  movsd  D216(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D69(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  D35(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  0x30(%rbp),%xmm2
  movsd  D69(%rip),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x40(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  movsd  D69(%rip),%xmm1
  movq   %rax,%xmm7
  subsd  %xmm1,%xmm7
  movapd %xmm7,%xmm0
  mulsd  -0x140(%rbp),%xmm0
  addsd  -0x138(%rbp),%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0x70(%rbp),%xmm0
  mulsd  -0x48(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  movq   D263(%rip),%xmm0
  movq   %rax,%xmm5
  xorpd  %xmm0,%xmm5
  movapd %xmm5,%xmm0
  movsd  -0xa8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0xd0(%rbp),%xmm0
  subsd  0x18(%rbp),%xmm0
  movsd  -0xd8(%rbp),%xmm1
  subsd  -0x120(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0xd0(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  0x18(%rbp),%xmm1
  movsd  D239(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x88(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  0x18(%rbp),%xmm2
  movsd  -0xd0(%rbp),%xmm1
  addsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D35(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movq   D263(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  divsd  -0x50(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  movsd  -0xa0(%rbp),%xmm1
  movsd  D69(%rip),%xmm2
  subsd  %xmm2,%xmm1
  divsd  -0x30(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  movsd  -0xd0(%rbp),%xmm0
  subsd  0x18(%rbp),%xmm0
  movq   D263(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  D35(%rip),%xmm1
  divsd  %xmm1,%xmm0
  divsd  -0x100(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0xd8(%rbp),%xmm0
  subsd  -0x120(%rbp),%xmm0
  mulsd  -0x100(%rbp),%xmm0
  movsd  D35(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movapd %xmm0,%xmm2
  addsd  -0x30(%rbp),%xmm2
  movsd  -0x68(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movapd %xmm0,%xmm3
  divsd  -0x30(%rbp),%xmm3
  movsd  -0xd0(%rbp),%xmm0
  subsd  0x18(%rbp),%xmm0
  mulsd  -0x88(%rbp),%xmm0
  movsd  D35(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  -0x88(%rbp),%xmm1
  movapd %xmm1,%xmm4
  mulsd  0x18(%rbp),%xmm4
  movsd  -0xd0(%rbp),%xmm1
  addsd  %xmm1,%xmm4
  movapd %xmm0,%xmm1
  divsd  %xmm4,%xmm1
  movsd  D69(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x88(%rbp),%xmm0
  movsd  D35(%rip),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  -0xd0(%rbp),%xmm0
  subsd  0x18(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x88(%rbp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  0x18(%rbp),%xmm2
  movsd  -0xd0(%rbp),%xmm1
  addsd  %xmm1,%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D35(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movq   D263(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x20(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x130(%rbp),%xmm1
  movsd  -0x18(%rbp),%xmm0
  mulsd  0x28(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x110(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x90(%rbp)
L63:
  movsd  -0xa8(%rbp),%xmm0
  mulsd  -0x98(%rbp),%xmm0
  movsd  -0x90(%rbp),%xmm1
  mulsd  -0xb0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0xb8(%rbp),%xmm1
  mulsd  -0x98(%rbp),%xmm1
  movsd  -0xa0(%rbp),%xmm2
  mulsd  -0xb0(%rbp),%xmm2
  subsd  %xmm2,%xmm1
  mov    -0xe8(%rbp),%rax
  add    $0x8,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0xa8(%rbp),%xmm0
  mulsd  -0xa0(%rbp),%xmm0
  movsd  -0x90(%rbp),%xmm1
  mulsd  -0xb8(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  -0xb0(%rbp),%xmm1
  mulsd  -0xa0(%rbp),%xmm1
  movsd  -0x98(%rbp),%xmm2
  mulsd  -0xb8(%rbp),%xmm2
  subsd  %xmm2,%xmm1
  mov    -0xe8(%rbp),%rax
  add    $0x10,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  D69(%rip),%xmm0
  divsd  -0x60(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x60(%rbp),%xmm1
  mov    -0xe8(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0xe8(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  D69(%rip),%xmm0
  divsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x58(%rbp),%xmm1
  mov    -0xe8(%rbp),%rax
  add    $0x10,%rax
  movsd  (%rax),%xmm0
  mov    -0xe8(%rbp),%rax
  add    $0x18,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L54
  call   <__stack_chk_fail@plt>
L54:
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D69
type: f64
name: extract_symbolic_data
arguments: data_tag: D331
type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D82
type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D216
type: f64
name: extract_symbolic_data
arguments: data_tag: D35
type: f64
name: extract_symbolic_data
arguments: data_tag: D263
type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D239
type: f64

tool: D69: 1

tool: {"D331": [0x00,0x00,0x00,0x00,0x00,0x40,0x8f,0x40]}

tool: {"D82": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: D216: 3

tool: D35: 2

tool: {"D263": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

tool: D239: 0.5

assistant: ```c
void linear_GRP_solver_LAG(double *direvative, double *mid, double rho_L, double rho_R, double s_rho_L, double s_rho_R, double u_L, double u_R, double s_u_L, double s_u_R, double p_L, double p_R, double s_p_L, double s_p_R, double gamma, double eps) {
  double dist;
  int CRW[2];
  double c_L, c_R, g_L, g_R;
  double rho_star, u_star, p_star, rho_star_L, rho_star_R, c_star_L, c_star_R, g_star_L, g_star_R, beta_star;
  double a_L, b_L, d_L, a_R, b_R, d_R, L_rho, L_u, L_p, sigma;
  double zeta = (gamma - 1.0) / (gamma + 1.0);
  c_L = sqrt(gamma * p_L / rho_L);
  c_R = sqrt(gamma * p_R / rho_R);
  g_L = rho_L * c_L;
  g_R = rho_R * c_R;
  Riemann_solver_exact(&u_star, &p_star, gamma, u_L, u_R, p_L, p_R, c_L, c_R, CRW, eps, 500);
  if (p_star > p_L)
    rho_star_L = rho_L * (p_star + zeta * p_L) / (p_L + zeta * p_star);
  else
    rho_star_L = rho_L * pow(p_star / p_L, 1.0 / gamma);
  if (p_star > p_R)
    rho_star_R = rho_R * (p_star + zeta * p_R) / (p_R + zeta * p_star);
  else
    rho_star_R = rho_R * pow(p_star / p_R, 1.0 / gamma);
  mid[1] = u_star;
  mid[2] = p_star;
  mid[0] = rho_star_L;
  mid[3] = rho_star_R;
  c_star_L = sqrt(gamma * p_star / rho_star_L);
  c_star_R = sqrt(gamma * p_star / rho_star_R);
  g_star_R = rho_star_R * c_star_R;
  g_star_L = rho_star_L * c_star_L;
  beta_star = g_star_L / g_L;
  dist = sqrt((u_L - u_R) * (u_L - u_R) + (p_L - p_R) * (p_L - p_R));
  if ((dist < eps * 1000)) {
    a_L = 1.0;
    b_L = 1.0 / g_star_L;
    d_L = -1.0 / g_star_L * g_L * (g_L * s_u_L + s_p_L);
    a_R = -1.0;
    b_R = 1.0 / g_star_R;
    d_R = -1.0 / g_star_R * g_R * (g_R * s_u_R - s_p_R);
  } else {
    a_L = 1.0;
    b_L = 1.0 / g_star_L;
    d_L = (s_u_L + s_p_L / g_L) + 1.0 / g_L / (3.0 * gamma - 1.0) * (c_L * c_L * s_rho_L - s_p_L) * (pow(beta_star, (3.0 * gamma - 1.0) / 2.0 / (gamma + 1.0)) - 1.0);
    d_L = (-1.0) * sqrt(g_L * g_star_L) * d_L;
    sigma = (p_star - p_R) / (u_star - u_R);
    a_R = 2.0 - 0.5 * (p_star - p_R) / (p_star + zeta * p_R);
    b_R = -1.0 * sigma / g_star_R / g_star_R - (a_R - 1) / sigma;
    L_rho = -1.0 * (p_star - p_R) / 2.0 / rho_R;
    L_u = sigma + rho_R * (u_star - u_R) / 2.0 + g_R * g_R / sigma * (1.0 + zeta * (p_star - p_R) / 2.0 / (p_star + zeta * p_R));
    L_p = -1.0 * (2 + zeta / 2.0 * (p_star - p_R) / (p_star + zeta * p_R));
    d_R = L_u * s_u_R + L_p * s_p_R + L_rho * s_rho_R;
  }
  direvative[1] = (d_L * b_R - d_R * b_L) / (a_L * b_R - a_R * b_L);
  direvative[2] = (d_L * a_R - d_R * a_L) / (b_L * a_R - b_R * a_L);
  direvative[0] = 1.0 / c_star_L / c_star_L * direvative[2];
  direvative[3] = 1.0 / c_star_R / c_star_R * direvative[2];
}
```