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

```asm
<bioD3_3d>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    0x90(%rsp),%rax
  mov    0xa8(%rsp),%r12d
  mov    %rdx,0x20(%rsp)
  mov    0xb0(%rsp),%r13d
  mov    0xb8(%rsp),%ebx
  mov    %rsi,0x18(%rsp)
  mov    %rax,0x38(%rsp)
  mov    0x98(%rsp),%rax
  mov    %rcx,0x28(%rsp)
  mov    0xa0(%rsp),%r14
  mov    %rax,0x40(%rsp)
  xor    %eax,%eax
  cmp    $0x1,%r12d
  mov    0xc0(%rsp),%r15d
  setg   %al
  xor    %edx,%edx
  cmp    $0x1,%r13d
  mov    %r8,0x30(%rsp)
  setg   %dl
  add    %edx,%eax
  xor    %edx,%edx
  cmp    $0x1,%ebx
  setg   %dl
  add    %edx,%eax
  cmp    $0x2,%eax
  je     L66
  cmp    $0x1,%eax
  jne    L47
  movsd  D63(%rip),%xmm0
  movss  D1(%rip),%xmm1
  jmp    L36
  nopl   0x0(%rax,%rax,1)
L47:
  test   %eax,%eax
  jne    L18
  movsd  D28(%rip),%xmm0
  movss  D68(%rip),%xmm1
L36:
  mov    $0x1,%r8d
  mov    %ebx,%ecx
  mov    %r13d,%edx
  mov    %r12d,%esi
  mov    %rbp,%rdi
  mov    $0x1,%eax
  mov    %r9,0x48(%rsp)
  movss  %xmm1,0x14(%rsp)
  movsd  %xmm0,0x8(%rsp)
  call   <bio_3d_premult@plt>
  test   %r15d,%r15d
  movsd  0x8(%rsp),%xmm0
  movss  0x14(%rsp),%xmm1
  je     L58
  mov    0x40(%rsp),%rax
  mov    %ebx,0xb8(%rsp)
  mov    %rbp,%rdi
  cvtss2sd %xmm1,%xmm1
  mov    %r13d,0xb0(%rsp)
  mov    0x48(%rsp),%r9
  mov    %rax,0x98(%rsp)
  mov    0x38(%rsp),%rax
  mov    %r12d,0xa8(%rsp)
  mov    0x30(%rsp),%r8
  mov    %r14,0xa0(%rsp)
  mov    0x28(%rsp),%rcx
  mov    %rax,0x90(%rsp)
  mov    0x20(%rsp),%rdx
  mov    $0x2,%eax
  mov    0x18(%rsp),%rsi
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1080 <bioD1__3d@plt>
  cs nopw 0x0(%rax,%rax,1)
L58:
  add    $0x58,%rsp
  mov    %ebx,%r8d
  mov    %r13d,%ecx
  mov    %r12d,%edx
  pop    %rbx
  mov    %r14,%rsi
  mov    %rbp,%rdi
  mov    $0x2,%eax
  pop    %rbp
  pop    %r12
  cvtss2sd %xmm1,%xmm1
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <bioD1_3dskip@plt>
  nopl   0x0(%rax,%rax,1)
L66:
  movsd  D63(%rip),%xmm0
  movss  D68(%rip),%xmm1
  jmp    L36
L18:
  movsd  D63(%rip),%xmm0
  movss  D16(%rip),%xmm1
  jmp    L36
```

tool_calls: name: extract_symbolic_data
arguments: label: D63
datatype: byte[8]
name: extract_symbolic_data
arguments: label: D1
datatype: byte[4]
name: extract_symbolic_data
arguments: label: D28
datatype: f64
name: extract_symbolic_data
arguments: label: D68
datatype: byte[4]
name: extract_symbolic_data
arguments: label: D16
datatype: byte[4]

tool: {"D63": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0xbf]}

tool: {"D1": [0x00,0x00,0x00,0xc0]}

tool: D28: 1

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

tool: {"D16": [0x00,0x00,0x00,0xbf]}

assistant: ```c
void bioD3_3d(invector, HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, xlength, ylength, zlength, ifnotAllskip) float *invector;
float *HLL, *HLH, *HHL, *HHH;
float *LLL, *LLH, *LHL, *LHH;
int xlength, ylength, zlength;
int ifnotAllskip;
{
  float x0 = -0.5;
  float X0 = x0;
  float X1 = x0;
  int pairity;
  float Norm;
  int dim;
  dim = (xlength > 1) + (ylength > 1) + (zlength > 1);
  ;
  if (dim == 3)
    Norm = x0;
  if (dim == 2)
    Norm = 1.0;
  if (dim == 1)
    Norm = 1.0 / x0;
  if (xlength > 1)
    Norm /= (1.0);
  if (ylength > 1)
    Norm /= (1.0);
  if (zlength > 1)
    Norm /= (1.0);
  if (dim == 0) {
    Norm = 1.0;
    X0 = 1.0, X1 = 1.0;
  }
  pairity = 1;
  bio_3d_premult(invector, xlength, ylength, zlength, X0, pairity);
  if (ifnotAllskip)
    bioD1__3d(invector, HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, xlength, ylength, zlength, X1, Norm);
  else
    bioD1_3dskip(invector, LLL, xlength, ylength, zlength, X1, Norm);
}
```