void sprspm(float sa[], unsigned long ija[], float sb[], unsigned long ijb[],
            float sc[], unsigned long ijc[]) {
  void nrerror(char error_text[]);
  unsigned long i, ijma, ijmb, j, m, ma, mb, mbb, mn;
  float sum;
  if (ija[1] != ijb[1] || ija[1] != ijc[1])
    nrerror("sprspm: sizes do not match");
  for (i = 1; i <= ijc[1] - 2; i++) {
    j = m = i;
    mn = ijc[i];
    sum = sa[i] * sb[i];
    for (;;) {
      mb = ijb[j];
      for (ma = ija[i]; ma <= ija[i + 1] - 1; ma++) {
        ijma = ija[ma];
        if (ijma == j)
          sum += sa[ma] * sb[j];
        else {
          while (mb < ijb[j + 1]) {
            ijmb = ijb[mb];
            if (ijmb == i) {
              sum += sa[i] * sb[mb++];
              continue;
            } else if (ijmb < ijma) {
              mb++;
              continue;
            } else if (ijmb == ijma) {
              sum += sa[ma] * sb[mb++];
              continue;
            }
            break;
          }
        }
      }
      for (mbb = mb; mbb <= ijb[j + 1] - 1; mbb++) {
        if (ijb[mbb] == i)
          sum += sa[i] * sb[mbb];
      }
      sc[m] = sum;
      sum = 0.0;
      if (mn >= ijc[i + 1])
        break;
      j = ijc[m = mn++];
    }
  }
}


