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

```asm
<neuralops_batchnorm2d_bwd_input>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %rsi,-0x60(%rbp)
  mov    %rdx,-0x68(%rbp)
  mov    %rcx,-0x70(%rbp)
  mov    %r8,-0x78(%rbp)
  mov    %r9,-0x80(%rbp)
  movss  %xmm0,-0x84(%rbp)
  movq   $0x0,-0x28(%rbp)
  movq   $0x0,-0x20(%rbp)
  jmp    L63
L94:
  movq   $0x0,-0x18(%rbp)
  jmp    L72
L17:
  mov    -0x18(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x44(%rbp)
  mov    -0x18(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    0x18(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x40(%rbp)
  mov    -0x18(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x3c(%rbp)
  mov    -0x60(%rbp),%rax
  imul   -0x68(%rbp),%rax
  imul   -0x58(%rbp),%rax
  test   %rax,%rax
  js     L88
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  jmp    L96
L88:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2ss %rdx,%xmm0
  addss  %xmm0,%xmm0
L96:
  movss  D13(%rip),%xmm1
  divss  %xmm0,%xmm1
  movss  %xmm1,-0x38(%rbp)
  movss  -0x40(%rbp),%xmm0
  addss  -0x84(%rbp),%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <sqrtf@plt>
  movaps %xmm0,%xmm1
  movss  D13(%rip),%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x34(%rbp)
  mov    -0x18(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    0x20(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x30(%rbp)
  mov    -0x60(%rbp),%rax
  imul   -0x68(%rbp),%rax
  imul   -0x58(%rbp),%rax
  sub    $0x1,%rax
  test   %rax,%rax
  js     L70
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  jmp    L79
L70:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2ss %rdx,%xmm0
  addss  %xmm0,%xmm0
L79:
  movss  D28(%rip),%xmm1
  divss  %xmm0,%xmm1
  movss  %xmm1,-0x2c(%rbp)
  movq   $0x0,-0x10(%rbp)
  jmp    L26
L53:
  movq   $0x0,-0x8(%rbp)
  jmp    L100
L49:
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    0x28(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x34(%rbp),%xmm1
  movss  -0x3c(%rbp),%xmm0
  mulss  -0x38(%rbp),%xmm0
  addss  %xmm0,%xmm1
  movss  -0x30(%rbp),%xmm0
  movaps %xmm0,%xmm2
  mulss  -0x2c(%rbp),%xmm2
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  subss  -0x44(%rbp),%xmm0
  mulss  %xmm2,%xmm0
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    0x30(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addq   $0x1,-0x28(%rbp)
  addq   $0x1,-0x8(%rbp)
L100:
  mov    -0x8(%rbp),%rax
  cmp    -0x60(%rbp),%rax
  jb     L49
  addq   $0x1,-0x10(%rbp)
L26:
  mov    -0x10(%rbp),%rax
  cmp    -0x68(%rbp),%rax
  jb     L53
  addq   $0x1,-0x18(%rbp)
L72:
  mov    -0x18(%rbp),%rax
  cmp    -0x70(%rbp),%rax
  jb     L17
  addq   $0x1,-0x20(%rbp)
L63:
  mov    -0x20(%rbp),%rax
  cmp    -0x58(%rbp),%rax
  jb     L94
  nop
  nop
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D13
type_guessed: f32
name: disassemble_data
arguments: tag: D28
type_guessed: byte[4]

tool: D13: 1

tool: {"D28": [0x00,0x00,0x00,0x40]}

assistant: ```c
void neuralops_batchnorm2d_bwd_input(size_t batch_sz, size_t width, size_t height, size_t chan, const float *in_buf, const float *mean, const float *mean_grad, const float *var, const float *var_grad, const float *out_grad, float *in_grad, float epsilon) {
  size_t p = 0;
  for (size_t idx = 0; idx < batch_sz; idx += 1) {
    for (size_t a = 0; a < chan; a += 1) {
      float m = mean[a];
      float v = var[a];
      float dm = mean_grad[a];
      float cm = 1.0f / ((float)(width * height * batch_sz));
      float rs = 1.0f / sqrtf(v + epsilon);
      float dv = var_grad[a];
      float cv = 2.0f / ((float)(width * height * batch_sz - 1));
      for (size_t y = 0; y < height; y += 1) {
        for (size_t x = 0; x < width; x += 1) {
          in_grad[p] = out_grad[p] * rs + dm * cm + dv * cv * (in_buf[p] - m);
          p += 1;
        }
      }
    }
  }
}
```