unsigned long long countErrorBank(int nShell, int *shellMap,
                                  double *Schwarz_Shell, double cutoff) {
  unsigned long long nErrorBank = 0;
  int pS, qS, iS, jS;
  int n;
  for (pS = nShell - 1; pS >= 0; pS--)
    for (qS = pS; qS >= 0; qS--)
      for (iS = pS; iS >= 0; iS--) {
        if (iS == pS)
          n = qS;
        else
          n = iS;
        for (jS = n; jS >= 0; jS--) {
          if (Schwarz_Shell[pS * nShell + qS] *
                  Schwarz_Shell[iS * nShell + jS] <
              cutoff)
            continue;
          nErrorBank++;
        }
      }
  return nErrorBank;
}


x86 assembly (clang, -O0):
<countErrorBank>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jl     120f <countErrorBank+0x10f>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x30(%rbp)
jl     11fc <countErrorBank+0xfc>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     11e9 <countErrorBank+0xe9>
mov    -0x34(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jne    1166 <countErrorBank+0x66>
mov    -0x30(%rbp),%eax
mov    %eax,-0x3c(%rbp)
jmp    116c <countErrorBank+0x6c>
mov    -0x34(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x38(%rbp)
jl     11d6 <countErrorBank+0xd6>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11bc <countErrorBank+0xbc>
jmp    11c8 <countErrorBank+0xc8>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
jmp    1172 <countErrorBank+0x72>
jmp    11db <countErrorBank+0xdb>
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1145 <countErrorBank+0x45>
jmp    11ee <countErrorBank+0xee>
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    1135 <countErrorBank+0x35>
jmp    1201 <countErrorBank+0x101>
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    1125 <countErrorBank+0x25>
mov    -0x28(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<countErrorBank>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
jle    12d6 <countErrorBank+0x1d6>
mov    %edi,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x8(%rdx),%rax
mov    %rax,-0x18(%rsp)
lea    -0x1(%rcx),%rax
imul   %rcx,%rax
mov    %rax,-0x30(%rsp)
lea    (%rdx,%rax,8),%rax
mov    %rax,-0x38(%rsp)
lea    0x0(,%rcx,8),%r13
neg    %r13
xor    %eax,%eax
mov    %rcx,%r10
jmp    1166 <countErrorBank+0x66>
nopl   (%rax)
sub    %rcx,-0x30(%rsp)
add    %r13,-0x38(%rsp)
cmpq   $0x2,-0x28(%rsp)
jl     12d8 <countErrorBank+0x1d8>
mov    %r10,-0x28(%rsp)
add    $0xffffffffffffffff,%r10
mov    %r10,%rsi
imul   %rcx,%rsi
mov    %rsi,-0x20(%rsp)
mov    %r10,%r15
mov    %r10,-0x10(%rsp)
jmp    119c <countErrorBank+0x9c>
cs nopw 0x0(%rax,%rax,1)
nop
lea    -0x1(%r15),%rsi
test   %r15,%r15
mov    %rsi,%r15
jle    1150 <countErrorBank+0x50>
mov    -0x20(%rsp),%rsi
lea    (%r15,%rsi,1),%r14
mov    -0x38(%rsp),%rsi
mov    -0x30(%rsp),%rdi
mov    %r10,%r12
mov    %r14,-0x8(%rsp)
jmp    11d3 <countErrorBank+0xd3>
nopl   0x0(%rax)
lea    -0x1(%r12),%rbx
sub    %rcx,%rdi
add    %r13,%rsi
test   %r12,%r12
mov    %rbx,%r12
jle    1190 <countErrorBank+0x90>
cmp    %r10,%r12
mov    %r12,%rbp
cmove  %r15,%rbp
test   $0x80000000,%ebp
jne    11c0 <countErrorBank+0xc0>
movslq %ebp,%r11
test   %r11,%r11
mov    $0x0,%ebp
cmovs  %r11,%rbp
movsd  (%rdx,%r14,8),%xmm2
neg    %rbp
lea    (%r11,%rbp,1),%r9
add    $0x1,%r9
cmp    $0x4,%r9
jb     129f <countErrorBank+0x19f>
mov    %r9,%r8
and    $0xfffffffffffffffc,%r8
movq   %rax,%xmm3
movapd %xmm2,%xmm4
unpcklpd %xmm2,%xmm4
mov    %r8,%rax
neg    %rax
lea    (%rdi,%r11,1),%rbp
sub    %r8,%r11
mov    -0x18(%rsp),%rbx
lea    (%rbx,%rbp,8),%r14
pxor   %xmm5,%xmm5
xor    %r10d,%r10d
movupd -0x10(%r14,%r10,8),%xmm6
movupd (%r14,%r10,8),%xmm7
shufpd $0x1,%xmm7,%xmm7
shufpd $0x1,%xmm6,%xmm6
mulpd  %xmm4,%xmm7
mulpd  %xmm4,%xmm6
cmpnltpd %xmm1,%xmm7
psubq  %xmm7,%xmm3
cmpnltpd %xmm1,%xmm6
psubq  %xmm6,%xmm5
add    $0xfffffffffffffffc,%r10
cmp    %r10,%rax
jne    1240 <countErrorBank+0x140>
paddq  %xmm3,%xmm5
pshufd $0xee,%xmm5,%xmm3
paddq  %xmm5,%xmm3
movq   %xmm3,%rax
cmp    %r8,%r9
mov    -0x10(%rsp),%r10
mov    -0x8(%rsp),%r14
je     11c0 <countErrorBank+0xc0>
add    $0x1,%r11
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x8(%rsi,%r11,8),%xmm3
mulsd  %xmm2,%xmm3
xor    %ebp,%ebp
ucomisd %xmm3,%xmm0
setbe  %bpl
add    %rbp,%rax
add    $0xffffffffffffffff,%r11
test   %r11,%r11
jg     12b0 <countErrorBank+0x1b0>
jmp    11c0 <countErrorBank+0xc0>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i64 @countErrorBank(i32 %0, i32* %1, double* %2, double %3) {
%5 = alloca i32
%6 = alloca i32*
%7 = alloca double*
%8 = alloca double
%9 = alloca i64
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
store i32 %0, i32* %5
store i32* %1, i32** %6
store double* %2, double** %7
store double %3, double* %8
store i64 0, i64* %9
%15 = load i32, i32* %5
%16 = sub nsw i32 %15, 1
store i32 %16, i32* %10
br label %17
17: ; preds = %81, %4
%18 = load i32, i32* %10
%19 = icmp sge i32 %18, 0
br i1 %19, label %20, label %84
20: ; preds = %17
%21 = load i32, i32* %10
store i32 %21, i32* %11
br label %22
22: ; preds = %77, %20
%23 = load i32, i32* %11
%24 = icmp sge i32 %23, 0
br i1 %24, label %25, label %80
25: ; preds = %22
%26 = load i32, i32* %10
store i32 %26, i32* %12
br label %27
27: ; preds = %73, %25
%28 = load i32, i32* %12
%29 = icmp sge i32 %28, 0
br i1 %29, label %30, label %76
30: ; preds = %27
%31 = load i32, i32* %12
%32 = load i32, i32* %10
%33 = icmp eq i32 %31, %32
br i1 %33, label %34, label %36
34: ; preds = %30
%35 = load i32, i32* %11
store i32 %35, i32* %14
br label %38
36: ; preds = %30
%37 = load i32, i32* %12
store i32 %37, i32* %14
br label %38
38: ; preds = %36, %34
%39 = load i32, i32* %14
store i32 %39, i32* %13
br label %40
40: ; preds = %69, %38
%41 = load i32, i32* %13
%42 = icmp sge i32 %41, 0
br i1 %42, label %43, label %72
43: ; preds = %40
%44 = load double*, double** %7
%45 = load i32, i32* %10
%46 = load i32, i32* %5
%47 = mul nsw i32 %45, %46
%48 = load i32, i32* %11
%49 = add nsw i32 %47, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %44, i64 %50
%52 = load double, double* %51
%53 = load double*, double** %7
%54 = load i32, i32* %12
%55 = load i32, i32* %5
%56 = mul nsw i32 %54, %55
%57 = load i32, i32* %13
%58 = add nsw i32 %56, %57
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %53, i64 %59
%61 = load double, double* %60
%62 = fmul double %52, %61
%63 = load double, double* %8
%64 = fcmp olt double %62, %63
br i1 %64, label %65, label %66
65: ; preds = %43
br label %69
66: ; preds = %43
%67 = load i64, i64* %9
%68 = add i64 %67, 1
store i64 %68, i64* %9
br label %69
69: ; preds = %66, %65
%70 = load i32, i32* %13
%71 = add nsw i32 %70, -1
store i32 %71, i32* %13
br label %40
72: ; preds = %40
br label %73
73: ; preds = %72
%74 = load i32, i32* %12
%75 = add nsw i32 %74, -1
store i32 %75, i32* %12
br label %27
76: ; preds = %27
br label %77
77: ; preds = %76
%78 = load i32, i32* %11
%79 = add nsw i32 %78, -1
store i32 %79, i32* %11
br label %22
80: ; preds = %22
br label %81
81: ; preds = %80
%82 = load i32, i32* %10
%83 = add nsw i32 %82, -1
store i32 %83, i32* %10
br label %17
84: ; preds = %17
%85 = load i64, i64* %9
ret i64 %85
}

LLVM IR (-O2):
define i64 @countErrorBank(i32 %0, i32* nocapture readnone %1, double* nocapture %2, double %3) {
%5 = icmp sgt i32 %0, 0
br i1 %5, label %6, label %102
6: ; preds = %4
%7 = zext i32 %0 to i64
%8 = insertelement <2 x double> poison, double %3, i64 0
%9 = shufflevector <2 x double> %8, <2 x double> poison, <2 x i32> zeroinitializer
%10 = insertelement <2 x double> poison, double %3, i64 0
%11 = shufflevector <2 x double> %10, <2 x double> poison, <2 x i32> zeroinitializer
br label %14
12: ; preds = %99
%13 = icmp sgt i64 %15, 1
br i1 %13, label %14, label %102
14: ; preds = %12, %6
%15 = phi i64 [ %7, %6 ], [ %19, %12 ]
%16 = phi i64 [ %7, %6 ], [ %18, %12 ]
%17 = phi i64 [ 0, %6 ], [ %96, %12 ]
%18 = add nsw i64 %16, -1
%19 = add nsw i64 %15, -1
%20 = mul nsw i64 %19, %7
br label %21
21: ; preds = %99, %14
%22 = phi i64 [ %18, %14 ], [ %100, %99 ]
%23 = phi i64 [ %17, %14 ], [ %96, %99 ]
%24 = add nsw i64 %22, %20
%25 = getelementptr inbounds double, double* %2, i64 %24
br label %26
26: ; preds = %21, %95
%27 = phi i64 [ %18, %21 ], [ %97, %95 ]
%28 = phi i64 [ %23, %21 ], [ %96, %95 ]
%29 = icmp eq i64 %27, %19
%30 = select i1 %29, i64 %22, i64 %27
%31 = and i64 %30, 2147483648
%32 = icmp eq i64 %31, 0
br i1 %32, label %33, label %95
33: ; preds = %26
%34 = load double, double* %25
%35 = mul nsw i64 %27, %7
%36 = shl i64 %30, 32
%37 = ashr exact i64 %36, 32
%38 = add nsw i64 %37, 1
%39 = call i64 @llvm.smin.i64(i64 %37, i64 0)
%40 = sub i64 %38, %39
%41 = icmp ult i64 %40, 4
br i1 %41, label %80, label %42
42: ; preds = %33
%43 = and i64 %40, -4
%44 = sub i64 %37, %43
%45 = insertelement <2 x i64> <i64 poison, i64 0>, i64 %28, i64 0
%46 = insertelement <2 x double> poison, double %34, i64 0
%47 = shufflevector <2 x double> %46, <2 x double> poison, <2 x i32> zeroinitializer
%48 = insertelement <2 x double> poison, double %34, i64 0
%49 = shufflevector <2 x double> %48, <2 x double> poison, <2 x i32> zeroinitializer
br label %50
50: ; preds = %50, %42
%51 = phi i64 [ 0, %42 ], [ %74, %50 ]
%52 = phi <2 x i64> [ %45, %42 ], [ %72, %50 ]
%53 = phi <2 x i64> [ zeroinitializer, %42 ], [ %73, %50 ]
%54 = sub i64 %37, %51
%55 = add nsw i64 %54, %35
%56 = getelementptr inbounds double, double* %2, i64 %55
%57 = getelementptr inbounds double, double* %56, i64 -1
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = shufflevector <2 x double> %59, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%61 = getelementptr inbounds double, double* %56, i64 -2
%62 = getelementptr inbounds double, double* %61, i64 -1
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = shufflevector <2 x double> %64, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%66 = fmul <2 x double> %47, %60
%67 = fmul <2 x double> %49, %65
%68 = fcmp uge <2 x double> %66, %9
%69 = fcmp uge <2 x double> %67, %11
%70 = zext <2 x i1> %68 to <2 x i64>
%71 = zext <2 x i1> %69 to <2 x i64>
%72 = add <2 x i64> %52, %70
%73 = add <2 x i64> %53, %71
%74 = add nuw i64 %51, 4
%75 = icmp eq i64 %74, %43
br i1 %75, label %76, label %50
76: ; preds = %50
%77 = add <2 x i64> %73, %72
%78 = call i64 @llvm.vector.reduce.add.v2i64(<2 x i64> %77)
%79 = icmp eq i64 %40, %43
br i1 %79, label %95, label %80
80: ; preds = %33, %76
%81 = phi i64 [ %37, %33 ], [ %44, %76 ]
%82 = phi i64 [ %28, %33 ], [ %78, %76 ]
br label %83
83: ; preds = %80, %83
%84 = phi i64 [ %93, %83 ], [ %81, %80 ]
%85 = phi i64 [ %92, %83 ], [ %82, %80 ]
%86 = add nsw i64 %84, %35
%87 = getelementptr inbounds double, double* %2, i64 %86
%88 = load double, double* %87
%89 = fmul double %34, %88
%90 = fcmp uge double %89, %3
%91 = zext i1 %90 to i64
%92 = add i64 %85, %91
%93 = add nsw i64 %84, -1
%94 = icmp sgt i64 %84, 0
br i1 %94, label %83, label %95
95: ; preds = %83, %76, %26
%96 = phi i64 [ %28, %26 ], [ %78, %76 ], [ %92, %83 ]
%97 = add nsw i64 %27, -1
%98 = icmp sgt i64 %27, 0
br i1 %98, label %26, label %99
99: ; preds = %95
%100 = add nsw i64 %22, -1
%101 = icmp sgt i64 %22, 0
br i1 %101, label %21, label %12
102: ; preds = %12, %4
%103 = phi i64 [ 0, %4 ], [ %96, %12 ]
ret i64 %103
}