#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__ calloc(int, int);
scalar_t__ dot_real(double *, double *, int);
double fabs(double);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ pflslv(int, double **, double *, double *);
int /*<<< orphan*/ printf(char *, int);
int pfefct(int neqns, double **penv, double *diag) {
  double *ptenv;
  int iband, i, jj, ifirst;
  double *work;
  for (i = 1; i < neqns; i++) {
    ptenv = penv[i];
    iband = penv[i + 1] - ptenv;
    work = (double *)calloc(iband, sizeof(double));
    if (iband > 0) {
      ifirst = i - iband;
      pflslv(iband, penv + ifirst, diag + ifirst, ptenv);
      for (jj = 0; jj < iband; jj++) {
        work[jj] = ptenv[jj];
        ptenv[jj] = ptenv[jj] / diag[i + jj - iband];
      }
      diag[i] -= dot_real(ptenv, work, iband);
    }
    free(work);
    if (fabs(diag[i]) < 1.0e-60) {
      printf("!!! pfefct(): diagonal %d is zero !!!\n", i);
      return (1);
    }
  }
  return (0);
}


x86 assembly (clang, -O0):
<pfefct>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12e9 <pfefct+0x199>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
sub    %rcx,%rax
sar    $0x3,%rax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%edi
mov    $0x8,%esi
call   1060 <calloc@plt>
mov    %rax,-0x38(%rbp)
cmpl   $0x0,-0x24(%rbp)
jle    1288 <pfefct+0x138>
mov    -0x28(%rbp),%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x24(%rbp),%edi
mov    -0x10(%rbp),%rsi
movslq -0x30(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x18(%rbp),%rdx
movslq -0x30(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdx
mov    -0x20(%rbp),%rcx
call   1040 <pflslv@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    125d <pfefct+0x10d>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x2c(%rbp),%ecx
sub    -0x24(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11fa <pfefct+0xaa>
mov    -0x20(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x24(%rbp),%edx
call   1070 <dot_real@plt>
cvtsi2sd %rax,%xmm1
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xd6b(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xd4f(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    12d6 <pfefct+0x186>
mov    -0x28(%rbp),%esi
lea    0xd5b(%rip),%rdi
call   1050 <printf@plt>
movl   $0x1,-0x4(%rbp)
jmp    12f0 <pfefct+0x1a0>
jmp    12db <pfefct+0x18b>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    116a <pfefct+0x1a>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pfefct>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rdx,%r13
mov    %rsi,%r14
cmp    $0x2,%edi
mov    $0x1,%eax
cmovge %edi,%eax
mov    %rax,0x10(%rsp)
mov    $0x1,%ebp
movl   $0x0,0xc(%rsp)
jmp    11e5 <pfefct+0x95>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r15,%rdi
mov    %rbx,%rsi
mov    %ebp,%edx
call   1070 <dot_real@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0x0(%r13,%r12,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%r12,8)
lea    0x1(%r12),%rbp
mov    %rbx,%rdi
call   1030 <free@plt>
movsd  0x0(%r13,%r12,8),%xmm0
andpd  0xe2d(%rip),%xmm0
movsd  0xe35(%rip),%xmm1
ucomisd %xmm0,%xmm1
ja     1366 <pfefct+0x216>
cmp    0x10(%rsp),%rbp
je     1362 <pfefct+0x212>
mov    %rbp,%r12
addl   $0x1,0xc(%rsp)
mov    (%r14,%rbp,8),%r15
mov    0x8(%r14,%rbp,8),%rbp
sub    %r15,%rbp
shr    $0x3,%rbp
mov    %ebp,%edi
mov    $0x8,%esi
call   1060 <calloc@plt>
mov    %rax,%rbx
test   %ebp,%ebp
jle    11b7 <pfefct+0x67>
mov    %r12d,%eax
sub    %ebp,%eax
cltq
lea    (%r14,%rax,8),%rsi
lea    0x0(,%rax,8),%rdx
add    %r13,%rdx
mov    %ebp,%edi
mov    %r15,%rcx
call   1040 <pflslv@plt>
mov    %ebp,%eax
cmp    $0x2,%rax
jb     1320 <pfefct+0x1d0>
lea    -0x1(%rax),%rcx
mov    0xc(%rsp),%edx
sub    %ebp,%edx
lea    (%rcx,%rdx,1),%esi
cmp    %edx,%esi
jl     1320 <pfefct+0x1d0>
shr    $0x20,%rcx
mov    $0x0,%ecx
jne    1322 <pfefct+0x1d2>
lea    (%rbx,%rax,8),%rsi
lea    (%r15,%rax,8),%r10
movslq %edx,%rcx
lea    0x0(,%rcx,8),%rdx
add    %r13,%rdx
add    %rax,%rcx
lea    0x0(,%rcx,8),%rdi
add    %r13,%rdi
cmp    %r10,%rbx
setb   %r11b
cmp    %rsi,%r15
setb   %r8b
cmp    %rdi,%rbx
setb   %cl
cmp    %rsi,%rdx
setb   %sil
cmp    %rdi,%r15
setb   %r9b
cmp    %r10,%rdx
setb   %dl
test   %r8b,%r11b
jne    1320 <pfefct+0x1d0>
and    %sil,%cl
jne    1320 <pfefct+0x1d0>
mov    $0x0,%ecx
and    %dl,%r9b
jne    1322 <pfefct+0x1d2>
mov    %ebp,%r8d
and    $0x1,%r8d
mov    %rax,%rcx
sub    %r8,%rcx
mov    %r12d,%esi
sub    %ebp,%esi
xor    %edi,%edi
nop
movupd (%r15,%rdi,8),%xmm0
movupd %xmm0,(%rbx,%rdi,8)
lea    (%rsi,%rdi,1),%edx
movslq %edx,%rdx
movupd 0x0(%r13,%rdx,8),%xmm1
divpd  %xmm1,%xmm0
movupd %xmm0,(%r15,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rcx
jne    12e0 <pfefct+0x190>
test   %r8,%r8
je     1190 <pfefct+0x40>
jmp    1322 <pfefct+0x1d2>
cs nopw 0x0(%rax,%rax,1)
xor    %ecx,%ecx
mov    %r12d,%edx
sub    %ebp,%edx
nopw   0x0(%rax,%rax,1)
movsd  (%r15,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
movsd  (%r15,%rcx,8),%xmm0
lea    (%rdx,%rcx,1),%esi
movslq %esi,%rsi
divsd  0x0(%r13,%rsi,8),%xmm0
movsd  %xmm0,(%r15,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1330 <pfefct+0x1e0>
jmp    1190 <pfefct+0x40>
xor    %eax,%eax
jmp    137a <pfefct+0x22a>
lea    0xcab(%rip),%rdi
mov    %r12d,%esi
call   1050 <printf@plt>
mov    $0x1,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [39 x i8] c"!!! pfefct(): diagonal %d is zero !!!\0A\00"
define i32 @pfefct(i32 %0, double** %1, double* %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca double**
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double*
store i32 %0, i32* %5
store double** %1, double*** %6
store double* %2, double** %7
store i32 1, i32* %10
br label %14
14: ; preds = %118, %3
%15 = load i32, i32* %10
%16 = load i32, i32* %5
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %121
18: ; preds = %14
%19 = load double**, double*** %6
%20 = load i32, i32* %10
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double*, double** %19, i64 %21
%23 = load double*, double** %22
store double* %23, double** %8
%24 = load double**, double*** %6
%25 = load i32, i32* %10
%26 = add nsw i32 %25, 1
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double*, double** %24, i64 %27
%29 = load double*, double** %28
%30 = load double*, double** %8
%31 = ptrtoint double* %29 to i64
%32 = ptrtoint double* %30 to i64
%33 = sub i64 %31, %32
%34 = sdiv exact i64 %33, 8
%35 = trunc i64 %34 to i32
store i32 %35, i32* %9
%36 = load i32, i32* %9
%37 = call i64 @calloc(i32 %36, i32 8)
%38 = inttoptr i64 %37 to double*
store double* %38, double** %13
%39 = load i32, i32* %9
%40 = icmp sgt i32 %39, 0
br i1 %40, label %41, label %104
41: ; preds = %18
%42 = load i32, i32* %10
%43 = load i32, i32* %9
%44 = sub nsw i32 %42, %43
store i32 %44, i32* %12
%45 = load i32, i32* %9
%46 = load double**, double*** %6
%47 = load i32, i32* %12
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double*, double** %46, i64 %48
%50 = load double*, double** %7
%51 = load i32, i32* %12
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %50, i64 %52
%54 = load double*, double** %8
%55 = call i32 @pflslv(i32 %45, double** %49, double* %53, double* %54)
store i32 0, i32* %11
br label %56
56: ; preds = %89, %41
%57 = load i32, i32* %11
%58 = load i32, i32* %9
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %92
60: ; preds = %56
%61 = load double*, double** %8
%62 = load i32, i32* %11
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %61, i64 %63
%65 = load double, double* %64
%66 = load double*, double** %13
%67 = load i32, i32* %11
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %66, i64 %68
store double %65, double* %69
%70 = load double*, double** %8
%71 = load i32, i32* %11
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %70, i64 %72
%74 = load double, double* %73
%75 = load double*, double** %7
%76 = load i32, i32* %10
%77 = load i32, i32* %11
%78 = add nsw i32 %76, %77
%79 = load i32, i32* %9
%80 = sub nsw i32 %78, %79
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %75, i64 %81
%83 = load double, double* %82
%84 = fdiv double %74, %83
%85 = load double*, double** %8
%86 = load i32, i32* %11
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %85, i64 %87
store double %84, double* %88
br label %89
89: ; preds = %60
%90 = load i32, i32* %11
%91 = add nsw i32 %90, 1
store i32 %91, i32* %11
br label %56
92: ; preds = %56
%93 = load double*, double** %8
%94 = load double*, double** %13
%95 = load i32, i32* %9
%96 = call i64 @dot_real(double* %93, double* %94, i32 %95)
%97 = sitofp i64 %96 to double
%98 = load double*, double** %7
%99 = load i32, i32* %10
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %98, i64 %100
%102 = load double, double* %101
%103 = fsub double %102, %97
store double %103, double* %101
br label %104
104: ; preds = %92, %18
%105 = load double*, double** %13
%106 = call i32 @free(double* %105)
%107 = load double*, double** %7
%108 = load i32, i32* %10
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %107, i64 %109
%111 = load double, double* %110
%112 = call double @llvm.fabs.f64(double %111)
%113 = fcmp olt double %112, 0x3379B604AAACA626
br i1 %113, label %114, label %117
114: ; preds = %104
%115 = load i32, i32* %10
%116 = call i32 @printf(i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str, i64 0, i64 0), i32 %115)
store i32 1, i32* %4
br label %122
117: ; preds = %104
br label %118
118: ; preds = %117
%119 = load i32, i32* %10
%120 = add nsw i32 %119, 1
store i32 %120, i32* %10
br label %14
121: ; preds = %14
store i32 0, i32* %4
br label %122
122: ; preds = %121, %114
%123 = load i32, i32* %4
ret i32 %123
}

LLVM IR (-O2):
@.str = [39 x i8] c"!!! pfefct(): diagonal %d is zero !!!\0A\00"
define i32 @pfefct(i32 %0, double** %1, double* %2) {
%4 = call i32 @llvm.smax.i32(i32 %0, i32 1)
%5 = zext i32 %4 to i64
br label %6
6: ; preds = %109, %3
%7 = phi i32 [ %115, %109 ], [ 0, %3 ]
%8 = phi i64 [ %15, %109 ], [ 1, %3 ]
%9 = add i32 %7, 1
%10 = add i32 %7, 1
%11 = icmp eq i64 %8, %5
br i1 %11, label %119, label %12
12: ; preds = %6
%13 = getelementptr inbounds double*, double** %1, i64 %8
%14 = load double*, double** %13
%15 = add nuw nsw i64 %8, 1
%16 = getelementptr inbounds double*, double** %1, i64 %15
%17 = load double*, double** %16
%18 = ptrtoint double* %17 to i64
%19 = ptrtoint double* %14 to i64
%20 = sub i64 %18, %19
%21 = lshr exact i64 %20, 3
%22 = trunc i64 %21 to i32
%23 = tail call i64 @calloc(i32 %22, i32 8)
%24 = inttoptr i64 %23 to double*
%25 = inttoptr i64 %23 to i8*
%26 = icmp sgt i32 %22, 0
br i1 %26, label %27, label %109
27: ; preds = %12
%28 = trunc i64 %8 to i32
%29 = sub i32 %28, %22
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double*, double** %1, i64 %30
%32 = getelementptr inbounds double, double* %2, i64 %30
%33 = tail call i32 @pflslv(i32 %22, double** %31, double* %32, double* %14)
%34 = and i64 %21, 4294967295
%35 = icmp ult i64 %34, 2
br i1 %35, label %87, label %36
36: ; preds = %27
%37 = add nsw i64 %34, -1
%38 = sub i32 %10, %22
%39 = trunc i64 %37 to i32
%40 = add i32 %38, %39
%41 = icmp slt i32 %40, %38
%42 = icmp ugt i64 %37, 4294967295
%43 = or i1 %41, %42
br i1 %43, label %87, label %44
44: ; preds = %36
%45 = getelementptr double, double* %24, i64 %34
%46 = getelementptr double, double* %14, i64 %34
%47 = bitcast double* %46 to i8*
%48 = sub i32 %9, %22
%49 = sext i32 %48 to i64
%50 = getelementptr double, double* %2, i64 %49
%51 = add nsw i64 %34, %49
%52 = getelementptr double, double* %2, i64 %51
%53 = bitcast double* %52 to i8*
%54 = icmp ult i8* %25, %47
%55 = icmp ult double* %14, %45
%56 = and i1 %54, %55
%57 = icmp ult i8* %25, %53
%58 = icmp ult double* %50, %45
%59 = and i1 %57, %58
%60 = or i1 %56, %59
%61 = icmp ult double* %14, %52
%62 = icmp ult double* %50, %46
%63 = and i1 %61, %62
%64 = or i1 %60, %63
br i1 %64, label %87, label %65
65: ; preds = %44
%66 = and i64 %21, 1
%67 = sub nsw i64 %34, %66
br label %68
68: ; preds = %68, %65
%69 = phi i64 [ 0, %65 ], [ %83, %68 ]
%70 = trunc i64 %69 to i32
%71 = getelementptr inbounds double, double* %14, i64 %69
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = getelementptr inbounds double, double* %24, i64 %69
%75 = bitcast double* %74 to <2 x double>*
store <2 x double> %73, <2 x double>* %75
%76 = add i32 %29, %70
%77 = sext i32 %76 to i64
%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> %73, %80
%82 = bitcast double* %71 to <2 x double>*
store <2 x double> %81, <2 x double>* %82
%83 = add nuw i64 %69, 2
%84 = icmp eq i64 %83, %67
br i1 %84, label %85, label %68
85: ; preds = %68
%86 = icmp eq i64 %66, 0
br i1 %86, label %103, label %87
87: ; preds = %44, %36, %27, %85
%88 = phi i64 [ 0, %44 ], [ 0, %36 ], [ 0, %27 ], [ %67, %85 ]
br label %89
89: ; preds = %87, %89
%90 = phi i64 [ %101, %89 ], [ %88, %87 ]
%91 = getelementptr inbounds double, double* %14, i64 %90
%92 = load double, double* %91
%93 = getelementptr inbounds double, double* %24, i64 %90
store double %92, double* %93
%94 = load double, double* %91
%95 = trunc i64 %90 to i32
%96 = add i32 %29, %95
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %2, i64 %97
%99 = load double, double* %98
%100 = fdiv double %94, %99
store double %100, double* %91
%101 = add nuw nsw i64 %90, 1
%102 = icmp eq i64 %101, %34
br i1 %102, label %103, label %89
103: ; preds = %89, %85
%104 = tail call i64 @dot_real(double* %14, double* %24, i32 %22)
%105 = sitofp i64 %104 to double
%106 = getelementptr inbounds double, double* %2, i64 %8
%107 = load double, double* %106
%108 = fsub double %107, %105
store double %108, double* %106
br label %109
109: ; preds = %103, %12
%110 = tail call i32 @free(double* %24)
%111 = getelementptr inbounds double, double* %2, i64 %8
%112 = load double, double* %111
%113 = tail call double @llvm.fabs.f64(double %112)
%114 = fcmp olt double %113, 0x3379B604AAACA626
%115 = add i32 %7, 1
br i1 %114, label %116, label %6
116: ; preds = %109
%117 = trunc i64 %8 to i32
%118 = tail call i32 @printf(i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str, i64 0, i64 0), i32 %117)
br label %119
119: ; preds = %6, %116
%120 = phi i32 [ 1, %116 ], [ 0, %6 ]
ret i32 %120
}