#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 */
int abs(int);
int /*<<< orphan*/ dallocateA(int, int, double **, int **, int **);
int /*<<< orphan*/ printf(char *);
void dband(int n, int b, int nonz, double **nzval, int **rowind, int **colptr) {
  int iseed[] = {1992, 1993, 1994, 1995};
  register int i, j, ub, lb, ilow, ihigh, lasta = 0;
  double *a;
  int *asub, *xa;
  double *val;
  int *row;
  extern double dlaran_();
  printf("A banded matrix.");
  dallocateA(n, nonz, nzval, rowind, colptr);
  a = *nzval;
  asub = *rowind;
  xa = *colptr;
  ub = lb = b;
  for (i = 0; i < 4; ++i)
    iseed[i] = abs(iseed[i]) % 4096;
  if (iseed[3] % 2 != 1)
    ++iseed[3];
  for (j = 0; j < n; ++j) {
    xa[j] = lasta;
    val = &a[lasta];
    row = &asub[lasta];
    ilow = ((0) > (j - ub) ? (0) : (j - ub));
    ihigh = ((n - 1) < (j + lb) ? (n - 1) : (j + lb));
    for (i = ilow; i <= ihigh; ++i) {
      val[i - ilow] = dlaran_(iseed);
      row[i - ilow] = i;
    }
    lasta += ihigh - ilow + 1;
  }
  xa[n] = lasta;
}


