user: What would the C source look like for the assembly code provided here?

```asm
<pet_c_oplem_dda_att>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %esi,%r12d
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    0xc0(%rsp),%rax
  mov    0xf8(%rsp),%ebx
  mov    %rdx,0x20(%rsp)
  xor    %edx,%edx
  mov    0x108(%rsp),%r14d
  mov    %rdi,0x18(%rsp)
  mov    %rax,0x30(%rsp)
  mov    0xd0(%rsp),%rax
  mov    0x110(%rsp),%rdi
  mov    %r8,0x28(%rsp)
  mov    %rax,0x38(%rsp)
  mov    0xe0(%rsp),%rax
  mov    %rdi,0x78(%rsp)
  mov    0x150(%rsp),%edi
  mov    %rax,0x40(%rsp)
  mov    0x100(%rsp),%eax
  mov    %edi,0x64(%rsp)
  mov    0xf0(%rsp),%rbp
  mov    %eax,%r13d
  mov    0x130(%rsp),%r15
  imul   %r14d,%r13d
  cmp    %ebx,%eax
  cmovl  %ebx,%edx
  cmp    %r14d,%edx
  mov    %r13d,0x50(%rsp)
  cmovl  %r14d,%edx
  imul   %r13d,%ebx
  movslq %edx,%rax
  mov    %edx,%ecx
  imul   $0x55555556,%rax,%rax
  sar    $0x1f,%ecx
  shr    $0x20,%rax
  sub    %ecx,%eax
  lea    (%rax,%rdx,1),%edi
  movslq %edi,%rdi
  shl    $0x2,%rdi
  call   <malloc@plt>
  mov    %ebx,0x70(%rsp)
  mov    %rax,0x10(%rsp)
  lea    0x0(,%rbx,4),%eax
  mov    %rax,%rdi
  mov    %rax,0x68(%rsp)
  call   <malloc@plt>
  mov    %rax,0x58(%rsp)
  test   %ebx,%ebx
  jle    L74
  mov    0x70(%rsp),%edi
  mov    0x78(%rsp),%rbx
  movss  D58(%rip),%xmm1
  lea    -0x1(%rdi),%edx
  mov    %rbx,%rax
  lea    0x4(%rbx,%rdx,4),%rdx
L97:
  movaps %xmm1,%xmm0
  add    $0x4,%rax
  divss  -0x4(%rax),%xmm0
  movss  %xmm0,-0x4(%rax)
  cmp    %rax,%rdx
  jne    L97
L74:
  mov    0x64(%rsp),%eax
  test   %eax,%eax
  jle    L82
  pxor   %xmm1,%xmm1
  pxor   %xmm0,%xmm0
  pxor   %xmm7,%xmm7
  xor    %eax,%eax
  cvtsi2ss %r12d,%xmm1
  mov    %eax,0x54(%rsp)
  cvtsi2ssl 0x64(%rsp),%xmm0
  cvtsi2ss %eax,%xmm7
  divss  %xmm0,%xmm1
  movss  %xmm7,0x60(%rsp)
  movss  %xmm1,0x74(%rsp)
  nopl   0x0(%rax,%rax,1)
L148:
  movss  0x60(%rsp),%xmm0
  movss  0x74(%rsp),%xmm6
  xor    %esi,%esi
  addl   $0x1,0x54(%rsp)
  mov    0x54(%rsp),%eax
  mulss  %xmm6,%xmm0
  mov    0x68(%rsp),%rdx
  mov    0x58(%rsp),%rdi
  addss  D19(%rip),%xmm0
  cvttss2si %xmm0,%r13d
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  movss  %xmm0,0x60(%rsp)
  mulss  %xmm6,%xmm0
  addss  D19(%rip),%xmm0
  cvttss2si %xmm0,%eax
  mov    %eax,%ebx
  mov    %eax,0xc(%rsp)
  call   <memset@plt>
  cmp    %ebx,%r13d
  pxor   %xmm3,%xmm3
  movss  D58(%rip),%xmm4
  movss  D5(%rip),%xmm2
  jge    L72
  mov    0x10(%rsp),%rax
  movslq %r13d,%r13
  add    $0x4,%rax
  mov    %rax,0x48(%rsp)
  nopl   0x0(%rax)
L94:
  mov    0x18(%rsp),%rax
  pxor   %xmm5,%xmm5
  pxor   %xmm1,%xmm1
  pxor   %xmm0,%xmm0
  movzwl (%rax,%r13,2),%ecx
  mov    0x20(%rsp),%rax
  movzwl (%rax,%r13,2),%esi
  mov    0x28(%rsp),%rax
  movzwl (%rax,%r13,2),%edi
  mov    0x30(%rsp),%rax
  movzwl (%rax,%r13,2),%r8d
  mov    0x38(%rsp),%rax
  movzwl (%rax,%r13,2),%r9d
  mov    0x40(%rsp),%rax
  sub    %ecx,%r8d
  movzwl (%rax,%r13,2),%edx
  mov    %r8d,%eax
  sub    %esi,%r9d
  sub    %edi,%edx
  neg    %eax
  cmovs  %r8d,%eax
  cvtsi2ss %eax,%xmm5
  mov    %edx,%eax
  neg    %eax
  cmovs  %edx,%eax
  cvtsi2ss %eax,%xmm1
  mov    %r9d,%eax
  neg    %eax
  cmovs  %r9d,%eax
  cvtsi2ss %eax,%xmm0
  cvttss2si %xmm0,%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  comiss %xmm0,%xmm5
  jbe    L152
  cvttss2si %xmm5,%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
L152:
  comiss %xmm0,%xmm1
  jbe    L78
  cvttss2si %xmm1,%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
L78:
  shl    $0x17,%ecx
  movaps %xmm4,%xmm1
  shl    $0x17,%esi
  divss  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  shl    $0x17,%edi
  cvtsi2ss %r8d,%xmm0
  mulss  %xmm1,%xmm0
  mulss  %xmm2,%xmm0
  cvttss2si %xmm0,%r8d
  pxor   %xmm0,%xmm0
  cvtsi2ss %r9d,%xmm0
  mulss  %xmm1,%xmm0
  mulss  %xmm2,%xmm0
  cvttss2si %xmm0,%r9d
  pxor   %xmm0,%xmm0
  cvtsi2ss %edx,%xmm0
  mulss  %xmm1,%xmm0
  mulss  %xmm2,%xmm0
  cvttss2si %xmm0,%r10d
  test   %eax,%eax
  je     L48
  mov    0x10(%rsp),%rdx
  mov    0x48(%rsp),%rbx
  sub    $0x1,%eax
  movaps %xmm3,%xmm0
  mov    0x50(%rsp),%r12d
  movaps %xmm3,%xmm1
  mov    %rdx,(%rsp)
  lea    (%rbx,%rax,4),%rbx
L40:
  mov    %edi,%eax
  mov    %esi,%r11d
  add    $0x4,%rdx
  add    %r9d,%esi
  sar    $0x17,%eax
  sar    $0x17,%r11d
  add    %r10d,%edi
  imul   %r12d,%eax
  imul   %r14d,%r11d
  add    %eax,%r11d
  mov    %ecx,%eax
  add    %r8d,%ecx
  sar    $0x17,%eax
  add    %r11d,%eax
  mov    %eax,-0x4(%rdx)
  cltq
  addss  0x0(%rbp,%rax,4),%xmm1
  subss  (%r15,%rax,4),%xmm0
  cmp    %rdx,%rbx
  jne    L40
  ucomiss %xmm3,%xmm1
  mov    (%rsp),%r12
  jp     L76
  jne    L76
L48:
  add    $0x1,%r13
  cmp    %r13d,0xc(%rsp)
  jg     L94
L72:
  movslq 0x70(%rsp),%rdx
  xor    %eax,%eax
  test   %edx,%edx
  jle    L122
  mov    0x58(%rsp),%rcx
  mov    0x78(%rsp),%rsi
  nopl   (%rax)
L23:
  movss  0x0(%rbp,%rax,4),%xmm0
  mulss  (%rcx,%rax,4),%xmm0
  mulss  (%rsi,%rax,4),%xmm0
  movss  %xmm0,0x0(%rbp,%rax,4)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L23
L122:
  mov    0x64(%rsp),%ebx
  cmp    %ebx,0x54(%rsp)
  jne    L148
L82:
  mov    0x10(%rsp),%rdi
  call   <free@plt>
  mov    0x58(%rsp),%rdi
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1080 <free@plt>
  nopw   0x0(%rax,%rax,1)
L76:
  mulss  D30(%rip),%xmm0
  movss  %xmm1,(%rsp)
  cvtss2sd %xmm0,%xmm0
  call   <exp@plt>
  movss  (%rsp),%xmm1
  movss  D58(%rip),%xmm4
  mov    0x58(%rsp),%rdx
  movss  D5(%rip),%xmm2
  pxor   %xmm3,%xmm3
  cvtss2sd %xmm1,%xmm1
  mulsd  %xmm0,%xmm1
  movaps %xmm4,%xmm5
  cvtsd2ss %xmm1,%xmm1
  divss  %xmm1,%xmm5
  nopw   0x0(%rax,%rax,1)
L146:
  movslq (%r12),%rax
  add    $0x4,%r12
  lea    (%rdx,%rax,4),%rax
  movss  (%rax),%xmm0
  addss  %xmm5,%xmm0
  movss  %xmm0,(%rax)
  cmp    %r12,%rbx
  jne    L146
  jmp    L48
```

