#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*/ dger_3l(int, int, double, double *, int, double *, int,
                           double *, int);
int /*<<< orphan*/ dscal_3l(int, double, double *);
int /*<<< orphan*/ dswap_3l(int, double *, int, double *, int);
int idamax_3l(int, double *);
void dgetf2_3l(int m, int n, double *A, int lda, int *ipiv, int *info) {
  if (m <= 0 || n <= 0)
    return;
  int i, j, jp;
  double Ajj;
  int size_min = (m < n ? m : n);
  for (j = 0; j < size_min; j++) {
    jp = j + idamax_3l(m - j, &A[j + lda * j]);
    ipiv[j] = jp;
    if (A[jp + lda * j] != 0) {
      if (jp != j) {
        dswap_3l(n, &A[j], lda, &A[jp], lda);
      }
      if (j < m - 1) {
        Ajj = A[j + lda * j];
        if ((Ajj > 0 ? Ajj : -Ajj) >= 2.22e-16) {
          dscal_3l(m - j - 1, 1.0 / Ajj, &A[j + 1 + lda * j]);
        } else {
          for (i = j + 1; i < m; i++) {
            A[i + lda * j] /= Ajj;
          }
        }
      }
    } else if (*info == 0) {
      *info = j + 1;
    }
    if (j < size_min) {
      dger_3l(m - j - 1, n - j - 1, -1.0, &A[j + 1 + lda * j], 1,
              &A[j + lda * (j + 1)], lda, &A[j + 1 + lda * (j + 1)], lda);
    }
  }
  return;
}


