void sprstx(float sa[], unsigned long ija[], float x[], float b[],
            unsigned long n) {
  void nrerror(char error_text[]);
  unsigned long i, j, k;
  if (ija[1] != n + 2)
    nrerror("mismatched vector and matrix in sprstx");
  for (i = 1; i <= n; i++)
    b[i] = sa[i] * x[i];
  for (i = 1; i <= n; i++) {
    for (k = ija[i]; k <= ija[i + 1] - 1; k++) {
      j = ija[k];
      b[j] += sa[k] * x[i];
    }
  }
}


x86 assembly (clang, -O0):
<sprstx>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x28(%rbp),%rcx
add    $0x2,%rcx
cmp    %rcx,%rax
je     1151 <sprstx+0x41>
lea    0xeb4(%rip),%rdi
call   1030 <nrerror@plt>
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     119f <sprstx+0x8f>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1159 <sprstx+0x49>
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     124e <sprstx+0x13e>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x30(%rbp),%rdx
mov    0x8(%rcx,%rdx,8),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
ja     1238 <sprstx+0x128>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11c5 <sprstx+0xb5>
jmp    123d <sprstx+0x12d>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11a7 <sprstx+0x97>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sprstx>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8,%r14
mov    %rcx,%rbx
mov    %rdx,%r15
mov    %rsi,%r12
mov    %rdi,%r13
lea    0x2(%r8),%rax
cmp    %rax,0x8(%rsi)
je     113e <sprstx+0x2e>
lea    0xec7(%rip),%rdi
call   1030 <nrerror@plt>
test   %r14,%r14
je     1376 <sprstx+0x266>
mov    $0x1,%eax
cmp    $0x8,%r14
jae    1213 <sprstx+0x103>
mov    %r14d,%edx
sub    %eax,%edx
add    $0x1,%edx
mov    %r14,%rcx
sub    %rax,%rcx
and    $0x3,%rdx
je     118c <sprstx+0x7c>
nopw   0x0(%rax,%rax,1)
movss  0x0(%r13,%rax,4),%xmm0
mulss  (%r15,%rax,4),%xmm0
movss  %xmm0,(%rbx,%rax,4)
add    $0x1,%rax
add    $0xffffffffffffffff,%rdx
jne    1170 <sprstx+0x60>
cmp    $0x3,%rcx
jb     1318 <sprstx+0x208>
mov    %r14,%rcx
sub    %rax,%rcx
lea    (%rbx,%rax,4),%rdx
add    $0xc,%rdx
lea    (%r15,%rax,4),%rsi
add    $0xc,%rsi
lea    0xc(,%rax,4),%rax
add    %r13,%rax
mov    $0xffffffffffffffff,%rdi
xchg   %ax,%ax
movss  -0x8(%rax,%rdi,4),%xmm0
mulss  -0x8(%rsi,%rdi,4),%xmm0
movss  %xmm0,-0x8(%rdx,%rdi,4)
movss  -0x4(%rax,%rdi,4),%xmm0
mulss  -0x4(%rsi,%rdi,4),%xmm0
movss  %xmm0,-0x4(%rdx,%rdi,4)
movss  (%rax,%rdi,4),%xmm0
mulss  (%rsi,%rdi,4),%xmm0
movss  %xmm0,(%rdx,%rdi,4)
movss  0x4(%rax,%rdi,4),%xmm0
mulss  0x4(%rsi,%rdi,4),%xmm0
movss  %xmm0,0x4(%rdx,%rdi,4)
add    $0x4,%rdi
cmp    %rdi,%rcx
jne    11c0 <sprstx+0xb0>
jmp    1318 <sprstx+0x208>
lea    0x4(%rbx),%rcx
lea    (%rbx,%r14,4),%rdx
add    $0x4,%rdx
lea    0x4(%r13),%r9
lea    0x4(,%r14,4),%rdi
add    %r13,%rdi
lea    0x4(%r15),%r8
lea    (%r15,%r14,4),%rsi
add    $0x4,%rsi
cmp    %rdi,%rcx
setb   %dil
cmp    %rdx,%r9
setb   %r9b
cmp    %rsi,%rcx
setb   %cl
cmp    %rdx,%r8
setb   %dl
test   %r9b,%dil
jne    1156 <sprstx+0x46>
and    %dl,%cl
jne    1156 <sprstx+0x46>
mov    %r14,%rax
and    $0xfffffffffffffff8,%rax
lea    -0x8(%rax),%rcx
mov    %rcx,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
test   %rcx,%rcx
je     1389 <sprstx+0x279>
mov    %rdx,%rsi
and    $0xfffffffffffffffe,%rsi
mov    $0x1,%ecx
movups 0x0(%r13,%rcx,4),%xmm0
movups 0x10(%r13,%rcx,4),%xmm1
movups (%r15,%rcx,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x10(%r15,%rcx,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,(%rbx,%rcx,4)
movups %xmm0,0x10(%rbx,%rcx,4)
movups 0x20(%r13,%rcx,4),%xmm0
movups 0x30(%r13,%rcx,4),%xmm1
movups 0x20(%r15,%rcx,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x30(%r15,%rcx,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,0x20(%rbx,%rcx,4)
movups %xmm0,0x30(%rbx,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1290 <sprstx+0x180>
test   $0x1,%dl
je     1313 <sprstx+0x203>
movups 0x0(%r13,%rcx,4),%xmm0
movups 0x10(%r13,%rcx,4),%xmm1
movups (%r15,%rcx,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x10(%r15,%rcx,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,(%rbx,%rcx,4)
movups %xmm0,0x10(%rbx,%rcx,4)
cmp    %r14,%rax
jne    1380 <sprstx+0x270>
test   %r14,%r14
je     1376 <sprstx+0x266>
mov    0x8(%r12),%rax
mov    $0x1,%ecx
jmp    133c <sprstx+0x22c>
nopl   0x0(%rax)
lea    0x1(%rcx),%rdx
cmp    %r14,%rcx
mov    %rdx,%rcx
je     1376 <sprstx+0x266>
mov    %rax,%rdx
mov    0x8(%r12,%rcx,8),%rax
lea    -0x1(%rax),%rsi
cmp    %rsi,%rdx
ja     1330 <sprstx+0x220>
nopl   (%rax)
mov    (%r12,%rdx,8),%rdi
movss  0x0(%r13,%rdx,4),%xmm0
mulss  (%r15,%rcx,4),%xmm0
addss  (%rbx,%rdi,4),%xmm0
movss  %xmm0,(%rbx,%rdi,4)
add    $0x1,%rdx
cmp    %rsi,%rdx
jbe    1350 <sprstx+0x240>
jmp    1330 <sprstx+0x220>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
or     $0x1,%rax
jmp    1156 <sprstx+0x46>
mov    $0x1,%ecx
test   $0x1,%dl
jne    12ed <sprstx+0x1dd>
jmp    1313 <sprstx+0x203>

LLVM IR (-O0):
@.str = [39 x i8] c"mismatched vector and matrix in sprstx\00"
define void @sprstx(float* %0, i64* %1, float* %2, float* %3, i64 %4) {
%6 = alloca float*
%7 = alloca i64*
%8 = alloca float*
%9 = alloca float*
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
store float* %0, float** %6
store i64* %1, i64** %7
store float* %2, float** %8
store float* %3, float** %9
store i64 %4, i64* %10
%14 = load i64*, i64** %7
%15 = getelementptr inbounds i64, i64* %14, i64 1
%16 = load i64, i64* %15
%17 = load i64, i64* %10
%18 = add i64 %17, 2
%19 = icmp ne i64 %16, %18
br i1 %19, label %20, label %21
20: ; preds = %5
br label %21
21: ; preds = %20, %5
store i64 1, i64* %11
br label %22
22: ; preds = %39, %21
%23 = load i64, i64* %11
%24 = load i64, i64* %10
%25 = icmp ule i64 %23, %24
br i1 %25, label %26, label %42
26: ; preds = %22
%27 = load float*, float** %6
%28 = load i64, i64* %11
%29 = getelementptr inbounds float, float* %27, i64 %28
%30 = load float, float* %29
%31 = load float*, float** %8
%32 = load i64, i64* %11
%33 = getelementptr inbounds float, float* %31, i64 %32
%34 = load float, float* %33
%35 = fmul float %30, %34
%36 = load float*, float** %9
%37 = load i64, i64* %11
%38 = getelementptr inbounds float, float* %36, i64 %37
store float %35, float* %38
br label %39
39: ; preds = %26
%40 = load i64, i64* %11
%41 = add i64 %40, 1
store i64 %41, i64* %11
br label %22
42: ; preds = %22
store i64 1, i64* %11
br label %43
43: ; preds = %83, %42
%44 = load i64, i64* %11
%45 = load i64, i64* %10
%46 = icmp ule i64 %44, %45
br i1 %46, label %47, label %86
47: ; preds = %43
%48 = load i64*, i64** %7
%49 = load i64, i64* %11
%50 = getelementptr inbounds i64, i64* %48, i64 %49
%51 = load i64, i64* %50
store i64 %51, i64* %13
br label %52
52: ; preds = %79, %47
%53 = load i64, i64* %13
%54 = load i64*, i64** %7
%55 = load i64, i64* %11
%56 = add i64 %55, 1
%57 = getelementptr inbounds i64, i64* %54, i64 %56
%58 = load i64, i64* %57
%59 = sub i64 %58, 1
%60 = icmp ule i64 %53, %59
br i1 %60, label %61, label %82
61: ; preds = %52
%62 = load i64*, i64** %7
%63 = load i64, i64* %13
%64 = getelementptr inbounds i64, i64* %62, i64 %63
%65 = load i64, i64* %64
store i64 %65, i64* %12
%66 = load float*, float** %6
%67 = load i64, i64* %13
%68 = getelementptr inbounds float, float* %66, i64 %67
%69 = load float, float* %68
%70 = load float*, float** %8
%71 = load i64, i64* %11
%72 = getelementptr inbounds float, float* %70, i64 %71
%73 = load float, float* %72
%74 = load float*, float** %9
%75 = load i64, i64* %12
%76 = getelementptr inbounds float, float* %74, i64 %75
%77 = load float, float* %76
%78 = call float @llvm.fmuladd.f32(float %69, float %73, float %77)
store float %78, float* %76
br label %79
79: ; preds = %61
%80 = load i64, i64* %13
%81 = add i64 %80, 1
store i64 %81, i64* %13
br label %52
82: ; preds = %52
br label %83
83: ; preds = %82
%84 = load i64, i64* %11
%85 = add i64 %84, 1
store i64 %85, i64* %11
br label %43
86: ; preds = %43
ret void
}

LLVM IR (-O2):
@.str = [39 x i8] c"mismatched vector and matrix in sprstx\00"
define void @sprstx(float* nocapture %0, i64* nocapture %1, float* nocapture %2, float* nocapture %3, i64 %4) {
%6 = getelementptr inbounds i64, i64* %1, i64 1
%7 = load i64, i64* %6
%8 = add i64 %4, 2
%9 = icmp eq i64 %7, %8
br i1 %9, label %11, label %10
10: ; preds = %5
tail call void @nrerror(i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str, i64 0, i64 0))
br label %11
11: ; preds = %10, %5
%12 = icmp eq i64 %4, 0
br i1 %12, label %190, label %13
13: ; preds = %11
%14 = icmp ult i64 %4, 8
br i1 %14, label %110, label %15
15: ; preds = %13
%16 = getelementptr float, float* %3, i64 1
%17 = add i64 %4, 1
%18 = getelementptr float, float* %3, i64 %17
%19 = getelementptr float, float* %0, i64 1
%20 = getelementptr float, float* %0, i64 %17
%21 = getelementptr float, float* %2, i64 1
%22 = getelementptr float, float* %2, i64 %17
%23 = icmp ult float* %16, %20
%24 = icmp ult float* %19, %18
%25 = and i1 %23, %24
%26 = icmp ult float* %16, %22
%27 = icmp ult float* %21, %18
%28 = and i1 %26, %27
%29 = or i1 %25, %28
br i1 %29, label %110, label %30
30: ; preds = %15
%31 = and i64 %4, -8
%32 = or i64 %31, 1
%33 = add i64 %31, -8
%34 = lshr exact i64 %33, 3
%35 = add nuw nsw i64 %34, 1
%36 = and i64 %35, 1
%37 = icmp eq i64 %33, 0
br i1 %37, label %86, label %38
38: ; preds = %30
%39 = and i64 %35, 4611686018427387902
br label %40
40: ; preds = %40, %38
%41 = phi i64 [ 0, %38 ], [ %81, %40 ]
%42 = phi i64 [ 0, %38 ], [ %82, %40 ]
%43 = or i64 %41, 1
%44 = getelementptr inbounds float, float* %0, i64 %43
%45 = bitcast float* %44 to <4 x float>*
%46 = load <4 x float>, <4 x float>* %45
%47 = getelementptr inbounds float, float* %44, i64 4
%48 = bitcast float* %47 to <4 x float>*
%49 = load <4 x float>, <4 x float>* %48
%50 = getelementptr inbounds float, float* %2, i64 %43
%51 = bitcast float* %50 to <4 x float>*
%52 = load <4 x float>, <4 x float>* %51
%53 = getelementptr inbounds float, float* %50, i64 4
%54 = bitcast float* %53 to <4 x float>*
%55 = load <4 x float>, <4 x float>* %54
%56 = fmul <4 x float> %46, %52
%57 = fmul <4 x float> %49, %55
%58 = getelementptr inbounds float, float* %3, i64 %43
%59 = bitcast float* %58 to <4 x float>*
store <4 x float> %56, <4 x float>* %59
%60 = getelementptr inbounds float, float* %58, i64 4
%61 = bitcast float* %60 to <4 x float>*
store <4 x float> %57, <4 x float>* %61
%62 = or i64 %41, 9
%63 = getelementptr inbounds float, float* %0, i64 %62
%64 = bitcast float* %63 to <4 x float>*
%65 = load <4 x float>, <4 x float>* %64
%66 = getelementptr inbounds float, float* %63, i64 4
%67 = bitcast float* %66 to <4 x float>*
%68 = load <4 x float>, <4 x float>* %67
%69 = getelementptr inbounds float, float* %2, i64 %62
%70 = bitcast float* %69 to <4 x float>*
%71 = load <4 x float>, <4 x float>* %70
%72 = getelementptr inbounds float, float* %69, i64 4
%73 = bitcast float* %72 to <4 x float>*
%74 = load <4 x float>, <4 x float>* %73
%75 = fmul <4 x float> %65, %71
%76 = fmul <4 x float> %68, %74
%77 = getelementptr inbounds float, float* %3, i64 %62
%78 = bitcast float* %77 to <4 x float>*
store <4 x float> %75, <4 x float>* %78
%79 = getelementptr inbounds float, float* %77, i64 4
%80 = bitcast float* %79 to <4 x float>*
store <4 x float> %76, <4 x float>* %80
%81 = add nuw i64 %41, 16
%82 = add i64 %42, 2
%83 = icmp eq i64 %82, %39
br i1 %83, label %84, label %40
84: ; preds = %40
%85 = or i64 %81, 1
br label %86
86: ; preds = %84, %30
%87 = phi i64 [ 1, %30 ], [ %85, %84 ]
%88 = icmp eq i64 %36, 0
br i1 %88, label %108, label %89
89: ; preds = %86
%90 = getelementptr inbounds float, float* %0, i64 %87
%91 = bitcast float* %90 to <4 x float>*
%92 = load <4 x float>, <4 x float>* %91
%93 = getelementptr inbounds float, float* %90, i64 4
%94 = bitcast float* %93 to <4 x float>*
%95 = load <4 x float>, <4 x float>* %94
%96 = getelementptr inbounds float, float* %2, i64 %87
%97 = bitcast float* %96 to <4 x float>*
%98 = load <4 x float>, <4 x float>* %97
%99 = getelementptr inbounds float, float* %96, i64 4
%100 = bitcast float* %99 to <4 x float>*
%101 = load <4 x float>, <4 x float>* %100
%102 = fmul <4 x float> %92, %98
%103 = fmul <4 x float> %95, %101
%104 = getelementptr inbounds float, float* %3, i64 %87
%105 = bitcast float* %104 to <4 x float>*
store <4 x float> %102, <4 x float>* %105
%106 = getelementptr inbounds float, float* %104, i64 4
%107 = bitcast float* %106 to <4 x float>*
store <4 x float> %103, <4 x float>* %107
br label %108
108: ; preds = %86, %89
%109 = icmp eq i64 %31, %4
br i1 %109, label %132, label %110
110: ; preds = %15, %13, %108
%111 = phi i64 [ 1, %15 ], [ 1, %13 ], [ %32, %108 ]
%112 = add i64 %4, 1
%113 = sub i64 %112, %111
%114 = sub i64 %4, %111
%115 = and i64 %113, 3
%116 = icmp eq i64 %115, 0
br i1 %116, label %129, label %117
117: ; preds = %110, %117
%118 = phi i64 [ %126, %117 ], [ %111, %110 ]
%119 = phi i64 [ %127, %117 ], [ 0, %110 ]
%120 = getelementptr inbounds float, float* %0, i64 %118
%121 = load float, float* %120
%122 = getelementptr inbounds float, float* %2, i64 %118
%123 = load float, float* %122
%124 = fmul float %121, %123
%125 = getelementptr inbounds float, float* %3, i64 %118
store float %124, float* %125
%126 = add nuw i64 %118, 1
%127 = add i64 %119, 1
%128 = icmp eq i64 %127, %115
br i1 %128, label %129, label %117
129: ; preds = %117, %110
%130 = phi i64 [ %111, %110 ], [ %126, %117 ]
%131 = icmp ult i64 %114, 3
br i1 %131, label %132, label %135
132: ; preds = %129, %135, %108
br i1 %12, label %190, label %133
133: ; preds = %132
%134 = load i64, i64* %6
br label %168
135: ; preds = %129, %135
%136 = phi i64 [ %164, %135 ], [ %130, %129 ]
%137 = getelementptr inbounds float, float* %0, i64 %136
%138 = load float, float* %137
%139 = getelementptr inbounds float, float* %2, i64 %136
%140 = load float, float* %139
%141 = fmul float %138, %140
%142 = getelementptr inbounds float, float* %3, i64 %136
store float %141, float* %142
%143 = add nuw i64 %136, 1
%144 = getelementptr inbounds float, float* %0, i64 %143
%145 = load float, float* %144
%146 = getelementptr inbounds float, float* %2, i64 %143
%147 = load float, float* %146
%148 = fmul float %145, %147
%149 = getelementptr inbounds float, float* %3, i64 %143
store float %148, float* %149
%150 = add nuw i64 %136, 2
%151 = getelementptr inbounds float, float* %0, i64 %150
%152 = load float, float* %151
%153 = getelementptr inbounds float, float* %2, i64 %150
%154 = load float, float* %153
%155 = fmul float %152, %154
%156 = getelementptr inbounds float, float* %3, i64 %150
store float %155, float* %156
%157 = add nuw i64 %136, 3
%158 = getelementptr inbounds float, float* %0, i64 %157
%159 = load float, float* %158
%160 = getelementptr inbounds float, float* %2, i64 %157
%161 = load float, float* %160
%162 = fmul float %159, %161
%163 = getelementptr inbounds float, float* %3, i64 %157
store float %162, float* %163
%164 = add nuw i64 %136, 4
%165 = icmp eq i64 %157, %4
br i1 %165, label %132, label %135
166: ; preds = %178, %168
%167 = icmp eq i64 %170, %4
br i1 %167, label %190, label %168
168: ; preds = %133, %166
%169 = phi i64 [ %173, %166 ], [ %134, %133 ]
%170 = phi i64 [ %171, %166 ], [ 1, %133 ]
%171 = add nuw i64 %170, 1
%172 = getelementptr inbounds i64, i64* %1, i64 %171
%173 = load i64, i64* %172
%174 = add i64 %173, -1
%175 = icmp ugt i64 %169, %174
br i1 %175, label %166, label %176
176: ; preds = %168
%177 = getelementptr inbounds float, float* %2, i64 %170
br label %178
178: ; preds = %176, %178
%179 = phi i64 [ %169, %176 ], [ %188, %178 ]
%180 = getelementptr inbounds i64, i64* %1, i64 %179
%181 = load i64, i64* %180
%182 = getelementptr inbounds float, float* %0, i64 %179
%183 = load float, float* %182
%184 = load float, float* %177
%185 = getelementptr inbounds float, float* %3, i64 %181
%186 = load float, float* %185
%187 = tail call float @llvm.fmuladd.f32(float %183, float %184, float %186)
store float %187, float* %185
%188 = add i64 %179, 1
%189 = icmp ugt i64 %188, %174
br i1 %189, label %166, label %178
190: ; preds = %166, %11, %132
ret void
}