user: What is the c source code of the assembly code below?

```asm
<ELossCL>:
  endbr64
  push   %r15
  pxor   %xmm3,%xmm3
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  comisd %xmm1,%xmm3
  movsd  %xmm0,0x18(%rsp)
  movsd  %xmm2,0x38(%rsp)
  ja     L65
  comisd %xmm2,%xmm3
  jae    L40
  movapd %xmm1,%xmm4
  movapd %xmm2,%xmm1
  mov    %edi,%ebx
  mov    %rsi,%rbp
  movsd  %xmm4,0x8(%rsp)
  mov    %rdx,%r12
  mov    %rcx,%r13
  mov    %r8d,%r15d
  mov    %r9d,%r14d
  call   <Range@plt>
  movsd  0x8(%rsp),%xmm4
  pxor   %xmm3,%xmm3
  movapd %xmm0,%xmm2
  subsd  %xmm4,%xmm2
  comisd %xmm2,%xmm3
  movsd  %xmm2,0x8(%rsp)
  jae    L40
  movsd  0x38(%rsp),%xmm7
  mov    %r14d,%r9d
  mov    %r15d,%r8d
  mov    %r13,%rcx
  movsd  D26(%rip),%xmm0
  mov    %r12,%rdx
  mov    %rbp,%rsi
  mov    %ebx,%edi
  movapd %xmm7,%xmm1
  divsd  D40(%rip),%xmm4
  movsd  %xmm4,0x20(%rsp)
  mulsd  %xmm7,%xmm0
  subsd  %xmm0,%xmm1
  movsd  0x18(%rsp),%xmm0
  movsd  %xmm1,0x10(%rsp)
  call   <Range@plt>
  movsd  0x8(%rsp),%xmm2
  pxor   %xmm3,%xmm3
  movsd  0x38(%rsp),%xmm7
  movsd  0x10(%rsp),%xmm1
  movsd  %xmm3,0x30(%rsp)
  subsd  %xmm2,%xmm0
  andpd  D48(%rip),%xmm0
  movsd  %xmm7,0x28(%rsp)
  comisd 0x20(%rsp),%xmm0
  jae    L43
  jmp    L83
  nopl   0x0(%rax)
L44:
  movapd %xmm1,%xmm3
  subsd  0x30(%rsp),%xmm3
  movsd  %xmm1,0x28(%rsp)
  mulsd  D26(%rip),%xmm3
  subsd  %xmm3,%xmm1
L56:
  subsd  %xmm2,%xmm0
  andpd  D48(%rip),%xmm0
  comisd 0x20(%rsp),%xmm0
  jb     L83
L43:
  movsd  0x18(%rsp),%xmm0
  mov    %r14d,%r9d
  mov    %r15d,%r8d
  mov    %r13,%rcx
  mov    %r12,%rdx
  mov    %rbp,%rsi
  mov    %ebx,%edi
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm1,0x8(%rsp)
  call   <Range@plt>
  movsd  0x10(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm1
  comisd %xmm2,%xmm0
  ja     L44
  comisd %xmm0,%xmm2
  jbe    L83
  movsd  0x28(%rsp),%xmm3
  movsd  %xmm1,0x30(%rsp)
  subsd  %xmm1,%xmm3
  mulsd  D26(%rip),%xmm3
  addsd  %xmm3,%xmm1
  jmp    L56
  nopw   0x0(%rax,%rax,1)
L40:
  movsd  0x38(%rsp),%xmm0
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L65:
  add    $0x48,%rsp
  movapd %xmm3,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L83:
  movsd  0x38(%rsp),%xmm0
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  subsd  %xmm1,%xmm0
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: unpack_data
arguments: tag: D26
data_type: f64
name: unpack_data
arguments: tag: D40
data_type: f32[2]
name: unpack_data
arguments: tag: D48
data_type: qword[2]

tool: D26: 0.681

tool: {"D40": [0.0,3.390625]}

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

assistant: ```c
double ELossCL(int ZP, double MP, int ZT[], double MT[], double CT[], int NTarEl, int IsGas, double TT, double E) {
  double EIn = 0.0, ELow = 0.0, EUp = 0.0, TTStep = 0.0, R = 0.0, DR = 0.0, FFF1 = 0.681;
  if (TT < 0)
    return 0;
  if (E <= 0)
    return E;
  if ((DR = Range(ZP, MP, ZT, MT, CT, NTarEl, IsGas, E) - TT) <= 0)
    return E;
  EIn = EUp = E;
  E = E - E * FFF1;
  TTStep = TT / 100;
  R = Range(ZP, MP, ZT, MT, CT, NTarEl, IsGas, E);
  while (fabs(R - DR) >= TTStep) {
    R = Range(ZP, MP, ZT, MT, CT, NTarEl, IsGas, E);
    if (R > DR) {
      EUp = E;
      E = EUp - (EUp - ELow) * FFF1;
    } else if (R < DR) {
      ELow = E;
      E = ELow + (EUp - ELow) * FFF1;
    } else
      return (EIn - E);
  };
  return (EIn - E);
}
```