#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 */
scalar_t__ Malloc(int);
int /*<<< orphan*/ free(double *);
int linpackDTRDI(double *, int, int);
double **matrixNew(int, int);
double **matrixLowerTriangularInverse(double **matrix, int n) {
  int i, j, k, info;
  double *t, **result;
  t = (double *)Malloc(n * n * sizeof(double));
  k = 0;
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      t[k] = matrix[j][i];
      k++;
    }
  }
  info = linpackDTRDI(t, n, n);
  result = matrixNew(n, n);
  k = 0;
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      result[j][i] = i > j ? 0.0 : t[k];
      k++;
    }
  }
  free(t);
  return (result);
}


x86 assembly (clang, -O0):
<matrixLowerTriangularInverse>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <Malloc@plt>
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x18(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11e0 <matrixLowerTriangularInverse+0xa0>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11cd <matrixLowerTriangularInverse+0x8d>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1188 <matrixLowerTriangularInverse+0x48>
jmp    11d2 <matrixLowerTriangularInverse+0x92>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1175 <matrixLowerTriangularInverse+0x35>
mov    -0x28(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0xc(%rbp),%edx
call   1050 <linpackDTRDI@plt>
mov    %eax,-0x1c(%rbp)
mov    -0xc(%rbp),%edi
mov    -0xc(%rbp),%esi
call   1060 <matrixNew@plt>
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x18(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    129d <matrixLowerTriangularInverse+0x15d>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    128a <matrixLowerTriangularInverse+0x14a>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jle    1247 <matrixLowerTriangularInverse+0x107>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    1259 <matrixLowerTriangularInverse+0x119>
mov    -0x28(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1222 <matrixLowerTriangularInverse+0xe2>
jmp    128f <matrixLowerTriangularInverse+0x14f>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    120f <matrixLowerTriangularInverse+0xcf>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x30(%rbp),%rax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matrixLowerTriangularInverse>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
mov    %esi,%r15d
mov    %rdi,%rbx
mov    %esi,%edi
imul   %esi,%edi
shl    $0x3,%edi
call   1040 <Malloc@plt>
mov    %rax,%r14
mov    %r15d,%r12d
test   %r15d,%r15d
jle    1240 <matrixLowerTriangularInverse+0x100>
lea    -0x1(%r12),%r9
mov    %r12d,%ecx
and    $0x3,%ecx
mov    %r12d,%r11d
and    $0xfffffffc,%r11d
lea    0x18(%r14),%r8
xor    %edi,%edi
xor    %r10d,%r10d
jmp    119d <matrixLowerTriangularInverse+0x5d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %r12,%rdi
je     1240 <matrixLowerTriangularInverse+0x100>
movslq %r10d,%r10
cmp    $0x3,%r9
jae    11b0 <matrixLowerTriangularInverse+0x70>
xor    %eax,%eax
test   %rcx,%rcx
jne    1210 <matrixLowerTriangularInverse+0xd0>
jmp    1190 <matrixLowerTriangularInverse+0x50>
nop
lea    (%r8,%r10,8),%rbp
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
mov    (%rbx,%rax,8),%rsi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,-0x18(%rbp,%rax,8)
mov    0x8(%rbx,%rax,8),%rsi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,-0x10(%rbp,%rax,8)
mov    0x10(%rbx,%rax,8),%rsi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,-0x8(%rbp,%rax,8)
mov    0x18(%rbx,%rax,8),%rsi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    $0x4,%rax
cmp    %rax,%r11
jne    11c0 <matrixLowerTriangularInverse+0x80>
add    %rax,%r10
test   %rcx,%rcx
je     1190 <matrixLowerTriangularInverse+0x50>
lea    (%r14,%r10,8),%rbp
lea    (%rbx,%rax,8),%rsi
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
mov    (%rsi,%rax,8),%rdx
movsd  (%rdx,%rdi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    $0x1,%rax
cmp    %rax,%rcx
jne    1220 <matrixLowerTriangularInverse+0xe0>
add    %rax,%r10
jmp    1190 <matrixLowerTriangularInverse+0x50>
mov    %r14,%rdi
mov    %r15d,%esi
mov    %r15d,%edx
call   1050 <linpackDTRDI@plt>
mov    %r15d,%edi
mov    %r15d,%esi
call   1060 <matrixNew@plt>
mov    %rax,%rbx
test   %r15d,%r15d
jle    1324 <matrixLowerTriangularInverse+0x1e4>
mov    %r12d,%eax
and    $0xfffffffe,%eax
mov    %r14,%r8
add    $0x8,%r8
xor    %edx,%edx
xor    %r9d,%r9d
jmp    129a <matrixLowerTriangularInverse+0x15a>
nopl   0x0(%rax)
mov    (%rbx,%rbp,8),%rcx
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%r9
add    $0x1,%rdx
cmp    %r12,%rdx
je     1324 <matrixLowerTriangularInverse+0x1e4>
movslq %r9d,%r9
cmp    $0x1,%r15d
jne    12b0 <matrixLowerTriangularInverse+0x170>
xor    %ebp,%ebp
test   $0x1,%r12b
je     128d <matrixLowerTriangularInverse+0x14d>
jmp    130d <matrixLowerTriangularInverse+0x1cd>
nopl   (%rax)
lea    (%r8,%r9,8),%rcx
xor    %edi,%edi
jmp    12d3 <matrixLowerTriangularInverse+0x193>
nopl   0x0(%rax,%rax,1)
mov    0x8(%rbx,%rbp,8),%rsi
movsd  %xmm0,(%rsi,%rdx,8)
add    $0x1,%rdi
cmp    %rax,%rdi
je     1300 <matrixLowerTriangularInverse+0x1c0>
mov    %rdi,%rbp
xorps  %xmm0,%xmm0
xorps  %xmm1,%xmm1
cmp    %rdi,%rdx
ja     12e7 <matrixLowerTriangularInverse+0x1a7>
movsd  -0x8(%rcx,%rbp,8),%xmm1
mov    (%rbx,%rbp,8),%rdi
movsd  %xmm1,(%rdi,%rdx,8)
lea    0x1(%rbp),%rdi
cmp    %rdi,%rdx
ja     12c0 <matrixLowerTriangularInverse+0x180>
movsd  (%rcx,%rbp,8),%xmm0
jmp    12c0 <matrixLowerTriangularInverse+0x180>
add    $0x2,%rbp
add    %rbp,%r9
test   $0x1,%r12b
je     128d <matrixLowerTriangularInverse+0x14d>
xorps  %xmm0,%xmm0
cmp    %rbp,%rdx
ja     1280 <matrixLowerTriangularInverse+0x140>
movsd  (%r14,%r9,8),%xmm0
jmp    1280 <matrixLowerTriangularInverse+0x140>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %rbx,%rax
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double** @matrixLowerTriangularInverse(double** %0, i32 %1) {
%3 = alloca double**
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca double**
store double** %0, double*** %3
store i32 %1, i32* %4
%11 = load i32, i32* %4
%12 = load i32, i32* %4
%13 = mul nsw i32 %11, %12
%14 = sext i32 %13 to i64
%15 = mul i64 %14, 8
%16 = trunc i64 %15 to i32
%17 = call i64 @Malloc(i32 %16)
%18 = inttoptr i64 %17 to double*
store double* %18, double** %9
store i32 0, i32* %7
store i32 0, i32* %5
br label %19
19: ; preds = %48, %2
%20 = load i32, i32* %5
%21 = load i32, i32* %4
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %51
23: ; preds = %19
store i32 0, i32* %6
br label %24
24: ; preds = %44, %23
%25 = load i32, i32* %6
%26 = load i32, i32* %4
%27 = icmp slt i32 %25, %26
br i1 %27, label %28, label %47
28: ; preds = %24
%29 = load double**, double*** %3
%30 = load i32, i32* %6
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double*, double** %29, i64 %31
%33 = load double*, double** %32
%34 = load i32, i32* %5
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %33, i64 %35
%37 = load double, double* %36
%38 = load double*, double** %9
%39 = load i32, i32* %7
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
store double %37, double* %41
%42 = load i32, i32* %7
%43 = add nsw i32 %42, 1
store i32 %43, i32* %7
br label %44
44: ; preds = %28
%45 = load i32, i32* %6
%46 = add nsw i32 %45, 1
store i32 %46, i32* %6
br label %24
47: ; preds = %24
br label %48
48: ; preds = %47
%49 = load i32, i32* %5
%50 = add nsw i32 %49, 1
store i32 %50, i32* %5
br label %19
51: ; preds = %19
%52 = load double*, double** %9
%53 = load i32, i32* %4
%54 = load i32, i32* %4
%55 = call i32 @linpackDTRDI(double* %52, i32 %53, i32 %54)
store i32 %55, i32* %8
%56 = load i32, i32* %4
%57 = load i32, i32* %4
%58 = call double** @matrixNew(i32 %56, i32 %57)
store double** %58, double*** %10
store i32 0, i32* %7
store i32 0, i32* %5
br label %59
59: ; preds = %95, %51
%60 = load i32, i32* %5
%61 = load i32, i32* %4
%62 = icmp slt i32 %60, %61
br i1 %62, label %63, label %98
63: ; preds = %59
store i32 0, i32* %6
br label %64
64: ; preds = %91, %63
%65 = load i32, i32* %6
%66 = load i32, i32* %4
%67 = icmp slt i32 %65, %66
br i1 %67, label %68, label %94
68: ; preds = %64
%69 = load i32, i32* %5
%70 = load i32, i32* %6
%71 = icmp sgt i32 %69, %70
br i1 %71, label %72, label %73
72: ; preds = %68
br label %79
73: ; preds = %68
%74 = load double*, double** %9
%75 = load i32, i32* %7
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
%78 = load double, double* %77
br label %79
79: ; preds = %73, %72
%80 = phi double [ 0.000000e+00, %72 ], [ %78, %73 ]
%81 = load double**, double*** %10
%82 = load i32, i32* %6
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds double*, double** %81, i64 %83
%85 = load double*, double** %84
%86 = load i32, i32* %5
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %85, i64 %87
store double %80, double* %88
%89 = load i32, i32* %7
%90 = add nsw i32 %89, 1
store i32 %90, i32* %7
br label %91
91: ; preds = %79
%92 = load i32, i32* %6
%93 = add nsw i32 %92, 1
store i32 %93, i32* %6
br label %64
94: ; preds = %64
br label %95
95: ; preds = %94
%96 = load i32, i32* %5
%97 = add nsw i32 %96, 1
store i32 %97, i32* %5
br label %59
98: ; preds = %59
%99 = load double*, double** %9
%100 = call i32 @free(double* %99)
%101 = load double**, double*** %10
ret double** %101
}

LLVM IR (-O2):
define double** @matrixLowerTriangularInverse(double** nocapture %0, i32 %1) {
%3 = shl i32 %1, 3
%4 = mul i32 %3, %1
%5 = tail call i64 @Malloc(i32 %4)
%6 = inttoptr i64 %5 to double*
%7 = icmp sgt i32 %1, 0
br i1 %7, label %8, label %75
8: ; preds = %2
%9 = zext i32 %1 to i64
%10 = add nsw i64 %9, -1
%11 = and i64 %9, 3
%12 = icmp ult i64 %10, 3
%13 = and i64 %9, 4294967292
%14 = icmp eq i64 %11, 0
br label %15
15: ; preds = %71, %8
%16 = phi i64 [ 0, %8 ], [ %73, %71 ]
%17 = phi i64 [ 0, %8 ], [ %72, %71 ]
%18 = shl i64 %17, 32
%19 = ashr exact i64 %18, 32
br i1 %12, label %54, label %20
20: ; preds = %15, %20
%21 = phi i64 [ %50, %20 ], [ %19, %15 ]
%22 = phi i64 [ %51, %20 ], [ 0, %15 ]
%23 = phi i64 [ %52, %20 ], [ 0, %15 ]
%24 = getelementptr inbounds double*, double** %0, i64 %22
%25 = load double*, double** %24
%26 = getelementptr inbounds double, double* %25, i64 %16
%27 = load double, double* %26
%28 = getelementptr inbounds double, double* %6, i64 %21
store double %27, double* %28
%29 = add nsw i64 %21, 1
%30 = or i64 %22, 1
%31 = getelementptr inbounds double*, double** %0, i64 %30
%32 = load double*, double** %31
%33 = getelementptr inbounds double, double* %32, i64 %16
%34 = load double, double* %33
%35 = getelementptr inbounds double, double* %6, i64 %29
store double %34, double* %35
%36 = add nsw i64 %21, 2
%37 = or i64 %22, 2
%38 = getelementptr inbounds double*, double** %0, i64 %37
%39 = load double*, double** %38
%40 = getelementptr inbounds double, double* %39, i64 %16
%41 = load double, double* %40
%42 = getelementptr inbounds double, double* %6, i64 %36
store double %41, double* %42
%43 = add nsw i64 %21, 3
%44 = or i64 %22, 3
%45 = getelementptr inbounds double*, double** %0, i64 %44
%46 = load double*, double** %45
%47 = getelementptr inbounds double, double* %46, i64 %16
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %6, i64 %43
store double %48, double* %49
%50 = add nsw i64 %21, 4
%51 = add nuw nsw i64 %22, 4
%52 = add i64 %23, 4
%53 = icmp eq i64 %52, %13
br i1 %53, label %54, label %20
54: ; preds = %20, %15
%55 = phi i64 [ undef, %15 ], [ %50, %20 ]
%56 = phi i64 [ %19, %15 ], [ %50, %20 ]
%57 = phi i64 [ 0, %15 ], [ %51, %20 ]
br i1 %14, label %71, label %58
58: ; preds = %54, %58
%59 = phi i64 [ %67, %58 ], [ %56, %54 ]
%60 = phi i64 [ %68, %58 ], [ %57, %54 ]
%61 = phi i64 [ %69, %58 ], [ 0, %54 ]
%62 = getelementptr inbounds double*, double** %0, i64 %60
%63 = load double*, double** %62
%64 = getelementptr inbounds double, double* %63, i64 %16
%65 = load double, double* %64
%66 = getelementptr inbounds double, double* %6, i64 %59
store double %65, double* %66
%67 = add nsw i64 %59, 1
%68 = add nuw nsw i64 %60, 1
%69 = add i64 %61, 1
%70 = icmp eq i64 %69, %11
br i1 %70, label %71, label %58
71: ; preds = %58, %54
%72 = phi i64 [ %55, %54 ], [ %67, %58 ]
%73 = add nuw nsw i64 %16, 1
%74 = icmp eq i64 %73, %9
br i1 %74, label %75, label %15
75: ; preds = %71, %2
%76 = tail call i32 @linpackDTRDI(double* %6, i32 %1, i32 %1)
%77 = tail call double** @matrixNew(i32 %1, i32 %1)
br i1 %7, label %78, label %136
78: ; preds = %75
%79 = zext i32 %1 to i64
%80 = and i64 %79, 1
%81 = icmp eq i32 %1, 1
%82 = and i64 %79, 4294967294
%83 = icmp eq i64 %80, 0
br label %84
84: ; preds = %132, %78
%85 = phi i64 [ 0, %78 ], [ %134, %132 ]
%86 = phi i64 [ 0, %78 ], [ %133, %132 ]
%87 = shl i64 %86, 32
%88 = ashr exact i64 %87, 32
br i1 %81, label %117, label %89
89: ; preds = %84, %108
%90 = phi i64 [ %113, %108 ], [ %88, %84 ]
%91 = phi i64 [ %114, %108 ], [ 0, %84 ]
%92 = phi i64 [ %115, %108 ], [ 0, %84 ]
%93 = icmp ugt i64 %85, %91
br i1 %93, label %97, label %94
94: ; preds = %89
%95 = getelementptr inbounds double, double* %6, i64 %90
%96 = load double, double* %95
br label %97
97: ; preds = %89, %94
%98 = phi double [ %96, %94 ], [ 0.000000e+00, %89 ]
%99 = getelementptr inbounds double*, double** %77, i64 %91
%100 = load double*, double** %99
%101 = getelementptr inbounds double, double* %100, i64 %85
store double %98, double* %101
%102 = or i64 %91, 1
%103 = icmp ugt i64 %85, %102
br i1 %103, label %108, label %104
104: ; preds = %97
%105 = add nsw i64 %90, 1
%106 = getelementptr inbounds double, double* %6, i64 %105
%107 = load double, double* %106
br label %108
108: ; preds = %104, %97
%109 = phi double [ %107, %104 ], [ 0.000000e+00, %97 ]
%110 = getelementptr inbounds double*, double** %77, i64 %102
%111 = load double*, double** %110
%112 = getelementptr inbounds double, double* %111, i64 %85
store double %109, double* %112
%113 = add nsw i64 %90, 2
%114 = add nuw nsw i64 %91, 2
%115 = add i64 %92, 2
%116 = icmp eq i64 %115, %82
br i1 %116, label %117, label %89
117: ; preds = %108, %84
%118 = phi i64 [ undef, %84 ], [ %113, %108 ]
%119 = phi i64 [ %88, %84 ], [ %113, %108 ]
%120 = phi i64 [ 0, %84 ], [ %114, %108 ]
br i1 %83, label %132, label %121
121: ; preds = %117
%122 = icmp ugt i64 %85, %120
br i1 %122, label %126, label %123
123: ; preds = %121
%124 = getelementptr inbounds double, double* %6, i64 %119
%125 = load double, double* %124
br label %126
126: ; preds = %123, %121
%127 = phi double [ %125, %123 ], [ 0.000000e+00, %121 ]
%128 = getelementptr inbounds double*, double** %77, i64 %120
%129 = load double*, double** %128
%130 = getelementptr inbounds double, double* %129, i64 %85
store double %127, double* %130
%131 = add nsw i64 %119, 1
br label %132
132: ; preds = %117, %126
%133 = phi i64 [ %118, %117 ], [ %131, %126 ]
%134 = add nuw nsw i64 %85, 1
%135 = icmp eq i64 %134, %79
br i1 %135, label %136, label %84
136: ; preds = %132, %75
%137 = tail call i32 @free(double* %6)
ret double** %77
}