user: What is the c source code of the assembly code below:

```asm
<bioR1_3dskip>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %edx,-0x54(%rbp)
  mov    %ecx,-0x58(%rbp)
  mov    %r8d,-0x5c(%rbp)
  movss  %xmm0,-0x60(%rbp)
  movss  %xmm1,-0x64(%rbp)
  mov    -0x5c(%rbp),%eax
  add    $0x1,%eax
  sar    %eax
  mov    %eax,%edx
  mov    -0x58(%rbp),%eax
  add    $0x1,%eax
  sar    %eax
  imul   %edx,%eax
  mov    -0x54(%rbp),%edx
  add    $0x1,%edx
  sar    %edx
  imul   %edx,%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x5c(%rbp),%eax
  imul   -0x58(%rbp),%eax
  imul   -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x28(%rbp)
  cmpl   $0x1,-0x54(%rbp)
  setg   %al
  movzbl %al,%edx
  cmpl   $0x1,-0x58(%rbp)
  setg   %al
  movzbl %al,%eax
  add    %eax,%edx
  cmpl   $0x1,-0x5c(%rbp)
  setg   %al
  movzbl %al,%eax
  add    %edx,%eax
  mov    %eax,-0x34(%rbp)
  cmpl   $0x3,-0x34(%rbp)
  jne    L395
  movss  -0x60(%rbp),%xmm0
  movaps %xmm0,%xmm1
  mulss  %xmm0,%xmm1
  movss  -0x64(%rbp),%xmm0
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x38(%rbp)
L395:
  cmpl   $0x2,-0x34(%rbp)
  jne    L247
  movss  -0x64(%rbp),%xmm0
  divss  -0x60(%rbp),%xmm0
  movss  %xmm0,-0x38(%rbp)
L247:
  cmpl   $0x1,-0x34(%rbp)
  jne    L348
  movss  -0x64(%rbp),%xmm0
  movss  %xmm0,-0x38(%rbp)
L348:
  cmpl   $0x0,-0x34(%rbp)
  jne    L133
  movss  -0x60(%rbp),%xmm0
  mulss  -0x64(%rbp),%xmm0
  movss  %xmm0,-0x38(%rbp)
L133:
  mov    -0x50(%rbp),%rax
  sub    $0x4,%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x5c(%rbp),%eax
  imul   -0x58(%rbp),%eax
  imul   -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x28(%rbp)
  mov    -0x28(%rbp),%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x28(%rbp),%rax
  mov    %rax,-0x20(%rbp)
  mov    -0x5c(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L154
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x18(%rbp)
  mov    -0x58(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L158
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  add    $0x1,%eax
  sar    %eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x28(%rbp),%rax
  sub    -0x8(%rbp),%rax
  mov    %eax,%edx
  mov    -0x8(%rbp),%rax
  add    $0x4,%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L129
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L129
L470:
  subq   $0x4,-0x28(%rbp)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L129:
  mov    -0x28(%rbp),%rax
  cmp    -0x8(%rbp),%rax
  ja     L470
  jmp    L294
L442:
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L294:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L442
L158:
  mov    -0x54(%rbp),%eax
  add    $0x1,%eax
  sar    %eax
  mov    %eax,%edx
  mov    -0x58(%rbp),%eax
  sar    %eax
  imul   %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x28(%rbp),%rax
  sub    -0x8(%rbp),%rax
  mov    %eax,%edx
  mov    -0x8(%rbp),%rax
  add    $0x4,%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  jmp    L38
L401:
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x28(%rbp)
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L203
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L203
L71:
  subq   $0x4,-0x28(%rbp)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L203:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L71
L38:
  mov    -0x20(%rbp),%rax
  cmp    -0x8(%rbp),%rax
  ja     L401
  jmp    L160
L156:
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x28(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%eax
  mov    %eax,%edx
  mov    -0x28(%rbp),%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L94
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L94
L23:
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L94:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L23
L160:
  mov    -0x20(%rbp),%rax
  cmp    -0x18(%rbp),%rax
  ja     L156
L154:
  mov    -0x54(%rbp),%eax
  add    $0x1,%eax
  sar    %eax
  mov    %eax,%edx
  mov    -0x58(%rbp),%eax
  add    $0x1,%eax
  sar    %eax
  imul   %edx,%eax
  mov    -0x5c(%rbp),%edx
  sar    %edx
  imul   %edx,%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x28(%rbp),%rax
  sub    -0x8(%rbp),%rax
  mov    %eax,%edx
  mov    -0x8(%rbp),%rax
  add    $0x4,%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  add    %rax,-0x8(%rbp)
  jmp    L55
L406:
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x18(%rbp)
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x28(%rbp)
  mov    -0x58(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L237
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L78
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L78
L204:
  subq   $0x4,-0x28(%rbp)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L78:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L204
  jmp    L237
L187:
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x28(%rbp)
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L420
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L420
L447:
  subq   $0x4,-0x28(%rbp)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L420:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L447
L237:
  mov    -0x20(%rbp),%rax
  cmp    -0x18(%rbp),%rax
  ja     L187
L55:
  mov    -0x18(%rbp),%rax
  cmp    -0x8(%rbp),%rax
  ja     L406
  jmp    L336
L112:
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x18(%rbp)
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x28(%rbp)
  mov    -0x54(%rbp),%eax
  imul   -0x58(%rbp),%eax
  cltq
  shl    $0x2,%eax
  mov    %eax,%edx
  mov    -0x28(%rbp),%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  mov    -0x58(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L234
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L359
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L359
L115:
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L359:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L115
L234:
  mov    -0x54(%rbp),%eax
  add    $0x1,%eax
  sar    %eax
  mov    %eax,%edx
  mov    -0x58(%rbp),%eax
  sar    %eax
  imul   %edx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x28(%rbp),%rax
  sub    -0x8(%rbp),%rax
  mov    %eax,%edx
  mov    -0x8(%rbp),%rax
  add    $0x4,%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  add    %rax,-0x8(%rbp)
  jmp    L392
L225:
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x28(%rbp)
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L293
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L293
L367:
  subq   $0x4,-0x28(%rbp)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L293:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L367
L392:
  mov    -0x20(%rbp),%rax
  cmp    -0x8(%rbp),%rax
  ja     L225
  jmp    L435
L45:
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x20(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%rax
  neg    %rax
  add    %rax,-0x28(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  shl    $0x2,%eax
  mov    %eax,%edx
  mov    -0x28(%rbp),%rax
  mov    $0x0,%esi
  mov    %rax,%rdi
  call   <memset@plt>
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L302
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L302
L40:
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x30(%rbp)
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  lea    -0x4(%rax),%rdx
  mov    %rdx,-0x28(%rbp)
  mulss  -0x38(%rbp),%xmm0
  movss  %xmm0,(%rax)
L302:
  mov    -0x28(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  ja     L40
L435:
  mov    -0x20(%rbp),%rax
  cmp    -0x18(%rbp),%rax
  ja     L45
L336:
  mov    -0x18(%rbp),%rax
  cmp    -0x10(%rbp),%rax
  ja     L112
  nop
  nop
  leave
  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--);
      }
    }
  }
}
```