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

```asm
<find_out_point>:
  endbr64
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %rsi,%r10
  mov    %rdx,%r11
  mov    %rcx,%rsi
  mov    0x80(%rsp),%rbx
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  movl   $0xffffffff,(%r9)
  movl   $0xffffffff,0x4(%r9)
  movl   $0xffffffff,0x8(%r9)
  mov    $0x0,%edx
L5:
  mov    %edx,%eax
  shr    $0x1f,%eax
  add    %edx,%eax
  sar    %eax
  cltq
  mov    %edx,%ebp
  shr    $0x1f,%ebp
  lea    0x0(%rbp,%rdx,1),%ecx
  and    $0x1,%ecx
  sub    %ebp,%ecx
  add    (%r11,%rax,4),%ecx
  pxor   %xmm0,%xmm0
  cvtsi2sd %ecx,%xmm0
  mulsd  (%r10,%rax,8),%xmm0
  addsd  (%rdi,%rax,8),%xmm0
  movsd  %xmm0,0x20(%rsp,%rdx,8)
  add    $0x1,%rdx
  cmp    $0x6,%rdx
  jne    L5
  mov    $0x0,%eax
  pxor   %xmm2,%xmm2
  movapd %xmm2,%xmm4
  lea    0x20(%rsp),%rdx
  movsd  D61(%rip),%xmm3
  jmp    L268
L161:
  comisd %xmm0,%xmm4
  ja     L51
L84:
  add    $0x8,%rax
  cmp    $0x18,%rax
  je     L73
L268:
  movsd  (%r8,%rax,1),%xmm0
  comisd %xmm2,%xmm0
  jbe    L161
  movsd  0x8(%rdx,%rax,2),%xmm1
  subsd  (%rsi,%rax,1),%xmm1
  addsd  %xmm3,%xmm0
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rsp,%rax,1)
  jmp    L84
L51:
  movsd  (%rsi,%rax,1),%xmm1
  subsd  (%rdx,%rax,2),%xmm1
  movsd  D61(%rip),%xmm5
  subsd  %xmm0,%xmm5
  divsd  %xmm5,%xmm1
  movsd  %xmm1,(%rsp,%rax,1)
  jmp    L84
L73:
  movsd  0x20(%rsp),%xmm4
  movsd  0x28(%rsp),%xmm5
  movsd  0x30(%rsp),%xmm6
  movsd  0x38(%rsp),%xmm7
  movsd  0x40(%rsp),%xmm9
  movsd  0x48(%rsp),%xmm10
  mov    $0x0,%eax
  movsd  D75(%rip),%xmm1
  mov    $0xffffffff,%edi
  movsd  D61(%rip),%xmm3
  jmp    L110
L11:
  add    $0x1,%rax
  cmp    $0x3,%rax
  je     L81
L110:
  mov    %eax,%ecx
  movsd  (%rsp,%rax,8),%xmm0
  comisd %xmm0,%xmm1
  jbe    L11
  comisd %xmm3,%xmm0
  jbe    L11
  mov    $0x0,%edx
L269:
  movapd %xmm0,%xmm2
  mulsd  (%r8,%rdx,1),%xmm2
  addsd  (%rsi,%rdx,1),%xmm2
  movsd  %xmm2,(%rbx,%rdx,1)
  add    $0x8,%rdx
  cmp    $0x18,%rdx
  jne    L269
  movsd  (%rbx),%xmm2
  movapd %xmm4,%xmm8
  subsd  D61(%rip),%xmm8
  comisd %xmm8,%xmm2
  jb     L11
  movapd %xmm5,%xmm8
  addsd  D61(%rip),%xmm8
  comisd %xmm2,%xmm8
  jb     L11
  movsd  0x8(%rbx),%xmm2
  movapd %xmm6,%xmm8
  subsd  D61(%rip),%xmm8
  comisd %xmm8,%xmm2
  jb     L11
  movapd %xmm7,%xmm8
  addsd  D61(%rip),%xmm8
  comisd %xmm2,%xmm8
  jb     L11
  movsd  0x10(%rbx),%xmm2
  movapd %xmm9,%xmm8
  subsd  D61(%rip),%xmm8
  comisd %xmm8,%xmm2
  jb     L11
  movapd %xmm10,%xmm11
  addsd  D61(%rip),%xmm11
  movapd %xmm2,%xmm8
  cmpnlesd %xmm11,%xmm8
  andpd  %xmm8,%xmm1
  andnpd %xmm0,%xmm8
  orpd   %xmm8,%xmm1
  ucomisd %xmm2,%xmm11
  cmovae %ecx,%edi
  jmp    L11
L81:
  cmp    $0x2,%edi
  ja     L116
  mov    $0x0,%eax
L70:
  movapd %xmm1,%xmm0
  mulsd  (%r8,%rax,1),%xmm0
  addsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,(%rbx,%rax,1)
  add    $0x8,%rax
  cmp    $0x18,%rax
  jne    L70
  mov    $0x1,%esi
  jmp    L21
L116:
  movsd  (%rsi),%xmm0
  movsd  %xmm0,(%rbx)
  movsd  0x8(%rsi),%xmm0
  movsd  %xmm0,0x8(%rbx)
  movsd  0x10(%rsi),%xmm0
  movsd  %xmm0,0x10(%rbx)
  mov    $0x0,%esi
L21:
  lea    0x20(%rsp),%rdx
  mov    $0x0,%eax
  movq   D112(%rip),%xmm3
  movsd  D61(%rip),%xmm2
  movapd %xmm2,%xmm4
  jmp    L137
L69:
  add    $0x1,%ecx
  subsd  0x8(%rdx),%xmm0
  andpd  D112(%rip),%xmm0
  comisd %xmm0,%xmm4
  jbe    L232
  mov    %ecx,(%r9,%rax,4)
L232:
  add    $0x1,%rax
  add    $0x10,%rdx
  cmp    $0x3,%rax
  je     L157
L137:
  movsd  (%rbx,%rax,8),%xmm0
  lea    (%rax,%rax,1),%ecx
  movapd %xmm0,%xmm1
  subsd  (%rdx),%xmm1
  andpd  %xmm3,%xmm1
  comisd %xmm1,%xmm2
  jbe    L69
  mov    %ecx,(%r9,%rax,4)
  jmp    L232
L157:
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L146
  mov    %esi,%eax
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  ret
L146:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D61
type: f64
name: extract_symbolic_data
arguments: data_tag: D75
type: byte[8]
name: extract_symbolic_data
arguments: data_tag: D112
type: qword[2]

tool: D61: 1e-08

tool: {"D75": [0x00,0xa0,0xd8,0x85,0x57,0x34,0x76,0x43]}

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

assistant: ```c
int find_out_point(double orig_xyz[3], double r_dxyz[3], int ijk[3], double in_point[3], double ray_direction[3], int face_sect[3], double out_point[3]) {
  int i, j, mincomp, intersection;
  double minmax[6], t[3], mint;
  int i1, j1;
  for (i = 0; i < 3; i++)
    face_sect[i] = -1;
  for (i = 0; i < 6; i++) {
    i1 = i / 2;
    j1 = i % 2;
    minmax[i] = orig_xyz[i1] + (ijk[i1] + j1) * r_dxyz[i1];
  }
  for (i = 0; i < 3; i++) {
    if (ray_direction[i] > 0)
      t[i] = (minmax[i * 2 + 1] - in_point[i]) / (ray_direction[i] + 0.00000001);
    else if (ray_direction[i] < 0)
      t[i] = (in_point[i] - minmax[i * 2]) / (-ray_direction[i] + 0.00000001);
  }
  mint = 1.0E+17;
  mincomp = -1;
  for (i = 0; i < 3; i++) {
    if ((mint > t[i]) && (t[i] > 0.00000001)) {
      for (j = 0; j < 3; j++)
        out_point[j] = in_point[j] + t[i] * ray_direction[j];
      if ((out_point[0] >= minmax[0] - 0.00000001) && (out_point[0] <= minmax[1] + 0.00000001) && (out_point[1] >= minmax[2] - 0.00000001) && (out_point[1] <= minmax[3] + 0.00000001) && (out_point[2] >= minmax[4] - 0.00000001) && (out_point[2] <= minmax[5] + 0.00000001)) {
        mint = t[i];
        mincomp = i;
      }
    }
  }
  if ((mincomp >= 0) && (mincomp <= 2)) {
    intersection = 1;
    for (i = 0; i < 3; i++)
      out_point[i] = in_point[i] + mint * ray_direction[i];
  } else {
    intersection = 0;
    for (i = 0; i < 3; i++)
      out_point[i] = in_point[i];
  }
  for (i = 0; i < 3; i++) {
    if (fabs(out_point[i] - minmax[i * 2]) < 0.00000001)
      face_sect[i] = i * 2;
    else if (fabs(out_point[i] - minmax[i * 2 + 1]) < 0.00000001)
      face_sect[i] = i * 2 + 1;
  }
  return (intersection);
}
```