x86 assembly (clang, -O0):
<dband>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    0xe99(%rip),%rax
mov    %rax,-0x40(%rbp)
mov    0xe96(%rip),%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x5c(%rbp)
lea    0xe8c(%rip),%rdi
call   1050 <printf@plt>
mov    -0x4(%rbp),%edi
mov    -0xc(%rbp),%esi
mov    -0x18(%rbp),%rdx
mov    -0x20(%rbp),%rcx
mov    -0x28(%rbp),%r8
call   1060 <dallocateA@plt>
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x68(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x70(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x78(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x44(%rbp)
cmpl   $0x4,-0x44(%rbp)
jge    1206 <dband+0xc6>
movslq -0x44(%rbp),%rax
mov    -0x40(%rbp,%rax,4),%edi
call   1040 <abs@plt>
mov    $0x1000,%ecx
cltd
idiv   %ecx
movslq -0x44(%rbp),%rax
mov    %edx,-0x40(%rbp,%rax,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11d1 <dband+0x91>
mov    -0x34(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x1,%edx
je     1223 <dband+0xe3>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1349 <dband+0x209>
mov    -0x5c(%rbp),%edx
mov    -0x78(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x68(%rbp),%rax
movslq -0x5c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x80(%rbp)
mov    -0x70(%rbp),%rax
movslq -0x5c(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x88(%rbp)
mov    -0x48(%rbp),%ecx
sub    -0x4c(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    128a <dband+0x14a>
xor    %eax,%eax
mov    %eax,-0x8c(%rbp)
jmp    1296 <dband+0x156>
mov    -0x48(%rbp),%eax
sub    -0x4c(%rbp),%eax
mov    %eax,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    -0x48(%rbp),%ecx
add    -0x50(%rbp),%ecx
cmp    %ecx,%eax
jge    12c4 <dband+0x184>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    12d0 <dband+0x190>
mov    -0x48(%rbp),%eax
add    -0x50(%rbp),%eax
mov    %eax,-0x90(%rbp)
mov    -0x90(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x58(%rbp),%eax
jg     132c <dband+0x1ec>
lea    -0x40(%rbp),%rdi
mov    $0x0,%al
call   1030 <dlaran_@plt>
mov    -0x80(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    -0x54(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%edx
mov    -0x88(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    -0x54(%rbp),%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12df <dband+0x19f>
mov    -0x58(%rbp),%eax
sub    -0x54(%rbp),%eax
add    $0x1,%eax
add    -0x5c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    122a <dband+0xea>
mov    -0x5c(%rbp),%edx
mov    -0x78(%rbp),%rax
movslq -0x4(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dband>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r9,%r15
mov    %r8,%r13
mov    %rcx,%rbp
mov    %edx,%r12d
mov    %esi,%ebx
mov    %edi,%r14d
lea    0xeba(%rip),%rdi
call   1040 <printf@plt>
mov    %r14d,%edi
mov    %r12d,%esi
mov    %rbp,%rdx
mov    %r13,%rcx
mov    %r15,%r8
call   1050 <dallocateA@plt>
mov    0x0(%rbp),%r12
mov    0x0(%r13),%rbp
mov    (%r15),%rax
mov    %rax,0x18(%rsp)
movaps 0xe7a(%rip),%xmm0
movaps %xmm0,0x50(%rsp)
test   %r14d,%r14d
jle    128a <dband+0x15a>
lea    -0x1(%r14),%eax
mov    %eax,0x10(%rsp)
mov    %ebx,%eax
neg    %eax
mov    %eax,0xc(%rsp)
mov    %r14d,%eax
mov    %rax,0x30(%rsp)
lea    0x50(%rsp),%r13
xor    %ecx,%ecx
xor    %edx,%edx
mov    %r14,0x40(%rsp)
mov    %rbx,0x38(%rsp)
jmp    1206 <dband+0xd6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x28(%rsp),%rdx
sub    0x14(%rsp),%edx
mov    0x20(%rsp),%rax
add    %eax,%edx
add    $0x1,%edx
mov    0x48(%rsp),%rcx
add    $0x1,%rcx
addl   $0x1,0xc(%rsp)
cmp    0x30(%rsp),%rcx
mov    0x40(%rsp),%r14
mov    0x38(%rsp),%rbx
je     128c <dband+0x15c>
mov    0x18(%rsp),%rax
mov    %rdx,0x28(%rsp)
mov    %edx,(%rax,%rcx,4)
mov    %ecx,%edx
sub    %ebx,%edx
mov    $0x0,%eax
cmovs  %eax,%edx
mov    %rcx,0x48(%rsp)
add    %ebx,%ecx
cmp    %r14d,%ecx
cmovge 0x10(%rsp),%ecx
mov    %edx,0x14(%rsp)
mov    %rcx,%rax
mov    %rcx,0x20(%rsp)
cmp    %ecx,%edx
jg     11d0 <dband+0xa0>
movslq 0x28(%rsp),%rbx
mov    0xc(%rsp),%eax
test   %eax,%eax
mov    $0x0,%r14d
cmovg  %eax,%r14d
movslq 0x20(%rsp),%r15
nopl   0x0(%rax,%rax,1)
mov    %r13,%rdi
xor    %eax,%eax
call   1030 <dlaran_@plt>
movsd  %xmm0,(%r12,%rbx,8)
mov    %r14d,0x0(%rbp,%rbx,4)
lea    0x1(%r14),%rax
add    $0x1,%rbx
cmp    %r15,%r14
mov    %rax,%r14
jl     1260 <dband+0x130>
jmp    11d0 <dband+0xa0>
xor    %edx,%edx
movslq %r14d,%rax
mov    0x18(%rsp),%rcx
mov    %edx,(%rcx,%rax,4)
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@__const.dband.iseed = [4 x i32] [i32 1992, i32 1993, i32 1994, i32 1995]
@.str = [17 x i8] c"A banded matrix.\00"
define void @dband(i32 %0, i32 %1, i32 %2, double** %3, i32** %4, i32** %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double**
%11 = alloca i32**
%12 = alloca i32**
%13 = alloca [4 x i32]
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca double*
%22 = alloca i32*
%23 = alloca i32*
%24 = alloca double*
%25 = alloca i32*
store i32 %0, i32* %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store double** %3, double*** %10
store i32** %4, i32*** %11
store i32** %5, i32*** %12
%26 = bitcast [4 x i32]* %13 to i8*
store i32 0, i32* %20
%27 = call i32 @printf(i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str, i64 0, i64 0))
%28 = load i32, i32* %7
%29 = load i32, i32* %9
%30 = load double**, double*** %10
%31 = load i32**, i32*** %11
%32 = load i32**, i32*** %12
%33 = call i32 @dallocateA(i32 %28, i32 %29, double** %30, i32** %31, i32** %32)
%34 = load double**, double*** %10
%35 = load double*, double** %34
store double* %35, double** %21
%36 = load i32**, i32*** %11
%37 = load i32*, i32** %36
store i32* %37, i32** %22
%38 = load i32**, i32*** %12
%39 = load i32*, i32** %38
store i32* %39, i32** %23
%40 = load i32, i32* %8
store i32 %40, i32* %17
store i32 %40, i32* %16
store i32 0, i32* %14
br label %41
41: ; preds = %54, %6
%42 = load i32, i32* %14
%43 = icmp slt i32 %42, 4
br i1 %43, label %44, label %57
44: ; preds = %41
%45 = load i32, i32* %14
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds [4 x i32], [4 x i32]* %13, i64 0, i64 %46
%48 = load i32, i32* %47
%49 = call i32 @abs(i32 %48)
%50 = srem i32 %49, 4096
%51 = load i32, i32* %14
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds [4 x i32], [4 x i32]* %13, i64 0, i64 %52
store i32 %50, i32* %53
br label %54
54: ; preds = %44
%55 = load i32, i32* %14
%56 = add nsw i32 %55, 1
store i32 %56, i32* %14
br label %41
57: ; preds = %41
%58 = getelementptr inbounds [4 x i32], [4 x i32]* %13, i64 0, i64 3
%59 = load i32, i32* %58
%60 = srem i32 %59, 2
%61 = icmp ne i32 %60, 1
br i1 %61, label %62, label %66
62: ; preds = %57
%63 = getelementptr inbounds [4 x i32], [4 x i32]* %13, i64 0, i64 3
%64 = load i32, i32* %63
%65 = add nsw i32 %64, 1
store i32 %65, i32* %63
br label %66
66: ; preds = %62, %57
store i32 0, i32* %15
br label %67
67: ; preds = %142, %66
%68 = load i32, i32* %15
%69 = load i32, i32* %7
%70 = icmp slt i32 %68, %69
br i1 %70, label %71, label %145
71: ; preds = %67
%72 = load i32, i32* %20
%73 = load i32*, i32** %23
%74 = load i32, i32* %15
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds i32, i32* %73, i64 %75
store i32 %72, i32* %76
%77 = load double*, double** %21
%78 = load i32, i32* %20
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %77, i64 %79
store double* %80, double** %24
%81 = load i32*, i32** %22
%82 = load i32, i32* %20
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds i32, i32* %81, i64 %83
store i32* %84, i32** %25
%85 = load i32, i32* %15
%86 = load i32, i32* %16
%87 = sub nsw i32 %85, %86
%88 = icmp sgt i32 0, %87
br i1 %88, label %89, label %90
89: ; preds = %71
br label %94
90: ; preds = %71
%91 = load i32, i32* %15
%92 = load i32, i32* %16
%93 = sub nsw i32 %91, %92
br label %94
94: ; preds = %90, %89
%95 = phi i32 [ 0, %89 ], [ %93, %90 ]
store i32 %95, i32* %18
%96 = load i32, i32* %7
%97 = sub nsw i32 %96, 1
%98 = load i32, i32* %15
%99 = load i32, i32* %17
%100 = add nsw i32 %98, %99
%101 = icmp slt i32 %97, %100
br i1 %101, label %102, label %105
102: ; preds = %94
%103 = load i32, i32* %7
%104 = sub nsw i32 %103, 1
br label %109
105: ; preds = %94
%106 = load i32, i32* %15
%107 = load i32, i32* %17
%108 = add nsw i32 %106, %107
br label %109
109: ; preds = %105, %102
%110 = phi i32 [ %104, %102 ], [ %108, %105 ]
store i32 %110, i32* %19
%111 = load i32, i32* %18
store i32 %111, i32* %14
br label %112
112: ; preds = %132, %109
%113 = load i32, i32* %14
%114 = load i32, i32* %19
%115 = icmp sle i32 %113, %114
br i1 %115, label %116, label %135
116: ; preds = %112
%117 = getelementptr inbounds [4 x i32], [4 x i32]* %13, i64 0, i64 0
%118 = call double (i32*, ...) bitcast (double (...)* @dlaran_ to double (i32*, ...)*)(i32* %117)
%119 = load double*, double** %24
%120 = load i32, i32* %14
%121 = load i32, i32* %18
%122 = sub nsw i32 %120, %121
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double, double* %119, i64 %123
store double %118, double* %124
%125 = load i32, i32* %14
%126 = load i32*, i32** %25
%127 = load i32, i32* %14
%128 = load i32, i32* %18
%129 = sub nsw i32 %127, %128
%130 = sext i32 %129 to i64
%131 = getelementptr inbounds i32, i32* %126, i64 %130
store i32 %125, i32* %131
br label %132
132: ; preds = %116
%133 = load i32, i32* %14
%134 = add nsw i32 %133, 1
store i32 %134, i32* %14
br label %112
135: ; preds = %112
%136 = load i32, i32* %19
%137 = load i32, i32* %18
%138 = sub nsw i32 %136, %137
%139 = add nsw i32 %138, 1
%140 = load i32, i32* %20
%141 = add nsw i32 %140, %139
store i32 %141, i32* %20
br label %142
142: ; preds = %135
%143 = load i32, i32* %15
%144 = add nsw i32 %143, 1
store i32 %144, i32* %15
br label %67
145: ; preds = %67
%146 = load i32, i32* %20
%147 = load i32*, i32** %23
%148 = load i32, i32* %7
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds i32, i32* %147, i64 %149
store i32 %146, i32* %150
ret void
}

LLVM IR (-O2):
@.str = [17 x i8] c"A banded matrix.\00"
define void @dband(i32 %0, i32 %1, i32 %2, double** %3, i32** %4, i32** %5) {
%7 = alloca <4 x i32>
%8 = bitcast <4 x i32>* %7 to i8*
%9 = tail call i32 @printf(i8* dereferenceable(1) getelementptr inbounds ([17 x i8], [17 x i8]* @.str, i64 0, i64 0))
%10 = tail call i32 @dallocateA(i32 %0, i32 %2, double** %3, i32** %4, i32** %5)
%11 = load double*, double** %3
%12 = load i32*, i32** %4
%13 = load i32*, i32** %5
%14 = getelementptr inbounds <4 x i32>, <4 x i32>* %7, i64 0, i64 0
store <4 x i32> <i32 1992, i32 1993, i32 1994, i32 1995>, <4 x i32>* %7
%15 = icmp sgt i32 %0, 0
br i1 %15, label %16, label %57
16: ; preds = %6
%17 = add nsw i32 %0, -1
%18 = sub i32 0, %1
%19 = zext i32 %0 to i64
br label %20
20: ; preds = %16, %50
%21 = phi i64 [ 0, %16 ], [ %54, %50 ]
%22 = phi i32 [ %18, %16 ], [ %55, %50 ]
%23 = phi i32 [ 0, %16 ], [ %53, %50 ]
%24 = getelementptr inbounds i32, i32* %13, i64 %21
store i32 %23, i32* %24
%25 = sext i32 %23 to i64
%26 = getelementptr inbounds double, double* %11, i64 %25
%27 = getelementptr inbounds i32, i32* %12, i64 %25
%28 = trunc i64 %21 to i32
%29 = sub nsw i32 %28, %1
%30 = icmp sgt i32 %29, 0
%31 = select i1 %30, i32 %29, i32 0
%32 = add nsw i32 %28, %1
%33 = icmp slt i32 %32, %0
%34 = select i1 %33, i32 %32, i32 %17
%35 = icmp sgt i32 %31, %34
br i1 %35, label %50, label %36
36: ; preds = %20
%37 = call i32 @llvm.smax.i32(i32 %22, i32 0)
%38 = zext i32 %37 to i64
%39 = zext i32 %31 to i64
%40 = sext i32 %34 to i64
br label %41
41: ; preds = %36, %41
%42 = phi i64 [ %38, %36 ], [ %48, %41 ]
%43 = call double (i32*, ...) bitcast (double (...)* @dlaran_ to double (i32*, ...)*)(i32* %14)
%44 = sub nuw nsw i64 %42, %39
%45 = getelementptr inbounds double, double* %26, i64 %44
store double %43, double* %45
%46 = getelementptr inbounds i32, i32* %27, i64 %44
%47 = trunc i64 %42 to i32
store i32 %47, i32* %46
%48 = add nuw nsw i64 %42, 1
%49 = icmp slt i64 %42, %40
br i1 %49, label %41, label %50
50: ; preds = %41, %20
%51 = add i32 %23, 1
%52 = sub i32 %51, %31
%53 = add i32 %52, %34
%54 = add nuw nsw i64 %21, 1
%55 = add i32 %22, 1
%56 = icmp eq i64 %54, %19
br i1 %56, label %57, label %20
57: ; preds = %50, %6
%58 = phi i32 [ 0, %6 ], [ %53, %50 ]
%59 = sext i32 %0 to i64
%60 = getelementptr inbounds i32, i32* %13, i64 %59
store i32 %58, i32* %60
ret void
}