void bioR5_3d(HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, outvector, xlength,
              ylength, zlength, ifnotAllskip) float *HLL,
    *HLH, *HHL, *HHH;
float *LLL, *LLH, *LHL, *LHH;
float *outvector;
int xlength, ylength;
int ifnotAllskip;
{
  float ix0 = 0.50;
  float ix1 = -0.25;
  float iX0, iX1, iX2;
  int pairity;
  float Norm;
  int dim;
  dim = (xlength > 1) + (ylength > 1) + (zlength > 1);
  if (ifnotAllskip) {
    iX2 = ix1;
    iX1 = 1.0 / (ix1 * ix0);
    iX0 = 1.0 / ix0;
  } else {
    iX2 = ix0;
    iX1 = 1.0;
    iX0 = 1.0 / ix0;
  }
  if (dim == 3)
    Norm = 1.0 / ix0;
  if (dim == 2)
    Norm = 1.0;
  if (dim == 1)
    Norm = ix0;
  if (xlength > 1)
    Norm *= (0.5);
  if (ylength > 1)
    Norm *= (0.5);
  if (zlength > 1)
    Norm *= (0.5);
  if (dim == 0) {
    Norm = 1.0;
    iX0 = 1.0;
    iX1 = 1.0;
    iX2 = 1.0;
  }
  if (ifnotAllskip) {
    bioR1__3d(HHH, HHL, HLH, HLL, LHH, LHL, LLH, LLL, outvector, xlength,
              ylength, zlength, iX2, Norm);
    pairity = 0;
    bio_3d_postmult(outvector, xlength, ylength, zlength, iX1, pairity);
  } else
    bioR1_3dskip(LLL, outvector, xlength, ylength, zlength, iX2, Norm);
  pairity = 1;
  bio_3d_postmult(outvector, xlength, ylength, zlength, iX0, pairity);
}