x86 assembly (clang, -O0):
<sprspm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x20(%rbp),%rcx
cmp    0x8(%rcx),%rax
jne    115f <sprspm+0x4f>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x30(%rbp),%rcx
cmp    0x8(%rcx),%rax
je     116b <sprspm+0x5b>
lea    0xe9a(%rip),%rdi
call   1030 <nrerror@plt>
movq   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    0x8(%rcx),%rcx
sub    $0x2,%rcx
cmp    %rcx,%rax
ja     144b <sprspm+0x33b>
mov    -0x38(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    %rax,-0x50(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x78(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x7c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x68(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x60(%rbp)
mov    -0x60(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x38(%rbp),%rdx
mov    0x8(%rcx,%rdx,8),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
ja     135d <sprspm+0x24d>
mov    -0x10(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x50(%rbp),%rax
jne    1254 <sprspm+0x144>
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x7c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x7c(%rbp)
jmp    1347 <sprspm+0x237>
jmp    1259 <sprspm+0x149>
mov    -0x68(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x50(%rbp),%rdx
cmp    0x8(%rcx,%rdx,8),%rax
jae    1342 <sprspm+0x232>
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x38(%rbp),%rax
jne    12ca <sprspm+0x1ba>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    %rcx,%rdx
add    $0x1,%rdx
mov    %rdx,-0x68(%rbp)
movss  (%rax,%rcx,4),%xmm2
movss  -0x7c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x7c(%rbp)
jmp    1259 <sprspm+0x149>
mov    -0x48(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    12e9 <sprspm+0x1d9>
mov    -0x68(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x68(%rbp)
jmp    1259 <sprspm+0x149>
mov    -0x48(%rbp),%rax
cmp    -0x40(%rbp),%rax
jne    1333 <sprspm+0x223>
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    %rcx,%rdx
add    $0x1,%rdx
mov    %rdx,-0x68(%rbp)
movss  (%rax,%rcx,4),%xmm2
movss  -0x7c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x7c(%rbp)
jmp    1259 <sprspm+0x149>
jmp    1338 <sprspm+0x228>
jmp    133d <sprspm+0x22d>
jmp    1342 <sprspm+0x232>
jmp    1347 <sprspm+0x237>
jmp    134c <sprspm+0x23c>
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    11e7 <sprspm+0xd7>
mov    -0x68(%rbp),%rax
mov    %rax,-0x70(%rbp)
mov    -0x70(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x50(%rbp),%rdx
mov    0x8(%rcx,%rdx,8),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
ja     13db <sprspm+0x2cb>
mov    -0x20(%rbp),%rax
mov    -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
cmp    -0x38(%rbp),%rax
jne    13c5 <sprspm+0x2b5>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x70(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x7c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x7c(%rbp)
jmp    13ca <sprspm+0x2ba>
mov    -0x70(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x70(%rbp)
jmp    1365 <sprspm+0x255>
movss  -0x7c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x7c(%rbp)
mov    -0x78(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%rdx
cmp    0x8(%rcx,%rdx,8),%rax
jb     1411 <sprspm+0x301>
jmp    1435 <sprspm+0x325>
mov    -0x30(%rbp),%rax
mov    -0x78(%rbp),%rcx
mov    %rcx,%rdx
add    $0x1,%rdx
mov    %rdx,-0x78(%rbp)
mov    %rcx,-0x58(%rbp)
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x50(%rbp)
jmp    11c7 <sprspm+0xb7>
jmp    143a <sprspm+0x32a>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1173 <sprspm+0x63>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sprspm>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r8,0x30(%rsp)
mov    %rcx,%rbx
mov    %rdx,%r12
mov    %rsi,%r13
mov    %rdi,%rbp
mov    0x8(%rsi),%rdx
cmp    0x8(%rcx),%rdx
mov    %r9,(%rsp)
jne    1143 <sprspm+0x33>
cmp    0x8(%r9),%rdx
je     1157 <sprspm+0x47>
lea    0xeb6(%rip),%rdi
call   1030 <nrerror@plt>
mov    (%rsp),%r9
mov    0x8(%r9),%rdx
cmp    $0x2,%rdx
jne    116c <sprspm+0x5c>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    -0x1(%rdx),%rax
cmp    $0x3,%rax
mov    $0x2,%ecx
cmovae %rax,%rcx
mov    %rcx,0x18(%rsp)
mov    0x8(%r13),%rax
mov    $0x1,%edi
jmp    11a1 <sprspm+0x91>
nopl   (%rax)
mov    0x20(%rsp),%rdi
cmp    0x18(%rsp),%rdi
mov    0x28(%rsp),%rax
je     115d <sprspm+0x4d>
mov    %rdx,0x10(%rsp)
mov    %rax,0x38(%rsp)
lea    0x0(,%rdi,4),%rax
add    %rbp,%rax
movss  0x0(%rbp,%rdi,4),%xmm0
mulss  (%r12,%rdi,4),%xmm0
lea    0x1(%rdi),%rcx
mov    %rcx,0x20(%rsp)
mov    0x8(%r13,%rdi,8),%rcx
mov    %rcx,0x28(%rsp)
add    $0xffffffffffffffff,%rcx
mov    0x8(%r9,%rdi,8),%rdx
mov    %rdx,0x40(%rsp)
mov    %rdi,%r9
mov    %rdi,0x8(%rsp)
nopl   0x0(%rax,%rax,1)
mov    (%rbx,%r9,8),%r8
mov    0x38(%rsp),%rsi
mov    %rsi,%rdx
cmp    %rcx,%rsi
jbe    1290 <sprspm+0x180>
mov    0x8(%rbx,%r9,8),%rdx
add    $0xffffffffffffffff,%rdx
cmp    %rdx,%r8
mov    (%rsp),%r9
jbe    1259 <sprspm+0x149>
mov    0x30(%rsp),%rdx
mov    0x8(%rsp),%rsi
movss  %xmm0,(%rdx,%rsi,4)
mov    0x40(%rsp),%rdx
mov    0x10(%rsp),%rsi
cmp    %rdx,%rsi
jae    1190 <sprspm+0x80>
mov    (%r9,%rsi,8),%r9
mov    %rsi,0x8(%rsp)
lea    0x1(%rsi),%rdx
xorps  %xmm0,%xmm0
mov    %rdx,0x10(%rsp)
jmp    11f0 <sprspm+0xe0>
add    $0x1,%r8
cmp    %rdx,%r8
ja     1217 <sprspm+0x107>
cmp    %rdi,(%rbx,%r8,8)
jne    1250 <sprspm+0x140>
movss  (%rax),%xmm1
mulss  (%r12,%r8,4),%xmm1
addss  %xmm1,%xmm0
jmp    1250 <sprspm+0x140>
nop
movss  0x0(%rbp,%rdx,4),%xmm1
mulss  (%r12,%r9,4),%xmm1
addss  %xmm0,%xmm1
add    $0x1,%rdx
movaps %xmm1,%xmm0
cmp    %rcx,%rdx
ja     1205 <sprspm+0xf5>
mov    0x0(%r13,%rdx,8),%r11
cmp    %r9,%r11
je     1270 <sprspm+0x160>
mov    0x8(%rbx,%r9,8),%r15
cmp    %r15,%r8
jae    1320 <sprspm+0x210>
lea    0x0(,%rdx,4),%r14
add    %rbp,%r14
movaps %xmm0,%xmm1
mov    %r8,%r10
nopl   0x0(%rax)
mov    (%rbx,%r10,8),%rsi
cmp    %rdi,%rsi
je     12e0 <sprspm+0x1d0>
cmp    %r11,%rsi
jae    12f0 <sprspm+0x1e0>
add    $0x1,%r10
cmp    %r15,%r10
jb     12c0 <sprspm+0x1b0>
jmp    1330 <sprspm+0x220>
nopl   0x0(%rax)
mov    %rax,%rsi
jmp    12f5 <sprspm+0x1e5>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r14,%rsi
jne    1340 <sprspm+0x230>
movss  (%rsi),%xmm0
mulss  (%r12,%r10,4),%xmm0
addss  %xmm1,%xmm0
lea    0x1(%r10),%r8
cmp    %r15,%r8
jb     12b3 <sprspm+0x1a3>
add    $0x1,%r10
mov    %r10,%r8
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movaps %xmm0,%xmm1
jmp    1280 <sprspm+0x170>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
cmp    %r8,%r15
cmova  %r15,%r8
jmp    1280 <sprspm+0x170>
mov    %r10,%r8
jmp    1280 <sprspm+0x170>

LLVM IR (-O0):
@.str = [27 x i8] c"sprspm: sizes do not match\00"
define void @sprspm(float* %0, i64* %1, float* %2, i64* %3, float* %4, i64* %5) {
%7 = alloca float*
%8 = alloca i64*
%9 = alloca float*
%10 = alloca i64*
%11 = alloca float*
%12 = alloca i64*
%13 = alloca i64
%14 = alloca i64
%15 = alloca i64
%16 = alloca i64
%17 = alloca i64
%18 = alloca i64
%19 = alloca i64
%20 = alloca i64
%21 = alloca i64
%22 = alloca float
store float* %0, float** %7
store i64* %1, i64** %8
store float* %2, float** %9
store i64* %3, i64** %10
store float* %4, float** %11
store i64* %5, i64** %12
%23 = load i64*, i64** %8
%24 = getelementptr inbounds i64, i64* %23, i64 1
%25 = load i64, i64* %24
%26 = load i64*, i64** %10
%27 = getelementptr inbounds i64, i64* %26, i64 1
%28 = load i64, i64* %27
%29 = icmp ne i64 %25, %28
br i1 %29, label %38, label %30
30: ; preds = %6
%31 = load i64*, i64** %8
%32 = getelementptr inbounds i64, i64* %31, i64 1
%33 = load i64, i64* %32
%34 = load i64*, i64** %12
%35 = getelementptr inbounds i64, i64* %34, i64 1
%36 = load i64, i64* %35
%37 = icmp ne i64 %33, %36
br i1 %37, label %38, label %39
38: ; preds = %30, %6
br label %39
39: ; preds = %38, %30
store i64 1, i64* %13
br label %40
40: ; preds = %212, %39
%41 = load i64, i64* %13
%42 = load i64*, i64** %12
%43 = getelementptr inbounds i64, i64* %42, i64 1
%44 = load i64, i64* %43
%45 = sub i64 %44, 2
%46 = icmp ule i64 %41, %45
br i1 %46, label %47, label %215
47: ; preds = %40
%48 = load i64, i64* %13
store i64 %48, i64* %17
store i64 %48, i64* %16
%49 = load i64*, i64** %12
%50 = load i64, i64* %13
%51 = getelementptr inbounds i64, i64* %49, i64 %50
%52 = load i64, i64* %51
store i64 %52, i64* %21
%53 = load float*, float** %7
%54 = load i64, i64* %13
%55 = getelementptr inbounds float, float* %53, i64 %54
%56 = load float, float* %55
%57 = load float*, float** %9
%58 = load i64, i64* %13
%59 = getelementptr inbounds float, float* %57, i64 %58
%60 = load float, float* %59
%61 = fmul float %56, %60
store float %61, float* %22
br label %62
62: ; preds = %205, %47
%63 = load i64*, i64** %10
%64 = load i64, i64* %16
%65 = getelementptr inbounds i64, i64* %63, i64 %64
%66 = load i64, i64* %65
store i64 %66, i64* %19
%67 = load i64*, i64** %8
%68 = load i64, i64* %13
%69 = getelementptr inbounds i64, i64* %67, i64 %68
%70 = load i64, i64* %69
store i64 %70, i64* %18
br label %71
71: ; preds = %156, %62
%72 = load i64, i64* %18
%73 = load i64*, i64** %8
%74 = load i64, i64* %13
%75 = add i64 %74, 1
%76 = getelementptr inbounds i64, i64* %73, i64 %75
%77 = load i64, i64* %76
%78 = sub i64 %77, 1
%79 = icmp ule i64 %72, %78
br i1 %79, label %80, label %159
80: ; preds = %71
%81 = load i64*, i64** %8
%82 = load i64, i64* %18
%83 = getelementptr inbounds i64, i64* %81, i64 %82
%84 = load i64, i64* %83
store i64 %84, i64* %14
%85 = load i64, i64* %14
%86 = load i64, i64* %16
%87 = icmp eq i64 %85, %86
br i1 %87, label %88, label %99
88: ; preds = %80
%89 = load float*, float** %7
%90 = load i64, i64* %18
%91 = getelementptr inbounds float, float* %89, i64 %90
%92 = load float, float* %91
%93 = load float*, float** %9
%94 = load i64, i64* %16
%95 = getelementptr inbounds float, float* %93, i64 %94
%96 = load float, float* %95
%97 = load float, float* %22
%98 = call float @llvm.fmuladd.f32(float %92, float %96, float %97)
store float %98, float* %22
br label %155
99: ; preds = %80
br label %100
100: ; preds = %139, %132, %116, %99
%101 = load i64, i64* %19
%102 = load i64*, i64** %10
%103 = load i64, i64* %16
%104 = add i64 %103, 1
%105 = getelementptr inbounds i64, i64* %102, i64 %104
%106 = load i64, i64* %105
%107 = icmp ult i64 %101, %106
br i1 %107, label %108, label %154
108: ; preds = %100
%109 = load i64*, i64** %10
%110 = load i64, i64* %19
%111 = getelementptr inbounds i64, i64* %109, i64 %110
%112 = load i64, i64* %111
store i64 %112, i64* %15
%113 = load i64, i64* %15
%114 = load i64, i64* %13
%115 = icmp eq i64 %113, %114
br i1 %115, label %116, label %128
116: ; preds = %108
%117 = load float*, float** %7
%118 = load i64, i64* %13
%119 = getelementptr inbounds float, float* %117, i64 %118
%120 = load float, float* %119
%121 = load float*, float** %9
%122 = load i64, i64* %19
%123 = add i64 %122, 1
store i64 %123, i64* %19
%124 = getelementptr inbounds float, float* %121, i64 %122
%125 = load float, float* %124
%126 = load float, float* %22
%127 = call float @llvm.fmuladd.f32(float %120, float %125, float %126)
store float %127, float* %22
br label %100
128: ; preds = %108
%129 = load i64, i64* %15
%130 = load i64, i64* %14
%131 = icmp ult i64 %129, %130
br i1 %131, label %132, label %135
132: ; preds = %128
%133 = load i64, i64* %19
%134 = add i64 %133, 1
store i64 %134, i64* %19
br label %100
135: ; preds = %128
%136 = load i64, i64* %15
%137 = load i64, i64* %14
%138 = icmp eq i64 %136, %137
br i1 %138, label %139, label %151
139: ; preds = %135
%140 = load float*, float** %7
%141 = load i64, i64* %18
%142 = getelementptr inbounds float, float* %140, i64 %141
%143 = load float, float* %142
%144 = load float*, float** %9
%145 = load i64, i64* %19
%146 = add i64 %145, 1
store i64 %146, i64* %19
%147 = getelementptr inbounds float, float* %144, i64 %145
%148 = load float, float* %147
%149 = load float, float* %22
%150 = call float @llvm.fmuladd.f32(float %143, float %148, float %149)
store float %150, float* %22
br label %100
151: ; preds = %135
br label %152
152: ; preds = %151
br label %153
153: ; preds = %152
br label %154
154: ; preds = %153, %100
br label %155
155: ; preds = %154, %88
br label %156
156: ; preds = %155
%157 = load i64, i64* %18
%158 = add i64 %157, 1
store i64 %158, i64* %18
br label %71
159: ; preds = %71
%160 = load i64, i64* %19
store i64 %160, i64* %20
br label %161
161: ; preds = %189, %159
%162 = load i64, i64* %20
%163 = load i64*, i64** %10
%164 = load i64, i64* %16
%165 = add i64 %164, 1
%166 = getelementptr inbounds i64, i64* %163, i64 %165
%167 = load i64, i64* %166
%168 = sub i64 %167, 1
%169 = icmp ule i64 %162, %168
br i1 %169, label %170, label %192
170: ; preds = %161
%171 = load i64*, i64** %10
%172 = load i64, i64* %20
%173 = getelementptr inbounds i64, i64* %171, i64 %172
%174 = load i64, i64* %173
%175 = load i64, i64* %13
%176 = icmp eq i64 %174, %175
br i1 %176, label %177, label %188
177: ; preds = %170
%178 = load float*, float** %7
%179 = load i64, i64* %13
%180 = getelementptr inbounds float, float* %178, i64 %179
%181 = load float, float* %180
%182 = load float*, float** %9
%183 = load i64, i64* %20
%184 = getelementptr inbounds float, float* %182, i64 %183
%185 = load float, float* %184
%186 = load float, float* %22
%187 = call float @llvm.fmuladd.f32(float %181, float %185, float %186)
store float %187, float* %22
br label %188
188: ; preds = %177, %170
br label %189
189: ; preds = %188
%190 = load i64, i64* %20
%191 = add i64 %190, 1
store i64 %191, i64* %20
br label %161
192: ; preds = %161
%193 = load float, float* %22
%194 = load float*, float** %11
%195 = load i64, i64* %17
%196 = getelementptr inbounds float, float* %194, i64 %195
store float %193, float* %196
store float 0.000000e+00, float* %22
%197 = load i64, i64* %21
%198 = load i64*, i64** %12
%199 = load i64, i64* %13
%200 = add i64 %199, 1
%201 = getelementptr inbounds i64, i64* %198, i64 %200
%202 = load i64, i64* %201
%203 = icmp uge i64 %197, %202
br i1 %203, label %204, label %205
204: ; preds = %192
br label %211
205: ; preds = %192
%206 = load i64*, i64** %12
%207 = load i64, i64* %21
%208 = add i64 %207, 1
store i64 %208, i64* %21
store i64 %207, i64* %17
%209 = getelementptr inbounds i64, i64* %206, i64 %207
%210 = load i64, i64* %209
store i64 %210, i64* %16
br label %62
211: ; preds = %204
br label %212
212: ; preds = %211
%213 = load i64, i64* %13
%214 = add i64 %213, 1
store i64 %214, i64* %13
br label %40
215: ; preds = %40
ret void
}

LLVM IR (-O2):
@.str = [27 x i8] c"sprspm: sizes do not match\00"
define void @sprspm(float* nocapture %0, i64* nocapture %1, float* nocapture %2, i64* nocapture %3, float* nocapture writeonly %4, i64* nocapture %5) {
%7 = getelementptr inbounds i64, i64* %1, i64 1
%8 = load i64, i64* %7
%9 = getelementptr inbounds i64, i64* %3, i64 1
%10 = load i64, i64* %9
%11 = icmp eq i64 %8, %10
br i1 %11, label %12, label %16
12: ; preds = %6
%13 = getelementptr inbounds i64, i64* %5, i64 1
%14 = load i64, i64* %13
%15 = icmp eq i64 %8, %14
br i1 %15, label %19, label %16
16: ; preds = %12, %6
tail call void @nrerror(i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str, i64 0, i64 0))
%17 = getelementptr inbounds i64, i64* %5, i64 1
%18 = load i64, i64* %17
br label %19
19: ; preds = %16, %12
%20 = phi i64 [ %18, %16 ], [ %8, %12 ]
%21 = icmp eq i64 %20, 2
br i1 %21, label %133, label %22
22: ; preds = %19
%23 = add i64 %20, -1
%24 = call i64 @llvm.umax.i64(i64 %23, i64 2)
%25 = load i64, i64* %7
br label %28
26: ; preds = %125
%27 = icmp eq i64 %37, %24
br i1 %27, label %133, label %28
28: ; preds = %22, %26
%29 = phi i64 [ %39, %26 ], [ %25, %22 ]
%30 = phi i64 [ %43, %26 ], [ %20, %22 ]
%31 = phi i64 [ %37, %26 ], [ 1, %22 ]
%32 = getelementptr inbounds float, float* %0, i64 %31
%33 = load float, float* %32
%34 = getelementptr inbounds float, float* %2, i64 %31
%35 = load float, float* %34
%36 = fmul float %33, %35
%37 = add nuw i64 %31, 1
%38 = getelementptr inbounds i64, i64* %1, i64 %37
%39 = load i64, i64* %38
%40 = add i64 %39, -1
%41 = icmp ugt i64 %29, %40
%42 = getelementptr inbounds i64, i64* %5, i64 %37
%43 = load i64, i64* %42
br label %44
44: ; preds = %129, %28
%45 = phi i64 [ %31, %28 ], [ %132, %129 ]
%46 = phi i64 [ %31, %28 ], [ %47, %129 ]
%47 = phi i64 [ %30, %28 ], [ %130, %129 ]
%48 = phi float [ %36, %28 ], [ 0.000000e+00, %129 ]
%49 = getelementptr inbounds i64, i64* %3, i64 %45
%50 = load i64, i64* %49
%51 = add i64 %45, 1
br i1 %41, label %55, label %52
52: ; preds = %44
%53 = getelementptr inbounds i64, i64* %3, i64 %51
%54 = getelementptr inbounds float, float* %2, i64 %45
br label %62
55: ; preds = %105, %44
%56 = phi i64 [ %50, %44 ], [ %106, %105 ]
%57 = phi float [ %48, %44 ], [ %107, %105 ]
%58 = getelementptr inbounds i64, i64* %3, i64 %51
%59 = load i64, i64* %58
%60 = add i64 %59, -1
%61 = icmp ugt i64 %56, %60
br i1 %61, label %125, label %110
62: ; preds = %52, %105
%63 = phi float [ %48, %52 ], [ %107, %105 ]
%64 = phi i64 [ %50, %52 ], [ %106, %105 ]
%65 = phi i64 [ %29, %52 ], [ %108, %105 ]
%66 = getelementptr inbounds i64, i64* %1, i64 %65
%67 = load i64, i64* %66
%68 = icmp eq i64 %67, %45
br i1 %68, label %74, label %69
69: ; preds = %62
%70 = load i64, i64* %53
%71 = icmp ult i64 %64, %70
br i1 %71, label %72, label %105
72: ; preds = %69
%73 = getelementptr inbounds float, float* %0, i64 %65
br label %92
74: ; preds = %62
%75 = getelementptr inbounds float, float* %0, i64 %65
%76 = load float, float* %75
%77 = load float, float* %54
%78 = tail call float @llvm.fmuladd.f32(float %76, float %77, float %63)
br label %105
79: ; preds = %92, %97
%80 = phi i64 [ %94, %92 ], [ %98, %97 ]
%81 = getelementptr inbounds i64, i64* %3, i64 %80
%82 = load i64, i64* %81
%83 = icmp eq i64 %82, %31
br i1 %83, label %84, label %95
84: ; preds = %79, %100
%85 = phi float* [ %73, %100 ], [ %32, %79 ]
%86 = load float, float* %85
%87 = getelementptr inbounds float, float* %2, i64 %80
%88 = load float, float* %87
%89 = tail call float @llvm.fmuladd.f32(float %86, float %88, float %93)
%90 = add nuw i64 %80, 1
%91 = icmp ult i64 %90, %70
br i1 %91, label %92, label %105
92: ; preds = %72, %84
%93 = phi float [ %63, %72 ], [ %89, %84 ]
%94 = phi i64 [ %64, %72 ], [ %90, %84 ]
br label %79
95: ; preds = %79
%96 = icmp ult i64 %82, %67
br i1 %96, label %97, label %100
97: ; preds = %95
%98 = add nuw i64 %80, 1
%99 = icmp ult i64 %98, %70
br i1 %99, label %79, label %102
100: ; preds = %95
%101 = icmp eq i64 %82, %67
br i1 %101, label %84, label %105
102: ; preds = %97
%103 = add nuw i64 %94, 1
%104 = call i64 @llvm.umax.i64(i64 %70, i64 %103)
br label %105
105: ; preds = %84, %100, %102, %69, %74
%106 = phi i64 [ %64, %74 ], [ %64, %69 ], [ %104, %102 ], [ %90, %84 ], [ %80, %100 ]
%107 = phi float [ %78, %74 ], [ %63, %69 ], [ %93, %102 ], [ %89, %84 ], [ %93, %100 ]
%108 = add i64 %65, 1
%109 = icmp ugt i64 %108, %40
br i1 %109, label %55, label %62
110: ; preds = %55, %121
%111 = phi float [ %122, %121 ], [ %57, %55 ]
%112 = phi i64 [ %123, %121 ], [ %56, %55 ]
%113 = getelementptr inbounds i64, i64* %3, i64 %112
%114 = load i64, i64* %113
%115 = icmp eq i64 %114, %31
br i1 %115, label %116, label %121
116: ; preds = %110
%117 = load float, float* %32
%118 = getelementptr inbounds float, float* %2, i64 %112
%119 = load float, float* %118
%120 = tail call float @llvm.fmuladd.f32(float %117, float %119, float %111)
br label %121
121: ; preds = %110, %116
%122 = phi float [ %120, %116 ], [ %111, %110 ]
%123 = add i64 %112, 1
%124 = icmp ugt i64 %123, %60
br i1 %124, label %125, label %110
125: ; preds = %121, %55
%126 = phi float [ %57, %55 ], [ %122, %121 ]
%127 = getelementptr inbounds float, float* %4, i64 %46
store float %126, float* %127
%128 = icmp ult i64 %47, %43
br i1 %128, label %129, label %26
129: ; preds = %125
%130 = add nuw i64 %47, 1
%131 = getelementptr inbounds i64, i64* %5, i64 %47
%132 = load i64, i64* %131
br label %44
133: ; preds = %26, %19
ret void
}