#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 *, int, int, int);
void dblockdiag(int nb, int bs, int nonz, double **nzval, int **rowind,
                int **colptr) {
  int iseed[] = {1992, 1993, 1994, 1995};
  register int i, j, b, n, lasta = 0, cstart, rstart;
  double *a;
  int *asub, *xa;
  double *val;
  int *row;
  extern double dlaran_();
  n = bs * nb;
  printf("A block diagonal matrix: nb %d, bs %d, n %d\n", nb, bs, n);
  dallocateA(n, nonz, nzval, rowind, colptr);
  a = *nzval;
  asub = *rowind;
  xa = *colptr;
  for (i = 0; i < 4; ++i)
    iseed[i] = abs(iseed[i]) % 4096;
  if (iseed[3] % 2 != 1)
    ++iseed[3];
  for (b = 0; b < nb; ++b) {
    cstart = b * bs;
    rstart = b * bs;
    for (j = cstart; j < cstart + bs; ++j) {
      xa[j] = lasta;
      val = &a[lasta];
      row = &asub[lasta];
      for (i = 0; i < bs; ++i) {
        val[i] = dlaran_(iseed);
        row[i] = i + rstart;
      }
      lasta += bs;
    }
  }
  xa[n] = lasta;
}


x86 assembly (clang, -O0):
<dblockdiag>:
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,-0x54(%rbp)
mov    -0x8(%rbp),%eax
imul   -0x4(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x4(%rbp),%esi
mov    -0x8(%rbp),%edx
mov    -0x50(%rbp),%ecx
lea    0xe79(%rip),%rdi
call   1050 <printf@plt>
mov    -0x50(%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)
movl   $0x0,-0x44(%rbp)
cmpl   $0x4,-0x44(%rbp)
jge    1210 <dblockdiag+0xd0>
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    11db <dblockdiag+0x9b>
mov    -0x34(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x1,%edx
je     122d <dblockdiag+0xed>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1319 <dblockdiag+0x1d9>
mov    -0x4c(%rbp),%eax
imul   -0x8(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x4c(%rbp),%eax
imul   -0x8(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    -0x58(%rbp),%ecx
add    -0x8(%rbp),%ecx
cmp    %ecx,%eax
jge    1306 <dblockdiag+0x1c6>
mov    -0x54(%rbp),%edx
mov    -0x78(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x68(%rbp),%rax
movslq -0x54(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x80(%rbp)
mov    -0x70(%rbp),%rax
movslq -0x54(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x88(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12ef <dblockdiag+0x1af>
lea    -0x40(%rbp),%rdi
mov    $0x0,%al
call   1030 <dlaran_@plt>
mov    -0x80(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%edx
add    -0x5c(%rbp),%edx
mov    -0x88(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12a9 <dblockdiag+0x169>
mov    -0x8(%rbp),%eax
add    -0x54(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    125a <dblockdiag+0x11a>
jmp    130b <dblockdiag+0x1cb>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1234 <dblockdiag+0xf4>
mov    -0x54(%rbp),%edx
mov    -0x78(%rbp),%rax
movslq -0x50(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dblockdiag>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9,%r14
mov    %r8,%r12
mov    %rcx,%r13
mov    %edx,(%rsp)
mov    %esi,%r15d
mov    %edi,%ebx
mov    %esi,%ebp
imul   %edi,%ebp
lea    0xeb5(%rip),%rdi
mov    %ebx,%esi
mov    %r15d,%edx
mov    %ebp,%ecx
call   1040 <printf@plt>
mov    %ebp,0x8(%rsp)
mov    %ebp,%edi
mov    (%rsp),%esi
mov    %r13,%rdx
mov    %r12,%rcx
mov    %r14,%r8
call   1050 <dallocateA@plt>
mov    0x0(%r13),%rax
mov    %rax,0x28(%rsp)
mov    (%r12),%rax
mov    %rax,0x20(%rsp)
mov    (%r14),%rax
mov    %rax,0x10(%rsp)
movaps 0xe61(%rip),%xmm0
movaps %xmm0,0x60(%rsp)
test   %ebx,%ebx
jle    12b4 <dblockdiag+0x184>
movslq %r15d,%rax
mov    %ebx,%ecx
mov    %rcx,0x30(%rsp)
mov    %r15d,%r12d
lea    0x0(,%rax,4),%rcx
mov    %rcx,0x48(%rsp)
mov    %rax,0x18(%rsp)
lea    0x0(,%rax,8),%rax
mov    %rax,0x40(%rsp)
lea    0x60(%rsp),%rbx
xor    %r14d,%r14d
xor    %eax,%eax
xor    %edx,%edx
mov    %r15d,0xc(%rsp)
jmp    1202 <dblockdiag+0xd2>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
add    %r12,%r14
cmp    0x30(%rsp),%rax
je     12b6 <dblockdiag+0x186>
test   %r15d,%r15d
jle    11f0 <dblockdiag+0xc0>
mov    %rax,0x38(%rsp)
mov    %rax,%rcx
mov    0x18(%rsp),%rax
imul   %rax,%rcx
add    %rax,%rcx
mov    %rcx,0x50(%rsp)
movslq %r14d,%rcx
movslq %edx,%rdx
mov    0x20(%rsp),%rax
lea    (%rax,%rdx,4),%r13
mov    0x28(%rsp),%rax
lea    (%rax,%rdx,8),%r15
nopl   0x0(%rax,%rax,1)
mov    0x10(%rsp),%rax
mov    %rdx,(%rsp)
mov    %rcx,0x58(%rsp)
mov    %edx,(%rax,%rcx,4)
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rbx,%rdi
xor    %eax,%eax
call   1030 <dlaran_@plt>
movsd  %xmm0,(%r15,%rbp,8)
lea    (%r14,%rbp,1),%eax
mov    %eax,0x0(%r13,%rbp,4)
add    $0x1,%rbp
cmp    %rbp,%r12
jne    1260 <dblockdiag+0x130>
mov    (%rsp),%rdx
add    0x18(%rsp),%rdx
mov    0x58(%rsp),%rcx
add    $0x1,%rcx
add    0x48(%rsp),%r13
add    0x40(%rsp),%r15
cmp    0x50(%rsp),%rcx
jl     1240 <dblockdiag+0x110>
mov    0xc(%rsp),%r15d
mov    0x38(%rsp),%rax
jmp    11f0 <dblockdiag+0xc0>
xor    %edx,%edx
movslq 0x8(%rsp),%rax
mov    0x10(%rsp),%rcx
mov    %edx,(%rcx,%rax,4)
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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