user: What would the C source look like for the assembly code provided here:

```asm
<psnd_paperplot_grid>:
  endbr64
  push   %r13
  movaps %xmm0,%xmm6
  push   %r12
  mov    %edi,%r12d
  push   %rbp
  push   %rbx
  mov    %edx,%ebx
  sub    $0x28,%rsp
  test   %esi,%esi
  jne    L77
  test   %ebx,%ebx
  jne    L14
L5:
  mulss  %xmm5,%xmm2
  movaps %xmm6,%xmm7
  mulss  %xmm4,%xmm7
  mulss  %xmm1,%xmm4
  mulss  %xmm3,%xmm5
  movss  %xmm2,0x18(%rsp)
  movaps %xmm7,%xmm0
  movss  %xmm7,0x10(%rsp)
  maxss  %xmm4,%xmm0
  movss  %xmm4,0x14(%rsp)
  movss  %xmm5,0x1c(%rsp)
  call   <ceil@plt>
  movss  0x10(%rsp),%xmm6
  minss  0x14(%rsp),%xmm6
  mov    %eax,%ebx
  movaps %xmm6,%xmm0
  call   <floor@plt>
  movss  0x1c(%rsp),%xmm3
  mov    %r12d,%edi
  movss  0x14(%rsp),%xmm2
  movss  0x18(%rsp),%xmm1
  movss  0x10(%rsp),%xmm0
  mov    %eax,%ebp
  call   <g_set_world@plt>
  cmp    %ebp,%ebx
  jl     L64
  cvttss2si 0x18(%rsp),%r13d
  cvttss2si 0x1c(%rsp),%r12d
  nopw   0x0(%rax,%rax,1)
L37:
  pxor   %xmm0,%xmm0
  mov    %r13d,%edi
  cvtsi2ss %ebx,%xmm0
  sub    $0x1,%ebx
  movss  %xmm0,0xc(%rsp)
  call   <g_moveto@plt>
  movss  0xc(%rsp),%xmm0
  mov    %r12d,%edi
  call   <g_lineto@plt>
  cmp    %ebx,%ebp
  jle    L37
L64:
  movss  0x18(%rsp),%xmm0
  maxss  0x1c(%rsp),%xmm0
  call   <ceil@plt>
  movss  0x18(%rsp),%xmm0
  minss  0x1c(%rsp),%xmm0
  mov    %eax,%ebx
  call   <floor@plt>
  mov    %eax,%ebp
  cmp    %eax,%ebx
  jl     L75
  nopl   (%rax)
L33:
  movss  0x10(%rsp),%xmm0
  mov    %ebx,%edi
  call   <g_moveto@plt>
  movss  0x14(%rsp),%xmm0
  mov    %ebx,%edi
  sub    $0x1,%ebx
  call   <g_lineto@plt>
  cmp    %ebx,%ebp
  jle    L33
L75:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
  nopl   0x0(%rax)
L77:
  movss  %xmm0,0x10(%rsp)
  subss  %xmm1,%xmm0
  movss  %xmm5,0x1c(%rsp)
  movss  %xmm3,0x18(%rsp)
  movss  %xmm2,0x14(%rsp)
  movss  %xmm1,0xc(%rsp)
  call   <fabs@plt>
  call   <get_long_ticks@plt>
  movss  0x1c(%rsp),%xmm5
  movss  0x18(%rsp),%xmm3
  movss  0x14(%rsp),%xmm2
  movss  0x10(%rsp),%xmm6
  movaps %xmm0,%xmm4
  movss  0xc(%rsp),%xmm1
  test   %ebx,%ebx
  je     L5
L14:
  movaps %xmm2,%xmm0
  movss  %xmm4,0x1c(%rsp)
  subss  %xmm3,%xmm0
  movss  %xmm1,0x18(%rsp)
  movss  %xmm6,0x14(%rsp)
  movss  %xmm2,0x10(%rsp)
  movss  %xmm3,0xc(%rsp)
  call   <fabs@plt>
  call   <get_long_ticks@plt>
  movss  0x1c(%rsp),%xmm4
  movss  0x18(%rsp),%xmm1
  movss  0x14(%rsp),%xmm6
  movss  0x10(%rsp),%xmm2
  movaps %xmm0,%xmm5
  movss  0xc(%rsp),%xmm3
  jmp    L5
```

assistant: ```c
void psnd_paperplot_grid(int vp, int autoscalex, int autoscaley, float xminp, float xmaxp, float yminp, float ymaxp, float xlines_ppm, float ylines_ppm) {
  float xtmp1, xtmp2, ytmp1, ytmp2;
  int n, start, stop;
  if (autoscalex) {
    float len = fabs(xminp - xmaxp);
    xlines_ppm = get_long_ticks(len);
  }
  if (autoscaley) {
    float len = fabs(yminp - ymaxp);
    ylines_ppm = get_long_ticks(len);
  }
  xtmp1 = xminp * xlines_ppm;
  xtmp2 = xmaxp * xlines_ppm;
  ytmp1 = yminp * ylines_ppm;
  ytmp2 = ymaxp * ylines_ppm;
  start = ceil((((xtmp1) > (xtmp2)) ? (xtmp1) : (xtmp2)));
  stop = floor((((xtmp1) < (xtmp2)) ? (xtmp1) : (xtmp2)));
  g_set_world(vp, xtmp1, ytmp1, xtmp2, ytmp2);
  for (n = start; n >= stop; n--) {
    g_moveto(n, ytmp1);
    g_lineto(n, ytmp2);
  }
  start = ceil((((ytmp1) > (ytmp2)) ? (ytmp1) : (ytmp2)));
  stop = floor((((ytmp1) < (ytmp2)) ? (ytmp1) : (ytmp2)));
  for (n = start; n >= stop; n--) {
    g_moveto(xtmp1, n);
    g_lineto(xtmp2, n);
  }
}
```