#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 /*<<< orphan*/ dpotrf_(char *, int *, double *, int *, int *);
int dpotrf_wrap(int n, const double *a, double *l) {
  char uplo = 'L';
  int i;
  for (i = 0; i < n * n; i++) {
    l[i] = a[i];
  }
  int info;
  dpotrf_(&uplo, &n, l, &n, &info);
  for (i = 0; i < n; i++) {
    int j;
    for (j = i + 1; j < n; j++) {
      l[j * n + i] = l[i * n + j];
      l[i * n + j] = 0.0;
    }
  }
  return (info);
}


x86 assembly (clang, -O0):
<dpotrf_wrap>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movb   $0x4c,-0x19(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x4(%rbp),%ecx
cmp    %ecx,%eax
jge    1168 <dpotrf_wrap+0x58>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112e <dpotrf_wrap+0x1e>
mov    -0x18(%rbp),%rdx
lea    -0x19(%rbp),%rdi
lea    -0x4(%rbp),%rcx
lea    -0x24(%rbp),%r8
mov    %rcx,%rsi
call   1030 <dpotrf_@plt>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    120e <dpotrf_wrap+0xfe>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11fb <dpotrf_wrap+0xeb>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    119c <dpotrf_wrap+0x8c>
jmp    1200 <dpotrf_wrap+0xf0>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1187 <dpotrf_wrap+0x77>
mov    -0x24(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dpotrf_wrap>:
push   %rbx
sub    $0x10,%rsp
mov    %rdx,%rbx
mov    %edi,0x8(%rsp)
movb   $0x4c,0x7(%rsp)
imul   %edi,%edi
test   %edi,%edi
je     1294 <dpotrf_wrap+0x184>
mov    %edi,%r9d
cmp    $0x4,%edi
jb     114a <dpotrf_wrap+0x3a>
lea    (%rsi,%r9,8),%rax
cmp    %rbx,%rax
jbe    11bc <dpotrf_wrap+0xac>
lea    (%rbx,%r9,8),%rax
cmp    %rsi,%rax
jbe    11bc <dpotrf_wrap+0xac>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %r9,%rdx
mov    %r9,%rdi
and    $0x3,%rdi
je     1174 <dpotrf_wrap+0x64>
xchg   %ax,%ax
movsd  (%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rdi
jne    1160 <dpotrf_wrap+0x50>
cmp    $0x3,%rdx
jb     1294 <dpotrf_wrap+0x184>
xchg   %ax,%ax
movsd  (%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
movsd  0x8(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rcx,8)
movsd  0x10(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rbx,%rcx,8)
movsd  0x18(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rbx,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%r9
jne    1180 <dpotrf_wrap+0x70>
jmp    1294 <dpotrf_wrap+0x184>
mov    %r9d,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rax
mov    %rax,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%r8d
and    $0x3,%r8d
cmp    $0xc,%rax
jae    11e2 <dpotrf_wrap+0xd2>
xor    %edx,%edx
jmp    1248 <dpotrf_wrap+0x138>
and    $0xfffffffffffffffc,%rdi
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movups (%rsi,%rdx,8),%xmm0
movups 0x10(%rsi,%rdx,8),%xmm1
movups %xmm0,(%rbx,%rdx,8)
movups %xmm1,0x10(%rbx,%rdx,8)
movups 0x20(%rsi,%rdx,8),%xmm0
movups 0x30(%rsi,%rdx,8),%xmm1
movups %xmm0,0x20(%rbx,%rdx,8)
movups %xmm1,0x30(%rbx,%rdx,8)
movups 0x40(%rsi,%rdx,8),%xmm0
movups 0x50(%rsi,%rdx,8),%xmm1
movups %xmm0,0x40(%rbx,%rdx,8)
movups %xmm1,0x50(%rbx,%rdx,8)
movups 0x60(%rsi,%rdx,8),%xmm0
movups 0x70(%rsi,%rdx,8),%xmm1
movups %xmm0,0x60(%rbx,%rdx,8)
movups %xmm1,0x70(%rbx,%rdx,8)
add    $0x10,%rdx
add    $0xfffffffffffffffc,%rdi
jne    11f0 <dpotrf_wrap+0xe0>
test   %r8,%r8
je     128b <dpotrf_wrap+0x17b>
lea    (%rbx,%rdx,8),%rdi
add    $0x10,%rdi
lea    (%rsi,%rdx,8),%rdx
add    $0x10,%rdx
shl    $0x5,%r8
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movups -0x10(%rdx,%rax,1),%xmm0
movups (%rdx,%rax,1),%xmm1
movups %xmm0,-0x10(%rdi,%rax,1)
movups %xmm1,(%rdi,%rax,1)
add    $0x20,%rax
cmp    %rax,%r8
jne    1270 <dpotrf_wrap+0x160>
cmp    %r9,%rcx
jne    114c <dpotrf_wrap+0x3c>
lea    0x7(%rsp),%rdi
lea    0x8(%rsp),%rcx
lea    0xc(%rsp),%r8
mov    %rcx,%rsi
mov    %rbx,%rdx
call   1030 <dpotrf_@plt>
mov    0x8(%rsp),%r8d
test   %r8d,%r8d
jle    131f <dpotrf_wrap+0x20f>
lea    (%rbx,%r8,8),%rdi
lea    0x0(,%r8,8),%rdx
lea    0x8(,%r8,8),%r9
xor    %r10d,%r10d
mov    %r8,%rax
jmp    12ef <dpotrf_wrap+0x1df>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    %r9,%rdi
add    %r9,%rbx
add    $0xffffffffffffffff,%rax
cmp    %r8,%r10
je     131f <dpotrf_wrap+0x20f>
add    $0x1,%r10
cmp    %r8,%r10
jae    12e0 <dpotrf_wrap+0x1d0>
mov    $0x1,%esi
mov    %rdi,%rcx
movsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,(%rcx)
movq   $0x0,(%rbx,%rsi,8)
add    %rdx,%rcx
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1300 <dpotrf_wrap+0x1f0>
jmp    12e0 <dpotrf_wrap+0x1d0>
mov    0xc(%rsp),%eax
add    $0x10,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define i32 @dpotrf_wrap(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca i8
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
store i8 76, i8* %7
store i32 0, i32* %8
br label %11
11: ; preds = %27, %3
%12 = load i32, i32* %8
%13 = load i32, i32* %4
%14 = load i32, i32* %4
%15 = mul nsw i32 %13, %14
%16 = icmp slt i32 %12, %15
br i1 %16, label %17, label %30
17: ; preds = %11
%18 = load double*, double** %5
%19 = load i32, i32* %8
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double, double* %18, i64 %20
%22 = load double, double* %21
%23 = load double*, double** %6
%24 = load i32, i32* %8
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double, double* %23, i64 %25
store double %22, double* %26
br label %27
27: ; preds = %17
%28 = load i32, i32* %8
%29 = add nsw i32 %28, 1
store i32 %29, i32* %8
br label %11
30: ; preds = %11
%31 = load double*, double** %6
%32 = call i32 @dpotrf_(i8* %7, i32* %4, double* %31, i32* %4, i32* %9)
store i32 0, i32* %8
br label %33
33: ; preds = %74, %30
%34 = load i32, i32* %8
%35 = load i32, i32* %4
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %77
37: ; preds = %33
%38 = load i32, i32* %8
%39 = add nsw i32 %38, 1
store i32 %39, i32* %10
br label %40
40: ; preds = %70, %37
%41 = load i32, i32* %10
%42 = load i32, i32* %4
%43 = icmp slt i32 %41, %42
br i1 %43, label %44, label %73
44: ; preds = %40
%45 = load double*, double** %6
%46 = load i32, i32* %8
%47 = load i32, i32* %4
%48 = mul nsw i32 %46, %47
%49 = load i32, i32* %10
%50 = add nsw i32 %48, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %45, i64 %51
%53 = load double, double* %52
%54 = load double*, double** %6
%55 = load i32, i32* %10
%56 = load i32, i32* %4
%57 = mul nsw i32 %55, %56
%58 = load i32, i32* %8
%59 = add nsw i32 %57, %58
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %54, i64 %60
store double %53, double* %61
%62 = load double*, double** %6
%63 = load i32, i32* %8
%64 = load i32, i32* %4
%65 = mul nsw i32 %63, %64
%66 = load i32, i32* %10
%67 = add nsw i32 %65, %66
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %62, i64 %68
store double 0.000000e+00, double* %69
br label %70
70: ; preds = %44
%71 = load i32, i32* %10
%72 = add nsw i32 %71, 1
store i32 %72, i32* %10
br label %40
73: ; preds = %40
br label %74
74: ; preds = %73
%75 = load i32, i32* %8
%76 = add nsw i32 %75, 1
store i32 %76, i32* %8
br label %33
77: ; preds = %33
%78 = load i32, i32* %9
ret i32 %78
}

LLVM IR (-O2):
define i32 @dpotrf_wrap(i32 %0, double* nocapture %1, double* %2) {
%4 = alloca i32
%5 = alloca i8
%6 = alloca i32
store i32 %0, i32* %4
store i8 76, i8* %5
%7 = mul i32 %0, %0
%8 = icmp eq i32 %7, 0
br i1 %8, label %134, label %9
9: ; preds = %3
%10 = zext i32 %7 to i64
%11 = icmp ult i32 %7, 4
br i1 %11, label %97, label %12
12: ; preds = %9
%13 = getelementptr double, double* %2, i64 %10
%14 = getelementptr double, double* %1, i64 %10
%15 = icmp ugt double* %14, %2
%16 = icmp ugt double* %13, %1
%17 = and i1 %15, %16
br i1 %17, label %97, label %18
18: ; preds = %12
%19 = and i64 %10, 4294967292
%20 = add nsw i64 %19, -4
%21 = lshr exact i64 %20, 2
%22 = add nuw nsw i64 %21, 1
%23 = and i64 %22, 3
%24 = icmp ult i64 %20, 12
br i1 %24, label %76, label %25
25: ; preds = %18
%26 = and i64 %22, 9223372036854775804
br label %27
27: ; preds = %27, %25
%28 = phi i64 [ 0, %25 ], [ %73, %27 ]
%29 = phi i64 [ 0, %25 ], [ %74, %27 ]
%30 = getelementptr inbounds double, double* %1, i64 %28
%31 = bitcast double* %30 to <2 x double>*
%32 = load <2 x double>, <2 x double>* %31
%33 = getelementptr inbounds double, double* %30, i64 2
%34 = bitcast double* %33 to <2 x double>*
%35 = load <2 x double>, <2 x double>* %34
%36 = getelementptr inbounds double, double* %2, i64 %28
%37 = bitcast double* %36 to <2 x double>*
store <2 x double> %32, <2 x double>* %37
%38 = getelementptr inbounds double, double* %36, i64 2
%39 = bitcast double* %38 to <2 x double>*
store <2 x double> %35, <2 x double>* %39
%40 = or i64 %28, 4
%41 = getelementptr inbounds double, double* %1, i64 %40
%42 = bitcast double* %41 to <2 x double>*
%43 = load <2 x double>, <2 x double>* %42
%44 = getelementptr inbounds double, double* %41, i64 2
%45 = bitcast double* %44 to <2 x double>*
%46 = load <2 x double>, <2 x double>* %45
%47 = getelementptr inbounds double, double* %2, i64 %40
%48 = bitcast double* %47 to <2 x double>*
store <2 x double> %43, <2 x double>* %48
%49 = getelementptr inbounds double, double* %47, i64 2
%50 = bitcast double* %49 to <2 x double>*
store <2 x double> %46, <2 x double>* %50
%51 = or i64 %28, 8
%52 = getelementptr inbounds double, double* %1, i64 %51
%53 = bitcast double* %52 to <2 x double>*
%54 = load <2 x double>, <2 x double>* %53
%55 = getelementptr inbounds double, double* %52, i64 2
%56 = bitcast double* %55 to <2 x double>*
%57 = load <2 x double>, <2 x double>* %56
%58 = getelementptr inbounds double, double* %2, i64 %51
%59 = bitcast double* %58 to <2 x double>*
store <2 x double> %54, <2 x double>* %59
%60 = getelementptr inbounds double, double* %58, i64 2
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> %57, <2 x double>* %61
%62 = or i64 %28, 12
%63 = getelementptr inbounds double, double* %1, i64 %62
%64 = bitcast double* %63 to <2 x double>*
%65 = load <2 x double>, <2 x double>* %64
%66 = getelementptr inbounds double, double* %63, i64 2
%67 = bitcast double* %66 to <2 x double>*
%68 = load <2 x double>, <2 x double>* %67
%69 = getelementptr inbounds double, double* %2, i64 %62
%70 = bitcast double* %69 to <2 x double>*
store <2 x double> %65, <2 x double>* %70
%71 = getelementptr inbounds double, double* %69, i64 2
%72 = bitcast double* %71 to <2 x double>*
store <2 x double> %68, <2 x double>* %72
%73 = add nuw i64 %28, 16
%74 = add nuw i64 %29, 4
%75 = icmp eq i64 %74, %26
br i1 %75, label %76, label %27
76: ; preds = %27, %18
%77 = phi i64 [ 0, %18 ], [ %73, %27 ]
%78 = icmp eq i64 %23, 0
br i1 %78, label %95, label %79
79: ; preds = %76, %79
%80 = phi i64 [ %92, %79 ], [ %77, %76 ]
%81 = phi i64 [ %93, %79 ], [ 0, %76 ]
%82 = getelementptr inbounds double, double* %1, i64 %80
%83 = bitcast double* %82 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = getelementptr inbounds double, double* %82, i64 2
%86 = bitcast double* %85 to <2 x double>*
%87 = load <2 x double>, <2 x double>* %86
%88 = getelementptr inbounds double, double* %2, i64 %80
%89 = bitcast double* %88 to <2 x double>*
store <2 x double> %84, <2 x double>* %89
%90 = getelementptr inbounds double, double* %88, i64 2
%91 = bitcast double* %90 to <2 x double>*
store <2 x double> %87, <2 x double>* %91
%92 = add nuw i64 %80, 4
%93 = add i64 %81, 1
%94 = icmp eq i64 %93, %23
br i1 %94, label %95, label %79
95: ; preds = %79, %76
%96 = icmp eq i64 %19, %10
br i1 %96, label %134, label %97
97: ; preds = %12, %9, %95
%98 = phi i64 [ 0, %12 ], [ 0, %9 ], [ %19, %95 ]
%99 = xor i64 %98, -1
%100 = add nsw i64 %99, %10
%101 = and i64 %10, 3
%102 = icmp eq i64 %101, 0
br i1 %102, label %112, label %103
103: ; preds = %97, %103
%104 = phi i64 [ %109, %103 ], [ %98, %97 ]
%105 = phi i64 [ %110, %103 ], [ 0, %97 ]
%106 = getelementptr inbounds double, double* %1, i64 %104
%107 = load double, double* %106
%108 = getelementptr inbounds double, double* %2, i64 %104
store double %107, double* %108
%109 = add nuw nsw i64 %104, 1
%110 = add i64 %105, 1
%111 = icmp eq i64 %110, %101
br i1 %111, label %112, label %103
112: ; preds = %103, %97
%113 = phi i64 [ %98, %97 ], [ %109, %103 ]
%114 = icmp ult i64 %100, 3
br i1 %114, label %134, label %115
115: ; preds = %112, %115
%116 = phi i64 [ %132, %115 ], [ %113, %112 ]
%117 = getelementptr inbounds double, double* %1, i64 %116
%118 = load double, double* %117
%119 = getelementptr inbounds double, double* %2, i64 %116
store double %118, double* %119
%120 = add nuw nsw i64 %116, 1
%121 = getelementptr inbounds double, double* %1, i64 %120
%122 = load double, double* %121
%123 = getelementptr inbounds double, double* %2, i64 %120
store double %122, double* %123
%124 = add nuw nsw i64 %116, 2
%125 = getelementptr inbounds double, double* %1, i64 %124
%126 = load double, double* %125
%127 = getelementptr inbounds double, double* %2, i64 %124
store double %126, double* %127
%128 = add nuw nsw i64 %116, 3
%129 = getelementptr inbounds double, double* %1, i64 %128
%130 = load double, double* %129
%131 = getelementptr inbounds double, double* %2, i64 %128
store double %130, double* %131
%132 = add nuw nsw i64 %116, 4
%133 = icmp eq i64 %132, %10
br i1 %133, label %134, label %115
134: ; preds = %112, %115, %95, %3
%135 = bitcast i32* %6 to i8*
%136 = call i32 @dpotrf_(i8* %5, i32* %4, double* %2, i32* %4, i32* %6)
%137 = load i32, i32* %4
%138 = icmp sgt i32 %137, 0
br i1 %138, label %139, label %162
139: ; preds = %134
%140 = zext i32 %137 to i64
%141 = zext i32 %137 to i64
br label %145
142: ; preds = %152, %145
%143 = add nuw nsw i64 %147, 1
%144 = icmp eq i64 %148, %141
br i1 %144, label %162, label %145
145: ; preds = %139, %142
%146 = phi i64 [ 0, %139 ], [ %148, %142 ]
%147 = phi i64 [ 1, %139 ], [ %143, %142 ]
%148 = add nuw nsw i64 %146, 1
%149 = icmp ult i64 %148, %140
br i1 %149, label %150, label %142
150: ; preds = %145
%151 = mul nsw i64 %146, %140
br label %152
152: ; preds = %150, %152
%153 = phi i64 [ %147, %150 ], [ %160, %152 ]
%154 = add nuw nsw i64 %153, %151
%155 = getelementptr inbounds double, double* %2, i64 %154
%156 = load double, double* %155
%157 = mul nsw i64 %153, %140
%158 = add nuw nsw i64 %157, %146
%159 = getelementptr inbounds double, double* %2, i64 %158
store double %156, double* %159
store double 0.000000e+00, double* %155
%160 = add nuw nsw i64 %153, 1
%161 = icmp eq i64 %160, %141
br i1 %161, label %142, label %152
162: ; preds = %142, %134
%163 = load i32, i32* %6
ret i32 %163
}