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

```asm
<dct_inverse>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rdi,%r15
  mov    %rsi,0x18(%rsp)
  movl   $0x0,0xc(%rsp)
  jmp    L41
L54:
  mov    $0x2,%edi
  call   <sqrt@plt>
  movss  D67(%rip),%xmm1
  divss  %xmm0,%xmm1
  movss  %xmm1,0x8(%rsp)
  jmp    L88
L37:
  mov    0xc(%rsp),%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  addss  %xmm0,%xmm0
  addss  D67(%rip),%xmm0
  mulss  D102(%rip),%xmm0
  mov    %r12d,%eax
  pxor   %xmm1,%xmm1
  cvtsi2ss %rax,%xmm1
  mulss  %xmm1,%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  D119(%rip),%xmm0
  cvtsd2ss %xmm0,%xmm0
  call   <cosf@plt>
  movss  %xmm0,0x4(%rsp)
  mov    %r14d,%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  addss  %xmm0,%xmm0
  addss  D67(%rip),%xmm0
  mulss  D102(%rip),%xmm0
  mov    %ebx,%eax
  pxor   %xmm1,%xmm1
  cvtsi2ss %rax,%xmm1
  mulss  %xmm1,%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  D119(%rip),%xmm0
  cvtsd2ss %xmm0,%xmm0
  call   <cosf@plt>
  movaps %xmm0,%xmm1
  lea    0x0(%r13,%rbx,1),%eax
  pxor   %xmm0,%xmm0
  cvtsi2ssl (%r15,%rax,4),%xmm0
  mulss  0x4(%rsp),%xmm0
  mulss  %xmm1,%xmm0
  mulss  0x8(%rsp),%xmm0
  mulss  (%rsp),%xmm0
  addss  0x0(%rbp),%xmm0
  movss  %xmm0,0x0(%rbp)
  add    $0x1,%ebx
  cmp    $0x8,%ebx
  je     L32
L35:
  movss  D67(%rip),%xmm2
  movss  %xmm2,(%rsp)
  test   %ebx,%ebx
  jne    L37
  mov    $0x2,%edi
  call   <sqrt@plt>
  movss  D67(%rip),%xmm6
  divss  %xmm0,%xmm6
  movss  %xmm6,(%rsp)
  jmp    L37
L32:
  add    $0x1,%r12d
  cmp    $0x8,%r12d
  je     L64
L98:
  movss  D67(%rip),%xmm7
  movss  %xmm7,0x8(%rsp)
  test   %r12d,%r12d
  je     L54
L88:
  lea    0x0(,%r12,8),%r13d
  mov    $0x0,%ebx
  jmp    L35
L64:
  mulss  D29(%rip),%xmm0
  addss  D90(%rip),%xmm0
  movss  %xmm0,0x0(%rbp)
  add    $0x1,%r14d
  addq   $0x4,0x10(%rsp)
  cmp    $0x8,%r14d
  je     L13
L8:
  mov    0x10(%rsp),%rax
  mov    %rax,%rbp
  movl   $0x0,(%rax)
  mov    $0x0,%r12d
  jmp    L98
L13:
  addl   $0x1,0xc(%rsp)
  mov    0xc(%rsp),%eax
  addq   $0x20,0x18(%rsp)
  cmp    $0x8,%eax
  je     L7
L41:
  mov    0x18(%rsp),%rax
  mov    %rax,0x10(%rsp)
  mov    $0x0,%r14d
  jmp    L8
L7:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: data_label: D67
datatype: byte[4]
name: parse_data
arguments: data_label: D102
datatype: byte[4]
name: parse_data
arguments: data_label: D119
datatype: f64
name: parse_data
arguments: data_label: D29
datatype: byte[4]
name: parse_data
arguments: data_label: D90
datatype: byte[4]

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

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

tool: D119: 3.141592653589793

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

tool: {"D90": [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;
    }
  }
}
```