user: Please show me the C implementation that matches the assembly code listed below:

```asm
<ProbabilisticModel>:
  endbr64
  push   %r15
  mov    %esi,%eax
  push   %r14
  push   %r13
  mov    %rcx,%r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  movss  D8(%rip),%xmm1
  mov    0x78(%rsp),%rsi
  mov    %rdx,0x20(%rsp)
  mov    %edi,0x14(%rsp)
  mov    0x80(%rsp),%r14
  mov    (%rsi),%rdx
  mov    0x70(%rsp),%rdi
  mov    %r8,0x28(%rsp)
  movss  %xmm1,(%rdx)
  test   %eax,%eax
  jle    L55
  sub    $0x1,%eax
  mov    %r9,%rcx
  lea    0xc(,%rax,8),%r10
  mov    $0x4,%eax
L61:
  movss  -0x4(%rcx,%rax,1),%xmm0
  mov    (%rsi,%rax,2),%r9
  mov    0x8(%rsi,%rax,2),%r8
  movss  %xmm0,(%rdx,%rax,1)
  movss  (%rcx,%rax,1),%xmm0
  movss  %xmm0,0x4(%rdx,%rax,1)
  addss  -0x4(%rcx,%rax,1),%xmm0
  movss  (%rdx),%xmm2
  subss  %xmm0,%xmm2
  movss  %xmm2,(%rdx)
  movss  -0x4(%rdi,%rax,1),%xmm0
  movss  %xmm0,(%r9,%rax,1)
  movss  (%rdi,%rax,1),%xmm0
  movss  %xmm0,0x4(%r8,%rax,1)
  movaps %xmm1,%xmm0
  movl   $0x0,0x4(%r9,%rax,1)
  movl   $0x0,(%r8,%rax,1)
  subss  -0x4(%rdi,%rax,1),%xmm0
  movss  %xmm0,(%r9)
  movaps %xmm1,%xmm0
  subss  (%rdi,%rax,1),%xmm0
  add    $0x8,%rax
  movss  %xmm0,(%r8)
  cmp    %r10,%rax
  jne    L61
L55:
  mov    0x14(%rsp),%edx
  test   %edx,%edx
  jle    L70
  movslq 0x14(%rsp),%rax
  mov    0x98(%rsp),%r15
  mov    %rsi,%r12
  xor    %ebp,%ebp
  mov    %rax,0x18(%rsp)
  shl    $0x2,%rax
  mov    %rax,0x8(%rsp)
L22:
  mov    0x20(%rsp),%rax
  xor    %ebx,%ebx
  movss  (%rax,%rbp,4),%xmm0
  call   <log@plt>
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  mov    0x90(%rsp),%rax
  movss  %xmm0,(%rax,%rbp,4)
L2:
  mov    (%r12),%rax
  movss  (%rax,%rbx,1),%xmm0
  call   <log@plt>
  pxor   %xmm0,%xmm0
  mov    (%r15),%rdx
  cvtsi2ss %rax,%xmm0
  movss  %xmm0,(%rdx,%rbx,1)
  add    $0x4,%rbx
  cmp    %rbx,0x8(%rsp)
  jne    L2
  add    $0x1,%rbp
  add    $0x8,%r12
  add    $0x8,%r15
  cmp    0x18(%rsp),%rbp
  jne    L22
L70:
  mov    0x14(%rsp),%eax
  mov    0x28(%rsp),%rbp
  xor    %r12d,%r12d
  mov    0x88(%rsp),%rbx
  sub    $0x1,%eax
  lea    0x4(,%rax,4),%rax
  mov    %rax,0x8(%rsp)
L26:
  mov    0x14(%rsp),%eax
  xor    %r15d,%r15d
  test   %eax,%eax
  jle    L34
  xchg   %ax,%ax
L38:
  movss  0x0(%r13,%r12,4),%xmm0
  call   <log@plt>
  pxor   %xmm0,%xmm0
  mov    (%r14,%r12,8),%rdx
  cvtsi2ss %rax,%xmm0
  movss  %xmm0,(%rdx,%r15,1)
  add    $0x4,%r15
  cmp    %r15,0x8(%rsp)
  jne    L38
L34:
  xor    %r15d,%r15d
  nopl   (%rax)
L85:
  mov    0x0(%rbp),%rax
  movss  (%rax,%r15,1),%xmm0
  call   <log@plt>
  pxor   %xmm0,%xmm0
  mov    (%rbx),%rdx
  cvtsi2ss %rax,%xmm0
  movss  %xmm0,(%rdx,%r15,1)
  add    $0x4,%r15
  cmp    $0x400,%r15
  jne    L85
  add    $0x1,%r12
  add    $0x8,%rbp
  add    $0x8,%rbx
  cmp    $0x100,%r12
  jne    L26
  add    $0x38,%rsp
  mov    $0x1,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: unpack_data
arguments: data_tag: D8
type_infered: byte[4]

tool: {"D8": [0x00,0x00,0x80,0x3f]}

assistant: ```c
int ProbabilisticModel(int NumMatrixTypes, int NumInsertStates, float *initDistribMat, float *emitSingle, float **emitPairs, float *gapOpen, float *gapExtend, float **transMat, float **insProb, float **matchProb, float *initialDistribution, float **transProb) {
  int i, j;
  transMat[0][0] = 1;
  for (i = 0; i < NumInsertStates; i++) {
    transMat[0][2 * i + 1] = gapOpen[2 * i];
    transMat[0][2 * i + 2] = gapOpen[2 * i + 1];
    transMat[0][0] -= (gapOpen[2 * i] + gapOpen[2 * i + 1]);
    transMat[2 * i + 1][2 * i + 1] = gapExtend[2 * i];
    transMat[2 * i + 2][2 * i + 2] = gapExtend[2 * i + 1];
    transMat[2 * i + 1][2 * i + 2] = 0;
    transMat[2 * i + 2][2 * i + 1] = 0;
    transMat[2 * i + 1][0] = 1 - gapExtend[2 * i];
    transMat[2 * i + 2][0] = 1 - gapExtend[2 * i + 1];
  }
  for (i = 0; i < NumMatrixTypes; i++) {
    initialDistribution[i] = (float)log((float)initDistribMat[i]);
    for (j = 0; j < NumMatrixTypes; j++) {
      transProb[i][j] = (float)log((float)transMat[i][j]);
    }
  }
  for (i = 0; i < 256; i++) {
    for (j = 0; j < NumMatrixTypes; j++) {
      insProb[i][j] = (float)log((float)emitSingle[i]);
    }
    for (j = 0; j < 256; j++) {
      matchProb[i][j] = (float)log((float)emitPairs[i][j]);
    }
  }
  return 1;
}
```