#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 */
double sqrt(double);
void OneFac(double *A, long n1, long n2) {
  long i, j, k;
  for (j = 0; j < n1; j++) {
    for (k = 0; k < j; k++)
      for (i = j; i < n1; i++)
        A[i + n2 * j] -= A[j + n2 * k] * A[i + n2 * k];
    A[k + n2 * k] = sqrt(A[k + n2 * k]);
    for (i = j + 1; i < n1; i++)
      A[i + n2 * k] /= A[k + n2 * k];
  }
}


x86 assembly (clang, -O0):
<OneFac>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    12af <OneFac+0x19f>
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    11fc <OneFac+0xec>
mov    -0x28(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11e6 <OneFac+0xd6>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    -0x18(%rbp),%rdx
imul   -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x18(%rbp),%rdx
imul   -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x18(%rbp),%rdx
imul   -0x28(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1158 <OneFac+0x48>
jmp    11eb <OneFac+0xdb>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1142 <OneFac+0x32>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x18(%rbp),%rdx
imul   -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <sqrt@plt>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x18(%rbp),%rdx
imul   -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1299 <OneFac+0x189>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x18(%rbp),%rdx
imul   -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x18(%rbp),%rdx
imul   -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    123f <OneFac+0x12f>
jmp    129e <OneFac+0x18e>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    112c <OneFac+0x1c>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<OneFac>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
test   %rsi,%rsi
jle    1439 <OneFac+0x329>
mov    %rsi,%r15
mov    %rdi,%r13
lea    0x10(%rdi),%rbx
lea    0x0(,%rdx,8),%rax
mov    %rax,0x40(%rsp)
lea    0x8(,%rdx,8),%rax
mov    %rax,0x48(%rsp)
mov    $0x1,%edi
xor    %eax,%eax
mov    %rax,0x28(%rsp)
xorpd  %xmm5,%xmm5
movapd 0xe9d(%rip),%xmm6
mov    %r13,%r12
mov    %rbx,0x20(%rsp)
xor    %r8d,%r8d
mov    %rdx,0x38(%rsp)
jmp    11a3 <OneFac+0x93>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rdi
addq   $0xffffffffffffffff,0x28(%rsp)
addq   $0x8,0x20(%rsp)
add    0x48(%rsp),%rbx
add    0x40(%rsp),%r12
cmp    %r15,%r8
je     1439 <OneFac+0x329>
mov    %r8,%r14
cmp    %rdi,%r15
mov    %rdi,0x30(%rsp)
mov    %rdi,%rax
cmovg  %r15,%rax
lea    0x1(%r8),%rcx
cmp    %rcx,%r15
mov    %rcx,0x10(%rsp)
cmovg  %r15,%rcx
mov    %rcx,0x8(%rsp)
test   %r8,%r8
je     1380 <OneFac+0x270>
add    0x28(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    0x8(%rsp),%r8
mov    %r8,%rdi
sub    %r14,%rdi
mov    %r14,%rbp
shl    $0x4,%rbp
add    %r13,%rbp
mov    %r14,%rcx
not    %rcx
add    %r8,%rcx
mov    %rbp,0x70(%rsp)
lea    0x8(,%rcx,8),%rcx
add    %rbp,%rcx
mov    %rcx,0x68(%rsp)
cmp    $0x1,%rdx
sete   %cl
cmp    $0x4,%rdi
setae  %r9b
and    %cl,%r9b
mov    %rdi,0x60(%rsp)
and    $0xfffffffffffffffc,%rdi
mov    %rdi,0x58(%rsp)
lea    (%r14,%rdi,1),%rcx
mov    %rcx,0x50(%rsp)
mov    %r13,%rbp
mov    0x20(%rsp),%rdi
xor    %r8d,%r8d
mov    %r14,0x18(%rsp)
jmp    126d <OneFac+0x15d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r8
mov    0x40(%rsp),%rcx
add    %rcx,%rdi
add    %rcx,%rbp
mov    0x18(%rsp),%r14
cmp    %r14,%r8
je     1380 <OneFac+0x270>
mov    %r8,%r10
imul   %rdx,%r10
add    %r14,%r10
test   %r9b,%r9b
je     1350 <OneFac+0x240>
mov    0x18(%rsp),%r14
lea    (%r14,%r8,1),%rcx
lea    0x0(,%rcx,8),%rcx
add    %r13,%rcx
mov    0x8(%rsp),%rdx
add    %r8,%rdx
lea    0x0(%r13,%rdx,8),%r11
mov    0x10(%rsp),%rdx
lea    (%rdx,%r8,1),%rsi
lea    0x0(,%rsi,8),%rsi
add    %r13,%rsi
mov    %r9d,%edx
mov    0x70(%rsp),%r9
cmp    %rsi,%r9
setb   %sil
cmp    0x68(%rsp),%rcx
setb   %cl
cmp    %r11,%r9
mov    %edx,%r9d
setb   %r11b
test   %cl,%sil
mov    0x38(%rsp),%rdx
jne    1350 <OneFac+0x240>
mov    0x18(%rsp),%r14
and    %cl,%r11b
jne    1350 <OneFac+0x240>
movsd  0x0(%r13,%r10,8),%xmm0
unpcklpd %xmm0,%xmm0
xorpd  %xmm6,%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movupd -0x10(%rdi,%rcx,8),%xmm1
movupd (%rdi,%rcx,8),%xmm2
movupd -0x10(%rbx,%rcx,8),%xmm3
movupd (%rbx,%rcx,8),%xmm4
mulpd  %xmm0,%xmm1
addpd  %xmm3,%xmm1
mulpd  %xmm0,%xmm2
addpd  %xmm4,%xmm2
movupd %xmm1,-0x10(%rbx,%rcx,8)
movupd %xmm2,(%rbx,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    1300 <OneFac+0x1f0>
mov    0x50(%rsp),%r14
mov    0x58(%rsp),%rcx
cmp    %rcx,0x60(%rsp)
je     1250 <OneFac+0x140>
nop
movsd  0x0(%r13,%r10,8),%xmm0
movsd  (%r12,%r14,8),%xmm1
mulsd  0x0(%rbp,%r14,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%r12,%r14,8)
add    $0x1,%r14
cmp    %r15,%r14
jl     1350 <OneFac+0x240>
jmp    1250 <OneFac+0x140>
nopl   0x0(%rax)
mov    %r14,%rax
imul   %rdx,%rax
lea    (%rax,%r14,1),%rbp
movsd  0x0(%r13,%rbp,8),%xmm0
ucomisd %xmm5,%xmm0
jb     13a0 <OneFac+0x290>
sqrtsd %xmm0,%xmm0
jmp    13c0 <OneFac+0x2b0>
xchg   %ax,%ax
mov    %rax,0x8(%rsp)
call   1030 <sqrt@plt>
mov    0x8(%rsp),%rax
movapd 0xc49(%rip),%xmm6
xorpd  %xmm5,%xmm5
mov    0x38(%rsp),%rdx
mov    0x30(%rsp),%rdi
mov    0x10(%rsp),%r8
movsd  %xmm0,0x0(%r13,%rbp,8)
cmp    %r15,%r8
jge    1180 <OneFac+0x70>
add    %r8,%rax
movsd  0x0(%r13,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rax,8)
lea    0x2(%r14),%rax
cmp    %r15,%rax
jge    1180 <OneFac+0x70>
mov    $0x2,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x10(%rbx,%rax,8),%xmm0
divsd  0x0(%r13,%rbp,8),%xmm0
movsd  %xmm0,-0x10(%rbx,%rax,8)
lea    (%r14,%rax,1),%rcx
add    $0x1,%rcx
add    $0x1,%rax
cmp    %r15,%rcx
jl     1410 <OneFac+0x300>
jmp    1180 <OneFac+0x70>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @OneFac(double* nocapture %0, i64 %1, i64 %2) {
%4 = icmp sgt i64 %1, 0
br i1 %4, label %5, label %124
5: ; preds = %3
%6 = icmp eq i64 %2, 1
br label %9
7: ; preds = %115, %108, %100
%8 = icmp eq i64 %106, %1
br i1 %8, label %124, label %9
9: ; preds = %5, %7
%10 = phi i64 [ %106, %7 ], [ 0, %5 ]
%11 = add nuw i64 %10, 1
%12 = call i64 @llvm.smax.i64(i64 %1, i64 %11)
%13 = sub i64 %12, %10
%14 = shl nuw i64 %10, 1
%15 = getelementptr double, double* %0, i64 %14
%16 = or i64 %14, 1
%17 = getelementptr double, double* %0, i64 %16
%18 = add nuw i64 %10, 1
%19 = call i64 @llvm.smax.i64(i64 %1, i64 %18)
%20 = xor i64 %10, -1
%21 = add i64 %19, %20
%22 = getelementptr double, double* %17, i64 %21
%23 = icmp eq i64 %10, 0
br i1 %23, label %100, label %24
24: ; preds = %9
%25 = mul nsw i64 %10, %2
%26 = icmp ugt i64 %13, 3
%27 = and i1 %26, %6
%28 = and i64 %13, -4
%29 = add i64 %10, %28
%30 = icmp eq i64 %13, %28
br label %31
31: ; preds = %97, %24
%32 = phi i64 [ 0, %24 ], [ %98, %97 ]
%33 = add nuw i64 %10, %32
%34 = getelementptr double, double* %0, i64 %33
%35 = mul nsw i64 %32, %2
%36 = add nsw i64 %35, %10
%37 = getelementptr inbounds double, double* %0, i64 %36
br i1 %27, label %38, label %82
38: ; preds = %31
%39 = add i64 %19, %32
%40 = getelementptr double, double* %0, i64 %39
%41 = add i64 %18, %32
%42 = getelementptr double, double* %0, i64 %41
%43 = icmp ult double* %15, %42
%44 = icmp ult double* %34, %22
%45 = and i1 %43, %44
%46 = icmp ult double* %15, %40
%47 = icmp ult double* %34, %22
%48 = and i1 %46, %47
%49 = or i1 %45, %48
br i1 %49, label %82, label %50
50: ; preds = %38
%51 = load double, double* %37
%52 = insertelement <2 x double> poison, double %51, i64 0
%53 = shufflevector <2 x double> %52, <2 x double> poison, <2 x i32> zeroinitializer
%54 = insertelement <2 x double> poison, double %51, i64 0
%55 = shufflevector <2 x double> %54, <2 x double> poison, <2 x i32> zeroinitializer
%56 = fneg <2 x double> %53
%57 = fneg <2 x double> %55
br label %58
58: ; preds = %58, %50
%59 = phi i64 [ 0, %50 ], [ %79, %58 ]
%60 = add i64 %10, %59
%61 = add nsw i64 %60, %35
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = getelementptr inbounds double, double* %62, i64 2
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = add nsw i64 %60, %25
%69 = getelementptr inbounds double, double* %0, i64 %68
%70 = bitcast double* %69 to <2 x double>*
%71 = load <2 x double>, <2 x double>* %70
%72 = getelementptr inbounds double, double* %69, i64 2
%73 = bitcast double* %72 to <2 x double>*
%74 = load <2 x double>, <2 x double>* %73
%75 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %56, <2 x double> %64, <2 x double> %71)
%76 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %57, <2 x double> %67, <2 x double> %74)
%77 = bitcast double* %69 to <2 x double>*
store <2 x double> %75, <2 x double>* %77
%78 = bitcast double* %72 to <2 x double>*
store <2 x double> %76, <2 x double>* %78
%79 = add nuw i64 %59, 4
%80 = icmp eq i64 %79, %28
br i1 %80, label %81, label %58
81: ; preds = %58
br i1 %30, label %97, label %82
82: ; preds = %38, %31, %81
%83 = phi i64 [ %10, %38 ], [ %10, %31 ], [ %29, %81 ]
br label %84
84: ; preds = %82, %84
%85 = phi i64 [ %95, %84 ], [ %83, %82 ]
%86 = load double, double* %37
%87 = add nsw i64 %85, %35
%88 = getelementptr inbounds double, double* %0, i64 %87
%89 = load double, double* %88
%90 = add nsw i64 %85, %25
%91 = getelementptr inbounds double, double* %0, i64 %90
%92 = load double, double* %91
%93 = fneg double %86
%94 = tail call double @llvm.fmuladd.f64(double %93, double %89, double %92)
store double %94, double* %91
%95 = add nuw nsw i64 %85, 1
%96 = icmp slt i64 %95, %1
br i1 %96, label %84, label %97
97: ; preds = %84, %81
%98 = add nuw nsw i64 %32, 1
%99 = icmp eq i64 %98, %10
br i1 %99, label %100, label %31
100: ; preds = %97, %9
%101 = mul nsw i64 %10, %2
%102 = add nsw i64 %101, %10
%103 = getelementptr inbounds double, double* %0, i64 %102
%104 = load double, double* %103
%105 = tail call double @sqrt(double %104)
store double %105, double* %103
%106 = add nuw nsw i64 %10, 1
%107 = icmp slt i64 %106, %1
br i1 %107, label %108, label %7
108: ; preds = %100
%109 = add nsw i64 %106, %101
%110 = getelementptr inbounds double, double* %0, i64 %109
%111 = load double, double* %110
%112 = fdiv double %111, %105
store double %112, double* %110
%113 = add nuw nsw i64 %10, 2
%114 = icmp slt i64 %113, %1
br i1 %114, label %115, label %7
115: ; preds = %108, %115
%116 = phi i64 [ %122, %115 ], [ %113, %108 ]
%117 = load double, double* %103
%118 = add nsw i64 %116, %101
%119 = getelementptr inbounds double, double* %0, i64 %118
%120 = load double, double* %119
%121 = fdiv double %120, %117
store double %121, double* %119
%122 = add nuw nsw i64 %116, 1
%123 = icmp slt i64 %122, %1
br i1 %123, label %115, label %7
124: ; preds = %7, %3
ret void
}