#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(double);
int /*<<< orphan*/ t_error(char *);
void sprsin(double **a, int n, float thresh, long nmax, double sa[],
            long ija[]) {
  int i, j;
  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)
          t_error("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
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
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
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
ucomiss -0x10(%rbp),%xmm0
jbe    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
jle    1218 <sprsin+0xf8>
lea    0xded(%rip),%rdi
call   1040 <t_error@plt>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
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
mov    %rdx,0x20(%rsp)
movss  %xmm0,0xc(%rsp)
mov    %esi,%r13d
test   %esi,%esi
jle    12b1 <sprsin+0x191>
mov    %rdi,%rbx
mov    %r13d,%edx
lea    -0x1(%rdx),%rdi
mov    %edx,%eax
and    $0x3,%eax
mov    $0x1,%ebp
cmp    $0x3,%rdi
jb     11c1 <sprsin+0xa1>
and    $0xfffffffc,%edx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0x8(%rbx,%rbp,8),%rdi
movsd  0x8(%rdi,%rbp,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rbp,8)
mov    0x10(%rbx,%rbp,8),%rdi
movsd  0x10(%rdi,%rbp,8),%xmm0
movsd  %xmm0,0x10(%rcx,%rbp,8)
mov    0x18(%rbx,%rbp,8),%rdi
movsd  0x18(%rdi,%rbp,8),%xmm0
movsd  %xmm0,0x18(%rcx,%rbp,8)
mov    0x20(%rbx,%rbp,8),%rdi
movsd  0x20(%rdi,%rbp,8),%xmm0
movsd  %xmm0,0x20(%rcx,%rbp,8)
add    $0x4,%rbp
cmp    %rbp,%rdx
jne    1170 <sprsin+0x50>
add    $0x1,%rbp
test   %rax,%rax
je     11e8 <sprsin+0xc8>
cs nopw 0x0(%rax,%rax,1)
mov    (%rbx,%rbp,8),%rdx
movsd  (%rdx,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%rbp,8)
add    $0x1,%rbp
add    $0xffffffffffffffff,%rax
jne    11d0 <sprsin+0xb0>
mov    %rcx,0x18(%rsp)
lea    0x2(%r13),%eax
cltq
mov    %rax,0x8(%r8)
test   %r13d,%r13d
jle    12bc <sprsin+0x19c>
add    $0x1,%r13d
movslq %r13d,%rbp
mov    $0x1,%r12d
xor    %r14d,%r14d
mov    %r8,0x10(%rsp)
jmp    123f <sprsin+0x11f>
nopw   0x0(%rax,%rax,1)
lea    0x1(%rbp),%rax
mov    0x10(%rsp),%r8
mov    %rax,0x8(%r8,%r12,8)
add    $0x1,%r12
add    $0xffffffffffffffff,%r14
cmp    %r13,%r12
je     12bc <sprsin+0x19c>
mov    $0x1,%r15d
jmp    127c <sprsin+0x15c>
nopw   0x0(%rax,%rax,1)
mov    (%rbx,%r12,8),%rax
movsd  (%rax,%r15,8),%xmm0
mov    0x18(%rsp),%rax
movsd  %xmm0,0x8(%rax,%rbp,8)
mov    0x10(%rsp),%rax
mov    %r15,0x8(%rax,%rbp,8)
add    $0x1,%rbp
add    $0x1,%r15
cmp    %r15,%r13
je     1220 <sprsin+0x100>
mov    (%rbx,%r12,8),%rax
movsd  (%rax,%r15,8),%xmm0
call   1030 <fabs@plt>
ucomiss 0xc(%rsp),%xmm0
jbe    1273 <sprsin+0x153>
lea    (%r14,%r15,1),%rax
cmp    $0x1,%rax
je     1273 <sprsin+0x153>
cmp    0x20(%rsp),%rbp
jl     1250 <sprsin+0x130>
lea    0xd56(%rip),%rdi
call   1040 <t_error@plt>
jmp    1250 <sprsin+0x130>
add    $0x2,%r13d
movslq %r13d,%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(double** %0, i32 %1, float %2, i64 %3, double* %4, i64* %5) {
%7 = alloca double**
%8 = alloca i32
%9 = alloca float
%10 = alloca i64
%11 = alloca double*
%12 = alloca i64*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i64
store double** %0, double*** %7
store i32 %1, i32* %8
store float %2, float* %9
store i64 %3, i64* %10
store double* %4, double** %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 double**, double*** %7
%22 = load i32, i32* %14
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double*, double** %21, i64 %23
%25 = load double*, double** %24
%26 = load i32, i32* %14
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = load double*, double** %11
%31 = load i32, i32* %14
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
store double %29, double* %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 = %109, %37
%47 = load i32, i32* %13
%48 = load i32, i32* %8
%49 = icmp sle i32 %47, %48
br i1 %49, label %50, label %112
50: ; preds = %46
store i32 1, i32* %14
br label %51
51: ; preds = %98, %50
%52 = load i32, i32* %14
%53 = load i32, i32* %8
%54 = icmp sle i32 %52, %53
br i1 %54, label %55, label %101
55: ; preds = %51
%56 = load double**, double*** %7
%57 = load i32, i32* %13
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double*, double** %56, i64 %58
%60 = load double*, double** %59
%61 = load i32, i32* %14
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
%64 = load double, double* %63
%65 = call float @fabs(double %64)
%66 = load float, float* %9
%67 = fcmp ogt float %65, %66
br i1 %67, label %68, label %97
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 %97
72: ; preds = %68
%73 = load i64, i64* %15
%74 = add nsw i64 %73, 1
store i64 %74, i64* %15
%75 = load i64, i64* %10
%76 = icmp sgt i64 %74, %75
br i1 %76, label %77, label %79
77: ; preds = %72
%78 = call i32 @t_error(i8* getelementptr inbounds ([23 x i8], [23 x i8]* @.str, i64 0, i64 0))
br label %79
79: ; preds = %77, %72
%80 = load double**, double*** %7
%81 = load i32, i32* %13
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds double*, double** %80, i64 %82
%84 = load double*, double** %83
%85 = load i32, i32* %14
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
%88 = load double, double* %87
%89 = load double*, double** %11
%90 = load i64, i64* %15
%91 = getelementptr inbounds double, double* %89, i64 %90
store double %88, double* %91
%92 = load i32, i32* %14
%93 = sext i32 %92 to i64
%94 = load i64*, i64** %12
%95 = load i64, i64* %15
%96 = getelementptr inbounds i64, i64* %94, i64 %95
store i64 %93, i64* %96
br label %97
97: ; preds = %79, %68, %55
br label %98
98: ; preds = %97
%99 = load i32, i32* %14
%100 = add nsw i32 %99, 1
store i32 %100, i32* %14
br label %51
101: ; preds = %51
%102 = load i64, i64* %15
%103 = add nsw i64 %102, 1
%104 = load i64*, i64** %12
%105 = load i32, i32* %13
%106 = add nsw i32 %105, 1
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds i64, i64* %104, i64 %107
store i64 %103, i64* %108
br label %109
109: ; preds = %101
%110 = load i32, i32* %13
%111 = add nsw i32 %110, 1
store i32 %111, i32* %13
br label %46
112: ; preds = %46
ret void
}

LLVM IR (-O2):
@.str = [23 x i8] c"sprsin: nmax too small\00"
define void @sprsin(double** nocapture %0, i32 %1, float %2, i64 %3, double* 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 %106
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 double*, double** %0, i64 %22
%25 = load double*, double** %24
%26 = getelementptr inbounds double, double* %25, i64 %22
%27 = load double, double* %26
%28 = getelementptr inbounds double, double* %4, i64 %22
store double %27, double* %28
%29 = add nuw nsw i64 %22, 1
%30 = getelementptr inbounds double*, double** %0, i64 %29
%31 = load double*, double** %30
%32 = getelementptr inbounds double, double* %31, i64 %29
%33 = load double, double* %32
%34 = getelementptr inbounds double, double* %4, i64 %29
store double %33, double* %34
%35 = add nuw nsw i64 %22, 2
%36 = getelementptr inbounds double*, double** %0, i64 %35
%37 = load double*, double** %36
%38 = getelementptr inbounds double, double* %37, i64 %35
%39 = load double, double* %38
%40 = getelementptr inbounds double, double* %4, i64 %35
store double %39, double* %40
%41 = add nuw nsw i64 %22, 3
%42 = getelementptr inbounds double*, double** %0, i64 %41
%43 = load double*, double** %42
%44 = getelementptr inbounds double, double* %43, i64 %41
%45 = load double, double* %44
%46 = getelementptr inbounds double, double* %4, i64 %41
store double %45, double* %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 double*, double** %0, i64 %54
%57 = load double*, double** %56
%58 = getelementptr inbounds double, double* %57, i64 %54
%59 = load double, double* %58
%60 = getelementptr inbounds double, double* %4, i64 %54
store double %59, double* %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 %106, 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 = %101, %68
%73 = phi i64 [ 1, %68 ], [ %103, %101 ]
%74 = phi i64 [ %70, %68 ], [ %98, %101 ]
%75 = getelementptr inbounds double*, double** %0, i64 %73
br label %76
76: ; preds = %72, %97
%77 = phi i64 [ 1, %72 ], [ %99, %97 ]
%78 = phi i64 [ %74, %72 ], [ %98, %97 ]
%79 = load double*, double** %75
%80 = getelementptr inbounds double, double* %79, i64 %77
%81 = load double, double* %80
%82 = tail call float @fabs(double %81)
%83 = fcmp ule float %82, %2
%84 = icmp eq i64 %73, %77
%85 = select i1 %83, i1 true, i1 %84
br i1 %85, label %97, label %86
86: ; preds = %76
%87 = add nsw i64 %78, 1
%88 = icmp slt i64 %78, %3
br i1 %88, label %91, label %89
89: ; preds = %86
%90 = tail call i32 @t_error(i8* getelementptr inbounds ([23 x i8], [23 x i8]* @.str, i64 0, i64 0))
br label %91
91: ; preds = %89, %86
%92 = load double*, double** %75
%93 = getelementptr inbounds double, double* %92, i64 %77
%94 = load double, double* %93
%95 = getelementptr inbounds double, double* %4, i64 %87
store double %94, double* %95
%96 = getelementptr inbounds i64, i64* %5, i64 %87
store i64 %77, i64* %96
br label %97
97: ; preds = %76, %91
%98 = phi i64 [ %87, %91 ], [ %78, %76 ]
%99 = add nuw nsw i64 %77, 1
%100 = icmp eq i64 %99, %71
br i1 %100, label %101, label %76
101: ; preds = %97
%102 = add nsw i64 %98, 1
%103 = add nuw nsw i64 %73, 1
%104 = getelementptr inbounds i64, i64* %5, i64 %103
store i64 %102, i64* %104
%105 = icmp eq i64 %103, %71
br i1 %105, label %106, label %72
106: ; preds = %101, %8, %64
ret void
}