user: Converte this assembly snippet into its equivalent C code:

```asm
<ibfstages>:
  endbr64
  push   %r15
  lea    (%r8,%r8,1),%rax
  mov    %rsi,%r10
  mov    %rcx,%rsi
  push   %r14
  mov    %r10d,%ecx
  mov    %rdi,%r14
  mov    %rdx,%r11
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %r8,%rbx
  shl    $0x4,%rbx
  sub    $0xf0,%rsp
  mov    %rax,0xa8(%rsp)
  lea    0x0(,%r8,8),%rax
  mov    %rax,0x98(%rsp)
  lea    0x1(%rax),%rdi
  mov    $0x1,%eax
  shl    %cl,%rax
  mov    %rdx,0xd0(%rsp)
  xor    %edx,%edx
  shr    $0x3,%rax
  mov    %rbx,0xc8(%rsp)
  mov    %r8,%rbx
  mov    %r8,0xa0(%rsp)
  mov    %rax,%r8
  div    %rbx
  mov    %r9,0xe8(%rsp)
  mov    %rax,0xb8(%rsp)
  imul   %rsi,%rax
  lea    (%rax,%rax,1),%rbx
  mov    %rax,0x90(%rsp)
  shl    $0x2,%rax
  mov    %rbx,0xc0(%rsp)
  mov    %rax,-0x8(%rsp)
  test   %r9,%r9
  jle    L15
  imul   %rsi,%r8
  lea    -0x2(%r10),%ecx
  mov    %rsi,%rax
  shl    %cl,%rax
  lea    (%r11,%rax,4),%rcx
  lea    0x0(,%r8,4),%rax
  mov    %rcx,0xb0(%rsp)
  mov    %r14,%r8
  mov    %rax,%rbx
  add    %rax,%r11
  sub    %rax,%rcx
  mov    %rax,0x80(%rsp)
  neg    %rbx
  mov    %r11,0xd8(%rsp)
  mov    %rbx,0x88(%rsp)
  mov    %rcx,0xe0(%rsp)
  nopl   (%rax)
L31:
  mov    0xd0(%rsp),%rax
  mov    0xa8(%rsp),%rsi
  movss  (%rax),%xmm3
  mov    0xb0(%rsp),%rax
  movss  %xmm3,-0x30(%rsp)
  movss  (%rax),%xmm3
  mov    0xd8(%rsp),%rax
  movss  %xmm3,-0x2c(%rsp)
  movss  (%rax),%xmm3
  mov    0xe0(%rsp),%rax
  movss  %xmm3,-0x74(%rsp)
  movss  (%rax),%xmm3
  lea    0x0(,%rsi,4),%rax
  shl    $0x3,%rsi
  lea    (%r8,%rax,1),%rbx
  mov    %rsi,0xa8(%rsp)
  lea    (%rbx,%rax,1),%rcx
  movss  %xmm3,-0x70(%rsp)
  add    %rcx,%rax
  mov    %rcx,0x28(%rsp)
  mov    %rax,0x30(%rsp)
  mov    0xa0(%rsp),%rax
  mov    %rax,-0x28(%rsp)
  test   %rax,%rax
  je     L37
  movss  -0x2c(%rsp),%xmm3
  movss  -0x30(%rsp),%xmm1
  xor    %r9d,%r9d
  lea    0x0(,%rdi,4),%r14
  mov    0x98(%rsp),%rax
  mov    0xa0(%rsp),%rdi
  movss  %xmm3,-0x78(%rsp)
  lea    0x0(,%rax,4),%r13
  mov    0xc8(%rsp),%rax
  movss  %xmm1,-0x6c(%rsp)
  movss  %xmm3,-0x34(%rsp)
  lea    0x0(,%rax,4),%r11
  mov    %rdi,%rax
  movss  %xmm1,0x38(%rsp)
  shr    $0x3f,%rax
  mov    %r11,0x50(%rsp)
  add    %rdi,%rax
  mov    0xd0(%rsp),%rdi
  sar    %rax
  mov    %rax,0x18(%rsp)
  mov    0xc0(%rsp),%rax
  mov    %rdi,0x8(%rsp)
  shl    $0x2,%rax
  mov    %rdi,(%rsp)
  mov    %rax,0x58(%rsp)
  neg    %rax
  mov    %rax,0x60(%rsp)
  mov    0x90(%rsp),%rax
  mov    %rdi,-0x20(%rsp)
  shl    $0x2,%rax
  mov    %rax,0x68(%rsp)
  neg    %rax
  mov    %rax,0x70(%rsp)
  mov    0xb8(%rsp),%rax
  sub    $0x1,%rax
  mov    %rax,%r10
  mov    %rax,0x78(%rsp)
  mov    0xb0(%rsp),%rax
  imul   %r11,%r10
  mov    %r8,%r11
  mov    %rbx,%r8
  mov    %rax,0x10(%rsp)
  mov    %rax,-0x10(%rsp)
  mov    %rax,-0x18(%rsp)
  mov    %r9,%rax
  mov    %r10,%r9
  mov    %rax,%r10
  nopl   0x0(%rax,%rax,1)
L26:
  mov    0x28(%rsp),%rbx
  mov    0x30(%rsp),%rax
  lea    (%r8,%r10,1),%rcx
  lea    (%r11,%r10,1),%rsi
  movss  (%r11,%r10,1),%xmm5
  movss  0x4(%r11,%r10,1),%xmm13
  lea    (%rbx,%r10,1),%rdx
  mov    0x78(%rsp),%rbx
  movss  (%r8,%r10,1),%xmm4
  add    %r10,%rax
  movss  0x4(%r8,%r10,1),%xmm2
  test   %rbx,%rbx
  je     L35
  mov    0x50(%rsp),%rbp
  movss  -0x30(%rsp),%xmm8
  mov    %rbx,%rdi
  movaps %xmm13,%xmm1
  movss  -0x2c(%rsp),%xmm9
  movss  0x38(%rsp),%xmm10
  movaps %xmm2,%xmm3
  movss  -0x34(%rsp),%xmm11
  nopl   0x0(%rax)
L19:
  movaps %xmm8,%xmm6
  movaps %xmm9,%xmm2
  movss  (%rax),%xmm14
  movss  (%rdx),%xmm0
  mulss  %xmm4,%xmm6
  movss  0x4(%rax),%xmm12
  movss  0x4(%rdx),%xmm7
  mulss  %xmm3,%xmm2
  movaps %xmm0,%xmm15
  addss  %xmm0,%xmm0
  mulss  %xmm9,%xmm4
  mulss  %xmm8,%xmm3
  addss  %xmm5,%xmm6
  addss  %xmm5,%xmm5
  subss  %xmm2,%xmm6
  movaps %xmm4,%xmm2
  movaps %xmm14,%xmm4
  mulss  %xmm8,%xmm4
  addss  %xmm1,%xmm2
  addss  %xmm1,%xmm1
  subss  %xmm6,%xmm5
  addss  %xmm3,%xmm2
  movss  (%rsi,%r13,1),%xmm3
  subss  %xmm4,%xmm15
  movaps %xmm12,%xmm4
  mulss  %xmm9,%xmm4
  subss  %xmm2,%xmm1
  mulss  %xmm8,%xmm12
  movaps %xmm1,%xmm13
  movss  (%rsi,%r14,1),%xmm1
  addss  %xmm15,%xmm4
  movaps %xmm14,%xmm15
  movaps %xmm7,%xmm14
  mulss  %xmm9,%xmm15
  addss  %xmm7,%xmm7
  subss  %xmm4,%xmm0
  subss  %xmm15,%xmm14
  movaps %xmm0,%xmm15
  mulss  %xmm11,%xmm0
  mulss  %xmm10,%xmm15
  subss  %xmm12,%xmm14
  subss  %xmm14,%xmm7
  addss  %xmm2,%xmm0
  addss  %xmm2,%xmm2
  addss  %xmm6,%xmm15
  addss  %xmm6,%xmm6
  movaps %xmm7,%xmm12
  mulss  %xmm10,%xmm7
  mulss  %xmm11,%xmm12
  addss  %xmm7,%xmm0
  movaps %xmm4,%xmm7
  mulss  %xmm11,%xmm7
  subss  %xmm12,%xmm15
  mulss  %xmm10,%xmm4
  subss  %xmm0,%xmm2
  subss  %xmm15,%xmm6
  movss  %xmm15,-0x38(%rsp)
  movss  (%rcx,%r14,1),%xmm15
  movss  %xmm2,-0x68(%rsp)
  mulss  %xmm9,%xmm15
  movaps %xmm7,%xmm2
  movaps %xmm14,%xmm7
  mulss  %xmm10,%xmm7
  addss  %xmm5,%xmm2
  movss  %xmm6,-0x58(%rsp)
  addss  %xmm5,%xmm5
  mulss  %xmm11,%xmm14
  movaps %xmm7,%xmm6
  movaps %xmm2,%xmm7
  movaps %xmm13,%xmm2
  addss  %xmm6,%xmm7
  subss  %xmm4,%xmm2
  movss  (%rcx,%r14,1),%xmm4
  movss  (%rdx,%r14,1),%xmm6
  addss  %xmm13,%xmm13
  mulss  %xmm8,%xmm4
  addss  %xmm6,%xmm6
  subss  %xmm7,%xmm5
  addss  %xmm14,%xmm2
  movss  %xmm7,-0x50(%rsp)
  movss  %xmm5,-0x60(%rsp)
  movss  (%rcx,%r13,1),%xmm5
  subss  %xmm2,%xmm13
  movss  %xmm2,-0x48(%rsp)
  mulss  %xmm8,%xmm5
  addss  %xmm3,%xmm5
  addss  %xmm3,%xmm3
  subss  %xmm15,%xmm5
  movss  (%rcx,%r13,1),%xmm15
  mulss  %xmm9,%xmm15
  subss  %xmm5,%xmm3
  movaps %xmm15,%xmm2
  movss  (%rdx,%r13,1),%xmm15
  addss  %xmm1,%xmm2
  addss  %xmm1,%xmm1
  addss  %xmm4,%xmm2
  movss  (%rax,%r13,1),%xmm4
  mulss  %xmm8,%xmm4
  subss  %xmm2,%xmm1
  subss  %xmm4,%xmm15
  movss  (%rax,%r14,1),%xmm4
  mulss  %xmm9,%xmm4
  movaps %xmm15,%xmm12
  movss  (%rdx,%r14,1),%xmm15
  addss  %xmm4,%xmm12
  movss  (%rax,%r13,1),%xmm4
  mulss  %xmm9,%xmm4
  subss  %xmm4,%xmm15
  movss  (%rax,%r14,1),%xmm4
  mulss  %xmm8,%xmm4
  movaps %xmm15,%xmm7
  subss  %xmm4,%xmm7
  movss  (%rdx,%r13,1),%xmm4
  addss  %xmm4,%xmm4
  subss  %xmm7,%xmm6
  subss  %xmm12,%xmm4
  movaps %xmm6,%xmm15
  mulss  %xmm11,%xmm15
  mulss  %xmm10,%xmm6
  movaps %xmm4,%xmm14
  mulss  %xmm10,%xmm14
  mulss  %xmm11,%xmm4
  addss  %xmm5,%xmm14
  addss  %xmm5,%xmm5
  addss  %xmm2,%xmm4
  addss  %xmm2,%xmm2
  subss  %xmm15,%xmm14
  movaps %xmm12,%xmm15
  mulss  %xmm11,%xmm15
  addss  %xmm6,%xmm4
  mulss  %xmm10,%xmm12
  subss  %xmm14,%xmm5
  subss  %xmm4,%xmm2
  movss  %xmm5,-0x40(%rsp)
  movaps %xmm7,%xmm5
  addss  %xmm3,%xmm15
  addss  %xmm3,%xmm3
  mulss  %xmm10,%xmm5
  mulss  %xmm11,%xmm7
  addss  %xmm5,%xmm15
  movaps %xmm1,%xmm5
  addss  %xmm1,%xmm1
  subss  %xmm12,%xmm5
  movss  -0x38(%rsp),%xmm12
  movaps %xmm12,%xmm6
  subss  %xmm15,%xmm3
  addss  %xmm5,%xmm7
  movss  -0x6c(%rsp),%xmm5
  mulss  %xmm14,%xmm5
  mulss  -0x78(%rsp),%xmm14
  subss  %xmm7,%xmm1
  subss  %xmm5,%xmm6
  movss  -0x78(%rsp),%xmm5
  mulss  %xmm4,%xmm5
  mulss  -0x6c(%rsp),%xmm4
  addss  %xmm6,%xmm5
  movaps %xmm0,%xmm6
  addss  %xmm0,%xmm0
  subss  %xmm14,%xmm6
  movaps %xmm12,%xmm14
  addss  %xmm12,%xmm14
  movss  -0x74(%rsp),%xmm12
  mulss  %xmm3,%xmm12
  subss  %xmm4,%xmm6
  movss  -0x60(%rsp),%xmm4
  subss  %xmm5,%xmm14
  subss  %xmm6,%xmm0
  subss  %xmm12,%xmm4
  movss  -0x70(%rsp),%xmm12
  movss  %xmm5,(%rsi,%r13,1)
  movss  %xmm6,(%rsi,%r14,1)
  mulss  %xmm1,%xmm12
  movss  %xmm14,(%rsi)
  movss  -0x58(%rsp),%xmm14
  mulss  -0x74(%rsp),%xmm1
  movss  %xmm0,0x4(%rsi)
  add    %rbp,%rsi
  movss  -0x78(%rsp),%xmm0
  movss  (%rsi),%xmm5
  movaps %xmm14,%xmm6
  addss  %xmm12,%xmm4
  movss  -0x70(%rsp),%xmm12
  mulss  %xmm3,%xmm12
  movaps %xmm13,%xmm3
  addss  %xmm13,%xmm13
  subss  %xmm12,%xmm3
  subss  %xmm1,%xmm3
  movss  -0x60(%rsp),%xmm1
  addss  %xmm1,%xmm1
  subss  %xmm3,%xmm13
  movaps %xmm1,%xmm12
  movss  0x4(%rsi),%xmm1
  movss  %xmm4,(%rcx,%r13,1)
  subss  %xmm4,%xmm12
  movss  %xmm3,(%rcx,%r14,1)
  movss  %xmm13,0x4(%rcx)
  movss  -0x6c(%rsp),%xmm13
  movss  %xmm12,(%rcx)
  movss  -0x40(%rsp),%xmm12
  add    %rbp,%rcx
  movss  (%rcx),%xmm4
  movss  0x4(%rcx),%xmm3
  mulss  %xmm12,%xmm0
  subss  %xmm0,%xmm6
  movaps %xmm13,%xmm0
  mulss  %xmm2,%xmm0
  mulss  -0x78(%rsp),%xmm2
  subss  %xmm0,%xmm6
  movaps %xmm12,%xmm0
  mulss  %xmm13,%xmm0
  movss  -0x70(%rsp),%xmm13
  addss  -0x68(%rsp),%xmm0
  mulss  %xmm15,%xmm13
  movss  %xmm6,(%rdx)
  mulss  -0x74(%rsp),%xmm15
  subss  %xmm2,%xmm0
  movaps %xmm13,%xmm12
  movss  -0x50(%rsp),%xmm13
  movss  %xmm0,0x4(%rdx)
  movaps %xmm13,%xmm2
  subss  %xmm12,%xmm2
  movss  -0x74(%rsp),%xmm12
  mulss  %xmm7,%xmm12
  mulss  -0x70(%rsp),%xmm7
  subss  %xmm12,%xmm2
  movss  -0x48(%rsp),%xmm12
  addss  %xmm12,%xmm15
  subss  %xmm7,%xmm15
  movaps %xmm14,%xmm7
  addss  %xmm14,%xmm7
  subss  %xmm6,%xmm7
  movss  -0x68(%rsp),%xmm6
  addss  %xmm6,%xmm6
  movss  %xmm7,(%rdx,%r13,1)
  subss  %xmm0,%xmm6
  movaps %xmm13,%xmm0
  addss  %xmm13,%xmm0
  movss  %xmm6,(%rdx,%r14,1)
  add    %rbp,%rdx
  subss  %xmm2,%xmm0
  movss  %xmm2,(%rax)
  movss  %xmm15,0x4(%rax)
  movss  %xmm0,(%rax,%r13,1)
  movaps %xmm12,%xmm0
  addss  %xmm12,%xmm0
  subss  %xmm15,%xmm0
  movss  %xmm0,(%rax,%r14,1)
  add    %rbp,%rax
  sub    $0x1,%rdi
  jne    L19
  mov    0x28(%rsp),%rax
  movaps %xmm1,%xmm13
  movaps %xmm3,%xmm2
  lea    (%r11,%r9,1),%rsi
  mov    %rbp,0x50(%rsp)
  lea    (%r8,%r9,1),%rcx
  lea    (%rax,%r9,1),%rdx
  mov    0x30(%rsp),%rax
  add    %r9,%rax
L35:
  lea    (%rsi,%r14,1),%rdi
  movss  (%rax),%xmm7
  movss  (%rdx),%xmm0
  lea    (%rdx,%r13,1),%r15
  movss  -0x30(%rsp),%xmm10
  mov    %rdi,-0x58(%rsp)
  lea    (%rcx,%r14,1),%r12
  lea    (%rsi,%r13,1),%rbx
  movss  0x4(%rax),%xmm8
  movss  (%r15),%xmm9
  mov    %r15,-0x48(%rsp)
  lea    (%rdx,%r14,1),%r15
  movss  -0x2c(%rsp),%xmm14
  movaps %xmm10,%xmm1
  movss  (%r15),%xmm15
  mov    %r15,-0x40(%rsp)
  mulss  %xmm4,%xmm1
  lea    (%rax,%r14,1),%r15
  movss  0x4(%rdx),%xmm6
  movss  -0x34(%rsp),%xmm12
  mulss  %xmm14,%xmm4
  movaps %xmm14,%xmm3
  mov    %r12,-0x68(%rsp)
  lea    (%rax,%r13,1),%r12
  mulss  %xmm2,%xmm3
  mov    %r12,-0x60(%rsp)
  mov    -0x8(%rsp),%r12
  mulss  %xmm10,%xmm2
  mov    %r15,0x20(%rsp)
  addss  %xmm5,%xmm1
  addss  %xmm5,%xmm5
  addss  %xmm13,%xmm4
  addss  %xmm13,%xmm13
  subss  %xmm3,%xmm1
  movaps %xmm10,%xmm3
  movss  (%r15),%xmm10
  mov    %r12,%r15
  addss  %xmm2,%xmm4
  movss  (%rdi),%xmm2
  lea    (%rcx,%r13,1),%rdi
  neg    %r15
  movss  (%rdi),%xmm11
  movss  %xmm10,0x44(%rsp)
  movaps %xmm6,%xmm10
  addss  %xmm6,%xmm6
  movss  %xmm2,0x3c(%rsp)
  movaps %xmm3,%xmm2
  subss  %xmm1,%xmm5
  mov    %r15,%rbp
  mulss  %xmm7,%xmm2
  subss  %xmm4,%xmm13
  mov    %rdi,-0x50(%rsp)
  mov    -0x28(%rsp),%rdi
  movss  %xmm11,0x40(%rsp)
  movaps %xmm3,%xmm11
  movaps %xmm0,%xmm3
  addss  %xmm0,%xmm0
  mulss  %xmm14,%xmm7
  cmp    %rdi,0x18(%rsp)
  cmovne %r12,%rbp
  subss  %xmm2,%xmm3
  movaps %xmm14,%xmm2
  mov    %rbp,-0x8(%rsp)
  lea    0x0(,%rbp,4),%r15
  mulss  %xmm8,%xmm2
  add    %r15,-0x20(%rsp)
  subss  %xmm7,%xmm10
  addss  %xmm2,%xmm3
  movaps %xmm8,%xmm2
  mulss  %xmm11,%xmm2
  subss  %xmm3,%xmm0
  subss  %xmm2,%xmm10
  subss  %xmm10,%xmm6
  movaps %xmm10,%xmm7
  movss  0x38(%rsp),%xmm10
  movaps %xmm10,%xmm8
  mulss  %xmm0,%xmm8
  movaps %xmm6,%xmm2
  movaps %xmm12,%xmm6
  mulss  %xmm12,%xmm0
  mulss  %xmm2,%xmm6
  mulss  %xmm10,%xmm2
  addss  %xmm1,%xmm8
  addss  %xmm1,%xmm1
  addss  %xmm4,%xmm0
  addss  %xmm4,%xmm4
  subss  %xmm6,%xmm8
  movaps %xmm10,%xmm6
  addss  %xmm2,%xmm0
  subss  %xmm8,%xmm1
  movaps %xmm8,%xmm14
  subss  %xmm0,%xmm4
  movaps %xmm0,%xmm10
  movaps %xmm12,%xmm0
  mulss  %xmm3,%xmm0
  movss  %xmm1,-0x38(%rsp)
  mulss  %xmm6,%xmm3
  movaps %xmm6,%xmm1
  mulss  %xmm7,%xmm1
  movss  %xmm4,0x48(%rsp)
  mulss  %xmm12,%xmm7
  addss  %xmm5,%xmm0
  addss  %xmm5,%xmm5
  addss  %xmm1,%xmm0
  movaps %xmm0,%xmm11
  movaps %xmm13,%xmm0
  addss  %xmm13,%xmm13
  subss  %xmm3,%xmm0
  addss  %xmm7,%xmm0
  movaps %xmm0,%xmm12
  movaps %xmm5,%xmm0
  subss  %xmm11,%xmm0
  subss  %xmm12,%xmm13
  movss  %xmm0,0x4c(%rsp)
  mov    -0x20(%rsp),%rbp
  mov    0x68(%rsp),%r12
  movss  0x40(%rsp),%xmm1
  movss  -0x30(%rsp),%xmm7
  movss  -0x2c(%rsp),%xmm5
  add    %r12,0x8(%rsp)
  movaps %xmm1,%xmm0
  mov    0x70(%rsp),%r12
  movss  (%rbx),%xmm6
  movaps %xmm7,%xmm4
  mulss  %xmm7,%xmm0
  add    %r12,0x10(%rsp)
  mov    -0x68(%rsp),%r12
  mulss  %xmm5,%xmm1
  movss  0x3c(%rsp),%xmm3
  mov    0x58(%rsp),%rdi
  movss  0x44(%rsp),%xmm8
  sub    %r15,-0x18(%rsp)
  add    %rdi,(%rsp)
  mov    -0x18(%rsp),%r15
  mov    0x60(%rsp),%rdi
  add    %rdi,-0x10(%rsp)
  addss  %xmm0,%xmm6
  movaps %xmm5,%xmm0
  mov    -0x10(%rsp),%rdi
  mulss  (%r12),%xmm0
  movaps %xmm1,%xmm2
  movss  (%rbx),%xmm1
  addss  %xmm3,%xmm2
  addss  %xmm1,%xmm1
  subss  %xmm0,%xmm6
  movss  (%r12),%xmm0
  mov    -0x60(%rsp),%r12
  mulss  %xmm7,%xmm0
  mulss  (%r12),%xmm4
  subss  %xmm6,%xmm1
  addss  %xmm0,%xmm2
  movaps %xmm3,%xmm0
  addss  %xmm3,%xmm0
  movaps %xmm9,%xmm3
  addss  %xmm9,%xmm9
  subss  %xmm4,%xmm3
  movaps %xmm8,%xmm4
  mulss  %xmm5,%xmm4
  mulss  (%r12),%xmm5
  subss  %xmm2,%xmm0
  addss  %xmm4,%xmm3
  movaps %xmm15,%xmm4
  addss  %xmm15,%xmm15
  subss  %xmm5,%xmm4
  movaps %xmm8,%xmm5
  mulss  %xmm7,%xmm5
  movaps %xmm15,%xmm7
  movss  0x38(%rsp),%xmm15
  subss  %xmm5,%xmm4
  movaps %xmm9,%xmm5
  movss  -0x34(%rsp),%xmm9
  subss  %xmm3,%xmm5
  subss  %xmm4,%xmm7
  movaps %xmm7,%xmm8
  movaps %xmm15,%xmm7
  mulss  %xmm5,%xmm7
  mulss  -0x34(%rsp),%xmm5
  mulss  %xmm8,%xmm9
  mulss  %xmm15,%xmm8
  addss  %xmm6,%xmm7
  addss  %xmm6,%xmm6
  addss  %xmm2,%xmm5
  addss  %xmm2,%xmm2
  subss  %xmm9,%xmm7
  addss  %xmm8,%xmm5
  subss  %xmm7,%xmm6
  subss  %xmm5,%xmm2
  movaps %xmm6,%xmm9
  movaps %xmm15,%xmm6
  mulss  %xmm4,%xmm6
  movaps %xmm2,%xmm8
  movss  -0x34(%rsp),%xmm2
  mulss  -0x34(%rsp),%xmm4
  mulss  %xmm3,%xmm2
  mulss  %xmm15,%xmm3
  addss  %xmm1,%xmm2
  addss  %xmm1,%xmm1
  addss  %xmm6,%xmm2
  movaps %xmm0,%xmm6
  addss  %xmm0,%xmm0
  subss  %xmm3,%xmm6
  subss  %xmm2,%xmm1
  movaps %xmm6,%xmm3
  movss  0x0(%rbp),%xmm6
  addss  %xmm4,%xmm3
  movss  (%r15),%xmm4
  mov    (%rsp),%r15
  movss  %xmm6,-0x30(%rsp)
  movss  %xmm4,-0x2c(%rsp)
  movss  (%r15),%xmm4
  subss  %xmm3,%xmm0
  movss  %xmm4,0x38(%rsp)
  movss  (%rdi),%xmm4
  mov    -0x28(%rsp),%rdi
  movss  %xmm4,-0x34(%rsp)
  cmp    %rdi,0x18(%rsp)
  jb     L69
  xorps  D0(%rip),%xmm6
  movss  %xmm6,-0x30(%rsp)
L69:
  movss  -0x6c(%rsp),%xmm15
  movaps %xmm14,%xmm6
  add    $0x8,%r10
  add    $0x8,%r9
  mov    0x10(%rsp),%rdi
  movaps %xmm15,%xmm4
  mulss  %xmm7,%xmm4
  mulss  -0x78(%rsp),%xmm7
  subss  %xmm4,%xmm6
  movss  -0x78(%rsp),%xmm4
  mulss  %xmm5,%xmm4
  mulss  %xmm15,%xmm5
  addss  %xmm4,%xmm6
  movaps %xmm10,%xmm4
  addss  %xmm10,%xmm10
  subss  %xmm7,%xmm4
  movss  -0x78(%rsp),%xmm7
  mulss  %xmm9,%xmm7
  movss  %xmm6,(%rbx)
  mov    -0x58(%rsp),%rbx
  subss  %xmm5,%xmm4
  movaps %xmm14,%xmm5
  addss  %xmm14,%xmm5
  movss  -0x6c(%rsp),%xmm14
  mulss  %xmm14,%xmm9
  subss  %xmm4,%xmm10
  movaps %xmm5,%xmm15
  movss  -0x38(%rsp),%xmm5
  subss  %xmm6,%xmm15
  movss  -0x38(%rsp),%xmm6
  subss  %xmm7,%xmm5
  movaps %xmm14,%xmm7
  movss  0x48(%rsp),%xmm14
  mulss  %xmm8,%xmm7
  addss  %xmm6,%xmm6
  mulss  -0x78(%rsp),%xmm8
  addss  %xmm14,%xmm9
  subss  %xmm7,%xmm5
  subss  %xmm8,%xmm9
  movss  %xmm5,(%rdx)
  subss  %xmm5,%xmm6
  movss  %xmm4,(%rbx)
  mov    0x8(%rsp),%rbx
  movaps %xmm14,%xmm4
  movss  %xmm9,0x4(%rdx)
  addss  %xmm14,%xmm4
  mov    -0x48(%rsp),%rdx
  movss  (%rbx),%xmm7
  movss  %xmm7,-0x6c(%rsp)
  movss  (%rdi),%xmm7
  subss  %xmm9,%xmm4
  movss  %xmm15,(%rsi)
  movss  %xmm7,-0x78(%rsp)
  movss  -0x74(%rsp),%xmm7
  movss  %xmm6,(%rdx)
  movss  -0x70(%rsp),%xmm6
  movss  %xmm10,0x4(%rsi)
  mov    -0x40(%rsp),%rsi
  movss  0x4c(%rsp),%xmm10
  movss  %xmm4,(%rsi)
  movaps %xmm7,%xmm4
  movaps %xmm10,%xmm5
  mov    -0x50(%rsp),%rsi
  mulss  %xmm1,%xmm4
  movaps %xmm10,%xmm14
  mulss  %xmm6,%xmm1
  addss  %xmm10,%xmm14
  subss  %xmm4,%xmm5
  movaps %xmm6,%xmm4
  mulss  %xmm0,%xmm4
  mulss  %xmm7,%xmm0
  addss  %xmm5,%xmm4
  movaps %xmm1,%xmm5
  movaps %xmm13,%xmm1
  subss  %xmm5,%xmm1
  movaps %xmm6,%xmm5
  addss  %xmm13,%xmm13
  mulss  %xmm2,%xmm5
  mulss  %xmm7,%xmm2
  movss  %xmm4,(%rsi)
  subss  %xmm4,%xmm14
  subss  %xmm0,%xmm1
  movaps %xmm11,%xmm0
  addss  %xmm11,%xmm11
  subss  %xmm5,%xmm0
  movaps %xmm7,%xmm5
  subss  %xmm1,%xmm13
  mulss  %xmm3,%xmm5
  addss  %xmm12,%xmm2
  addss  %xmm12,%xmm12
  mulss  %xmm6,%xmm3
  subss  %xmm5,%xmm0
  subss  %xmm3,%xmm2
  movss  %xmm0,(%rax)
  mov    -0x68(%rsp),%rsi
  subss  %xmm0,%xmm11
  subq   $0x1,-0x28(%rsp)
  subss  %xmm2,%xmm12
  movss  %xmm1,(%rsi)
  movss  %xmm2,0x4(%rax)
  mov    0x80(%rsp),%rax
  movss  (%rbx,%rax,1),%xmm3
  mov    0x88(%rsp),%rax
  movss  %xmm3,-0x74(%rsp)
  movss  (%rdi,%rax,1),%xmm3
  mov    -0x60(%rsp),%rax
  movss  %xmm14,(%rcx)
  movss  %xmm11,(%rax)
  mov    0x20(%rsp),%rax
  movss  %xmm13,0x4(%rcx)
  movss  %xmm3,-0x70(%rsp)
  movss  %xmm12,(%rax)
  jne    L26
  mov    %r11,%r8
L37:
  shrq   $0x3,0x90(%rsp)
  mov    0x90(%rsp),%rax
  shlq   $0x3,0x98(%rsp)
  shrq   $0x3,0xb8(%rsp)
  shl    $0x2,%rax
  shrq   $0x3,0xc0(%rsp)
  shlq   $0x3,0xa0(%rsp)
  shlq   $0x3,0xc8(%rsp)
  subq   $0x1,0xe8(%rsp)
  mov    %rax,-0x8(%rsp)
  mov    0x98(%rsp),%rax
  lea    0x1(%rax),%rdi
  jne    L31
L15:
  add    $0xf0,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_tag: D0
type: byte[16]

tool: {"D0": [0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

assistant: ```c
void ibfstages(float *ioptr, long M, float *Utbl, long Ustride, long NDiffU, long StageCnt) {
  unsigned long pos;
  unsigned long posi;
  unsigned long pinc;
  unsigned long pnext;
  unsigned long NSameU;
  unsigned long Uinc;
  unsigned long Uinc2;
  unsigned long Uinc4;
  unsigned long DiffUCnt;
  unsigned long SameUCnt;
  unsigned long U2toU3;
  float *pstrt;
  float *p0r, *p1r, *p2r, *p3r;
  float *u0r, *u0i, *u1r, *u1i, *u2r, *u2i;
  float w0r, w0i, w1r, w1i, w2r, w2i, w3r, w3i;
  float f0r, f0i, f1r, f1i, f2r, f2i, f3r, f3i;
  float f4r, f4i, f5r, f5i, f6r, f6i, f7r, f7i;
  float t0r, t0i, t1r, t1i;
  const float Two = 2.0;
  pinc = NDiffU * 2;
  pnext = pinc * 8;
  pos = pinc * 4;
  posi = pos + 1;
  NSameU = ((unsigned long)1 << (M)) / 8 / NDiffU;
  Uinc = NSameU * Ustride;
  Uinc2 = Uinc * 2;
  Uinc4 = Uinc * 4;
  U2toU3 = (((unsigned long)1 << (M)) / 8) * Ustride;
  for (; StageCnt > 0; StageCnt--) {
    u0r = &Utbl[0];
    u0i = &Utbl[((unsigned long)1 << (M - 2)) * Ustride];
    u1r = u0r;
    u1i = u0i;
    u2r = u0r;
    u2i = u0i;
    w0r = *u0r;
    w0i = *u0i;
    w1r = *u1r;
    w1i = *u1i;
    w2r = *u2r;
    w2i = *u2i;
    w3r = *(u2r + U2toU3);
    w3i = *(u2i - U2toU3);
    pstrt = ioptr;
    p0r = pstrt;
    p1r = pstrt + pinc;
    p2r = p1r + pinc;
    p3r = p2r + pinc;
    for (DiffUCnt = NDiffU; DiffUCnt > 0; DiffUCnt--) {
      f0r = *p0r;
      f0i = *(p0r + 1);
      f1r = *p1r;
      f1i = *(p1r + 1);
      for (SameUCnt = NSameU - 1; SameUCnt > 0; SameUCnt--) {
        f2r = *p2r;
        f2i = *(p2r + 1);
        f3r = *p3r;
        f3i = *(p3r + 1);
        t0r = f0r + f1r * w0r - f1i * w0i;
        t0i = f0i + f1r * w0i + f1i * w0r;
        f1r = f0r * Two - t0r;
        f1i = f0i * Two - t0i;
        f4r = *(p0r + pos);
        f4i = *(p0r + posi);
        f5r = *(p1r + pos);
        f5i = *(p1r + posi);
        f6r = *(p2r + pos);
        f6i = *(p2r + posi);
        f7r = *(p3r + pos);
        f7i = *(p3r + posi);
        t1r = f2r - f3r * w0r + f3i * w0i;
        t1i = f2i - f3r * w0i - f3i * w0r;
        f2r = f2r * Two - t1r;
        f2i = f2i * Two - t1i;
        f0r = t0r + f2r * w1r - f2i * w1i;
        f0i = t0i + f2r * w1i + f2i * w1r;
        f2r = t0r * Two - f0r;
        f2i = t0i * Two - f0i;
        f3r = f1r + t1r * w1i + t1i * w1r;
        f3i = f1i - t1r * w1r + t1i * w1i;
        f1r = f1r * Two - f3r;
        f1i = f1i * Two - f3i;
        t0r = f4r + f5r * w0r - f5i * w0i;
        t0i = f4i + f5r * w0i + f5i * w0r;
        f5r = f4r * Two - t0r;
        f5i = f4i * Two - t0i;
        t1r = f6r - f7r * w0r + f7i * w0i;
        t1i = f6i - f7r * w0i - f7i * w0r;
        f6r = f6r * Two - t1r;
        f6i = f6i * Two - t1i;
        f4r = t0r + f6r * w1r - f6i * w1i;
        f4i = t0i + f6r * w1i + f6i * w1r;
        f6r = t0r * Two - f4r;
        f6i = t0i * Two - f4i;
        f7r = f5r + t1r * w1i + t1i * w1r;
        f7i = f5i - t1r * w1r + t1i * w1i;
        f5r = f5r * Two - f7r;
        f5i = f5i * Two - f7i;
        t0r = f0r - f4r * w2r + f4i * w2i;
        t0i = f0i - f4r * w2i - f4i * w2r;
        f0r = f0r * Two - t0r;
        f0i = f0i * Two - t0i;
        t1r = f1r - f5r * w3r + f5i * w3i;
        t1i = f1i - f5r * w3i - f5i * w3r;
        f1r = f1r * Two - t1r;
        f1i = f1i * Two - t1i;
        *(p0r + pos) = t0r;
        *(p0r + posi) = t0i;
        *p0r = f0r;
        *(p0r + 1) = f0i;
        p0r += pnext;
        f0r = *p0r;
        f0i = *(p0r + 1);
        *(p1r + pos) = t1r;
        *(p1r + posi) = t1i;
        *p1r = f1r;
        *(p1r + 1) = f1i;
        p1r += pnext;
        f1r = *p1r;
        f1i = *(p1r + 1);
        f4r = f2r - f6r * w2i - f6i * w2r;
        f4i = f2i + f6r * w2r - f6i * w2i;
        f6r = f2r * Two - f4r;
        f6i = f2i * Two - f4i;
        f5r = f3r - f7r * w3i - f7i * w3r;
        f5i = f3i + f7r * w3r - f7i * w3i;
        f7r = f3r * Two - f5r;
        f7i = f3i * Two - f5i;
        *p2r = f4r;
        *(p2r + 1) = f4i;
        *(p2r + pos) = f6r;
        *(p2r + posi) = f6i;
        p2r += pnext;
        *p3r = f5r;
        *(p3r + 1) = f5i;
        *(p3r + pos) = f7r;
        *(p3r + posi) = f7i;
        p3r += pnext;
      }
      f2r = *p2r;
      f2i = *(p2r + 1);
      f3r = *p3r;
      f3i = *(p3r + 1);
      t0r = f0r + f1r * w0r - f1i * w0i;
      t0i = f0i + f1r * w0i + f1i * w0r;
      f1r = f0r * Two - t0r;
      f1i = f0i * Two - t0i;
      f4r = *(p0r + pos);
      f4i = *(p0r + posi);
      f5r = *(p1r + pos);
      f5i = *(p1r + posi);
      f6r = *(p2r + pos);
      f6i = *(p2r + posi);
      f7r = *(p3r + pos);
      f7i = *(p3r + posi);
      t1r = f2r - f3r * w0r + f3i * w0i;
      t1i = f2i - f3r * w0i - f3i * w0r;
      f2r = f2r * Two - t1r;
      f2i = f2i * Two - t1i;
      f0r = t0r + f2r * w1r - f2i * w1i;
      f0i = t0i + f2r * w1i + f2i * w1r;
      f2r = t0r * Two - f0r;
      f2i = t0i * Two - f0i;
      f3r = f1r + t1r * w1i + t1i * w1r;
      f3i = f1i - t1r * w1r + t1i * w1i;
      f1r = f1r * Two - f3r;
      f1i = f1i * Two - f3i;
      if (DiffUCnt == NDiffU / 2)
        Uinc4 = -Uinc4;
      u0r += Uinc4;
      u0i -= Uinc4;
      u1r += Uinc2;
      u1i -= Uinc2;
      u2r += Uinc;
      u2i -= Uinc;
      pstrt += 2;
      t0r = f4r + f5r * w0r - f5i * w0i;
      t0i = f4i + f5r * w0i + f5i * w0r;
      f5r = f4r * Two - t0r;
      f5i = f4i * Two - t0i;
      t1r = f6r - f7r * w0r + f7i * w0i;
      t1i = f6i - f7r * w0i - f7i * w0r;
      f6r = f6r * Two - t1r;
      f6i = f6i * Two - t1i;
      f4r = t0r + f6r * w1r - f6i * w1i;
      f4i = t0i + f6r * w1i + f6i * w1r;
      f6r = t0r * Two - f4r;
      f6i = t0i * Two - f4i;
      f7r = f5r + t1r * w1i + t1i * w1r;
      f7i = f5i - t1r * w1r + t1i * w1i;
      f5r = f5r * Two - f7r;
      f5i = f5i * Two - f7i;
      w0r = *u0r;
      w0i = *u0i;
      w1r = *u1r;
      w1i = *u1i;
      if (DiffUCnt <= NDiffU / 2)
        w0r = -w0r;
      t0r = f0r - f4r * w2r + f4i * w2i;
      t0i = f0i - f4r * w2i - f4i * w2r;
      f0r = f0r * Two - t0r;
      f0i = f0i * Two - t0i;
      f4r = f2r - f6r * w2i - f6i * w2r;
      f4i = f2i + f6r * w2r - f6i * w2i;
      f6r = f2r * Two - f4r;
      f6i = f2i * Two - f4i;
      *(p0r + pos) = t0r;
      *p2r = f4r;
      *(p0r + posi) = t0i;
      *(p2r + 1) = f4i;
      w2r = *u2r;
      w2i = *u2i;
      *p0r = f0r;
      *(p2r + pos) = f6r;
      *(p0r + 1) = f0i;
      *(p2r + posi) = f6i;
      p0r = pstrt;
      p2r = pstrt + pinc + pinc;
      t1r = f1r - f5r * w3r + f5i * w3i;
      t1i = f1i - f5r * w3i - f5i * w3r;
      f1r = f1r * Two - t1r;
      f1i = f1i * Two - t1i;
      f5r = f3r - f7r * w3i - f7i * w3r;
      f5i = f3i + f7r * w3r - f7i * w3i;
      f7r = f3r * Two - f5r;
      f7i = f3i * Two - f5i;
      *(p1r + pos) = t1r;
      *p3r = f5r;
      *(p1r + posi) = t1i;
      *(p3r + 1) = f5i;
      w3r = *(u2r + U2toU3);
      w3i = *(u2i - U2toU3);
      *p1r = f1r;
      *(p3r + pos) = f7r;
      *(p1r + 1) = f1i;
      *(p3r + posi) = f7i;
      p1r = pstrt + pinc;
      p3r = p2r + pinc;
    }
    NSameU /= 8;
    Uinc /= 8;
    Uinc2 /= 8;
    Uinc4 = Uinc * 4;
    NDiffU *= 8;
    pinc *= 8;
    pnext *= 8;
    pos *= 8;
    posi = pos + 1;
  }
}
```