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

```asm
<bioR1_3dskip>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rsi,0x28(%rsp)
  mov    %edx,%ebx
  mov    %edx,0x1c(%rsp)
  mov    %ecx,0x4c(%rsp)
  mov    %r8d,0x38(%rsp)
  movss  %xmm1,0x18(%rsp)
  lea    0x1(%rcx),%eax
  sar    %eax
  mov    %eax,%edx
  mov    %eax,0x48(%rsp)
  lea    0x1(%rbx),%eax
  sar    %eax
  mov    %eax,%r10d
  mov    %eax,0x20(%rsp)
  lea    0x1(%r8),%eax
  sar    %eax
  imul   %edx,%eax
  imul   %r10d,%eax
  cltq
  lea    -0x4(%rdi,%rax,4),%r12
  mov    %r8d,%eax
  imul   %ecx,%eax
  imul   %ebx,%eax
  cltq
  lea    -0x4(%rsi,%rax,4),%rbp
  cmp    $0x1,%ebx
  setg   %al
  movzbl %al,%eax
  cmp    $0x1,%ecx
  setg   %dl
  movzbl %dl,%edx
  add    %edx,%eax
  cmp    $0x1,%r8d
  setg   %dl
  movzbl %dl,%edx
  add    %edx,%eax
  cmp    $0x3,%eax
  je     L135
  cmp    $0x2,%eax
  jne    L651
  movss  0x18(%rsp),%xmm3
  divss  %xmm0,%xmm3
  movss  %xmm3,0x18(%rsp)
  jmp    L104
L135:
  movaps %xmm0,%xmm1
  mulss  %xmm0,%xmm1
  movss  0x18(%rsp),%xmm2
  divss  %xmm1,%xmm2
L306:
  test   %eax,%eax
  jne    L620
  mulss  0x18(%rsp),%xmm0
  movss  %xmm0,0x18(%rsp)
L104:
  mov    0x28(%rsp),%rax
  sub    $0x4,%rax
  mov    %rax,0x30(%rsp)
  testb  $0x1,0x38(%rsp)
  je     L516
  mov    0x4c(%rsp),%edi
  mov    %edi,%eax
  imul   0x1c(%rsp),%eax
  cltq
  shl    $0x2,%rax
  mov    %rbp,%r13
  sub    %rax,%r13
  mov    %rbp,%rbx
  test   $0x1,%dil
  jne    L658
L3:
  mov    0x4c(%rsp),%eax
  sar    %eax
  imul   0x20(%rsp),%eax
  cltq
  lea    0x0(%r13,%rax,4),%r14
  lea    0x4(%r14),%rdi
  mov    %ebp,%edx
  sub    %r14d,%edx
  mov    $0x0,%esi
  call   <memset@plt>
  cmp    %r14,%rbx
  jbe    L465
  mov    0x1c(%rsp),%eax
  lea    (%rax,%rax,1),%ecx
  movslq %ecx,%rcx
  shl    $0x2,%rcx
  neg    %rcx
  movslq %eax,%rsi
  shl    $0x2,%rsi
  neg    %rsi
  and    $0x1,%eax
  mov    %eax,%edi
  jmp    L71
L651:
  cmp    $0x1,%eax
  je     L104
  jmp    L306
L620:
  movss  %xmm2,0x18(%rsp)
  jmp    L104
L658:
  mov    0x1c(%rsp),%r15d
  movslq %r15d,%rax
  shl    $0x2,%rax
  sub    %rax,%rbx
  movslq 0x20(%rsp),%rax
  lea    (%rbx,%rax,4),%r14
  lea    0x4(%r14),%rdi
  mov    %ebp,%edx
  sub    %r14d,%edx
  mov    $0x0,%esi
  call   <memset@plt>
  test   $0x1,%r15b
  je     L246
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
L246:
  cmp    %rbp,%r14
  jae    L570
  mov    %rbp,%rax
  mov    %r12,%rdx
L405:
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%r14
  jb     L405
  mov    %rbp,%rax
  sub    %r14,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
L570:
  cmp    %rbp,%rbx
  jae    L3
  mov    %rbp,%rax
  mov    %r12,%rdx
L62:
  movl   $0x0,(%rax)
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L62
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
  jmp    L3
L183:
  cmp    %rbp,%rbx
  jae    L150
  mov    %rbp,%rax
  mov    %r12,%rdx
L521:
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L521
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
L150:
  cmp    %rbx,%r14
  jae    L465
L71:
  add    %rcx,%rbx
  add    %rsi,%rbp
  test   %edi,%edi
  je     L183
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
  jmp    L183
L465:
  cmp    %rbx,%r13
  jae    L496
  mov    0x1c(%rsp),%edi
  lea    (%rdi,%rdi,1),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  mov    %rax,0x8(%rsp)
  movslq %edi,%rax
  shl    $0x2,%rax
  neg    %rax
  mov    %rax,0x10(%rsp)
  lea    0x0(,%rdi,4),%r15d
  and    $0x1,%edi
  mov    %edi,%r14d
  jmp    L346
L408:
  cmp    %rbp,%rbx
  jae    L197
  mov    %rbp,%rax
  mov    %r12,%rdx
L78:
  movl   $0x0,(%rax)
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L78
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
L197:
  cmp    %rbx,%r13
  jae    L496
L346:
  add    0x8(%rsp),%rbx
  add    0x10(%rsp),%rbp
  mov    %r15d,%edx
  mov    $0x0,%esi
  mov    %rbp,%rdi
  call   <memset@plt>
  test   %r14d,%r14d
  je     L408
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
  jmp    L408
L516:
  mov    %rbp,%r13
  mov    %rbp,%rbx
L496:
  mov    0x48(%rsp),%eax
  imul   0x20(%rsp),%eax
  mov    0x38(%rsp),%edx
  sar    %edx
  imul   %edx,%eax
  cltq
  shl    $0x2,%rax
  mov    0x28(%rsp),%rdi
  lea    -0x4(%rdi,%rax,1),%r14
  add    %rax,%rdi
  mov    %ebp,%edx
  sub    %r14d,%edx
  mov    $0x0,%esi
  call   <memset@plt>
  mov    0x4c(%rsp),%edx
  mov    %edx,%r10d
  mov    0x1c(%rsp),%edi
  imul   %edi,%r10d
  lea    (%r10,%r10,1),%eax
  cltq
  lea    0x0(,%rax,4),%r11
  mov    %r11,0x28(%rsp)
  add    %r11,%r14
  cmp    %r14,%r13
  jbe    L87
  movslq %r10d,%r8
  shl    $0x2,%r8
  neg    %r8
  movslq %edi,%rcx
  shl    $0x2,%rcx
  neg    %rcx
  lea    (%rdi,%rdi,1),%esi
  movslq %esi,%rsi
  shl    $0x2,%rsi
  neg    %rsi
  and    $0x1,%edx
  mov    %edx,%r9d
  and    $0x1,%edi
  jmp    L88
L81:
  add    %rcx,%rbx
  test   %edi,%edi
  je     L324
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
L324:
  cmp    %rbp,%rbx
  jae    L450
  mov    %rbp,%rax
  mov    %r12,%rdx
L489:
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L489
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
  jmp    L450
L6:
  cmp    %rbp,%rbx
  jae    L38
  mov    %rbp,%rax
  mov    %r12,%rdx
L241:
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L241
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
L38:
  cmp    %rbx,%r13
  jae    L482
L349:
  add    %rsi,%rbx
  add    %rcx,%rbp
  test   %edi,%edi
  je     L6
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
  jmp    L6
L482:
  cmp    %r13,%r14
  jae    L87
L88:
  sub    %r11,%r13
  add    %r8,%rbx
  add    %r8,%rbp
  test   %r9d,%r9d
  jne    L81
L450:
  cmp    %rbx,%r13
  jb     L349
  jmp    L482
L87:
  cmp    %r13,0x30(%rsp)
  jae    L488
  movslq %r10d,%rax
  shl    $0x2,%rax
  neg    %rax
  mov    %rax,0x38(%rsp)
  lea    0x0(,%r10,4),%eax
  mov    %eax,0x48(%rsp)
  mov    0x1c(%rsp),%edi
  movslq %edi,%rax
  shl    $0x2,%rax
  neg    %rax
  mov    %rax,0x10(%rsp)
  lea    (%rdi,%rdi,1),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,0x8(%rsp)
  lea    0x0(,%rdi,4),%r15d
  mov    0x28(%rsp),%rcx
  mov    %rcx,%rax
  neg    %rax
  mov    %rax,0x40(%rsp)
  mov    0x4c(%rsp),%edi
  mov    %edi,%eax
  sar    %eax
  imul   0x20(%rsp),%eax
  cltq
  lea    0x0(%r13,%rax,4),%rax
  sub    %rcx,%rax
  mov    %rax,0x20(%rsp)
  and    $0x1,%edi
  mov    %edi,0x4c(%rsp)
  jmp    L360
L159:
  add    0x10(%rsp),%rbx
  testb  $0x1,0x1c(%rsp)
  je     L649
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
L649:
  cmp    %rbp,%rbx
  jae    L111
  mov    %rbp,%rax
  mov    %r12,%rdx
L565:
  movl   $0x0,(%rax)
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L565
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
  jmp    L111
L125:
  cmp    %rbp,%rbx
  jae    L601
  mov    %rbp,%rax
  mov    %r12,%rdx
L523:
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L523
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
L601:
  cmp    %rcx,%rbx
  jbe    L312
L537:
  sub    0x8(%rsp),%rbx
  add    0x10(%rsp),%rbp
  test   %esi,%esi
  je     L125
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
  jmp    L125
L312:
  cmp    %rbx,%r13
  jae    L409
  mov    0x1c(%rsp),%r14d
  and    $0x1,%r14d
  jmp    L506
L255:
  cmp    %rbp,%rbx
  jae    L590
  mov    %rbp,%rax
  mov    %r12,%rdx
L14:
  movl   $0x0,(%rax)
  sub    $0x4,%rdx
  sub    $0x8,%rax
  movss  0x18(%rsp),%xmm0
  mulss  0x4(%rdx),%xmm0
  movss  %xmm0,0x4(%rax)
  cmp    %rax,%rbx
  jb     L14
  mov    %rbp,%rax
  sub    %rbx,%rax
  sub    $0x1,%rax
  shr    $0x3,%rax
  mov    %rax,%rdx
  not    %rdx
  lea    (%r12,%rdx,4),%r12
  lea    0x0(%rbp,%rdx,8),%rbp
L590:
  cmp    %rbx,%r13
  jae    L409
L506:
  sub    0x8(%rsp),%rbx
  add    0x10(%rsp),%rbp
  mov    %r15d,%edx
  mov    $0x0,%esi
  mov    %rbp,%rdi
  call   <memset@plt>
  test   %r14d,%r14d
  je     L255
  movss  0x18(%rsp),%xmm0
  mulss  (%r12),%xmm0
  movss  %xmm0,0x0(%rbp)
  lea    -0x4(%rbp),%rbp
  lea    -0x4(%r12),%r12
  jmp    L255
L409:
  mov    0x28(%rsp),%rdi
  sub    %rdi,0x20(%rsp)
  cmp    %r13,0x30(%rsp)
  jae    L488
L360:
  add    0x40(%rsp),%r13
  mov    0x38(%rsp),%rax
  add    %rax,%rbx
  add    %rax,%rbp
  mov    0x48(%rsp),%edx
  mov    $0x0,%esi
  mov    %rbp,%rdi
  call   <memset@plt>
  cmpl   $0x0,0x4c(%rsp)
  jne    L159
L111:
  mov    0x20(%rsp),%r14
  lea    0x4(%r14),%rdi
  mov    %ebp,%edx
  sub    %r14d,%edx
  mov    $0x0,%esi
  call   <memset@plt>
  mov    0x8(%rsp),%rcx
  add    %r14,%rcx
  cmp    %rcx,%rbx
  jbe    L312
  mov    0x1c(%rsp),%esi
  and    $0x1,%esi
  jmp    L537
L488:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void bioR1_3dskip(float *LLL, float *vector, int xlength, int ylength, int zlength, float X, float norm) {
  float *LLLptr;
  float *vptr, *vptrx0, *vptry0, *vptrz0, *vptrmem;
  float x_2;
  int dim;
  LLLptr = LLL + ((zlength + 1) >> 1) * ((ylength + 1) >> 1) * ((xlength + 1) >> 1) - 1;
  vptr = vector + zlength * ylength * xlength - 1;
  dim = (xlength > 1) + (ylength > 1) + (zlength > 1);
  if (dim == 3) {
    x_2 = norm / (X * X);
  }
  if (dim == 2) {
    x_2 = norm / X;
  }
  if (dim == 1) {
    x_2 = norm;
  }
  if (dim == 0) {
    x_2 = X * norm;
  }
  vptrz0 = vector - 1;
  vptr = vptrz0 + zlength * ylength * xlength;
  vptry0 = vptr;
  vptrx0 = vptr;
  if (1 & zlength) {
    vptry0 -= xlength * ylength;
    if (1 & ylength) {
      vptrx0 -= xlength;
      vptrmem = vptrx0 + ((xlength + 1) >> 1);
      memset(vptrmem + 1, 0, (vptr - vptrmem) * sizeof(float));
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrmem) {
        vptr--;
        *(vptr--) = x_2 * *(LLLptr--);
      }
      while (vptr > vptrx0) {
        *vptr-- = 0;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
    vptrmem = vptry0 + ((xlength + 1) >> 1) * (ylength >> 1);
    memset(vptrmem + 1, 0, (vptr - vptrmem) * sizeof(float));
    while (vptrx0 > vptrmem) {
      vptrx0 -= 2 * xlength;
      vptr -= xlength;
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrx0) {
        vptr--;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
    while (vptrx0 > vptry0) {
      vptrx0 -= 2 * xlength;
      vptr -= xlength;
      memset(vptr, 0, xlength * sizeof(float));
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrx0) {
        *vptr-- = 0;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
  }
  vptrmem = vector + ((xlength + 1) >> 1) * ((ylength + 1) >> 1) * (zlength >> 1) - 1;
  memset(vptrmem + 1, 0, (vptr - vptrmem) * sizeof(float));
  vptrmem += 2 * xlength * ylength;
  while (vptry0 > vptrmem) {
    vptry0 -= 2 * xlength * ylength;
    vptrx0 -= xlength * ylength;
    vptr -= xlength * ylength;
    if (1 & ylength) {
      vptrx0 -= xlength;
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrx0) {
        vptr--;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
    while (vptrx0 > vptry0) {
      vptrx0 -= 2 * xlength;
      vptr -= xlength;
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrx0) {
        vptr--;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
  }
  while (vptry0 > vptrz0) {
    vptry0 -= 2 * xlength * ylength;
    vptrx0 -= xlength * ylength;
    vptr -= xlength * ylength;
    memset(vptr, 0, xlength * ylength * sizeof(float));
    if (1 & ylength) {
      vptrx0 -= xlength;
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrx0) {
        *vptr-- = 0;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
    vptrmem = vptry0 + ((xlength + 1) >> 1) * (ylength >> 1);
    memset(vptrmem + 1, 0, (vptr - vptrmem) * sizeof(float));
    vptrmem += 2 * xlength;
    while (vptrx0 > vptrmem) {
      vptrx0 -= 2 * xlength;
      vptr -= xlength;
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrx0) {
        vptr--;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
    while (vptrx0 > vptry0) {
      vptrx0 -= 2 * xlength;
      vptr -= xlength;
      memset(vptr, 0, xlength * sizeof(float));
      if (1 & xlength)
        *(vptr--) = x_2 * *(LLLptr--);
      while (vptr > vptrx0) {
        *vptr-- = 0;
        *(vptr--) = x_2 * *(LLLptr--);
      }
    }
  }
}
```