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

```asm
<dct_inverse>:
  endbr64
  push   %r15
  push   %r14
  mov    %rdi,%r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rsi,0x28(%rsp)
  movl   $0x0,0x24(%rsp)
L73:
  pxor   %xmm7,%xmm7
  mov    0x28(%rsp),%rbx
  xor    %r13d,%r13d
  cvtsi2ssl 0x24(%rsp),%xmm7
  movss  %xmm7,0x20(%rsp)
  nopw   0x0(%rax,%rax,1)
L8:
  pxor   %xmm7,%xmm7
  movl   $0x0,(%rbx)
  mov    %r14,%rbp
  xor    %r12d,%r12d
  cvtsi2ss %r13d,%xmm7
  movss  %xmm7,0x1c(%rsp)
  nopl   0x0(%rax,%rax,1)
L22:
  movl   $0x3f800000,0x14(%rsp)
  test   %r12d,%r12d
  je     L29
L45:
  pxor   %xmm7,%xmm7
  xor    %r15d,%r15d
  cvtsi2ss %r12d,%xmm7
  movss  %xmm7,0x18(%rsp)
  jmp    L88
  nopl   0x0(%rax,%rax,1)
L47:
  movss  0x20(%rsp),%xmm3
  movss  %xmm1,0x10(%rsp)
  movaps %xmm3,%xmm0
  addss  %xmm3,%xmm0
  addss  D94(%rip),%xmm0
  mulss  D39(%rip),%xmm0
  mulss  0x18(%rsp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  D31(%rip),%xmm0
  cvtsd2ss %xmm0,%xmm0
  call   <cosf@plt>
  movss  0x1c(%rsp),%xmm5
  pxor   %xmm2,%xmm2
  movss  %xmm0,0xc(%rsp)
  movss  D94(%rip),%xmm0
  cvtsi2ss %r15d,%xmm2
  movaps %xmm5,%xmm4
  addss  %xmm5,%xmm4
  addss  %xmm4,%xmm0
  mulss  D39(%rip),%xmm0
  mulss  %xmm2,%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  D31(%rip),%xmm0
  cvtsd2ss %xmm0,%xmm0
  call   <cosf@plt>
  movss  0x10(%rsp),%xmm1
  movaps %xmm0,%xmm2
  pxor   %xmm0,%xmm0
  cvtsi2ssl 0x0(%rbp,%r15,4),%xmm0
  mulss  0xc(%rsp),%xmm0
  add    $0x1,%r15
  mulss  %xmm2,%xmm0
  mulss  0x14(%rsp),%xmm0
  mulss  %xmm1,%xmm0
  addss  (%rbx),%xmm0
  movss  %xmm0,(%rbx)
  cmp    $0x8,%r15
  je     L34
L88:
  movss  D94(%rip),%xmm1
  test   %r15,%r15
  jne    L47
  mov    $0x2,%edi
  call   <sqrt@plt>
  movss  D94(%rip),%xmm1
  divss  %xmm0,%xmm1
  jmp    L47
  nopl   (%rax)
L34:
  add    $0x1,%r12d
  add    $0x20,%rbp
  cmp    $0x8,%r12d
  jne    L22
  mulss  D75(%rip),%xmm0
  add    $0x1,%r13d
  add    $0x4,%rbx
  addss  D8(%rip),%xmm0
  movss  %xmm0,-0x4(%rbx)
  cmp    $0x8,%r13d
  jne    L8
  addl   $0x1,0x24(%rsp)
  mov    0x24(%rsp),%eax
  addq   $0x20,0x28(%rsp)
  cmp    $0x8,%eax
  jne    L73
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L29:
  mov    $0x2,%edi
  call   <sqrt@plt>
  movss  D94(%rip),%xmm7
  divss  %xmm0,%xmm7
  movss  %xmm7,0x14(%rsp)
  jmp    L45
```

tool_calls: name: unpack_data
arguments: label: D94
datatype: byte[4]
name: unpack_data
arguments: label: D39
datatype: byte[4]
name: unpack_data
arguments: label: D31
datatype: f64
name: unpack_data
arguments: label: D75
datatype: byte[4]
name: unpack_data
arguments: label: D8
datatype: byte[4]

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

tool: {"D39": [0x00,0x00,0x80,0x3d]}

tool: D31: 3.141592653589793

tool: {"D75": [0x00,0x00,0x80,0x3e]}

tool: {"D8": [0x00,0x00,0x00,0x43]}

assistant: ```c
void dct_inverse(int chunk[(8 * 8)], float pixels[8][8]) {
  unsigned int x, y;
  for (x = 0; x < 8; x++) {
    for (y = 0; y < 8; y++) {
      unsigned int u, v;
      pixels[x][y] = 0.0f;
      for (u = 0; u < 8; u++) {
        float u_scaling_factor = 1.0f;
        if (u == 0) {
          u_scaling_factor = 1.0f / sqrt(2);
        }
        for (v = 0; v < 8; v++) {
          float v_scaling_factor = 1.0f;
          if (v == 0) {
            v_scaling_factor = 1.0f / sqrt(2);
          }
          pixels[x][y] += (float)chunk[u * 8 + v] * cosf((2.0f * (float)x + 1.0f) / 16.0f * (float)u * 3.14159265358979323846264338327) * cosf((2.0f * (float)y + 1.0f) / 16.0f * (float)v * 3.14159265358979323846264338327) * u_scaling_factor * v_scaling_factor;
        }
      }
      pixels[x][y] *= 0.25f;
      pixels[x][y] += 128.0f;
    }
  }
}
```