#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*/ blasDAXPY(int, double, double *, int, double *, int);
int /*<<< orphan*/ blasDSCAL(int, double, double *, int);
int linpackDTRDI(double t[], int ldt, int n) {
  int j, k, info;
  double temp;
  info = 0;
  for (k = n; 1 <= k; k--) {
    if (t[k - 1 + (k - 1) * ldt] == 0.0) {
      info = k;
      break;
    }
    t[k - 1 + (k - 1) * ldt] = 1.0 / t[k - 1 + (k - 1) * ldt];
    temp = -t[k - 1 + (k - 1) * ldt];
    if (k != n) {
      blasDSCAL(n - k, temp, t + k + (k - 1) * ldt, 1);
    }
    for (j = 1; j <= k - 1; j++) {
      temp = t[k - 1 + (j - 1) * ldt];
      t[k - 1 + (j - 1) * ldt] = 0.0;
      blasDAXPY(n - k + 1, temp, t + k - 1 + (k - 1) * ldt, 1,
                t + k - 1 + (j - 1) * ldt, 1);
    }
  }
  return (info);
}


x86 assembly (clang, -O0):
<linpackDTRDI>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    $0x1,%eax
cmp    -0x18(%rbp),%eax
jg     132f <linpackDTRDI+0x20f>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1189 <linpackDTRDI+0x69>
jp     1189 <linpackDTRDI+0x69>
mov    -0x18(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    132f <linpackDTRDI+0x20f>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xe56(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
je     124a <linpackDTRDI+0x12a>
mov    -0x10(%rbp),%edi
sub    -0x18(%rbp),%edi
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rsi
movslq -0x18(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x18(%rbp),%eax
sub    $0x1,%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    $0x1,%edx
call   1040 <blasDSCAL@plt>
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     131c <linpackDTRDI+0x1fc>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x14(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x14(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%edi
sub    -0x18(%rbp),%edi
add    $0x1,%edi
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rsi
movslq -0x18(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
add    $0xfffffffffffffff8,%rsi
mov    -0x18(%rbp),%eax
sub    $0x1,%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x8(%rbp),%rcx
movslq -0x18(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rcx
add    $0xfffffffffffffff8,%rcx
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rcx
mov    $0x1,%r8d
mov    %r8d,%edx
call   1030 <blasDAXPY@plt>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1251 <linpackDTRDI+0x131>
jmp    1321 <linpackDTRDI+0x201>
mov    -0x18(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x18(%rbp)
jmp    113f <linpackDTRDI+0x1f>
mov    -0x1c(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<linpackDTRDI>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
xor    %eax,%eax
test   %edx,%edx
jle    1289 <linpackDTRDI+0x169>
mov    %esi,%r14d
mov    %rdi,%rbp
lea    0x1(%rdx),%eax
mov    %eax,0x8(%rsp)
movslq %esi,%rax
mov    %rdx,0x10(%rsp)
mov    %edx,%r13d
lea    -0x8(%rdi),%rcx
mov    %rcx,0x18(%rsp)
lea    -0x8(%rdi,%r13,8),%rcx
mov    %rax,0x30(%rsp)
lea    0x0(,%rax,8),%r15
mov    %r13,0x28(%rsp)
mov    %esi,0xc(%rsp)
mov    %rdi,0x20(%rsp)
nopl   0x0(%rax)
lea    -0x1(%r13),%rbx
mov    %ebx,%eax
imul   %r14d,%eax
cltq
add    %r13,%rax
add    $0xffffffffffffffff,%rax
movsd  0x0(%rbp,%rax,8),%xmm1
ucomisd 0xe67(%rip),%xmm1
jne    11a9 <linpackDTRDI+0x89>
jnp    1286 <linpackDTRDI+0x166>
mov    %rbx,0x40(%rsp)
imul   0x30(%rsp),%rbx
movsd  0xe44(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
cmp    0x28(%rsp),%r13
je     11fe <linpackDTRDI+0xde>
xorpd  0xe3b(%rip),%xmm0
lea    0x0(,%r13,8),%rax
add    %rbp,%rax
lea    (%rax,%rbx,8),%rsi
mov    0x10(%rsp),%rax
mov    %eax,%edi
sub    %r13d,%edi
mov    $0x1,%edx
mov    %rcx,%r12
call   1040 <blasDSCAL@plt>
mov    %r12,%rcx
cmp    $0x2,%r13
jb     1282 <linpackDTRDI+0x162>
mov    0x8(%rsp),%ebp
sub    %r13d,%ebp
mov    0x18(%rsp),%rax
lea    (%rax,%r13,8),%rax
lea    (%rax,%rbx,8),%rbx
mov    $0x1,%r12d
mov    %rcx,0x38(%rsp)
mov    %rcx,%r14
cs nopw 0x0(%rax,%rax,1)
movsd  (%r14),%xmm0
movq   $0x0,(%r14)
mov    %ebp,%edi
mov    %rbx,%rsi
mov    $0x1,%edx
mov    %r14,%rcx
mov    $0x1,%r8d
call   1030 <blasDAXPY@plt>
add    $0x1,%r12
add    %r15,%r14
cmp    %r12,%r13
jne    1230 <linpackDTRDI+0x110>
mov    0x38(%rsp),%rcx
add    $0xfffffffffffffff8,%rcx
cmp    $0x2,%r13
mov    0x40(%rsp),%r13
mov    0xc(%rsp),%r14d
mov    0x20(%rsp),%rbp
jae    1180 <linpackDTRDI+0x60>
xor    %eax,%eax
jmp    1289 <linpackDTRDI+0x169>
mov    %r13d,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @linpackDTRDI(double* %0, i32 %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %53
5: ; preds = %3
%6 = add nuw i32 %2, 1
%7 = sext i32 %1 to i64
%8 = zext i32 %2 to i64
%9 = zext i32 %2 to i64
%10 = getelementptr inbounds double, double* %0, i64 -1
br label %12
11: ; preds = %40
br i1 %35, label %12, label %53
12: ; preds = %5, %11
%13 = phi i64 [ %8, %5 ], [ %15, %11 ]
%14 = trunc i64 %13 to i32
%15 = add nsw i64 %13, -1
%16 = trunc i64 %15 to i32
%17 = mul nsw i64 %15, %7
%18 = mul nsw i32 %16, %1
%19 = sext i32 %18 to i64
%20 = add nsw i64 %15, %19
%21 = getelementptr inbounds double, double* %0, i64 %20
%22 = load double, double* %21
%23 = fcmp oeq double %22, 0.000000e+00
br i1 %23, label %51, label %24
24: ; preds = %12
%25 = fdiv double 1.000000e+00, %22
store double %25, double* %21
%26 = icmp eq i64 %13, %9
br i1 %26, label %34, label %27
27: ; preds = %24
%28 = fneg double %25
%29 = getelementptr inbounds double, double* %0, i64 %13
%30 = getelementptr inbounds double, double* %29, i64 %17
%31 = trunc i64 %13 to i32
%32 = sub i32 %2, %31
%33 = tail call i32 @blasDSCAL(i32 %32, double %28, double* %30, i32 1)
br label %34
34: ; preds = %27, %24
%35 = icmp ugt i64 %13, 1
br i1 %35, label %36, label %53
36: ; preds = %34
%37 = sub i32 %6, %14
%38 = getelementptr inbounds double, double* %10, i64 %13
%39 = getelementptr inbounds double, double* %38, i64 %17
br label %40
40: ; preds = %36, %40
%41 = phi i64 [ 1, %36 ], [ %49, %40 ]
%42 = add nsw i64 %41, -1
%43 = mul nsw i64 %42, %7
%44 = add nsw i64 %43, %15
%45 = getelementptr inbounds double, double* %0, i64 %44
%46 = load double, double* %45
store double 0.000000e+00, double* %45
%47 = getelementptr inbounds double, double* %38, i64 %43
%48 = tail call i32 @blasDAXPY(i32 %37, double %46, double* %39, i32 1, double* %47, i32 1)
%49 = add nuw nsw i64 %41, 1
%50 = icmp eq i64 %49, %13
br i1 %50, label %11, label %40
51: ; preds = %12
%52 = trunc i64 %13 to i32
br label %53
53: ; preds = %34, %11, %51, %3
%54 = phi i32 [ 0, %3 ], [ %52, %51 ], [ 0, %11 ], [ 0, %34 ]
ret i32 %54
}