x86 assembly (clang, -O0):
<bioR5_3d>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %rbx
sub    $0x88,%rsp
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
mov    %rcx,-0x38(%rbp)
mov    %r8,-0x40(%rbp)
mov    %r9,-0x48(%rbp)
movss  0xe8c(%rip),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  0xe7b(%rip),%xmm0
movss  %xmm0,-0x50(%rbp)
cmpl   $0x1,0x28(%rbp)
setg   %al
and    $0x1,%al
movzbl %al,%eax
cmpl   $0x1,0x30(%rbp)
setg   %cl
and    $0x1,%cl
movzbl %cl,%ecx
add    %ecx,%eax
cmpl   $0x1,0x38(%rbp)
setg   %cl
and    $0x1,%cl
movzbl %cl,%ecx
add    %ecx,%eax
mov    %eax,-0x68(%rbp)
cmpl   $0x0,0x40(%rbp)
je     1211 <bioR5_3d+0xe1>
movss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
movss  -0x50(%rbp),%xmm0
mulss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe2f(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe11(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
jmp    1246 <bioR5_3d+0x116>
movss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
movss  0xde5(%rip),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdd7(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
cmpl   $0x3,-0x68(%rbp)
jne    126e <bioR5_3d+0x13e>
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdaf(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
cmpl   $0x2,-0x68(%rbp)
jne    1285 <bioR5_3d+0x155>
movss  0xd88(%rip),%xmm0
movss  %xmm0,-0x64(%rbp)
cmpl   $0x1,-0x68(%rbp)
jne    1299 <bioR5_3d+0x169>
movss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x64(%rbp)
cmpl   $0x1,0x28(%rbp)
jle    12c1 <bioR5_3d+0x191>
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd64(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
cmpl   $0x1,0x30(%rbp)
jle    12e9 <bioR5_3d+0x1b9>
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd3c(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
cmpl   $0x1,0x38(%rbp)
jle    1311 <bioR5_3d+0x1e1>
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd14(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
cmpl   $0x0,-0x68(%rbp)
jne    134f <bioR5_3d+0x21f>
movss  0xce5(%rip),%xmm0
movss  %xmm0,-0x64(%rbp)
movss  0xcd8(%rip),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  0xccb(%rip),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  0xcbe(%rip),%xmm0
movss  %xmm0,-0x5c(%rbp)
cmpl   $0x0,0x40(%rbp)
je     13ea <bioR5_3d+0x2ba>
mov    -0x20(%rbp),%rdi
mov    -0x28(%rbp),%rsi
mov    -0x30(%rbp),%rdx
mov    -0x38(%rbp),%rcx
mov    -0x40(%rbp),%r8
mov    -0x48(%rbp),%r9
mov    0x10(%rbp),%r15
mov    0x18(%rbp),%r14
mov    0x20(%rbp),%rbx
mov    0x28(%rbp),%r11d
mov    0x30(%rbp),%r10d
mov    0x38(%rbp),%eax
movss  -0x5c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x64(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    %r15,(%rsp)
mov    %r14,0x8(%rsp)
mov    %rbx,0x10(%rsp)
mov    %r11d,0x18(%rsp)
mov    %r10d,0x20(%rsp)
mov    %eax,0x28(%rsp)
mov    $0x2,%al
call   1030 <bioR1__3d@plt>
movl   $0x0,-0x60(%rbp)
mov    0x20(%rbp),%rdi
mov    0x28(%rbp),%esi
mov    0x30(%rbp),%edx
mov    0x38(%rbp),%ecx
movss  -0x58(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x60(%rbp),%r8d
mov    $0x1,%al
call   1040 <bio_3d_postmult@plt>
jmp    1415 <bioR5_3d+0x2e5>
mov    0x18(%rbp),%rdi
mov    0x20(%rbp),%rsi
mov    0x28(%rbp),%edx
mov    0x30(%rbp),%ecx
mov    0x38(%rbp),%r8d
movss  -0x5c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x64(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mov    $0x2,%al
call   1050 <bioR1_3dskip@plt>
movl   $0x1,-0x60(%rbp)
mov    0x20(%rbp),%rdi
mov    0x28(%rbp),%esi
mov    0x30(%rbp),%edx
mov    0x38(%rbp),%ecx
movss  -0x54(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x60(%rbp),%r8d
mov    $0x1,%al
call   1040 <bio_3d_postmult@plt>
add    $0x88,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

x86 assembly (clang, -O2):
<bioR5_3d>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    0x68(%rsp),%r14d
mov    0x60(%rsp),%r15d
mov    0x58(%rsp),%r12d
xor    %ebp,%ebp
cmp    $0x2,%r12d
setge  %bpl
xor    %eax,%eax
cmp    $0x2,%r15d
setge  %al
add    %ebp,%eax
xor    %r10d,%r10d
cmp    $0x2,%r14d
setge  %r10b
lea    (%rax,%r10,1),%ebp
xor    %ebx,%ebx
cmp    $0x2,%ebp
sete   %r11b
movss  0xea8(%rip),%xmm0
movaps %xmm0,%xmm1
cmp    $0x1,%ebp
jne    11e0 <bioR5_3d+0xb0>
cmp    $0x2,%r12d
jge    11f6 <bioR5_3d+0xc6>
cmp    $0x2,%r15d
jl     1198 <bioR5_3d+0x68>
mulss  %xmm0,%xmm1
mov    0x70(%rsp),%r11d
cmp    $0x2,%r14d
jl     11a7 <bioR5_3d+0x77>
mulss  %xmm0,%xmm1
xor    %ebx,%ebx
test   %r11d,%r11d
sete   %bpl
add    %r10d,%eax
je     1202 <bioR5_3d+0xd2>
mov    %bpl,%bl
lea    0xe49(%rip),%rax
movss  (%rax,%rbx,4),%xmm3
lea    0xe45(%rip),%rax
movss  (%rax,%rbx,4),%xmm0
movsd  0xe40(%rip),%xmm2
movsd  %xmm2,0x8(%rsp)
jmp    121e <bioR5_3d+0xee>
mov    %r11b,%bl
lea    0xe16(%rip),%rbp
movss  0x0(%rbp,%rbx,4),%xmm1
cmp    $0x2,%r12d
jl     118e <bioR5_3d+0x5e>
mulss  %xmm0,%xmm1
cmp    $0x2,%r15d
jge    1194 <bioR5_3d+0x64>
jmp    1198 <bioR5_3d+0x68>
movsd  0xe16(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
movss  0xe14(%rip),%xmm3
movaps %xmm3,%xmm0
movaps %xmm3,%xmm1
mov    0x50(%rsp),%rbx
mov    0x48(%rsp),%r10
cvtss2sd %xmm0,%xmm0
cvtss2sd %xmm1,%xmm1
test   %r11d,%r11d
je     1275 <bioR5_3d+0x145>
mov    $0x2,%al
push   %r14
push   %r15
push   %r12
push   %rbx
push   %r10
push   0x68(%rsp)
movss  %xmm3,0x34(%rsp)
call   1030 <bioR1__3d@plt>
add    $0x30,%rsp
movss  0x4(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %rbx,%rdi
mov    %r12d,%esi
mov    %r15d,%edx
mov    %r14d,%ecx
xor    %r8d,%r8d
mov    $0x1,%al
call   1040 <bio_3d_postmult@plt>
jmp    128b <bioR5_3d+0x15b>
mov    %r10,%rdi
mov    %rbx,%rsi
mov    %r12d,%edx
mov    %r15d,%ecx
mov    %r14d,%r8d
mov    $0x2,%al
call   1050 <bioR1_3dskip@plt>
mov    %rbx,%rdi
mov    %r12d,%esi
mov    %r15d,%edx
mov    %r14d,%ecx
movsd  0x8(%rsp),%xmm0
mov    $0x1,%r8d
mov    $0x1,%al
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <bio_3d_postmult@plt>

LLVM IR (-O0):
define void @bioR5_3d(float* %0, float* %1, float* %2, float* %3, float* %4, float* %5, float* %6, float* %7, float* %8, i32 %9, i32 %10, i32 %11, i32 %12) {
%14 = alloca float*
%15 = alloca float*
%16 = alloca float*
%17 = alloca float*
%18 = alloca float*
%19 = alloca float*
%20 = alloca float*
%21 = alloca float*
%22 = alloca float*
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca float
%28 = alloca float
%29 = alloca float
%30 = alloca float
%31 = alloca float
%32 = alloca i32
%33 = alloca float
%34 = alloca i32
store float* %0, float** %14
store float* %1, float** %15
store float* %2, float** %16
store float* %3, float** %17
store float* %4, float** %18
store float* %5, float** %19
store float* %6, float** %20
store float* %7, float** %21
store float* %8, float** %22
store i32 %9, i32* %23
store i32 %10, i32* %24
store i32 %11, i32* %25
store i32 %12, i32* %26
store float 5.000000e-01, float* %27
store float -2.500000e-01, float* %28
%35 = load i32, i32* %23
%36 = icmp sgt i32 %35, 1
%37 = zext i1 %36 to i32
%38 = load i32, i32* %24
%39 = icmp sgt i32 %38, 1
%40 = zext i1 %39 to i32
%41 = add nsw i32 %37, %40
%42 = load i32, i32* %25
%43 = icmp sgt i32 %42, 1
%44 = zext i1 %43 to i32
%45 = add nsw i32 %41, %44
store i32 %45, i32* %34
%46 = load i32, i32* %26
%47 = icmp ne i32 %46, 0
br i1 %47, label %48, label %60
48: ; preds = %13
%49 = load float, float* %28
store float %49, float* %31
%50 = load float, float* %28
%51 = load float, float* %27
%52 = fmul float %50, %51
%53 = fpext float %52 to double
%54 = fdiv double 1.000000e+00, %53
%55 = fptrunc double %54 to float
store float %55, float* %30
%56 = load float, float* %27
%57 = fpext float %56 to double
%58 = fdiv double 1.000000e+00, %57
%59 = fptrunc double %58 to float
store float %59, float* %29
br label %66
60: ; preds = %13
%61 = load float, float* %27
store float %61, float* %31
store float 1.000000e+00, float* %30
%62 = load float, float* %27
%63 = fpext float %62 to double
%64 = fdiv double 1.000000e+00, %63
%65 = fptrunc double %64 to float
store float %65, float* %29
br label %66
66: ; preds = %60, %48
%67 = load i32, i32* %34
%68 = icmp eq i32 %67, 3
br i1 %68, label %69, label %74
69: ; preds = %66
%70 = load float, float* %27
%71 = fpext float %70 to double
%72 = fdiv double 1.000000e+00, %71
%73 = fptrunc double %72 to float
store float %73, float* %33
br label %74
74: ; preds = %69, %66
%75 = load i32, i32* %34
%76 = icmp eq i32 %75, 2
br i1 %76, label %77, label %78
77: ; preds = %74
store float 1.000000e+00, float* %33
br label %78
78: ; preds = %77, %74
%79 = load i32, i32* %34
%80 = icmp eq i32 %79, 1
br i1 %80, label %81, label %83
81: ; preds = %78
%82 = load float, float* %27
store float %82, float* %33
br label %83
83: ; preds = %81, %78
%84 = load i32, i32* %23
%85 = icmp sgt i32 %84, 1
br i1 %85, label %86, label %91
86: ; preds = %83
%87 = load float, float* %33
%88 = fpext float %87 to double
%89 = fmul double %88, 5.000000e-01
%90 = fptrunc double %89 to float
store float %90, float* %33
br label %91
91: ; preds = %86, %83
%92 = load i32, i32* %24
%93 = icmp sgt i32 %92, 1
br i1 %93, label %94, label %99
94: ; preds = %91
%95 = load float, float* %33
%96 = fpext float %95 to double
%97 = fmul double %96, 5.000000e-01
%98 = fptrunc double %97 to float
store float %98, float* %33
br label %99
99: ; preds = %94, %91
%100 = load i32, i32* %25
%101 = icmp sgt i32 %100, 1
br i1 %101, label %102, label %107
102: ; preds = %99
%103 = load float, float* %33
%104 = fpext float %103 to double
%105 = fmul double %104, 5.000000e-01
%106 = fptrunc double %105 to float
store float %106, float* %33
br label %107
107: ; preds = %102, %99
%108 = load i32, i32* %34
%109 = icmp eq i32 %108, 0
br i1 %109, label %110, label %111
110: ; preds = %107
store float 1.000000e+00, float* %33
store float 1.000000e+00, float* %29
store float 1.000000e+00, float* %30
store float 1.000000e+00, float* %31
br label %111
111: ; preds = %110, %107
%112 = load i32, i32* %26
%113 = icmp ne i32 %112, 0
br i1 %113, label %114, label %140
114: ; preds = %111
%115 = load float*, float** %14
%116 = load float*, float** %15
%117 = load float*, float** %16
%118 = load float*, float** %17
%119 = load float*, float** %18
%120 = load float*, float** %19
%121 = load float*, float** %20
%122 = load float*, float** %21
%123 = load float*, float** %22
%124 = load i32, i32* %23
%125 = load i32, i32* %24
%126 = load i32, i32* %25
%127 = load float, float* %31
%128 = fpext float %127 to double
%129 = load float, float* %33
%130 = fpext float %129 to double
%131 = call i32 (float*, float*, float*, float*, float*, float*, float*, float*, float*, i32, i32, i32, double, double, ...) bitcast (i32 (...)* @bioR1__3d to i32 (float*, float*, float*, float*, float*, float*, float*, float*, float*, i32, i32, i32, double, double, ...)*)(float* %115, float* %116, float* %117, float* %118, float* %119, float* %120, float* %121, float* %122, float* %123, i32 %124, i32 %125, i32 %126, double %128, double %130)
store i32 0, i32* %32
%132 = load float*, float** %22
%133 = load i32, i32* %23
%134 = load i32, i32* %24
%135 = load i32, i32* %25
%136 = load float, float* %30
%137 = fpext float %136 to double
%138 = load i32, i32* %32
%139 = call i32 (float*, i32, i32, i32, double, i32, ...) bitcast (i32 (...)* @bio_3d_postmult to i32 (float*, i32, i32, i32, double, i32, ...)*)(float* %132, i32 %133, i32 %134, i32 %135, double %137, i32 %138)
br label %151
140: ; preds = %111
%141 = load float*, float** %21
%142 = load float*, float** %22
%143 = load i32, i32* %23
%144 = load i32, i32* %24
%145 = load i32, i32* %25
%146 = load float, float* %31
%147 = fpext float %146 to double
%148 = load float, float* %33
%149 = fpext float %148 to double
%150 = call i32 (float*, float*, i32, i32, i32, double, double, ...) bitcast (i32 (...)* @bioR1_3dskip to i32 (float*, float*, i32, i32, i32, double, double, ...)*)(float* %141, float* %142, i32 %143, i32 %144, i32 %145, double %147, double %149)
br label %151
151: ; preds = %140, %114
store i32 1, i32* %32
%152 = load float*, float** %22
%153 = load i32, i32* %23
%154 = load i32, i32* %24
%155 = load i32, i32* %25
%156 = load float, float* %29
%157 = fpext float %156 to double
%158 = load i32, i32* %32
%159 = call i32 (float*, i32, i32, i32, double, i32, ...) bitcast (i32 (...)* @bio_3d_postmult to i32 (float*, i32, i32, i32, double, i32, ...)*)(float* %152, i32 %153, i32 %154, i32 %155, double %157, i32 %158)
ret void
}

LLVM IR (-O2):
define void @bioR5_3d(float* %0, float* %1, float* %2, float* %3, float* %4, float* %5, float* %6, float* %7, float* %8, i32 %9, i32 %10, i32 %11, i32 %12) {
%14 = icmp sgt i32 %9, 1
%15 = zext i1 %14 to i32
%16 = icmp sgt i32 %10, 1
%17 = zext i1 %16 to i32
%18 = add nuw nsw i32 %17, %15
%19 = icmp sgt i32 %11, 1
%20 = zext i1 %19 to i32
%21 = add nuw nsw i32 %18, %20
%22 = icmp eq i32 %12, 0
%23 = select i1 %22, float 1.000000e+00, float -8.000000e+00
%24 = select i1 %22, float 5.000000e-01, float -2.500000e-01
%25 = icmp eq i32 %21, 2
%26 = select i1 %25, float 1.000000e+00, float 2.000000e+00
%27 = icmp eq i32 %21, 1
%28 = select i1 %27, float 5.000000e-01, float %26
%29 = fmul float %28, 5.000000e-01
%30 = select i1 %14, float %29, float %28
%31 = fmul float %30, 5.000000e-01
%32 = select i1 %16, float %31, float %30
%33 = fmul float %32, 5.000000e-01
%34 = select i1 %19, float %33, float %32
%35 = icmp eq i32 %21, 0
br i1 %35, label %36, label %37
36: ; preds = %13
br label %37
37: ; preds = %36, %13
%38 = phi float [ 1.000000e+00, %36 ], [ %23, %13 ]
%39 = phi float [ 1.000000e+00, %36 ], [ %24, %13 ]
%40 = phi float [ 1.000000e+00, %36 ], [ %34, %13 ]
%41 = phi double [ 1.000000e+00, %36 ], [ 2.000000e+00, %13 ]
%42 = fpext float %39 to double
%43 = fpext float %40 to double
br i1 %22, label %48, label %44
44: ; preds = %37
%45 = tail call i32 (float*, float*, float*, float*, float*, float*, float*, float*, float*, i32, i32, i32, double, double, ...) bitcast (i32 (...)* @bioR1__3d to i32 (float*, float*, float*, float*, float*, float*, float*, float*, float*, i32, i32, i32, double, double, ...)*)(float* %0, float* %1, float* %2, float* %3, float* %4, float* %5, float* %6, float* %7, float* %8, i32 %9, i32 %10, i32 %11, double %42, double %43)
%46 = fpext float %38 to double
%47 = tail call i32 (float*, i32, i32, i32, double, i32, ...) bitcast (i32 (...)* @bio_3d_postmult to i32 (float*, i32, i32, i32, double, i32, ...)*)(float* %8, i32 %9, i32 %10, i32 %11, double %46, i32 0)
br label %50
48: ; preds = %37
%49 = tail call i32 (float*, float*, i32, i32, i32, double, double, ...) bitcast (i32 (...)* @bioR1_3dskip to i32 (float*, float*, i32, i32, i32, double, double, ...)*)(float* %7, float* %8, i32 %9, i32 %10, i32 %11, double %42, double %43)
br label %50
50: ; preds = %48, %44
%51 = tail call i32 (float*, i32, i32, i32, double, i32, ...) bitcast (i32 (...)* @bio_3d_postmult to i32 (float*, i32, i32, i32, double, i32, ...)*)(float* %8, i32 %9, i32 %10, i32 %11, double %41, i32 1)
ret void
}