tool_calls: name: parse_data
arguments: data_tag: D58
type_infered: byte[4]
name: parse_data
arguments: data_tag: D19
type_infered: byte[4]
name: parse_data
arguments: data_tag: D5
type_infered: byte[4]
name: parse_data
arguments: data_tag: D30
type_infered: byte[4]

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

tool: {"D19": [0x00,0x00,0x00,0x3f]}

tool: {"D5": [0x00,0x00,0x00,0x4b]}

tool: {"D30": [0x4c,0xa6,0xaa,0x3e]}

assistant: ```c
void pet_c_oplem_dda_att(unsigned short int *X1, int nx1, unsigned short int *Y1, int ny1, unsigned short int *Z1, int nz1, unsigned short int *X2, int nx2, unsigned short int *Y2, int ny2, unsigned short int *Z2, int nz2, float *im, int nim1, int nim2, int nim3, float *NM, int nm1, int nm2, int nm3, float *AM, int am1, int am2, int am3, int nsub) {
  int length, lengthy, lengthz, ilor, n;
  float flength, Qi, Ai;
  float lx, ly, lz;
  int fxinc, fyinc, fzinc, fx, fy, fz;
  int x1, y1, z1, x2, y2, z2, diffx, diffy, diffz;
  int step = nim2 * nim3;
  int ndata, ind;
  if (nim1 > nim2) {
    ndata = nim1;
  }
  if (nim3 > ndata) {
    ndata = nim3;
  }
  ndata += ndata / 3;
  int *buf = (int *)malloc(ndata * sizeof(int));
  step = nim2 * nim3;
  int lor_start, lor_stop, nlor, isub;
  int nvox = step * nim1;
  unsigned int mem_size_F = nvox * sizeof(float);
  float *F = (float *)malloc(mem_size_F);
  for (n = 0; n < nvox; ++n) {
    NM[n] = 1.0f / NM[n];
  }
  for (isub = 0; isub < nsub; ++isub) {
    lor_start = (int)((float)nx1 / nsub * isub + 0.5f);
    lor_stop = (int)((float)nx1 / nsub * (isub + 1) + 0.5f);
    nlor = lor_stop - lor_start;
    memset(F, 0, mem_size_F);
    for (ilor = lor_start; ilor < lor_stop; ++ilor) {
      Qi = 0.0f;
      Ai = 0.0f;
      x1 = X1[ilor];
      x2 = X2[ilor];
      y1 = Y1[ilor];
      y2 = Y2[ilor];
      z1 = Z1[ilor];
      z2 = Z2[ilor];
      diffx = x2 - x1;
      diffy = y2 - y1;
      diffz = z2 - z1;
      lx = abs(diffx);
      ly = abs(diffy);
      lz = abs(diffz);
      length = ly;
      if (lx > length) {
        length = lx;
      }
      if (lz > length) {
        length = lz;
      }
      flength = 1.0f / (float)length;
      fxinc = ((int)diffx * flength * (int)(pow(2, 23)));
      fyinc = ((int)diffy * flength * (int)(pow(2, 23)));
      fzinc = ((int)diffz * flength * (int)(pow(2, 23)));
      fx = ((int)x1 * (int)(pow(2, 23)));
      fy = ((int)y1 * (int)(pow(2, 23)));
      fz = ((int)z1 * (int)(pow(2, 23)));
      for (n = 0; n < length; ++n) {
        ind = (fz >> 23) * step + (fy >> 23) * nim3 + (fx >> 23);
        buf[n] = ind;
        Qi += im[ind];
        Ai -= AM[ind];
        fx = fx + fxinc;
        fy = fy + fyinc;
        fz = fz + fzinc;
      }
      if (Qi == 0.0f) {
        continue;
      }
      Qi = Qi * exp(Ai * 0.3333f);
      Qi = 1.0f / Qi;
      for (n = 0; n < length; ++n) {
        F[buf[n]] += Qi;
      }
    }
    for (n = 0; n < nvox; ++n) {
      im[n] = im[n] * F[n] * NM[n];
    }
  }
  free(buf);
  free(F);
}
```