#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
float fabs(float);
void sprsin(float **a, int n, float thresh, unsigned long nmax, float sa[],
            unsigned long ija[]) {
  void nrerror(char error_text[]);
  int i, j;
  unsigned long k;
  for (j = 1; j <= n; j++)
    sa[j] = a[j][j];
  ija[1] = n + 2;
  k = n + 1;
  for (i = 1; i <= n; i++) {
    for (j = 1; j <= n; j++) {
      if (fabs(a[i][j]) >= thresh && i != j) {
        if (++k > nmax)
          nrerror("sprsin: nmax too small");
        sa[k] = a[i][j];
        ija[k] = j;
      }
    }
    ija[i + 1] = k + 1;
  }
}


x86 assembly (clang, -O0):
<sprsin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1183 <sprsin+0x63>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1147 <sprsin+0x27>
mov    -0xc(%rbp),%eax
add    $0x2,%eax
movslq %eax,%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,0x8(%rax)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
cltq
mov    %rax,-0x38(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1284 <sprsin+0x164>
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     125d <sprsin+0x13d>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
call   1040 <fabs@plt>
ucomiss -0x10(%rbp),%xmm0
jb     124a <sprsin+0x12a>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
je     124a <sprsin+0x12a>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmp    -0x18(%rbp),%rax
jbe    1218 <sprsin+0xf8>
lea    0xded(%rip),%rdi
call   1030 <nrerror@plt>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movslq -0x30(%rbp),%rdx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
jmp    124f <sprsin+0x12f>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11ba <sprsin+0x9a>
mov    -0x38(%rbp),%rdx
add    $0x1,%rdx
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11a7 <sprsin+0x87>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sprsin>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
test   %esi,%esi
jle    12bd <sprsin+0x1ad>
mov    %rdx,%r15
mov    %rdi,%r13
mov    %esi,%ebx
lea    -0x1(%rbx),%rbp
mov    %ebx,%eax
and    $0x3,%eax
mov    $0x1,%edi
cmp    $0x3,%rbp
jb     11a1 <sprsin+0x91>
and    $0xfffffffc,%ebx
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
mov    0x8(%r13,%rdi,8),%rbp
movss  0x4(%rbp,%rdi,4),%xmm1
movss  %xmm1,0x4(%rcx,%rdi,4)
mov    0x10(%r13,%rdi,8),%rbp
movss  0x8(%rbp,%rdi,4),%xmm1
movss  %xmm1,0x8(%rcx,%rdi,4)
mov    0x18(%r13,%rdi,8),%rbp
movss  0xc(%rbp,%rdi,4),%xmm1
movss  %xmm1,0xc(%rcx,%rdi,4)
mov    0x20(%r13,%rdi,8),%rbp
movss  0x10(%rbp,%rdi,4),%xmm1
movss  %xmm1,0x10(%rcx,%rdi,4)
add    $0x4,%rdi
cmp    %rdi,%rbx
jne    1150 <sprsin+0x40>
add    $0x1,%rdi
test   %rax,%rax
je     11ca <sprsin+0xba>
cs nopw 0x0(%rax,%rax,1)
mov    0x0(%r13,%rdi,8),%rbp
movss  0x0(%rbp,%rdi,4),%xmm1
movss  %xmm1,(%rcx,%rdi,4)
add    $0x1,%rdi
add    $0xffffffffffffffff,%rax
jne    11b0 <sprsin+0xa0>
lea    0x2(%rsi),%eax
cltq
mov    %rax,0x8(%r8)
test   %esi,%esi
jle    12c7 <sprsin+0x1b7>
add    $0x1,%esi
movslq %esi,%rbx
mov    $0x1,%r10d
xor    %ebp,%ebp
movaps 0xe10(%rip),%xmm2
movss  %xmm0,0x4(%rsp)
mov    %rsi,0x10(%rsp)
mov    %r8,0x8(%rsp)
jmp    122a <sprsin+0x11a>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    0x1(%rbx),%rax
mov    %rax,0x8(%r8,%r10,8)
add    $0x1,%r10
add    $0xffffffffffffffff,%rbp
cmp    %rsi,%r10
je     12c7 <sprsin+0x1b7>
mov    0x0(%r13,%r10,8),%rax
mov    $0x1,%r14d
mov    %r10,0x18(%rsp)
jmp    125d <sprsin+0x14d>
nopl   0x0(%rax)
movss  (%rax,%r14,4),%xmm1
movss  %xmm1,0x4(%rcx,%rbx,4)
mov    %r14,0x8(%r8,%rbx,8)
mov    %r12,%rbx
add    $0x1,%r14
cmp    %r14,%rsi
je     1210 <sprsin+0x100>
movss  (%rax,%r14,4),%xmm1
andps  %xmm2,%xmm1
ucomiss %xmm0,%xmm1
jb     1254 <sprsin+0x144>
lea    (%r14,%rbp,1),%rdi
cmp    $0x1,%rdi
je     1254 <sprsin+0x144>
lea    0x1(%rbx),%r12
cmp    %r15,%r12
jbe    1240 <sprsin+0x130>
lea    0xd8b(%rip),%rdi
mov    %rbx,0x20(%rsp)
mov    %rcx,%rbx
call   1030 <nrerror@plt>
movaps 0xd67(%rip),%xmm2
mov    0x18(%rsp),%r10
movss  0x4(%rsp),%xmm0
mov    0x10(%rsp),%rsi
mov    %rbx,%rcx
mov    0x20(%rsp),%rbx
mov    0x8(%rsp),%r8
mov    0x0(%r13,%r10,8),%rax
jmp    1240 <sprsin+0x130>
add    $0x2,%esi
movslq %esi,%rax
mov    %rax,0x8(%r8)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [23 x i8] c"sprsin: nmax too small\00"
define void @sprsin(float** %0, i32 %1, float %2, i64 %3, float* %4, i64* %5) {
%7 = alloca float**
%8 = alloca i32
%9 = alloca float
%10 = alloca i64
%11 = alloca float*
%12 = alloca i64*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i64
store float** %0, float*** %7
store i32 %1, i32* %8
store float %2, float* %9
store i64 %3, i64* %10
store float* %4, float** %11
store i64* %5, i64** %12
store i32 1, i32* %14
br label %16
16: ; preds = %34, %6
%17 = load i32, i32* %14
%18 = load i32, i32* %8
%19 = icmp sle i32 %17, %18
br i1 %19, label %20, label %37
20: ; preds = %16
%21 = load float**, float*** %7
%22 = load i32, i32* %14
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds float*, float** %21, i64 %23
%25 = load float*, float** %24
%26 = load i32, i32* %14
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds float, float* %25, i64 %27
%29 = load float, float* %28
%30 = load float*, float** %11
%31 = load i32, i32* %14
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %30, i64 %32
store float %29, float* %33
br label %34
34: ; preds = %20
%35 = load i32, i32* %14
%36 = add nsw i32 %35, 1
store i32 %36, i32* %14
br label %16
37: ; preds = %16
%38 = load i32, i32* %8
%39 = add nsw i32 %38, 2
%40 = sext i32 %39 to i64
%41 = load i64*, i64** %12
%42 = getelementptr inbounds i64, i64* %41, i64 1
store i64 %40, i64* %42
%43 = load i32, i32* %8
%44 = add nsw i32 %43, 1
%45 = sext i32 %44 to i64
store i64 %45, i64* %15
store i32 1, i32* %13
br label %46
46: ; preds = %108, %37
%47 = load i32, i32* %13
%48 = load i32, i32* %8
%49 = icmp sle i32 %47, %48
br i1 %49, label %50, label %111
50: ; preds = %46
store i32 1, i32* %14
br label %51
51: ; preds = %97, %50
%52 = load i32, i32* %14
%53 = load i32, i32* %8
%54 = icmp sle i32 %52, %53
br i1 %54, label %55, label %100
55: ; preds = %51
%56 = load float**, float*** %7
%57 = load i32, i32* %13
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float*, float** %56, i64 %58
%60 = load float*, float** %59
%61 = load i32, i32* %14
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %60, i64 %62
%64 = load float, float* %63
%65 = call float @fabs(float %64)
%66 = load float, float* %9
%67 = fcmp oge float %65, %66
br i1 %67, label %68, label %96
68: ; preds = %55
%69 = load i32, i32* %13
%70 = load i32, i32* %14
%71 = icmp ne i32 %69, %70
br i1 %71, label %72, label %96
72: ; preds = %68
%73 = load i64, i64* %15
%74 = add i64 %73, 1
store i64 %74, i64* %15
%75 = load i64, i64* %10
%76 = icmp ugt i64 %74, %75
br i1 %76, label %77, label %78
77: ; preds = %72
br label %78
78: ; preds = %77, %72
%79 = load float**, float*** %7
%80 = load i32, i32* %13
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float*, float** %79, i64 %81
%83 = load float*, float** %82
%84 = load i32, i32* %14
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds float, float* %83, i64 %85
%87 = load float, float* %86
%88 = load float*, float** %11
%89 = load i64, i64* %15
%90 = getelementptr inbounds float, float* %88, i64 %89
store float %87, float* %90
%91 = load i32, i32* %14
%92 = sext i32 %91 to i64
%93 = load i64*, i64** %12
%94 = load i64, i64* %15
%95 = getelementptr inbounds i64, i64* %93, i64 %94
store i64 %92, i64* %95
br label %96
96: ; preds = %78, %68, %55
br label %97
97: ; preds = %96
%98 = load i32, i32* %14
%99 = add nsw i32 %98, 1
store i32 %99, i32* %14
br label %51
100: ; preds = %51
%101 = load i64, i64* %15
%102 = add i64 %101, 1
%103 = load i64*, i64** %12
%104 = load i32, i32* %13
%105 = add nsw i32 %104, 1
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds i64, i64* %103, i64 %106
store i64 %102, i64* %107
br label %108
108: ; preds = %100
%109 = load i32, i32* %13
%110 = add nsw i32 %109, 1
store i32 %110, i32* %13
br label %46
111: ; preds = %46
ret void
}

LLVM IR (-O2):
@.str = [23 x i8] c"sprsin: nmax too small\00"
define void @sprsin(float** nocapture %0, i32 %1, float %2, i64 %3, float* nocapture writeonly %4, i64* nocapture writeonly %5) {
%7 = icmp slt i32 %1, 1
br i1 %7, label %8, label %12
8: ; preds = %6
%9 = add nsw i32 %1, 2
%10 = sext i32 %9 to i64
%11 = getelementptr inbounds i64, i64* %5, i64 1
store i64 %10, i64* %11
br label %110
12: ; preds = %6
%13 = add nuw i32 %1, 1
%14 = zext i32 %13 to i64
%15 = add nsw i64 %14, -1
%16 = add nsw i64 %14, -2
%17 = and i64 %15, 3
%18 = icmp ult i64 %16, 3
br i1 %18, label %50, label %19
19: ; preds = %12
%20 = and i64 %15, -4
br label %21
21: ; preds = %21, %19
%22 = phi i64 [ 1, %19 ], [ %47, %21 ]
%23 = phi i64 [ 0, %19 ], [ %48, %21 ]
%24 = getelementptr inbounds float*, float** %0, i64 %22
%25 = load float*, float** %24
%26 = getelementptr inbounds float, float* %25, i64 %22
%27 = load float, float* %26
%28 = getelementptr inbounds float, float* %4, i64 %22
store float %27, float* %28
%29 = add nuw nsw i64 %22, 1
%30 = getelementptr inbounds float*, float** %0, i64 %29
%31 = load float*, float** %30
%32 = getelementptr inbounds float, float* %31, i64 %29
%33 = load float, float* %32
%34 = getelementptr inbounds float, float* %4, i64 %29
store float %33, float* %34
%35 = add nuw nsw i64 %22, 2
%36 = getelementptr inbounds float*, float** %0, i64 %35
%37 = load float*, float** %36
%38 = getelementptr inbounds float, float* %37, i64 %35
%39 = load float, float* %38
%40 = getelementptr inbounds float, float* %4, i64 %35
store float %39, float* %40
%41 = add nuw nsw i64 %22, 3
%42 = getelementptr inbounds float*, float** %0, i64 %41
%43 = load float*, float** %42
%44 = getelementptr inbounds float, float* %43, i64 %41
%45 = load float, float* %44
%46 = getelementptr inbounds float, float* %4, i64 %41
store float %45, float* %46
%47 = add nuw nsw i64 %22, 4
%48 = add i64 %23, 4
%49 = icmp eq i64 %48, %20
br i1 %49, label %50, label %21
50: ; preds = %21, %12
%51 = phi i64 [ 1, %12 ], [ %47, %21 ]
%52 = icmp eq i64 %17, 0
br i1 %52, label %64, label %53
53: ; preds = %50, %53
%54 = phi i64 [ %61, %53 ], [ %51, %50 ]
%55 = phi i64 [ %62, %53 ], [ 0, %50 ]
%56 = getelementptr inbounds float*, float** %0, i64 %54
%57 = load float*, float** %56
%58 = getelementptr inbounds float, float* %57, i64 %54
%59 = load float, float* %58
%60 = getelementptr inbounds float, float* %4, i64 %54
store float %59, float* %60
%61 = add nuw nsw i64 %54, 1
%62 = add i64 %55, 1
%63 = icmp eq i64 %62, %17
br i1 %63, label %64, label %53
64: ; preds = %53, %50
%65 = add nsw i32 %1, 2
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds i64, i64* %5, i64 1
store i64 %66, i64* %67
br i1 %7, label %110, label %68
68: ; preds = %64
%69 = add nuw i32 %1, 1
%70 = sext i32 %69 to i64
%71 = zext i32 %69 to i64
br label %72
72: ; preds = %105, %68
%73 = phi i64 [ 1, %68 ], [ %107, %105 ]
%74 = phi i64 [ %70, %68 ], [ %102, %105 ]
%75 = getelementptr inbounds float*, float** %0, i64 %73
%76 = load float*, float** %75
br label %77
77: ; preds = %72, %99
%78 = phi float* [ %76, %72 ], [ %100, %99 ]
%79 = phi float* [ %76, %72 ], [ %101, %99 ]
%80 = phi i64 [ 1, %72 ], [ %103, %99 ]
%81 = phi i64 [ %74, %72 ], [ %102, %99 ]
%82 = getelementptr inbounds float, float* %79, i64 %80
%83 = load float, float* %82
%84 = tail call float @llvm.fabs.f32(float %83)
%85 = fcmp ult float %84, %2
%86 = icmp eq i64 %73, %80
%87 = select i1 %85, i1 true, i1 %86
br i1 %87, label %99, label %88
88: ; preds = %77
%89 = add i64 %81, 1
%90 = icmp ugt i64 %89, %3
br i1 %90, label %91, label %93
91: ; preds = %88
tail call void @nrerror(i8* getelementptr inbounds ([23 x i8], [23 x i8]* @.str, i64 0, i64 0))
%92 = load float*, float** %75
br label %93
93: ; preds = %91, %88
%94 = phi float* [ %92, %91 ], [ %78, %88 ]
%95 = getelementptr inbounds float, float* %94, i64 %80
%96 = load float, float* %95
%97 = getelementptr inbounds float, float* %4, i64 %89
store float %96, float* %97
%98 = getelementptr inbounds i64, i64* %5, i64 %89
store i64 %80, i64* %98
br label %99
99: ; preds = %77, %93
%100 = phi float* [ %94, %93 ], [ %78, %77 ]
%101 = phi float* [ %94, %93 ], [ %79, %77 ]
%102 = phi i64 [ %89, %93 ], [ %81, %77 ]
%103 = add nuw nsw i64 %80, 1
%104 = icmp eq i64 %103, %71
br i1 %104, label %105, label %77
105: ; preds = %99
%106 = add i64 %102, 1
%107 = add nuw nsw i64 %73, 1
%108 = getelementptr inbounds i64, i64* %5, i64 %107
store i64 %106, i64* %108
%109 = icmp eq i64 %107, %71
br i1 %109, label %110, label %72
110: ; preds = %105, %8, %64
ret void
}