user: 请从给定的汇编指令推导出C源代码:

```asm
<neuralops_batchnorm2d_bwd_input>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %r9,0x20(%rsp)
  mov    0xb8(%rsp),%r13
  mov    %rdi,0x38(%rsp)
  mov    0xc0(%rsp),%r9
  mov    %rcx,0x18(%rsp)
  movss  %xmm0,0x40(%rsp)
  test   %rdi,%rdi
  je     L84
  mov    %rsi,%rax
  mov    %rsi,%rbp
  mov    %rdx,%r14
  imul   %rdx,%rax
  imul   %rdi,%rax
  lea    -0x1(%rax),%rsi
  mov    %rsi,0x10(%rsp)
  test   %rcx,%rcx
  je     L84
  mov    %r8,%r12
  test   %rax,%rax
  js     L65
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
L143:
  lea    -0x1(%rbp),%rax
  mov    %rbp,%r10
  mov    %r14,%r8
  xor    %ebx,%ebx
  movss  D20(%rip),%xmm9
  mov    %rax,(%rsp)
  mov    %rbp,%rax
  and    $0xfffffffffffffffc,%r10
  movq   $0x0,0x30(%rsp)
  shr    $0x2,%rax
  mov    %r9,%r15
  divss  %xmm0,%xmm9
  shl    $0x4,%rax
  mov    %rax,%rcx
  mov    0x10(%rsp),%rax
  mov    %rcx,%r14
  and    $0x1,%eax
  mov    %rax,0x28(%rsp)
L188:
  movq   $0x0,0x8(%rsp)
  nopl   (%rax)
L171:
  mov    0x8(%rsp),%rax
  mov    0x20(%rsp),%rsi
  pxor   %xmm5,%xmm5
  movss  0x40(%rsp),%xmm0
  movss  (%rsi,%rax,4),%xmm3
  mov    0xa0(%rsp),%rsi
  movss  (%rsi,%rax,4),%xmm7
  mov    0xa8(%rsp),%rsi
  addss  (%rsi,%rax,4),%xmm0
  ucomiss %xmm0,%xmm5
  ja     L51
  sqrtss %xmm0,%xmm0
  mov    %rax,%rsi
L15:
  mov    0xb0(%rsp),%rax
  movss  D20(%rip),%xmm10
  movss  (%rax,%rsi,4),%xmm8
  mov    0x10(%rsp),%rax
  divss  %xmm0,%xmm10
  test   %rax,%rax
  js     L90
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
L120:
  movss  D13(%rip),%xmm13
  xor    %esi,%esi
  divss  %xmm0,%xmm13
  test   %r8,%r8
  je     L127
  movaps %xmm10,%xmm12
  movaps %xmm3,%xmm11
  movaps %xmm10,%xmm6
  mov    %rbp,%r9
  movaps %xmm3,%xmm5
  sub    %r10,%r9
  shufps $0xe0,%xmm12,%xmm12
  shufps $0xe0,%xmm11,%xmm11
  shufps $0x0,%xmm6,%xmm6
  shufps $0x0,%xmm5,%xmm5
  cs nopw 0x0(%rax,%rax,1)
L8:
  test   %rbp,%rbp
  je     L172
  lea    0x0(,%rbx,4),%rax
  movaps %xmm7,%xmm1
  movaps %xmm8,%xmm2
  lea    (%r15,%rax,1),%rdx
  lea    0x4(%rax),%rcx
  mulss  %xmm9,%xmm1
  lea    (%r12,%rcx,1),%r11
  mulss  %xmm13,%xmm2
  mov    %rdx,%rdi
  sub    %r11,%rdi
  mov    %rdx,%r11
  cmp    $0x8,%rdi
  seta   %dil
  add    %r13,%rcx
  sub    %rcx,%r11
  cmp    $0x8,%r11
  seta   %cl
  test   %cl,%dil
  je     L116
  cmp    $0x1,%rbp
  je     L116
  cmpq   $0x2,(%rsp)
  jbe    L50
  lea    0x0(%r13,%rax,1),%rdi
  movaps %xmm1,%xmm14
  movaps %xmm2,%xmm4
  add    %r12,%rax
  xor    %ecx,%ecx
  shufps $0x0,%xmm14,%xmm14
  shufps $0x0,%xmm4,%xmm4
  nopl   0x0(%rax)
L211:
  movups (%rdi,%rcx,1),%xmm0
  mulps  %xmm6,%xmm0
  addps  %xmm14,%xmm0
  movaps %xmm0,%xmm15
  movups (%rax,%rcx,1),%xmm0
  subps  %xmm5,%xmm0
  mulps  %xmm4,%xmm0
  addps  %xmm15,%xmm0
  movups %xmm0,(%rdx,%rcx,1)
  add    $0x10,%rcx
  cmp    %rcx,%r14
  jne    L211
  lea    (%r10,%rbx,1),%rax
  lea    0x0(%rbp,%rbx,1),%rdx
  cmp    %rbp,%r10
  je     L135
  mov    %r9,%rdi
  cmp    $0x1,%r9
  je     L169
  mov    %r10,%rdx
L58:
  add    %rbx,%rdx
  movaps %xmm1,%xmm0
  movaps %xmm2,%xmm14
  mov    %rdi,%rcx
  movq   0x0(%r13,%rdx,4),%xmm4
  shufps $0xe0,%xmm0,%xmm0
  shufps $0xe0,%xmm14,%xmm14
  and    $0xfffffffffffffffe,%rcx
  add    %rcx,%rax
  mulps  %xmm12,%xmm4
  addps  %xmm0,%xmm4
  movq   (%r12,%rdx,4),%xmm0
  subps  %xmm11,%xmm0
  mulps  %xmm14,%xmm0
  addps  %xmm0,%xmm4
  movlps %xmm4,(%r15,%rdx,4)
  lea    0x0(%rbp,%rbx,1),%rdx
  cmp    %rcx,%rdi
  je     L135
L169:
  movss  0x0(%r13,%rax,4),%xmm0
  lea    0x0(%rbp,%rbx,1),%rdx
  mulss  %xmm10,%xmm0
  addss  %xmm1,%xmm0
  movss  (%r12,%rax,4),%xmm1
  subss  %xmm3,%xmm1
  mulss  %xmm2,%xmm1
  addss  %xmm0,%xmm1
  movss  %xmm1,(%r15,%rax,4)
L135:
  mov    %rdx,%rbx
L172:
  add    $0x1,%rsi
  cmp    %rsi,%r8
  jne    L8
L127:
  addq   $0x1,0x8(%rsp)
  mov    0x8(%rsp),%rax
  cmp    %rax,0x18(%rsp)
  jne    L171
  addq   $0x1,0x30(%rsp)
  mov    0x30(%rsp),%rax
  cmp    %rax,0x38(%rsp)
  jne    L188
L84:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L65:
  mov    %rax,%rdx
  and    $0x1,%eax
  pxor   %xmm0,%xmm0
  shr    %rdx
  or     %rax,%rdx
  cvtsi2ss %rdx,%xmm0
  addss  %xmm0,%xmm0
  jmp    L143
  nopl   0x0(%rax)
L116:
  lea    0x0(%rbp,%rbx,1),%rdx
  nopl   (%rax)
L45:
  movss  (%r12,%rbx,4),%xmm0
  movss  0x0(%r13,%rbx,4),%xmm4
  subss  %xmm3,%xmm0
  mulss  %xmm10,%xmm4
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm4
  addss  %xmm4,%xmm0
  movss  %xmm0,(%r15,%rbx,4)
  add    $0x1,%rbx
  cmp    %rdx,%rbx
  jne    L45
  jmp    L135
L50:
  mov    %rbp,%rdi
  mov    %rbx,%rax
  xor    %edx,%edx
  jmp    L58
L90:
  shr    %rax
  pxor   %xmm0,%xmm0
  or     0x28(%rsp),%rax
  cvtsi2ss %rax,%xmm0
  addss  %xmm0,%xmm0
  jmp    L120
L51:
  mov    %r8,0x58(%rsp)
  mov    %r10,0x48(%rsp)
  movss  %xmm9,0x54(%rsp)
  movss  %xmm7,0x50(%rsp)
  movss  %xmm3,0x44(%rsp)
  call   <sqrtf@plt>
  mov    0x58(%rsp),%r8
  mov    0x48(%rsp),%r10
  movss  0x54(%rsp),%xmm9
  movss  0x50(%rsp),%xmm7
  movss  0x44(%rsp),%xmm3
  mov    0x8(%rsp),%rsi
  jmp    L15
```

tool_calls: name: read_data
arguments: data_label: D20
datatype: f32
name: read_data
arguments: data_label: D13
datatype: byte[4]

tool: D20: 1

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