int SLS_asolve(const unsigned long n, const double b[], double x[],
               const int itrnsp, const double sa[], const unsigned long ija[]) {
#pragma omp parallel for
  for (long i = 1; i <= n; i++)
    x[i] = (sa[i] != 0.0 ? b[i] / sa[i] : b[i]);
  return 0;
}


x86 assembly (clang, -O0):
<SLS_asolve>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movq   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
ja     11af <SLS_asolve+0xaf>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1156 <SLS_asolve+0x56>
jp     1156 <SLS_asolve+0x56>
jmp    117a <SLS_asolve+0x7a>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    118c <SLS_asolve+0x8c>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1123 <SLS_asolve+0x23>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<SLS_asolve>:
push   %rbp
push   %rbx
test   %rdi,%rdi
je     1307 <SLS_asolve+0x207>
mov    %rdi,%r9
cmp    $0x1,%rdi
jne    11e3 <SLS_asolve+0xe3>
lea    0x1(%rdi),%rax
mov    %eax,%ebp
sub    %r9d,%ebp
mov    %r9,%rcx
test   $0x1,%bpl
je     1161 <SLS_asolve+0x61>
movsd  (%r8,%r9,8),%xmm0
movsd  (%rsi,%r9,8),%xmm1
xorpd  %xmm2,%xmm2
cmpneqsd %xmm0,%xmm2
andpd  %xmm2,%xmm0
movsd  0xec5(%rip),%xmm3
andnpd %xmm3,%xmm2
orpd   %xmm0,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rdx,%r9,8)
lea    0x1(%r9),%rcx
cmp    %rdi,%r9
je     1307 <SLS_asolve+0x207>
xorpd  %xmm0,%xmm0
movsd  0xe9a(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
movsd  (%r8,%rcx,8),%xmm2
movsd  (%rsi,%rcx,8),%xmm3
movapd %xmm2,%xmm4
cmpneqsd %xmm0,%xmm4
andpd  %xmm4,%xmm2
andnpd %xmm1,%xmm4
orpd   %xmm2,%xmm4
divsd  %xmm4,%xmm3
movsd  %xmm3,(%rdx,%rcx,8)
movsd  0x8(%r8,%rcx,8),%xmm2
movsd  0x8(%rsi,%rcx,8),%xmm3
movapd %xmm2,%xmm4
cmpneqsd %xmm0,%xmm4
andpd  %xmm4,%xmm2
andnpd %xmm1,%xmm4
orpd   %xmm2,%xmm4
divsd  %xmm4,%xmm3
movsd  %xmm3,0x8(%rdx,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rax
jne    1180 <SLS_asolve+0x80>
jmp    1307 <SLS_asolve+0x207>
lea    0x8(%rdx),%rax
lea    (%rdx,%rdi,8),%rcx
add    $0x8,%rcx
lea    0x8(%r8),%r9
lea    (%r8,%rdi,8),%r10
add    $0x8,%r10
lea    0x8(%rsi),%r11
lea    (%rsi,%rdi,8),%rbx
add    $0x8,%rbx
cmp    %r10,%rax
setb   %r10b
cmp    %rcx,%r9
setb   %bpl
cmp    %rbx,%rax
setb   %al
cmp    %rcx,%r11
setb   %cl
mov    $0x1,%r9d
test   %bpl,%r10b
jne    1118 <SLS_asolve+0x18>
and    %cl,%al
jne    1118 <SLS_asolve+0x18>
mov    %rdi,%r9
and    $0xfffffffffffffffe,%r9
lea    -0x2(%r9),%rax
mov    %rax,%r10
shr    %r10
add    $0x1,%r10
test   %rax,%rax
je     1318 <SLS_asolve+0x218>
mov    %r10,%rax
and    $0xfffffffffffffffe,%rax
mov    $0x1,%ecx
xorpd  %xmm0,%xmm0
movapd 0xd92(%rip),%xmm1
xchg   %ax,%ax
movupd (%r8,%rcx,8),%xmm2
movapd %xmm2,%xmm3
cmpneqpd %xmm0,%xmm3
movupd (%rsi,%rcx,8),%xmm4
andpd  %xmm3,%xmm2
andnpd %xmm1,%xmm3
orpd   %xmm2,%xmm3
divpd  %xmm3,%xmm4
movupd %xmm4,(%rdx,%rcx,8)
movupd 0x10(%r8,%rcx,8),%xmm2
movapd %xmm2,%xmm3
cmpneqpd %xmm0,%xmm3
movupd 0x10(%rsi,%rcx,8),%xmm4
andpd  %xmm3,%xmm2
andnpd %xmm1,%xmm3
orpd   %xmm2,%xmm3
divpd  %xmm3,%xmm4
movupd %xmm4,0x10(%rdx,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rax
jne    1270 <SLS_asolve+0x170>
test   $0x1,%r10b
je     1302 <SLS_asolve+0x202>
movupd (%r8,%rcx,8),%xmm0
xorpd  %xmm1,%xmm1
cmpneqpd %xmm0,%xmm1
andpd  %xmm1,%xmm0
andnpd 0xd10(%rip),%xmm1
movupd (%rsi,%rcx,8),%xmm2
orpd   %xmm0,%xmm1
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdx,%rcx,8)
cmp    %rdi,%r9
jne    130c <SLS_asolve+0x20c>
xor    %eax,%eax
pop    %rbx
pop    %rbp
ret
mov    %rdi,%r9
or     $0x1,%r9
jmp    1118 <SLS_asolve+0x18>
mov    $0x1,%ecx
test   $0x1,%r10b
jne    12d5 <SLS_asolve+0x1d5>
jmp    1302 <SLS_asolve+0x202>

LLVM IR (-O0):
define i32 @SLS_asolve(i64 %0, double* %1, double* %2, i32 %3, double* %4, i64* %5) {
%7 = alloca i64
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca double*
%12 = alloca i64*
%13 = alloca i64
store i64 %0, i64* %7
store double* %1, double** %8
store double* %2, double** %9
store i32 %3, i32* %10
store double* %4, double** %11
store i64* %5, i64** %12
store i64 1, i64* %13
br label %14
14: ; preds = %44, %6
%15 = load i64, i64* %13
%16 = load i64, i64* %7
%17 = icmp ule i64 %15, %16
br i1 %17, label %18, label %47
18: ; preds = %14
%19 = load double*, double** %11
%20 = load i64, i64* %13
%21 = getelementptr inbounds double, double* %19, i64 %20
%22 = load double, double* %21
%23 = fcmp une double %22, 0.000000e+00
br i1 %23, label %24, label %34
24: ; preds = %18
%25 = load double*, double** %8
%26 = load i64, i64* %13
%27 = getelementptr inbounds double, double* %25, i64 %26
%28 = load double, double* %27
%29 = load double*, double** %11
%30 = load i64, i64* %13
%31 = getelementptr inbounds double, double* %29, i64 %30
%32 = load double, double* %31
%33 = fdiv double %28, %32
br label %39
34: ; preds = %18
%35 = load double*, double** %8
%36 = load i64, i64* %13
%37 = getelementptr inbounds double, double* %35, i64 %36
%38 = load double, double* %37
br label %39
39: ; preds = %34, %24
%40 = phi double [ %33, %24 ], [ %38, %34 ]
%41 = load double*, double** %9
%42 = load i64, i64* %13
%43 = getelementptr inbounds double, double* %41, i64 %42
store double %40, double* %43
br label %44
44: ; preds = %39
%45 = load i64, i64* %13
%46 = add nsw i64 %45, 1
store i64 %46, i64* %13
br label %14
47: ; preds = %14
ret i32 0
}

LLVM IR (-O2):
define i32 @SLS_asolve(i64 %0, double* nocapture %1, double* nocapture writeonly %2, i32 %3, double* nocapture %4, i64* nocapture readnone %5) {
switch i64 %0, label %7 [
i64 0, label %100
i64 1, label %81
]
7: ; preds = %6
%8 = getelementptr double, double* %2, i64 1
%9 = add i64 %0, 1
%10 = getelementptr double, double* %2, i64 %9
%11 = getelementptr double, double* %4, i64 1
%12 = getelementptr double, double* %4, i64 %9
%13 = getelementptr double, double* %1, i64 1
%14 = getelementptr double, double* %1, i64 %9
%15 = icmp ult double* %8, %12
%16 = icmp ult double* %11, %10
%17 = and i1 %15, %16
%18 = icmp ult double* %8, %14
%19 = icmp ult double* %13, %10
%20 = and i1 %18, %19
%21 = or i1 %17, %20
br i1 %21, label %81, label %22
22: ; preds = %7
%23 = and i64 %0, -2
%24 = or i64 %0, 1
%25 = add i64 %23, -2
%26 = lshr exact i64 %25, 1
%27 = add nuw i64 %26, 1
%28 = and i64 %27, 1
%29 = icmp eq i64 %25, 0
br i1 %29, label %64, label %30
30: ; preds = %22
%31 = and i64 %27, -2
br label %32
32: ; preds = %32, %30
%33 = phi i64 [ 0, %30 ], [ %59, %32 ]
%34 = phi i64 [ 0, %30 ], [ %60, %32 ]
%35 = or i64 %33, 1
%36 = getelementptr inbounds double, double* %4, i64 %35
%37 = bitcast double* %36 to <2 x double>*
%38 = load <2 x double>, <2 x double>* %37
%39 = fcmp une <2 x double> %38, zeroinitializer
%40 = getelementptr inbounds double, double* %1, i64 %35
%41 = bitcast double* %40 to <2 x double>*
%42 = load <2 x double>, <2 x double>* %41
%43 = select <2 x i1> %39, <2 x double> %38, <2 x double> <double 1.000000e+00, double 1.000000e+00>
%44 = fdiv <2 x double> %42, %43
%45 = getelementptr inbounds double, double* %2, i64 %35
%46 = bitcast double* %45 to <2 x double>*
store <2 x double> %44, <2 x double>* %46
%47 = or i64 %33, 3
%48 = getelementptr inbounds double, double* %4, i64 %47
%49 = bitcast double* %48 to <2 x double>*
%50 = load <2 x double>, <2 x double>* %49
%51 = fcmp une <2 x double> %50, zeroinitializer
%52 = getelementptr inbounds double, double* %1, i64 %47
%53 = bitcast double* %52 to <2 x double>*
%54 = load <2 x double>, <2 x double>* %53
%55 = select <2 x i1> %51, <2 x double> %50, <2 x double> <double 1.000000e+00, double 1.000000e+00>
%56 = fdiv <2 x double> %54, %55
%57 = getelementptr inbounds double, double* %2, i64 %47
%58 = bitcast double* %57 to <2 x double>*
store <2 x double> %56, <2 x double>* %58
%59 = add nuw i64 %33, 4
%60 = add nuw i64 %34, 2
%61 = icmp eq i64 %60, %31
br i1 %61, label %62, label %32
62: ; preds = %32
%63 = or i64 %59, 1
br label %64
64: ; preds = %62, %22
%65 = phi i64 [ 1, %22 ], [ %63, %62 ]
%66 = icmp eq i64 %28, 0
br i1 %66, label %79, label %67
67: ; preds = %64
%68 = getelementptr inbounds double, double* %4, i64 %65
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = fcmp une <2 x double> %70, zeroinitializer
%72 = getelementptr inbounds double, double* %1, i64 %65
%73 = bitcast double* %72 to <2 x double>*
%74 = load <2 x double>, <2 x double>* %73
%75 = select <2 x i1> %71, <2 x double> %70, <2 x double> <double 1.000000e+00, double 1.000000e+00>
%76 = fdiv <2 x double> %74, %75
%77 = getelementptr inbounds double, double* %2, i64 %65
%78 = bitcast double* %77 to <2 x double>*
store <2 x double> %76, <2 x double>* %78
br label %79
79: ; preds = %64, %67
%80 = icmp eq i64 %23, %0
br i1 %80, label %100, label %81
81: ; preds = %6, %7, %79
%82 = phi i64 [ %0, %6 ], [ 1, %7 ], [ %24, %79 ]
%83 = add i64 %0, 1
%84 = sub i64 %83, %82
%85 = and i64 %84, 1
%86 = icmp eq i64 %85, 0
br i1 %86, label %97, label %87
87: ; preds = %81
%88 = getelementptr inbounds double, double* %4, i64 %82
%89 = load double, double* %88
%90 = fcmp une double %89, 0.000000e+00
%91 = getelementptr inbounds double, double* %1, i64 %82
%92 = load double, double* %91
%93 = select i1 %90, double %89, double 1.000000e+00
%94 = fdiv double %92, %93
%95 = getelementptr inbounds double, double* %2, i64 %82
store double %94, double* %95
%96 = add nuw i64 %82, 1
br label %97
97: ; preds = %87, %81
%98 = phi i64 [ %82, %81 ], [ %96, %87 ]
%99 = icmp eq i64 %82, %0
br i1 %99, label %100, label %101
100: ; preds = %97, %101, %79, %6
ret i32 0
101: ; preds = %97, %101
%102 = phi i64 [ %120, %101 ], [ %98, %97 ]
%103 = getelementptr inbounds double, double* %4, i64 %102
%104 = load double, double* %103
%105 = fcmp une double %104, 0.000000e+00
%106 = getelementptr inbounds double, double* %1, i64 %102
%107 = load double, double* %106
%108 = select i1 %105, double %104, double 1.000000e+00
%109 = fdiv double %107, %108
%110 = getelementptr inbounds double, double* %2, i64 %102
store double %109, double* %110
%111 = add nuw i64 %102, 1
%112 = getelementptr inbounds double, double* %4, i64 %111
%113 = load double, double* %112
%114 = fcmp une double %113, 0.000000e+00
%115 = getelementptr inbounds double, double* %1, i64 %111
%116 = load double, double* %115
%117 = select i1 %114, double %113, double 1.000000e+00
%118 = fdiv double %116, %117
%119 = getelementptr inbounds double, double* %2, i64 %111
store double %118, double* %119
%120 = add nuw i64 %102, 2
%121 = icmp eq i64 %111, %0
br i1 %121, label %100, label %101
}