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

```asm
<Roe_Goundov_solver>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x100,%rsp
  mov    %rdi,%rbp
  movsd  %xmm0,0x28(%rsp)
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm5,(%rsp)
  movsd  %xmm6,0x20(%rsp)
  mov    %rsi,%r12
  mov    %fs:0x28,%rax
  mov    %rax,0xf8(%rsp)
  xor    %eax,%eax
  mulsd  %xmm1,%xmm0
  divsd  %xmm2,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L221
L168:
  movsd  0x28(%rsp),%xmm0
  mulsd  0x18(%rsp),%xmm0
  divsd  (%rsp),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L137
L6:
  movsd  D45(%rip),%xmm0
  movsd  0x28(%rsp),%xmm1
  movapd %xmm1,%xmm6
  subsd  %xmm0,%xmm6
  movsd  %xmm6,0x28(%rsp)
  divsd  %xmm6,%xmm1
  movsd  D61(%rip),%xmm4
  movsd  0x10(%rsp),%xmm5
  movapd %xmm5,%xmm2
  mulsd  %xmm5,%xmm2
  mulsd  %xmm4,%xmm2
  movsd  0x30(%rsp),%xmm9
  movapd %xmm9,%xmm3
  mulsd  %xmm1,%xmm3
  movsd  0x8(%rsp),%xmm7
  divsd  %xmm7,%xmm3
  addsd  %xmm2,%xmm3
  movsd  %xmm3,0x40(%rsp)
  mulsd  0x18(%rsp),%xmm1
  movsd  (%rsp),%xmm11
  divsd  %xmm11,%xmm1
  movsd  0x20(%rsp),%xmm10
  movapd %xmm10,%xmm3
  mulsd  %xmm10,%xmm3
  mulsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm1
  movsd  %xmm1,0x48(%rsp)
  movsd  %xmm7,0x50(%rsp)
  movapd %xmm7,%xmm1
  mulsd  %xmm5,%xmm1
  movsd  %xmm1,0x58(%rsp)
  divsd  %xmm6,%xmm0
  mulsd  %xmm9,%xmm0
  divsd  %xmm7,%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm7,%xmm0
  movsd  %xmm0,0x60(%rsp)
  movapd %xmm7,%xmm0
  mulsd  %xmm11,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L25
  sqrtsd %xmm0,%xmm0
  movq   %xmm0,%r13
L218:
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm5
  ucomisd %xmm5,%xmm0
  ja     L244
  sqrtsd %xmm5,%xmm5
  movapd %xmm5,%xmm0
L12:
  mulsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x38(%rsp)
  pxor   %xmm0,%xmm0
  movsd  (%rsp),%xmm6
  ucomisd %xmm6,%xmm0
  ja     L99
  sqrtsd %xmm6,%xmm6
  movapd %xmm6,%xmm0
L164:
  mulsd  0x20(%rsp),%xmm0
  addsd  0x38(%rsp),%xmm0
  movq   %xmm0,%rbx
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm3
  ucomisd %xmm3,%xmm0
  ja     L97
  sqrtsd %xmm3,%xmm3
  movsd  %xmm3,0x38(%rsp)
L214:
  pxor   %xmm0,%xmm0
  movsd  (%rsp),%xmm4
  ucomisd %xmm4,%xmm0
  ja     L234
  sqrtsd %xmm4,%xmm4
  movapd %xmm4,%xmm0
L21:
  addsd  0x38(%rsp),%xmm0
  movq   %rbx,%xmm4
  divsd  %xmm0,%xmm4
  movq   %xmm4,%rbx
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm5
  ucomisd %xmm5,%xmm0
  ja     L258
  sqrtsd %xmm5,%xmm5
  movapd %xmm5,%xmm0
L205:
  mulsd  0x40(%rsp),%xmm0
  movsd  %xmm0,0x38(%rsp)
  pxor   %xmm0,%xmm0
  movsd  (%rsp),%xmm7
  ucomisd %xmm7,%xmm0
  ja     L193
  sqrtsd %xmm7,%xmm7
  movapd %xmm7,%xmm0
L106:
  mulsd  0x48(%rsp),%xmm0
  addsd  0x38(%rsp),%xmm0
  movq   %xmm0,%r14
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm6
  ucomisd %xmm6,%xmm0
  ja     L91
  sqrtsd %xmm6,%xmm6
  movsd  %xmm6,0x38(%rsp)
L147:
  pxor   %xmm0,%xmm0
  movsd  (%rsp),%xmm3
  ucomisd %xmm3,%xmm0
  ja     L188
  sqrtsd %xmm3,%xmm3
  movapd %xmm3,%xmm0
L57:
  addsd  0x38(%rsp),%xmm0
  movq   %r14,%xmm4
  divsd  %xmm0,%xmm4
  movsd  %xmm4,0x38(%rsp)
  movq   %rbx,%xmm1
  mulsd  D61(%rip),%xmm1
  movq   %rbx,%xmm5
  mulsd  %xmm5,%xmm1
  subsd  %xmm1,%xmm4
  movapd %xmm4,%xmm0
  mulsd  0x28(%rsp),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L199
  sqrtsd %xmm0,%xmm0
L265:
  movsd  D45(%rip),%xmm1
  movsd  %xmm1,0xb0(%rsp)
  movsd  %xmm1,0xb8(%rsp)
  movsd  %xmm1,0xc0(%rsp)
  movq   %rbx,%xmm7
  subsd  %xmm0,%xmm7
  movsd  %xmm7,0xc8(%rsp)
  mov    %rbx,0xd0(%rsp)
  movq   %rbx,%xmm6
  addsd  %xmm0,%xmm6
  movsd  %xmm6,0xd8(%rsp)
  movq   %rbx,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  0x38(%rsp),%xmm4
  movapd %xmm4,%xmm2
  subsd  %xmm1,%xmm2
  movsd  %xmm2,0xe0(%rsp)
  movq   %rbx,%xmm2
  mulsd  %xmm2,%xmm2
  movsd  D61(%rip),%xmm5
  mulsd  %xmm5,%xmm2
  movsd  %xmm2,0xe8(%rsp)
  addsd  %xmm4,%xmm1
  movsd  %xmm1,0xf0(%rsp)
  movsd  0x18(%rsp),%xmm2
  subsd  0x30(%rsp),%xmm2
  movq   %r13,%xmm3
  mulsd  %xmm0,%xmm3
  movsd  0x20(%rsp),%xmm1
  subsd  0x10(%rsp),%xmm1
  mulsd  %xmm1,%xmm3
  movapd %xmm0,%xmm4
  mulsd  %xmm0,%xmm4
  movapd %xmm2,%xmm1
  subsd  %xmm3,%xmm1
  mulsd  %xmm5,%xmm1
  divsd  %xmm4,%xmm1
  movsd  %xmm1,0x90(%rsp)
  movsd  (%rsp),%xmm1
  subsd  0x8(%rsp),%xmm1
  movapd %xmm2,%xmm8
  divsd  %xmm4,%xmm8
  subsd  %xmm8,%xmm1
  movsd  %xmm1,0x98(%rsp)
  addsd  %xmm3,%xmm2
  mulsd  %xmm5,%xmm2
  divsd  %xmm4,%xmm2
  movsd  %xmm2,0xa0(%rsp)
  movsd  %xmm7,0x70(%rsp)
  mov    %rbx,0x78(%rsp)
  movsd  %xmm6,0x80(%rsp)
  movq   %rbx,%xmm1
  andpd  D37(%rip),%xmm1
  addsd  %xmm0,%xmm1
  movsd  %xmm1,(%r12)
  mov    $0x0,%eax
  lea    0x70(%rsp),%rsi
  pxor   %xmm1,%xmm1
  lea    0xb0(%rsp),%r8
  lea    0x68(%rsp),%rdi
L242:
  comisd (%rax,%rsi,1),%xmm1
  jae    L165
L246:
  add    $0x8,%rax
  cmp    $0x18,%rax
  jne    L242
  movsd  0x50(%rsp),%xmm2
  movsd  %xmm2,0x0(%rbp)
  movsd  0x58(%rsp),%xmm1
  divsd  %xmm2,%xmm1
  movsd  %xmm1,0x8(%rbp)
  movsd  0x60(%rsp),%xmm0
  divsd  %xmm2,%xmm0
  movapd %xmm1,%xmm3
  mulsd  D61(%rip),%xmm3
  mulsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  mulsd  0x28(%rsp),%xmm0
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,0x10(%rbp)
  mov    0xf8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L54
  add    $0x100,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L221:
  call   <sqrt@plt>
  jmp    L168
L137:
  call   <sqrt@plt>
  jmp    L6
L25:
  call   <sqrt@plt>
  movq   %xmm0,%r13
  jmp    L218
L244:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L12
L99:
  movsd  (%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L164
L97:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  %xmm0,0x38(%rsp)
  jmp    L214
L234:
  movsd  (%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L21
L258:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L205
L193:
  movsd  (%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L106
L91:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  %xmm0,0x38(%rsp)
  jmp    L147
L188:
  movsd  (%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L57
L199:
  call   <sqrt@plt>
  jmp    L265
L165:
  movsd  0x90(%rsp,%rax,1),%xmm2
  lea    0x50(%rsp),%rdx
  lea    (%r8,%rax,1),%rcx
L127:
  movapd %xmm2,%xmm0
  mulsd  (%rcx),%xmm0
  addsd  (%rdx),%xmm0
  movsd  %xmm0,(%rdx)
  add    $0x8,%rdx
  add    $0x18,%rcx
  cmp    %rdx,%rdi
  jne    L127
  jmp    L246
L54:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D45
datatype: f64
name: read_data
arguments: data_label: D61
datatype: f64
name: read_data
arguments: data_label: D37
datatype: qword[2]

tool: D45: 1

tool: D61: 0.5

tool: {"D37": [0x7fffffffffffffff,0x0000000000000000]}

assistant: ```c
void Roe_Goundov_solver(double *F, double gamma, double P_L, double RHO_L, double U_L, double P_R, double RHO_R, double U_R, double *lambda_max, double delta) {
  double C_L, C_R;
  C_L = sqrt(gamma * P_L / RHO_L);
  C_R = sqrt(gamma * P_R / RHO_R);
  double z = 0.5 * (gamma - 1.0) / gamma;
  double H_L, H_R;
  H_L = gamma / (gamma - 1.0) * P_L / RHO_L + 0.5 * (U_L * U_L);
  H_R = gamma / (gamma - 1.0) * P_R / RHO_R + 0.5 * (U_R * U_R);
  double E_L, E_R;
  E_L = 1.0 / (gamma - 1.0) * P_L / RHO_L + 0.5 * (U_L * U_L);
  E_R = 1.0 / (gamma - 1.0) * P_R / RHO_R + 0.5 * (U_R * U_R);
  double U[3];
  U[0] = RHO_L;
  U[1] = RHO_L * U_L;
  U[2] = RHO_L * E_L;
  double RHO_S, U_S, H_S, C_S;
  RHO_S = sqrt(RHO_L * RHO_R);
  U_S = (U_L * sqrt(RHO_L) + U_R * sqrt(RHO_R)) / (sqrt(RHO_L) + sqrt(RHO_R));
  H_S = (H_L * sqrt(RHO_L) + H_R * sqrt(RHO_R)) / (sqrt(RHO_L) + sqrt(RHO_R));
  C_S = sqrt((gamma - 1.0) * (H_S - 0.5 * U_S * U_S));
  double R[3][3];
  double lambda[3], W[3];
  R[0][0] = 1.0;
  R[0][1] = 1.0;
  R[0][2] = 1.0;
  R[1][0] = U_S - C_S;
  R[1][1] = U_S;
  R[1][2] = U_S + C_S;
  R[2][0] = H_S - U_S * C_S;
  R[2][1] = 0.5 * (U_S * U_S);
  R[2][2] = H_S + U_S * C_S;
  int i, j;
  W[0] = 0.5 * ((P_R - P_L) - RHO_S * C_S * (U_R - U_L)) / (C_S * C_S);
  W[1] = (RHO_R - RHO_L) - (P_R - P_L) / (C_S * C_S);
  W[2] = 0.5 * ((P_R - P_L) + RHO_S * C_S * (U_R - U_L)) / (C_S * C_S);
  lambda[0] = U_S - C_S;
  lambda[1] = U_S;
  lambda[2] = U_S + C_S;
  *lambda_max = fabs(U_S) + C_S;
  for (j = 0; j < 3; j++) {
    if (lambda[j] <= 0) {
      for (i = 0; i < 3; i++) {
        U[i] += W[j] * R[i][j];
      }
    }
  }
  F[0] = U[0];
  F[1] = U[1] / U[0];
  F[2] = (U[2] / U[0] - 0.5 * F[1] * F[1]) * (gamma - 1.0) * F[0];
}
```