void dsprstx(double sa[], unsigned long ija[], double x[], double 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 dsprstx");
  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):
<dsprstx>:
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 <dsprstx+0x41>
lea    0xeb4(%rip),%rdi
call   1030 <nrerror@plt>
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     119f <dsprstx+0x8f>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1159 <dsprstx+0x49>
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     124e <dsprstx+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 <dsprstx+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
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11c5 <dsprstx+0xb5>
jmp    123d <dsprstx+0x12d>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11a7 <dsprstx+0x97>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dsprstx>:
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 <dsprstx+0x2e>
lea    0xec7(%rip),%rdi
call   1030 <nrerror@plt>
test   %r14,%r14
je     1396 <dsprstx+0x286>
mov    $0x1,%eax
cmp    $0x4,%r14
jae    1213 <dsprstx+0x103>
mov    %r14d,%edx
sub    %eax,%edx
add    $0x1,%edx
mov    %r14,%rcx
sub    %rax,%rcx
and    $0x3,%rdx
je     118c <dsprstx+0x7c>
nopw   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rax,8),%xmm0
mulsd  (%r15,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rdx
jne    1170 <dsprstx+0x60>
cmp    $0x3,%rcx
jb     1330 <dsprstx+0x220>
mov    %r14,%rcx
sub    %rax,%rcx
lea    (%rbx,%rax,8),%rdx
add    $0x18,%rdx
lea    (%r15,%rax,8),%rsi
add    $0x18,%rsi
lea    0x18(,%rax,8),%rax
add    %r13,%rax
mov    $0xffffffffffffffff,%rdi
xchg   %ax,%ax
movsd  -0x10(%rax,%rdi,8),%xmm0
mulsd  -0x10(%rsi,%rdi,8),%xmm0
movsd  %xmm0,-0x10(%rdx,%rdi,8)
movsd  -0x8(%rax,%rdi,8),%xmm0
mulsd  -0x8(%rsi,%rdi,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%rdi,8)
movsd  (%rax,%rdi,8),%xmm0
mulsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rdi,8)
movsd  0x8(%rax,%rdi,8),%xmm0
mulsd  0x8(%rsi,%rdi,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rdi,8)
add    $0x4,%rdi
cmp    %rdi,%rcx
jne    11c0 <dsprstx+0xb0>
jmp    1330 <dsprstx+0x220>
lea    0x8(%rbx),%rcx
lea    (%rbx,%r14,8),%rdx
add    $0x8,%rdx
lea    0x8(%r13),%r9
lea    0x8(,%r14,8),%rdi
add    %r13,%rdi
lea    0x8(%r15),%r8
lea    (%r15,%r14,8),%rsi
add    $0x8,%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 <dsprstx+0x46>
and    %dl,%cl
jne    1156 <dsprstx+0x46>
mov    %r14,%rax
and    $0xfffffffffffffffc,%rax
lea    -0x4(%rax),%rcx
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
test   %rcx,%rcx
je     13a9 <dsprstx+0x299>
mov    %rdx,%rsi
and    $0xfffffffffffffffe,%rsi
mov    $0x1,%ecx
movupd 0x0(%r13,%rcx,8),%xmm0
movupd 0x10(%r13,%rcx,8),%xmm1
movupd (%r15,%rcx,8),%xmm2
mulpd  %xmm0,%xmm2
movupd 0x10(%r15,%rcx,8),%xmm0
mulpd  %xmm1,%xmm0
movupd %xmm2,(%rbx,%rcx,8)
movupd %xmm0,0x10(%rbx,%rcx,8)
movupd 0x20(%r13,%rcx,8),%xmm0
movupd 0x30(%r13,%rcx,8),%xmm1
movupd 0x20(%r15,%rcx,8),%xmm2
mulpd  %xmm0,%xmm2
movupd 0x30(%r15,%rcx,8),%xmm0
mulpd  %xmm1,%xmm0
movupd %xmm2,0x20(%rbx,%rcx,8)
movupd %xmm0,0x30(%rbx,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1290 <dsprstx+0x180>
test   $0x1,%dl
je     132b <dsprstx+0x21b>
movupd 0x0(%r13,%rcx,8),%xmm0
movupd 0x10(%r13,%rcx,8),%xmm1
movupd (%r15,%rcx,8),%xmm2
mulpd  %xmm0,%xmm2
movupd 0x10(%r15,%rcx,8),%xmm0
mulpd  %xmm1,%xmm0
movupd %xmm2,(%rbx,%rcx,8)
movupd %xmm0,0x10(%rbx,%rcx,8)
cmp    %r14,%rax
jne    13a0 <dsprstx+0x290>
test   %r14,%r14
je     1396 <dsprstx+0x286>
mov    0x8(%r12),%rax
mov    $0x1,%ecx
jmp    135c <dsprstx+0x24c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0x1(%rcx),%rdx
cmp    %r14,%rcx
mov    %rdx,%rcx
je     1396 <dsprstx+0x286>
mov    %rax,%rdx
mov    0x8(%r12,%rcx,8),%rax
lea    -0x1(%rax),%rsi
cmp    %rsi,%rdx
ja     1350 <dsprstx+0x240>
nopl   (%rax)
mov    (%r12,%rdx,8),%rdi
movsd  0x0(%r13,%rdx,8),%xmm0
mulsd  (%r15,%rcx,8),%xmm0
addsd  (%rbx,%rdi,8),%xmm0
movsd  %xmm0,(%rbx,%rdi,8)
add    $0x1,%rdx
cmp    %rsi,%rdx
jbe    1370 <dsprstx+0x260>
jmp    1350 <dsprstx+0x240>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
or     $0x1,%rax
jmp    1156 <dsprstx+0x46>
mov    $0x1,%ecx
test   $0x1,%dl
jne    12fd <dsprstx+0x1ed>
jmp    132b <dsprstx+0x21b>

LLVM IR (-O0):
@.str = [40 x i8] c"mismatched vector and matrix in dsprstx\00"
define void @dsprstx(double* %0, i64* %1, double* %2, double* %3, i64 %4) {
%6 = alloca double*
%7 = alloca i64*
%8 = alloca double*
%9 = alloca double*
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
store double* %0, double** %6
store i64* %1, i64** %7
store double* %2, double** %8
store double* %3, double** %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 double*, double** %6
%28 = load i64, i64* %11
%29 = getelementptr inbounds double, double* %27, i64 %28
%30 = load double, double* %29
%31 = load double*, double** %8
%32 = load i64, i64* %11
%33 = getelementptr inbounds double, double* %31, i64 %32
%34 = load double, double* %33
%35 = fmul double %30, %34
%36 = load double*, double** %9
%37 = load i64, i64* %11
%38 = getelementptr inbounds double, double* %36, i64 %37
store double %35, double* %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 double*, double** %6
%67 = load i64, i64* %13
%68 = getelementptr inbounds double, double* %66, i64 %67
%69 = load double, double* %68
%70 = load double*, double** %8
%71 = load i64, i64* %11
%72 = getelementptr inbounds double, double* %70, i64 %71
%73 = load double, double* %72
%74 = load double*, double** %9
%75 = load i64, i64* %12
%76 = getelementptr inbounds double, double* %74, i64 %75
%77 = load double, double* %76
%78 = call double @llvm.fmuladd.f64(double %69, double %73, double %77)
store double %78, double* %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 = [40 x i8] c"mismatched vector and matrix in dsprstx\00"
define void @dsprstx(double* nocapture %0, i64* nocapture %1, double* nocapture %2, double* 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 ([40 x i8], [40 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, 4
br i1 %14, label %110, label %15
15: ; preds = %13
%16 = getelementptr double, double* %3, i64 1
%17 = add i64 %4, 1
%18 = getelementptr double, double* %3, i64 %17
%19 = getelementptr double, double* %0, i64 1
%20 = getelementptr double, double* %0, i64 %17
%21 = getelementptr double, double* %2, i64 1
%22 = getelementptr double, double* %2, i64 %17
%23 = icmp ult double* %16, %20
%24 = icmp ult double* %19, %18
%25 = and i1 %23, %24
%26 = icmp ult double* %16, %22
%27 = icmp ult double* %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, -4
%32 = or i64 %31, 1
%33 = add i64 %31, -4
%34 = lshr exact i64 %33, 2
%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, 9223372036854775806
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 double, double* %0, i64 %43
%45 = bitcast double* %44 to <2 x double>*
%46 = load <2 x double>, <2 x double>* %45
%47 = getelementptr inbounds double, double* %44, i64 2
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = getelementptr inbounds double, double* %2, i64 %43
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = getelementptr inbounds double, double* %50, i64 2
%54 = bitcast double* %53 to <2 x double>*
%55 = load <2 x double>, <2 x double>* %54
%56 = fmul <2 x double> %46, %52
%57 = fmul <2 x double> %49, %55
%58 = getelementptr inbounds double, double* %3, i64 %43
%59 = bitcast double* %58 to <2 x double>*
store <2 x double> %56, <2 x double>* %59
%60 = getelementptr inbounds double, double* %58, i64 2
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> %57, <2 x double>* %61
%62 = or i64 %41, 5
%63 = getelementptr inbounds double, double* %0, i64 %62
%64 = bitcast double* %63 to <2 x double>*
%65 = load <2 x double>, <2 x double>* %64
%66 = getelementptr inbounds double, double* %63, i64 2
%67 = bitcast double* %66 to <2 x double>*
%68 = load <2 x double>, <2 x double>* %67
%69 = getelementptr inbounds double, double* %2, i64 %62
%70 = bitcast double* %69 to <2 x double>*
%71 = load <2 x double>, <2 x double>* %70
%72 = getelementptr inbounds double, double* %69, i64 2
%73 = bitcast double* %72 to <2 x double>*
%74 = load <2 x double>, <2 x double>* %73
%75 = fmul <2 x double> %65, %71
%76 = fmul <2 x double> %68, %74
%77 = getelementptr inbounds double, double* %3, i64 %62
%78 = bitcast double* %77 to <2 x double>*
store <2 x double> %75, <2 x double>* %78
%79 = getelementptr inbounds double, double* %77, i64 2
%80 = bitcast double* %79 to <2 x double>*
store <2 x double> %76, <2 x double>* %80
%81 = add nuw i64 %41, 8
%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 double, double* %0, i64 %87
%91 = bitcast double* %90 to <2 x double>*
%92 = load <2 x double>, <2 x double>* %91
%93 = getelementptr inbounds double, double* %90, i64 2
%94 = bitcast double* %93 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = getelementptr inbounds double, double* %2, i64 %87
%97 = bitcast double* %96 to <2 x double>*
%98 = load <2 x double>, <2 x double>* %97
%99 = getelementptr inbounds double, double* %96, i64 2
%100 = bitcast double* %99 to <2 x double>*
%101 = load <2 x double>, <2 x double>* %100
%102 = fmul <2 x double> %92, %98
%103 = fmul <2 x double> %95, %101
%104 = getelementptr inbounds double, double* %3, i64 %87
%105 = bitcast double* %104 to <2 x double>*
store <2 x double> %102, <2 x double>* %105
%106 = getelementptr inbounds double, double* %104, i64 2
%107 = bitcast double* %106 to <2 x double>*
store <2 x double> %103, <2 x double>* %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 double, double* %0, i64 %118
%121 = load double, double* %120
%122 = getelementptr inbounds double, double* %2, i64 %118
%123 = load double, double* %122
%124 = fmul double %121, %123
%125 = getelementptr inbounds double, double* %3, i64 %118
store double %124, double* %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 double, double* %0, i64 %136
%138 = load double, double* %137
%139 = getelementptr inbounds double, double* %2, i64 %136
%140 = load double, double* %139
%141 = fmul double %138, %140
%142 = getelementptr inbounds double, double* %3, i64 %136
store double %141, double* %142
%143 = add nuw i64 %136, 1
%144 = getelementptr inbounds double, double* %0, i64 %143
%145 = load double, double* %144
%146 = getelementptr inbounds double, double* %2, i64 %143
%147 = load double, double* %146
%148 = fmul double %145, %147
%149 = getelementptr inbounds double, double* %3, i64 %143
store double %148, double* %149
%150 = add nuw i64 %136, 2
%151 = getelementptr inbounds double, double* %0, i64 %150
%152 = load double, double* %151
%153 = getelementptr inbounds double, double* %2, i64 %150
%154 = load double, double* %153
%155 = fmul double %152, %154
%156 = getelementptr inbounds double, double* %3, i64 %150
store double %155, double* %156
%157 = add nuw i64 %136, 3
%158 = getelementptr inbounds double, double* %0, i64 %157
%159 = load double, double* %158
%160 = getelementptr inbounds double, double* %2, i64 %157
%161 = load double, double* %160
%162 = fmul double %159, %161
%163 = getelementptr inbounds double, double* %3, i64 %157
store double %162, double* %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 double, double* %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 double, double* %0, i64 %179
%183 = load double, double* %182
%184 = load double, double* %177
%185 = getelementptr inbounds double, double* %3, i64 %181
%186 = load double, double* %185
%187 = tail call double @llvm.fmuladd.f64(double %183, double %184, double %186)
store double %187, double* %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
}