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

```asm
<evaluateGTRGAMMAINVAR_VECTOR>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %r9,%r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %r8,%rbx
  sub    $0x88,%rsp
  movsd  D47(%rip),%xmm2
  mov    %rdx,0x48(%rsp)
  mov    0xd0(%rsp),%rdx
  mov    0xd8(%rsp),%edi
  subsd  %xmm0,%xmm2
  mov    %rsi,0x18(%rsp)
  mulsd  D35(%rip),%xmm2
  mov    %rdx,0x30(%rsp)
  mov    0xc0(%rsp),%rax
  mov    0xe0(%rsp),%edx
  mov    %rcx,0x20(%rsp)
  mov    %edi,0x40(%rsp)
  mov    %rax,0x28(%rsp)
  mov    0xc8(%rsp),%rax
  mov    %edx,0x44(%rsp)
  mov    %fs:0x28,%rdx
  mov    %rdx,0x78(%rsp)
  xor    %edx,%edx
  movsd  %xmm2,0x10(%rsp)
  movsd  (%rax),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm2,0x50(%rsp)
  movsd  0x8(%rax),%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm2,0x58(%rsp)
  movsd  0x10(%rax),%xmm2
  mulsd  %xmm0,%xmm2
  mulsd  0x18(%rax),%xmm0
  movsd  %xmm2,0x60(%rsp)
  movsd  %xmm0,0x68(%rsp)
  movsd  D22(%rip),%xmm0
  maxsd  %xmm1,%xmm0
  call   <log@plt>
  pxor   %xmm6,%xmm6
  mov    $0x80,%edi
  cvtsi2sd %eax,%xmm6
  movsd  %xmm6,0x8(%rsp)
  call   <malloc@plt>
  mov    %rax,%r15
  mov    %rax,%r14
  lea    0x80(%rax),%r13
L68:
  movsd  (%r12),%xmm1
  mov    D47(%rip),%rax
  add    $0x20,%r14
  add    $0x8,%r12
  mov    %rax,-0x20(%r14)
  movapd %xmm1,%xmm0
  mulsd  (%rbx),%xmm0
  movsd  %xmm1,0x38(%rsp)
  mulsd  0x8(%rsp),%xmm0
  call   <exp@plt>
  movsd  0x38(%rsp),%xmm1
  movsd  %xmm0,-0x18(%r14)
  movsd  0x8(%rbx),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  0x8(%rsp),%xmm0
  call   <exp@plt>
  movsd  0x38(%rsp),%xmm1
  movsd  %xmm0,-0x10(%r14)
  mulsd  0x10(%rbx),%xmm1
  mulsd  0x8(%rsp),%xmm1
  movapd %xmm1,%xmm0
  call   <exp@plt>
  movsd  %xmm0,-0x8(%r14)
  cmp    %r14,%r13
  jne    L68
  mov    0x44(%rsp),%edi
  cmp    %edi,0x40(%rsp)
  jge    L41
  movslq 0x40(%rsp),%r13
  mov    0x44(%rsp),%ebx
  mov    0x48(%rsp),%r14
  mov    %r13,%rdx
  sub    $0x1,%ebx
  lea    0x1(%r13),%r12
  shl    $0x7,%rdx
  sub    %r13d,%ebx
  add    %rdx,%r14
  add    %r12,%rbx
  jmp    L12
  xchg   %ax,%ax
L58:
  addsd  0x50(%rsp,%rax,8),%xmm1
  mov    D11(%rip),%rax
  mov    0x0(%rbp,%r13,4),%edi
  movq   %rax,%xmm0
  movsd  %xmm1,0x8(%rsp)
  call   <pow@plt>
  mulsd  0x8(%rsp),%xmm0
  call   <log@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
L47:
  mov    0x20(%rsp),%rax
  sub    $0xffffffffffffff80,%r14
  movsd  %xmm0,(%rax,%r13,8)
  mov    %r12,%r13
  cmp    %r12,%rbx
  je     L41
  add    $0x1,%r12
L12:
  mov    0x30(%rsp),%rax
  mov    0x28(%rsp),%rcx
  movsd  (%r14),%xmm0
  movsd  0x8(%r14),%xmm1
  movsbl (%rax,%r13,1),%eax
  movsd  0x10(%r14),%xmm5
  movsd  0x18(%r14),%xmm4
  shl    $0x2,%eax
  cltq
  lea    (%rcx,%rax,8),%rax
  movsd  (%rax),%xmm3
  movsd  0x8(%rax),%xmm2
  mulsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  (%r15),%xmm0
  mulsd  0x8(%r15),%xmm1
  addsd  %xmm0,%xmm1
  movsd  0x10(%rax),%xmm0
  mulsd  %xmm0,%xmm5
  mulsd  0x10(%r15),%xmm5
  addsd  %xmm1,%xmm5
  movsd  0x18(%rax),%xmm1
  mulsd  %xmm1,%xmm4
  mulsd  0x18(%r15),%xmm4
  addsd  %xmm5,%xmm4
  movsd  0x20(%r14),%xmm5
  mulsd  %xmm3,%xmm5
  mulsd  0x20(%r15),%xmm5
  addsd  %xmm4,%xmm5
  movsd  0x28(%r14),%xmm4
  mulsd  %xmm2,%xmm4
  mulsd  0x28(%r15),%xmm4
  addsd  %xmm5,%xmm4
  movsd  0x30(%r14),%xmm5
  mulsd  %xmm0,%xmm5
  mulsd  0x30(%r15),%xmm5
  addsd  %xmm4,%xmm5
  movsd  0x38(%r14),%xmm4
  mulsd  %xmm1,%xmm4
  mulsd  0x38(%r15),%xmm4
  addsd  %xmm5,%xmm4
  movsd  0x40(%r14),%xmm5
  mulsd  %xmm3,%xmm5
  mulsd  0x40(%r15),%xmm5
  mulsd  0x60(%r14),%xmm3
  addsd  %xmm4,%xmm5
  movsd  0x48(%r14),%xmm4
  mulsd  %xmm2,%xmm4
  mulsd  0x48(%r15),%xmm4
  addsd  %xmm5,%xmm4
  movsd  0x50(%r14),%xmm5
  mulsd  %xmm0,%xmm5
  mulsd  0x50(%r15),%xmm5
  addsd  %xmm4,%xmm5
  movsd  0x58(%r14),%xmm4
  mulsd  %xmm1,%xmm4
  mulsd  0x58(%r15),%xmm4
  mulsd  0x60(%r15),%xmm3
  mov    0x18(%rsp),%rax
  mulsd  0x68(%r14),%xmm2
  mulsd  0x68(%r15),%xmm2
  movslq (%rax,%r13,4),%rax
  mulsd  0x70(%r14),%xmm0
  addsd  %xmm5,%xmm4
  mulsd  0x70(%r15),%xmm0
  mulsd  0x78(%r14),%xmm1
  mulsd  0x78(%r15),%xmm1
  addsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  mulsd  0x10(%rsp),%xmm1
  cmp    $0x3,%eax
  jle    L58
  movapd %xmm1,%xmm0
  call   <log@plt>
  mov    D11(%rip),%rdx
  mov    0x0(%rbp,%r13,4),%esi
  mov    %eax,0x8(%rsp)
  movq   %rdx,%xmm0
  mov    %esi,0x38(%rsp)
  call   <log@plt>
  pxor   %xmm0,%xmm0
  imul   0x38(%rsp),%eax
  add    0x8(%rsp),%eax
  cvtsi2sd %eax,%xmm0
  jmp    L47
  xchg   %ax,%ax
L41:
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L56
  add    $0x88,%rsp
  mov    %r15,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10a0 <free@plt>
L56:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: tag: D47
type: f64
name: retrieve_data
arguments: tag: D35
type: f64
name: retrieve_data
arguments: tag: D22
type: byte[8]
name: retrieve_data
arguments: tag: D11
type: byte[8]

tool: D47: 1

tool: D35: 0.25

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

tool: {"D11": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x2f]}

assistant: ```c
void evaluateGTRGAMMAINVAR_VECTOR(int *ex2, int *iptr, double *x2_start, double *v, double *EIGN, double *gammaRates, double *tipVector, double *tFreqs, double invariants, double pz, char *tipX1, int lower, int n) {
  double z, lz, term, ki;
  int i;
  double *diagptable;
  double freqs[4];
  double *x1, *x2;
  double scaler = 0.25 * (1.0 - invariants);
  freqs[0] = tFreqs[0] * invariants;
  freqs[1] = tFreqs[1] * invariants;
  freqs[2] = tFreqs[2] * invariants;
  freqs[3] = tFreqs[3] * invariants;
  z = pz;
  if (z < 1.0E-15)
    z = 1.0E-15;
  lz = log(z);
  diagptable = (double *)malloc(sizeof(double) * 16);
  for (i = 0; i < 4; i++) {
    ki = gammaRates[i];
    diagptable[i * 4] = 1.0;
    diagptable[i * 4 + 1] = exp(EIGN[0] * ki * lz);
    diagptable[i * 4 + 2] = exp(EIGN[1] * ki * lz);
    diagptable[i * 4 + 3] = exp(EIGN[2] * ki * lz);
  }
  for (i = lower; i < n; i++) {
    x1 = &(tipVector[4 * tipX1[i]]);
    x2 = &x2_start[16 * i];
    term = x1[0] * x2[0] * diagptable[0];
    term += x1[1] * x2[1] * diagptable[1];
    term += x1[2] * x2[2] * diagptable[2];
    term += x1[3] * x2[3] * diagptable[3];
    term += x1[0] * x2[4] * diagptable[4];
    term += x1[1] * x2[5] * diagptable[5];
    term += x1[2] * x2[6] * diagptable[6];
    term += x1[3] * x2[7] * diagptable[7];
    term += x1[0] * x2[8] * diagptable[8];
    term += x1[1] * x2[9] * diagptable[9];
    term += x1[2] * x2[10] * diagptable[10];
    term += x1[3] * x2[11] * diagptable[11];
    term += x1[0] * x2[12] * diagptable[12];
    term += x1[1] * x2[13] * diagptable[13];
    term += x1[2] * x2[14] * diagptable[14];
    term += x1[3] * x2[15] * diagptable[15];
    if (iptr[i] < 4)
      term = log((scaler * term + freqs[iptr[i]]) * pow((1.0 / 115792089237316195423570985008687907853269984665640564039457584007913129639936.0), ex2[i]));
    else
      term = log(scaler * term) + ex2[i] * log((1.0 / 115792089237316195423570985008687907853269984665640564039457584007913129639936.0));
    v[i] = term;
  }
  free(diagptable);
}
```