x86 assembly (clang, -O0):
<dgetf2_3l>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
cmpl   $0x0,-0x4(%rbp)
jle    1171 <dgetf2_3l+0x31>
cmpl   $0x0,-0x8(%rbp)
jg     1176 <dgetf2_3l+0x36>
jmp    1448 <dgetf2_3l+0x308>
mov    -0x4(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    118d <dgetf2_3l+0x4d>
mov    -0x4(%rbp),%eax
mov    %eax,-0x48(%rbp)
jmp    1193 <dgetf2_3l+0x53>
mov    -0x8(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x44(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x44(%rbp),%eax
jge    1443 <dgetf2_3l+0x303>
mov    -0x30(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4(%rbp),%edi
sub    -0x30(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x30(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x30(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
call   1040 <idamax_3l@plt>
mov    %eax,%ecx
mov    -0x4c(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%edx
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x14(%rbp),%edx
imul   -0x30(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    121e <dgetf2_3l+0xde>
jp     121e <dgetf2_3l+0xde>
jmp    1372 <dgetf2_3l+0x232>
mov    -0x34(%rbp),%eax
cmp    -0x30(%rbp),%eax
je     1257 <dgetf2_3l+0x117>
mov    -0x8(%rbp),%edi
mov    -0x10(%rbp),%rsi
movslq -0x30(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x14(%rbp),%edx
mov    -0x10(%rbp),%rcx
movslq -0x34(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rcx
mov    -0x14(%rbp),%r8d
call   1030 <dswap_3l@plt>
mov    -0x30(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    136d <dgetf2_3l+0x22d>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x14(%rbp),%edx
imul   -0x30(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12a6 <dgetf2_3l+0x166>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    12c7 <dgetf2_3l+0x187>
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  0xd2c(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     131a <dgetf2_3l+0x1da>
mov    -0x4(%rbp),%edi
sub    -0x30(%rbp),%edi
sub    $0x1,%edi
movsd  0xd19(%rip),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rsi
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    -0x14(%rbp),%ecx
imul   -0x30(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
call   1060 <dscal_3l@plt>
jmp    1368 <dgetf2_3l+0x228>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1363 <dgetf2_3l+0x223>
movsd  -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x14(%rbp),%edx
imul   -0x30(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1323 <dgetf2_3l+0x1e3>
jmp    1368 <dgetf2_3l+0x228>
jmp    136d <dgetf2_3l+0x22d>
jmp    1390 <dgetf2_3l+0x250>
mov    -0x28(%rbp),%rax
cmpl   $0x0,(%rax)
jne    138b <dgetf2_3l+0x24b>
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
jmp    1390 <dgetf2_3l+0x250>
mov    -0x30(%rbp),%eax
cmp    -0x44(%rbp),%eax
jge    1430 <dgetf2_3l+0x2f0>
mov    -0x4(%rbp),%edi
sub    -0x30(%rbp),%edi
sub    $0x1,%edi
mov    -0x8(%rbp),%esi
sub    -0x30(%rbp),%esi
sub    $0x1,%esi
mov    -0x10(%rbp),%rdx
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    -0x14(%rbp),%ecx
imul   -0x30(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
mov    -0x10(%rbp),%r8
mov    -0x30(%rbp),%eax
mov    -0x14(%rbp),%ecx
mov    -0x30(%rbp),%r9d
add    $0x1,%r9d
imul   %r9d,%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%r8
mov    -0x14(%rbp),%r9d
mov    -0x10(%rbp),%r10
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    -0x14(%rbp),%ecx
mov    -0x30(%rbp),%r11d
add    $0x1,%r11d
imul   %r11d,%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%r10
mov    -0x14(%rbp),%eax
movsd  0xbf2(%rip),%xmm0
mov    $0x1,%ecx
mov    %r10,(%rsp)
mov    %eax,0x8(%rsp)
call   1050 <dger_3l@plt>
jmp    1435 <dgetf2_3l+0x2f5>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11a0 <dgetf2_3l+0x60>
jmp    1448 <dgetf2_3l+0x308>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dgetf2_3l>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r9,0x10(%rsp)
mov    %r8,0x60(%rsp)
mov    %esi,0x4(%rsp)
test   %edi,%edi
jle    144e <dgetf2_3l+0x30e>
cmpl   $0x0,0x4(%rsp)
jle    144e <dgetf2_3l+0x30e>
mov    %ecx,%ebx
mov    %rdx,%r12
mov    %rdi,%rcx
mov    0x4(%rsp),%eax
cmp    %eax,%ecx
cmovb  %ecx,%eax
mov    %rax,0x58(%rsp)
mov    %ecx,%ebp
add    $0xffffffff,%ecx
mov    %rcx,0x38(%rsp)
movslq %ebx,%rax
lea    0x18(%rdx),%rcx
lea    0x0(,%rax,8),%rdx
mov    %rdx,0x48(%rsp)
mov    %rax,0x50(%rsp)
lea    0x8(,%rax,8),%rax
mov    %rax,0x40(%rsp)
lea    -0x3(%rbp),%rax
mov    $0x1,%edx
mov    %rdx,0x8(%rsp)
xor    %r15d,%r15d
mov    %r12,%r14
mov    %rbx,0x18(%rsp)
jmp    128d <dgetf2_3l+0x14d>
movsd  0xe3d(%rip),%xmm0
divsd  %xmm1,%xmm0
lea    (%r15,%rbx,1),%rax
lea    (%r12,%rax,8),%rsi
add    $0x8,%rsi
mov    0x30(%rsp),%rax
lea    -0x1(%rax),%edi
call   1060 <dscal_3l@plt>
nopl   0x0(%rax,%rax,1)
mov    %r15d,%esi
not    %esi
add    0x4(%rsp),%esi
mov    %rbp,%r13
lea    0x1(%r15),%rbp
lea    (%r15,%rbx,1),%rax
add    $0x1,%rax
lea    (%r12,%rax,8),%rdx
mov    %ebp,%eax
mov    0x18(%rsp),%rbx
imul   %ebx,%eax
cltq
lea    (%r15,%rax,1),%rcx
lea    (%r12,%rcx,8),%r8
add    %r15,%rax
add    $0x1,%rax
lea    (%r12,%rax,8),%rax
mov    0x30(%rsp),%rdi
add    $0xffffffff,%edi
movsd  0xdd5(%rip),%xmm0
mov    $0x1,%ecx
mov    %ebx,%r9d
push   %rbx
push   %rax
call   1050 <dger_3l@plt>
add    $0x10,%rsp
addq   $0x1,0x8(%rsp)
mov    0x28(%rsp),%rcx
add    0x40(%rsp),%rcx
mov    0x20(%rsp),%rax
add    $0xffffffffffffffff,%rax
add    0x48(%rsp),%r14
mov    %rbp,%r15
cmp    0x58(%rsp),%rbp
mov    %r13,%rbp
je     144e <dgetf2_3l+0x30e>
mov    %rax,0x20(%rsp)
mov    %rcx,0x28(%rsp)
mov    %rbp,%rdi
sub    %r15,%rdi
mov    %r15,%rcx
imul   0x50(%rsp),%rcx
mov    %r15d,%eax
imul   %ebx,%eax
mov    %rcx,%rbx
cltq
add    %r15,%rax
lea    (%r12,%rax,8),%r13
mov    %rdi,0x30(%rsp)
mov    %r13,%rsi
call   1040 <idamax_3l@plt>
lea    (%rax,%r15,1),%ecx
mov    0x60(%rsp),%rdx
mov    %ecx,(%rdx,%r15,4)
movslq %ecx,%rcx
lea    (%rbx,%rcx,1),%rdx
movsd  (%r12,%rdx,8),%xmm0
ucomisd 0xd41(%rip),%xmm0
jne    12ef <dgetf2_3l+0x1af>
jnp    13a0 <dgetf2_3l+0x260>
test   %eax,%eax
je     130c <dgetf2_3l+0x1cc>
lea    (%r12,%r15,8),%rsi
lea    (%r12,%rcx,8),%rcx
mov    0x4(%rsp),%edi
mov    0x18(%rsp),%rdx
mov    %edx,%r8d
call   1030 <dswap_3l@plt>
cmp    0x38(%rsp),%r15
jae    1200 <dgetf2_3l+0xc0>
movsd  0x0(%r13),%xmm1
movapd %xmm1,%xmm0
xorpd  0xcd7(%rip),%xmm0
xorpd  %xmm2,%xmm2
cmpltsd %xmm1,%xmm2
movapd %xmm2,%xmm3
andpd  %xmm1,%xmm3
andnpd %xmm0,%xmm2
orpd   %xmm3,%xmm2
ucomisd 0xcc6(%rip),%xmm2
jae    11d3 <dgetf2_3l+0x93>
lea    0x1(%r15),%rax
cmp    %rbp,%rax
jae    1200 <dgetf2_3l+0xc0>
mov    %r15,%rcx
not    %rcx
lea    (%rcx,%rbp,1),%rax
mov    0x8(%rsp),%rdx
cmp    $0x2,%rax
jb     1430 <dgetf2_3l+0x2f0>
lea    (%rcx,%rbp,1),%rdx
add    $0xfffffffffffffffe,%rdx
mov    %rdx,%rcx
shr    %rcx
add    $0x1,%rcx
movapd %xmm1,%xmm0
unpcklpd %xmm1,%xmm0
cmp    $0x2,%rdx
mov    0x28(%rsp),%rdi
jae    13be <dgetf2_3l+0x27e>
xor    %esi,%esi
jmp    13f8 <dgetf2_3l+0x2b8>
nop
mov    0x10(%rsp),%rax
cmpl   $0x0,(%rax)
jne    1200 <dgetf2_3l+0xc0>
lea    0x1(%r15),%eax
mov    0x10(%rsp),%rcx
mov    %eax,(%rcx)
jmp    1200 <dgetf2_3l+0xc0>
mov    0x20(%rsp),%rdx
shr    %rdx
add    $0x1,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %esi,%esi
movupd -0x10(%rdi,%rsi,8),%xmm2
movupd (%rdi,%rsi,8),%xmm3
divpd  %xmm0,%xmm2
movupd %xmm2,-0x10(%rdi,%rsi,8)
divpd  %xmm0,%xmm3
movupd %xmm3,(%rdi,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdx
jne    13d0 <dgetf2_3l+0x290>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
test   $0x1,%cl
je     141c <dgetf2_3l+0x2dc>
add    0x8(%rsp),%rsi
add    %rbx,%rsi
movupd (%r12,%rsi,8),%xmm2
divpd  %xmm0,%xmm2
movupd %xmm2,(%r12,%rsi,8)
cmp    %rdx,%rax
je     1200 <dgetf2_3l+0xc0>
add    0x8(%rsp),%rdx
nopw   0x0(%rax,%rax,1)
movsd  (%r14,%rdx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r14,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rbp
jne    1430 <dgetf2_3l+0x2f0>
jmp    1200 <dgetf2_3l+0xc0>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @dgetf2_3l(i32 %0, i32 %1, double* %2, i32 %3, i32* %4, i32* %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32*
%12 = alloca i32*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca i32
store i32 %0, i32* %7
store i32 %1, i32* %8
store double* %2, double** %9
store i32 %3, i32* %10
store i32* %4, i32** %11
store i32* %5, i32** %12
%18 = load i32, i32* %7
%19 = icmp sle i32 %18, 0
br i1 %19, label %23, label %20
20: ; preds = %6
%21 = load i32, i32* %8
%22 = icmp sle i32 %21, 0
br i1 %22, label %23, label %24
23: ; preds = %20, %6
br label %210
24: ; preds = %20
%25 = load i32, i32* %7
%26 = load i32, i32* %8
%27 = icmp slt i32 %25, %26
br i1 %27, label %28, label %30
28: ; preds = %24
%29 = load i32, i32* %7
br label %32
30: ; preds = %24
%31 = load i32, i32* %8
br label %32
32: ; preds = %30, %28
%33 = phi i32 [ %29, %28 ], [ %31, %30 ]
store i32 %33, i32* %17
store i32 0, i32* %14
br label %34
34: ; preds = %206, %32
%35 = load i32, i32* %14
%36 = load i32, i32* %17
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %209
38: ; preds = %34
%39 = load i32, i32* %14
%40 = load i32, i32* %7
%41 = load i32, i32* %14
%42 = sub nsw i32 %40, %41
%43 = load double*, double** %9
%44 = load i32, i32* %14
%45 = load i32, i32* %10
%46 = load i32, i32* %14
%47 = mul nsw i32 %45, %46
%48 = add nsw i32 %44, %47
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %43, i64 %49
%51 = call i32 @idamax_3l(i32 %42, double* %50)
%52 = add nsw i32 %39, %51
store i32 %52, i32* %15
%53 = load i32, i32* %15
%54 = load i32*, i32** %11
%55 = load i32, i32* %14
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds i32, i32* %54, i64 %56
store i32 %53, i32* %57
%58 = load double*, double** %9
%59 = load i32, i32* %15
%60 = load i32, i32* %10
%61 = load i32, i32* %14
%62 = mul nsw i32 %60, %61
%63 = add nsw i32 %59, %62
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %58, i64 %64
%66 = load double, double* %65
%67 = fcmp une double %66, 0.000000e+00
br i1 %67, label %68, label %152
68: ; preds = %38
%69 = load i32, i32* %15
%70 = load i32, i32* %14
%71 = icmp ne i32 %69, %70
br i1 %71, label %72, label %85
72: ; preds = %68
%73 = load i32, i32* %8
%74 = load double*, double** %9
%75 = load i32, i32* %14
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
%78 = load i32, i32* %10
%79 = load double*, double** %9
%80 = load i32, i32* %15
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %79, i64 %81
%83 = load i32, i32* %10
%84 = call i32 @dswap_3l(i32 %73, double* %77, i32 %78, double* %82, i32 %83)
br label %85
85: ; preds = %72, %68
%86 = load i32, i32* %14
%87 = load i32, i32* %7
%88 = sub nsw i32 %87, 1
%89 = icmp slt i32 %86, %88
br i1 %89, label %90, label %151
90: ; preds = %85
%91 = load double*, double** %9
%92 = load i32, i32* %14
%93 = load i32, i32* %10
%94 = load i32, i32* %14
%95 = mul nsw i32 %93, %94
%96 = add nsw i32 %92, %95
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %91, i64 %97
%99 = load double, double* %98
store double %99, double* %16
%100 = load double, double* %16
%101 = fcmp ogt double %100, 0.000000e+00
br i1 %101, label %102, label %104
102: ; preds = %90
%103 = load double, double* %16
br label %107
104: ; preds = %90
%105 = load double, double* %16
%106 = fneg double %105
br label %107
107: ; preds = %104, %102
%108 = phi double [ %103, %102 ], [ %106, %104 ]
%109 = fcmp oge double %108, 2.220000e-16
br i1 %109, label %110, label %127
110: ; preds = %107
%111 = load i32, i32* %7
%112 = load i32, i32* %14
%113 = sub nsw i32 %111, %112
%114 = sub nsw i32 %113, 1
%115 = load double, double* %16
%116 = fdiv double 1.000000e+00, %115
%117 = load double*, double** %9
%118 = load i32, i32* %14
%119 = add nsw i32 %118, 1
%120 = load i32, i32* %10
%121 = load i32, i32* %14
%122 = mul nsw i32 %120, %121
%123 = add nsw i32 %119, %122
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %117, i64 %124
%126 = call i32 @dscal_3l(i32 %114, double %116, double* %125)
br label %150
127: ; preds = %107
%128 = load i32, i32* %14
%129 = add nsw i32 %128, 1
store i32 %129, i32* %13
br label %130
130: ; preds = %146, %127
%131 = load i32, i32* %13
%132 = load i32, i32* %7
%133 = icmp slt i32 %131, %132
br i1 %133, label %134, label %149
134: ; preds = %130
%135 = load double, double* %16
%136 = load double*, double** %9
%137 = load i32, i32* %13
%138 = load i32, i32* %10
%139 = load i32, i32* %14
%140 = mul nsw i32 %138, %139
%141 = add nsw i32 %137, %140
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %136, i64 %142
%144 = load double, double* %143
%145 = fdiv double %144, %135
store double %145, double* %143
br label %146
146: ; preds = %134
%147 = load i32, i32* %13
%148 = add nsw i32 %147, 1
store i32 %148, i32* %13
br label %130
149: ; preds = %130
br label %150
150: ; preds = %149, %110
br label %151
151: ; preds = %150, %85
br label %161
152: ; preds = %38
%153 = load i32*, i32** %12
%154 = load i32, i32* %153
%155 = icmp eq i32 %154, 0
br i1 %155, label %156, label %160
156: ; preds = %152
%157 = load i32, i32* %14
%158 = add nsw i32 %157, 1
%159 = load i32*, i32** %12
store i32 %158, i32* %159
br label %160
160: ; preds = %156, %152
br label %161
161: ; preds = %160, %151
%162 = load i32, i32* %14
%163 = load i32, i32* %17
%164 = icmp slt i32 %162, %163
br i1 %164, label %165, label %205
165: ; preds = %161
%166 = load i32, i32* %7
%167 = load i32, i32* %14
%168 = sub nsw i32 %166, %167
%169 = sub nsw i32 %168, 1
%170 = load i32, i32* %8
%171 = load i32, i32* %14
%172 = sub nsw i32 %170, %171
%173 = sub nsw i32 %172, 1
%174 = load double*, double** %9
%175 = load i32, i32* %14
%176 = add nsw i32 %175, 1
%177 = load i32, i32* %10
%178 = load i32, i32* %14
%179 = mul nsw i32 %177, %178
%180 = add nsw i32 %176, %179
%181 = sext i32 %180 to i64
%182 = getelementptr inbounds double, double* %174, i64 %181
%183 = load double*, double** %9
%184 = load i32, i32* %14
%185 = load i32, i32* %10
%186 = load i32, i32* %14
%187 = add nsw i32 %186, 1
%188 = mul nsw i32 %185, %187
%189 = add nsw i32 %184, %188
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %183, i64 %190
%192 = load i32, i32* %10
%193 = load double*, double** %9
%194 = load i32, i32* %14
%195 = add nsw i32 %194, 1
%196 = load i32, i32* %10
%197 = load i32, i32* %14
%198 = add nsw i32 %197, 1
%199 = mul nsw i32 %196, %198
%200 = add nsw i32 %195, %199
%201 = sext i32 %200 to i64
%202 = getelementptr inbounds double, double* %193, i64 %201
%203 = load i32, i32* %10
%204 = call i32 @dger_3l(i32 %169, i32 %173, double -1.000000e+00, double* %182, i32 1, double* %191, i32 %192, double* %202, i32 %203)
br label %205
205: ; preds = %165, %161
br label %206
206: ; preds = %205
%207 = load i32, i32* %14
%208 = add nsw i32 %207, 1
store i32 %208, i32* %14
br label %34
209: ; preds = %34
br label %210
210: ; preds = %209, %23
ret void
}

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